Showing 1–50 of 100 results for author: Nara, Y

Enhanced Dilepton production near the color superconducting phase and the QCD critical point

Authors: Toru Nishimura, Yasushi Nara, Jan Steinheimer

Abstract: The dilepton production yields in relativistic heavy ion collisions are investigated along isentropic trajectories in the quark (Wigner) phase within the two-flavor Nambu-Jona-Lasinio model. An enhancement of the ultra-low energy dilepton yield in the vicinity of the color superconducting (CSC) phase and the QCD critical point (QCD-CP) is found, compared to the free quark gas. Furthermore, we have… ▽ More The dilepton production yields in relativistic heavy ion collisions are investigated along isentropic trajectories in the quark (Wigner) phase within the two-flavor Nambu-Jona-Lasinio model. An enhancement of the ultra-low energy dilepton yield in the vicinity of the color superconducting (CSC) phase and the QCD critical point (QCD-CP) is found, compared to the free quark gas. Furthermore, we have found a nontrivial structure in the beam energy dependence of the ultra-low energy dilepton yield. A local maximum and minimum of the dilepton yield as a function of entropy per baryon emerge when the trajectories are close to both locations of the CSC phase transition line and the QCD-CP. Only the maximum appears in the scenario without the CSC phase but with the QCD-CP. On the other hand, when only the CSC phase is considered, the dilepton yield monotonically increases as the beam energy decreases. These distinctive patterns could potentially serve as signals of the CSC phase and QCD-CP. In addition, it is found that the dilepton production yield and the location of the minimum strongly depend on the value of diquark coupling, suggesting the possibility that the value of the diquark coupling may be extracted from experimental data. △ Less

Submitted 23 November, 2023; originally announced November 2023.

Study of bulk properties of the system formed in U+U collisions at $\sqrt{s_{\mathrm NN}}$ =~2.12~GeV using JAM model

Authors: Aswini Kumar Sahoo, Xionghong He, Yasushi Nara, Subhash Singha

Abstract: The Lanzhou Cooling-Storage-Ring facility is set to conduct experiments involving Uranium-Uranium collisions at the center of mass energies ranging from 2.12 to 2.4 GeV. Our investigation is focused on various bulk observables, which include charged particle multiplicity ($N_{\text{ch}}$), average transverse momentum ($\langle p_{\text{T}}\rangle$), initial eccentricity ($ε_{n}$), and flow harmoni… ▽ More The Lanzhou Cooling-Storage-Ring facility is set to conduct experiments involving Uranium-Uranium collisions at the center of mass energies ranging from 2.12 to 2.4 GeV. Our investigation is focused on various bulk observables, which include charged particle multiplicity ($N_{\text{ch}}$), average transverse momentum ($\langle p_{\text{T}}\rangle$), initial eccentricity ($ε_{n}$), and flow harmonics ($v_{n}$), for different orientations of U+U collisions within the range of $0^{\circ} < θ< 120 ^{\circ}$ at $\sqrt{s_{\mathrm NN}} = 2.12$ GeV ($p_{\mathrm lab}$ = 500 MeV). Among the various collision configurations at this energy, the tip-tip scenario emerged with the highest average charged particle multiplicity, denoted as $\langle N_{\text{ch}} \rangle$. Notably, both the second and third-order eccentricities, $ε_{2,3}$, revealed intricate patterns as they varied with impact parameter across distinct configurations. The tip-tip configuration displayed the most pronounced magnitude of rapidity-odd directed flow ($v_{1}$), whereas the body-body configuration exhibited the least pronounced magnitude. Concerning elliptic flow ($v_{2}$) near mid-rapidity ($η< 1.0$), a negative sign is observed for all configurations except for the side-side exhibited a distinctly positive sign. Within the spectrum of configurations, the body-body scenario displayed the highest magnitude of $v_{2}$. For reaction plane correlated triangular flow ($v_{3}$), the body-body configuration emerged with the largest magnitude while the side-side exhibited the smallest magnitude. Our study seeks to establish a fundamental understanding of various U+U collision configurations in preparation for the forthcoming CEE experiment. △ Less

Submitted 1 November, 2023; originally announced November 2023.

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)

Repulsive $Λ$ potentials in dense neutron star matter and binding energy of $Λ$ in hypernuclei

Authors: Asanosuke Jinno, Koichi Murase, Yasushi Nara, Akira Ohnishi

Abstract: The repulsive three-body force between the lambda ($Λ$) hyperon and medium nucleons is a key element in solving the hyperon puzzle in neutron stars. We investigate the binding energies of $Λ$ hyperon in hypernuclei to verify the repulsive $Λ$ potentials from the chiral effective field theory ($χ$EFT) employing the Skyrme Hartree-Fock method. We find that the $χ$EFT $Λ$ potential with the $ΛNN$ thr… ▽ More The repulsive three-body force between the lambda ($Λ$) hyperon and medium nucleons is a key element in solving the hyperon puzzle in neutron stars. We investigate the binding energies of $Λ$ hyperon in hypernuclei to verify the repulsive $Λ$ potentials from the chiral effective field theory ($χ$EFT) employing the Skyrme Hartree-Fock method. We find that the $χ$EFT $Λ$ potential with the $ΛNN$ three-body forces reproduces the existing hypernuclear binding energy data, whereas the $Λ$ binding energies are overestimated without the $ΛNN$ three-body force. Additionally, we search for the parameter space of the $Λ$ potentials by varying the Taylor coefficients of the $Λ$ potential and the effective mass of $Λ$ at the saturation density. Our analysis demonstrates that the parameter region consistent with the $Λ$ binding energy data spans a wide range of the parameter space, including even more repulsive potentials than the $χ$EFT prediction. We confirm that these strong repulsive $Λ$ potentials suppress the presence of $Λ$ in the neutron star matter. We found that the $Λ$ potentials repulsive at high densities are favored when the depth of the $Λ$ potential at the saturation density, $U_Λ(ρ_0)=J_Λ$, is $J_Λ\gtrsim-29~\text{MeV}$, while attractive ones are favored when $J_Λ\lesssim -31~\text{MeV}$. This suggests that the future high-resolution data of hypernuclei could rule out the scenario in which $Λ$s appear through the precise determination of $J_Λ$ within the accuracy of $1~\text{MeV}$. △ Less

Submitted 13 December, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: 15 pages, 11 figures, 3 tables, figures updated and extended, (Published in Physical Review C)

Report number: YITP-23-78, KUNS-2968

Journal ref: Phys. Rev. C 108, 065803 (2023)

A Poincaré covariant cascade method for high-energy nuclear collisions

Authors: Yasushi Nara, Asanosuke Jinno, Tomoyuki Maruyama, Koichi Murase, Akira Ohnishi

Abstract: We present a Poincaré covariant cascade algorithm based on the constrained Hamiltonian dynamics in an $8N$-dimensional phase space to simulate the Boltzmann-type two-body collision term. We compare this covariant cascade algorithm with traditional $6N$-dimensional phase-space cascade algorithms. To validate the covariant cascade algorithm, we perform box calculations. We examine the frame dependen… ▽ More We present a Poincaré covariant cascade algorithm based on the constrained Hamiltonian dynamics in an $8N$-dimensional phase space to simulate the Boltzmann-type two-body collision term. We compare this covariant cascade algorithm with traditional $6N$-dimensional phase-space cascade algorithms. To validate the covariant cascade algorithm, we perform box calculations. We examine the frame dependence of the algorithm in a one-dimensionally expanding system as well as the compression stages of colliding two nuclei. We confirm that our covariant cascade method is reliable to simulate high-energy nuclear collisions. Furthermore, we present Lorentz-covariant equations of motion for the $N$-body system interacting via potentials, which can be efficiently solved numerically. △ Less

Submitted 8 September, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 13 pages, 10 figures, 1 table, typos fixed

Report number: YITP-23-77, KUNS-2969

Journal ref: Phys.Rev.C 108 (2023) 2, 024910

Directed flow of $Λ$ from heavy-ion collisions and hyperon puzzle of neutron stars

Authors: Akira Ohnishi, Asanosuke Jinno, Koichi Murase, Yasushi Nara

Abstract: We examine the $Λ$ potential from the chiral effective field theory ($χ$EFT) via the $Λ$ directed flow from heavy-ion collisions. We implement the $Λ$ potential obtained from the $χ$EFT in a vector potential version of relativistic quantum molecular dynamics. We find that the $Λ$ potentials obtained from the $χ$EFT assuming weak momentum dependence reproduce the $Λ$ directed flow measured by the S… ▽ More We examine the $Λ$ potential from the chiral effective field theory ($χ$EFT) via the $Λ$ directed flow from heavy-ion collisions. We implement the $Λ$ potential obtained from the $χ$EFT in a vector potential version of relativistic quantum molecular dynamics. We find that the $Λ$ potentials obtained from the $χ$EFT assuming weak momentum dependence reproduce the $Λ$ directed flow measured by the STAR collaboration in the Beam Energy Scan program. While the $Λ$ directed flow is not very sensitive to the density dependence of the potential, the directed flow at large rapidities is susceptible to the momentum dependence. Thus understanding the directed flow of hyperons in a wide range of beam energy and rapidity is helpful in understanding hyperon potentials in dense matter. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: Talk given at HYP2022 (the 14th International Conference on Hypernuclear and Strange Particle Physics), June 27-July 1, 2022, Prague, Czech Republic. 6 pages, 2 figures

Report number: YITP-22-122

The high-density equation of state in heavy-ion collisions: Constraints from proton flow

Authors: Jan Steinheimer, Anton Motornenko, Agnieszka Sorensen, Yasushi Nara, Volker Koch, Marcus Bleicher

Abstract: A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $m_T$. We also… ▽ More A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $m_T$. We also show that an increase in $v_2$ is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $n_{\mathrm{CE}}<4 n_0$, show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $4n_0$. △ Less

Submitted 14 October, 2022; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures, version accepted for publication

Journal ref: Eur. Phys. J. C 82, 911 (2022)

Directed flow of $Λ$ in high-energy heavy-ion collisions and $Λ$ potential in dense nuclear matter

Authors: Yasushi Nara, Asanosuke Jinno, Koichi Murase, Akira Ohnishi

Abstract: We investigate the sensitivity of the $Λ$ directed flow to the $Λ$ potential in mid-central Au + Au collisions at $\sqrt{s_{NN}}\approx3.0$--$30$ GeV. The $Λ$ potential obtained from the chiral effective field theory ($χ$EFT) is used in a microscopic transport model, a vector version of relativistic quantum molecular dynamics (RQMDv). We find that the density-dependent $Λ$ potentials, obtained fro… ▽ More We investigate the sensitivity of the $Λ$ directed flow to the $Λ$ potential in mid-central Au + Au collisions at $\sqrt{s_{NN}}\approx3.0$--$30$ GeV. The $Λ$ potential obtained from the chiral effective field theory ($χ$EFT) is used in a microscopic transport model, a vector version of relativistic quantum molecular dynamics (RQMDv). We find that the density-dependent $Λ$ potentials, obtained from the $χ$EFT assuming weak momentum dependence of the potential, reproduce the rapidity and the beam-energy dependence of the $Λ$ directed flow measured by the STAR collaboration in the Beam Energy Scan program. Although the $Λ$ directed flow is insensitive to the density dependence of the potential, it is susceptible to the momentum dependence. We also show that a hydrodynamics picture based on the blast-wave model predicts a similarity of the proton, $Λ$, and $Ξ$ directed flows, but the directed flow of $Ω$ baryons slightly deviates from other baryons. We also show that the quark coalescence predicts different rapidity dependence of the directed flows for hyperons. These investigations suggest that measurements of a wide range of the rapidity dependence of the directed flow of hyperons may provide important information about the properties of hot and dense matter created in high-energy heavy-ion collisions. △ Less

Submitted 29 October, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 11 pages, 9 figures. version accepted for publication

Report number: YITP-22-81, KUNS-2935

Journal ref: Phys. Rev. C106 (2022)4,044902

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

A Chiral Mean-Field Equation-of-State in UrQMD: Effects on the Heavy Ion Compression Stage

Authors: Manjunath Omana Kuttan, Anton Motornenko, Jan Steinheimer, Horst Stoecker, Yasushi Nara, Marcus Bleicher

Abstract: It is shown that the initial compression in central heavy ion collisions at beam energies of $E_\mathrm{lab}=1-10A$~GeV depends dominantly on the underlying equation of state and only marginally on the model used for the dynamical description. To do so, a procedure to incorporate any equation of state in the UrQMD transport model is introduced. In particular we compare the baryon density, temperat… ▽ More It is shown that the initial compression in central heavy ion collisions at beam energies of $E_\mathrm{lab}=1-10A$~GeV depends dominantly on the underlying equation of state and only marginally on the model used for the dynamical description. To do so, a procedure to incorporate any equation of state in the UrQMD transport model is introduced. In particular we compare the baryon density, temperature and pressure evolution as well as produced entropy in a relativistic ideal hydrodynamics approach and the UrQMD transport model, where the same equation of state is used in both approaches. Not only is the compression similar if the same equation of state is used in either dynamical model, but it also strongly depends on the actual equation of state. These results indicate that the equation of state can be studied with observables which are sensitive to the initial compression phase and maximum compression achieved in heavy ion collisions at these beam energies. △ Less

Submitted 12 May, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

Comments: 11 pages, 7 figures, version accepted by: The European Physical Journal C

Perturbative structure of two- and four-point functions of color charge in a non-Gaussian small-$x$ action

Authors: Andre V. Giannini, Yasushi Nara

Abstract: We compute the perturbative expansion of the two- and four-point functions of color charges in the Color Glass Condensate framework considering the quartic correction to the McLerran-Venugopalan (MV) model of Gaussian color charge fluctuations. Expressions for these correlators in the perturbative expansion for small and large non-Gaussian color charge fluctuations are derived for arbitrary orders… ▽ More We compute the perturbative expansion of the two- and four-point functions of color charges in the Color Glass Condensate framework considering the quartic correction to the McLerran-Venugopalan (MV) model of Gaussian color charge fluctuations. Expressions for these correlators in the perturbative expansion for small and large non-Gaussian color charge fluctuations are derived for arbitrary orders in perturbation theory. We explicitly show that the perturbative series does not converge at higher orders as expected. We apply the Borel-Padé resummation method to our problem to construct a convergent series. It is shown that the fully non-perturbative solution can be described by the Borel-Padé approximants constructed from the first few terms of the perturbative series for small non-Gaussian fluctuations. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 17 pages, 6 figures

Mean-field update in the JAM microscopic model: Mean-field effects on collective flow in high-energy heavy-ion collisions at $\sqrt{s_{NN}}=2-20$ GeV energies

Authors: Yasushi Nara, Akira Ohnishi

Abstract: The beam energy dependence of the directed flow is a sensitive probe for the properties of strongly interacting matter. We consider different implementations of momentum-dependent hadronic mean fields in the relativistic quantum molecular dynamics (RQMD) framework. First, Lorentz scalar implementation of a Skyrme type potential is examined. Then, full implementation of the Skyrme type potential as… ▽ More The beam energy dependence of the directed flow is a sensitive probe for the properties of strongly interacting matter. We consider different implementations of momentum-dependent hadronic mean fields in the relativistic quantum molecular dynamics (RQMD) framework. First, Lorentz scalar implementation of a Skyrme type potential is examined. Then, full implementation of the Skyrme type potential as a Lorentz vector in the RQMD approach is proposed. We find that scalar implementation of the Skyrme force is too weak to generate repulsion explaining observed data of sideward flows at $\sqrt{s_{NN}}<10$ GeV, while vector implementation gives collective flows compatible with the data for a wide range of beam energies $2.7 <\sqrt{s_{NN}}<20$ GeV. We show that our approach reproduces the negative proton directed flow at $\sqrt{s_{NN}}>10$ GeV discovered by experiments. We discuss the dynamical generation mechanisms of the directed flow within a conventional hadronic mean field. A positive slope of proton directed flow is generated predominantly during compression stages of heavy-ion collisions by the strong repulsive interaction due to high baryon densities. In contrast, in the expansion stages of the collision, the negative directed flow is generated more strongly than the positive one by the tilted expansion and shadowing by the spectator matter. At lower collision energies $\sqrt{s_{NN}}<10$ GeV, the positive flow wins against the negative flow because of a long compression time. On the other hand, at higher energies $\sqrt{s_{NN}}>10$ GeV, negative flow wins because of shorter compression time and longer expansion time. A transition beam energy from positive to negative flow is highly sensitive to the strength of the interaction. △ Less

Submitted 9 March, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 20 pages, 12 figures, 1 table, typos fixed

Report number: YITP-21-97

Journal ref: Phys. Rev. C105, 014911 (2022)

Momentum-dependent potential and collective flows within the relativistic quantum molecular dynamics approach based on relativistic mean-field theory

Authors: Yasushi Nara, Tomoyuki Maruyama, Horst Stoecker

Abstract: Relativistic quantum molecular dynamics based on the relativistic mean field theory (RQMD.RMF) is extended by including momentum-dependent potential. The equation of state (EoS) dependence of the directed and the elliptic flow of protons in the beam energy range of $2.3 < \sqrt{s_{NN}}< 20$ GeV is examined. It is found that the directed flow depends strongly on the optical potential at high energi… ▽ More Relativistic quantum molecular dynamics based on the relativistic mean field theory (RQMD.RMF) is extended by including momentum-dependent potential. The equation of state (EoS) dependence of the directed and the elliptic flow of protons in the beam energy range of $2.3 < \sqrt{s_{NN}}< 20$ GeV is examined. It is found that the directed flow depends strongly on the optical potential at high energies,$\sqrt{s_{NN}} > 3 $ GeV, where no information is available experimentally. The correlation between effective mass at saturation density and the optical potential is found: smaller values of effective mass require smaller strengths of the optical potential to describe the directed flow data.This correlation can also be seen in the beam energy dependence of the elliptic flow at $\sqrt{s_{NN}}>3$ GeV, although its effect is rather weak. On the other hand, stiff EoS is required to describe the elliptic flow at lower energies.Experimental constraints on the optical potential from $pA$ collisions will provide important information on the EoS at high energies.The proton directed and the elliptic flow are well described in the RQMD.RMF model from $\sqrt{s_{NN}}=2.3$ to 8.8 GeV. In contrast,to reproduce the collapse of the directed flow above 10 GeV, pressure has to be reduced, which indicates a softening of the EoS around $\sqrt{s_{NN}} =10 $ GeV. △ Less

Submitted 30 August, 2020; v1 submitted 12 April, 2020; originally announced April 2020.

Comments: 11 pages, 7 figures, published in PRC

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

Non-perturbative renormalization of the average color charge and multi-point correlators of color charge from a non-Gaussian small-$x$ action

Authors: Andre V. Giannini, Yasushi Nara

Abstract: The McLerran-Venugopalan (MV) model is a Gaussian effective theory of color charge fluctuations at small-$x$ in the limit of large valence charge density, {\it i}.{\it e}., a large nucleus made of uncorrelated color charges. In this work, we explore the effects of the first non-trivial (even C-parity) non-Gaussian correction on the color charge density to the MV model ("quartic" term) in SU(2) and… ▽ More The McLerran-Venugopalan (MV) model is a Gaussian effective theory of color charge fluctuations at small-$x$ in the limit of large valence charge density, {\it i}.{\it e}., a large nucleus made of uncorrelated color charges. In this work, we explore the effects of the first non-trivial (even C-parity) non-Gaussian correction on the color charge density to the MV model ("quartic" term) in SU(2) and SU(3) color group in the non-perturbative regime. We compare our (numerical) non-perturbative results to (analytical) perturbative ones in the limit of small or large non-Gaussian fluctuations. The couplings in the non-Gaussian action, $\barμ$ for the quadratic and $κ_4$ for the quartic term, need to be renormalized in order to match the two-point function in the Gaussian theory. We investigate three different choices for the renormalization of these couplings: i) $κ_{4}$ is proportional to a power of $\barμ$; ii) $κ_4$ is kept constant and iii) $\barμ$ is kept constant. We find that the first two choices lead to a scenario where the small-$x$ action evolves towards a theory dominated by large non-Gaussian fluctuations, regardless of the system size, while the last one allows for controlling the deviations from the MV model. △ Less

Submitted 10 March, 2021; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: v4: The right panel of figure 5 now shows the leading order result in the Z->0 region, in agreement with Eqs. (52) - (54) (now written for an arbitrary number of colors, Nc). Matches the published version

Effects of Pauli blocking on pion production in central collisions of neutron-rich nuclei

Authors: Natsumi Ikeno, Akira Ono, Yasushi Nara, Akira Ohnishi

Abstract: Pauli blocking is carefully investigated for the processes of $NN \rightarrow N Δ$ and $Δ\rightarrow N π$ in heavy-ion collisions, aiming at a more precise prediction of the $π^-/ π^+$ ratio which is an important observable to constrain the high-density symmetry energy. We use the AMD+JAM approach, which combines the antisymmetrized molecular dynamics for the time evolution of nucleons and the JAM… ▽ More Pauli blocking is carefully investigated for the processes of $NN \rightarrow N Δ$ and $Δ\rightarrow N π$ in heavy-ion collisions, aiming at a more precise prediction of the $π^-/ π^+$ ratio which is an important observable to constrain the high-density symmetry energy. We use the AMD+JAM approach, which combines the antisymmetrized molecular dynamics for the time evolution of nucleons and the JAM model to treat processes for $Δ$ resonances and pions. As is known in general transport-code simulations, it is difficult to treat Pauli blocking very precisely due to unphysical fluctuations and additional smearing of the phase-space distribution function, when Pauli blocking is treated in the standard method of JAM. We propose an improved method in AMD+JAM to use the Wigner function precisely calculated in AMD as the blocking probability. Different Pauli blocking methods are compared in heavy-ion collisions of neutron-rich nuclei, ${}^{132}\mathrm{Sn}+{}^{124}\mathrm{Sn}$, at 270 MeV/nucleon. With the more accurate method, we find that Pauli blocking is stronger, in particular for the neutron in the final state in $NN \rightarrow N Δ$ and $ Δ\to Nπ$, compared to the case with a proton in the final state. Consequently, the $π^-/π^+$ ratio becomes higher when the Pauli blocking is improved, the effect of which is found to be comparable to the sensitivity to the high-density symmetry energy. △ Less

Submitted 13 March, 2020; v1 submitted 21 October, 2019; originally announced October 2019.

Comments: 10 pages, 5 figures, 1 table: The title is changed. All figures are updated. Table 1 is added

Journal ref: Phys. Rev. C 101, 034607 (2020)

Sensitivity of the excitation functions of collective flow to relativistic scalar and vector meson interactions in the relativistic quantum molecular dynamics model RQMD.RMF

Authors: Yasushi Nara, Horst Stoecker

Abstract: Relativistic quantum molecular dynamics with scalar and vector interactions based on the relativistic mean meson field theory, RQMD.RMF, is developed.It is implemented into the microscopic transport code JAM.The sensitivity of the directed and of the elliptic proton flow in high energy heavy-ion collisions on the stiffness of the RMF equation of state (EoS) is examined. These new calculations are… ▽ More Relativistic quantum molecular dynamics with scalar and vector interactions based on the relativistic mean meson field theory, RQMD.RMF, is developed.It is implemented into the microscopic transport code JAM.The sensitivity of the directed and of the elliptic proton flow in high energy heavy-ion collisions on the stiffness of the RMF equation of state (EoS) is examined. These new calculations are compared to experimental data at mid-central Au + Au collisions in the beam energy range $2.5 < \sqrt{s_{NN}} < 20$ GeV. This new RQMD model with the relativistic mean field scalar and vector meson interactions does describe consistently, with one RMF parameter set,the beam energy dependence of both the directed flow and the elliptic flow,from SIS18 to AGS and RHIC BES-II energies, $\sqrt{s_{NN}}< 10$ GeV.There are different sensitivities of these different kinds of flow to the EoS: elliptic flow is most sensitive to the nuclear incompressibility constant,at the moderate beam energies $\sqrt{s_{NN}}<3$ GeV,whereas the directed flow is most sensitive to the effective baryon mass at saturation density at $3< \sqrt{s_{NN}}<5 $ GeV. Matters abruptly change in the next higher energy range,$\sqrt{s_{NN}}\gtrsim 10-20$ GeV:the directed flow data show a double change of sign of the slope of $v_1$, inverting twice in this energy range,in sudden contradiction to the RQMD.RMF calculation for a monotonous, stiff EoS. This surprising oscillating behavior,a double change of sign of the $v_1$ slope, points to the appearance of a hitherto unknown first-order phase transition in excited QCD matter at high baryon densities in mid-central Au + Au collisions. △ Less

Submitted 4 November, 2019; v1 submitted 8 June, 2019; originally announced June 2019.

Comments: 9 pages, 5 figures,accepted for publication in PRC

Journal ref: Phys. Rev. C 100, 054902 (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)

Dynamically Integrated Transport Approach for High-Energy Nuclear Collisions at High Baryon Density

Authors: Koichi Murase, Yukinao Akamatsu, Masayuki Asakawa, Tetsufumi Hirano, Masakiyo Kitazawa, Kenji Morita, Yasushi Nara, Chiho Nonaka, Akira Ohnishi

Abstract: To explore the structure of the QCD phase diagram in high baryon density domain, several high-energy nuclear collision experiments in a wide range of beam energies are currently performed or planned using many accelerator facilities. In these experiments search for a first-order phase transition and the QCD critical point is one of the most important topics. To find the signature of the phase tran… ▽ More To explore the structure of the QCD phase diagram in high baryon density domain, several high-energy nuclear collision experiments in a wide range of beam energies are currently performed or planned using many accelerator facilities. In these experiments search for a first-order phase transition and the QCD critical point is one of the most important topics. To find the signature of the phase transition, experimental data should be compared to appropriate dynamical models which quantitatively describe the process of the collisions. In this study we develop a new dynamical model on the basis of the non-equilibrium hadronic transport model JAM and 3+1D hydrodynamics. We show that the new model reproduce well the experimental beam-energy dependence of hadron yields and particle ratio by the partial thermalization of the system in our core-corona approach. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 4 pages, 3 figures; contribution to the proceedings of the 8th International Conference on Quarks and Nuclear Physics (QNP2018), Tsukuba, November 13-17, 2018

Centrality and energy dependence of charged particles in p+A and A+A collisions from running coupling $k_{T}$-factorization

Authors: Adrian Dumitru, Andre V. Giannini, Matthew Luzum, Yasushi Nara

Abstract: We extend the numerical analysis of the energy and centrality dependence of particle multiplicities at midrapidity in high-energy p+A and A+A collisions from a running coupling $k_T$-factorization formula made in~\cite{Dumitru:2018gjm} by considering two unintegrated gluon distributions that were left out. While a good agreement with the experimental data in A+A collisions is achieved, improving t… ▽ More We extend the numerical analysis of the energy and centrality dependence of particle multiplicities at midrapidity in high-energy p+A and A+A collisions from a running coupling $k_T$-factorization formula made in~\cite{Dumitru:2018gjm} by considering two unintegrated gluon distributions that were left out. While a good agreement with the experimental data in A+A collisions is achieved, improving the description of those observables in p+A collisions calls for a better understanding of the proton unintegrated gluon distribution at larger values of $x$ and also the use of a realistic impact parameter dependence. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 6 pages, 4 figures. Proceedings of the Diffraction and Low-x 2018 conference. To be published in Acta Physica Polonica B Proceeding Supplement

The enhancement of $v_4$ in nuclear collisions at the highest densities signals a first order phase transition

Authors: Yasushi Nara, Jan Steinheimer, Horst Stoecker

Abstract: The beam energy dependence of $v_4$ (the quadrupole moment of the transverse radial flow) is sensitive to the nuclear equation of state (EoS) in mid-central Au + Au collisions at the energy range of $3 < \sqrt{s_{NN}} < 30$ GeV, which is investigated within the hadronic transport model JAM. Different equations of state, namely, a free hadron gas, a first-order phase transition and a crossover are… ▽ More The beam energy dependence of $v_4$ (the quadrupole moment of the transverse radial flow) is sensitive to the nuclear equation of state (EoS) in mid-central Au + Au collisions at the energy range of $3 < \sqrt{s_{NN}} < 30$ GeV, which is investigated within the hadronic transport model JAM. Different equations of state, namely, a free hadron gas, a first-order phase transition and a crossover are compared. An enhancement of $v_4$ at $\sqrt{s_{NN}}\approx 6$ GeV is predicted for an EoS with a first-order phase transition. This enhanced $v_4$ flow is driven by both the enhancement of $v_2$ as well as the positive contribution to $v_4$ from the squeeze-out of spectator particles which turn into participants due to the admixture of the strong collective flow in the shocked, compressed nuclear matter. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 10pages, 3 figures

The cumulants of the baryon number from central Au+Au collision at $E_{lab}$= 1.23 GeV$/$nucleon reveal the nuclear mean-field potentials

Authors: Yunxiao Ye, Yongjia Wang, Jan Steinheimer, Yasushi Nara, Hao-jie Xu, Pengcheng Li, Dinghui Lu, Qingfeng Li, Horst Stoecker

Abstract: Fluctuations of the baryon number in relativistic heavy-ion collisions are a promising observable to explore the structure of the QCD phase diagram. The cumulant ratios in heavy ion collisions at intermediate energies ($\sqrt{s_{\textrm{NN}}} < 7$ GeV) have not been studied to date. We investigate the effects of mean field potential and clustering on the cumulant ratios of baryon and proton number… ▽ More Fluctuations of the baryon number in relativistic heavy-ion collisions are a promising observable to explore the structure of the QCD phase diagram. The cumulant ratios in heavy ion collisions at intermediate energies ($\sqrt{s_{\textrm{NN}}} < 7$ GeV) have not been studied to date. We investigate the effects of mean field potential and clustering on the cumulant ratios of baryon and proton number distributions in Au+Au collisions at beam energy of 1.23 GeV$/$nucleon as measured by the HADES Collaboration at GSI. Ultrarelativistic Quantum Molecular Dynamics (UrQMD) and the JAM model are used to calculate the cumulants with different mean field potentials. It is found that the cumulant ratios are strongly time dependent. At the early stage, the effects of the potentials on the fluctuations of the particle multiplicity in momentum space are relatively weak. The mean fields enhance the fluctuations during the expansion stage, especially for small rapidity acceptance windows. The enhancement of cumulant ratios for free protons is strongly suppressed as compared to that for all baryons. The mean field potentials and the clustering play an important role for the measured cumulant ratios at intermediate energy. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 10 pages, 11 figures

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

Dynamically integrated transport approach for heavy-ion collisions at high baryon density

Authors: Yukinao Akamatsu, Masayuki Asakawa, Tetsufumi Hirano, Masakiyo Kitazawa, Kenji Morita, Koichi Murase, Yasushi Nara, Chiho Nonaka, Akira Ohnishi

Abstract: We develop a new dynamical model for high energy heavy-ion collisions in the beam energy region of the highest net-baryon densities on the basis of non-equilibrium microscopic transport model JAM and macroscopic 3+1D hydrodynamics by utilizing a dynamical initialization method. In this model,dynamical fluidization of a system is controlled by the source terms of the hydrodynamic fields. In additio… ▽ More We develop a new dynamical model for high energy heavy-ion collisions in the beam energy region of the highest net-baryon densities on the basis of non-equilibrium microscopic transport model JAM and macroscopic 3+1D hydrodynamics by utilizing a dynamical initialization method. In this model,dynamical fluidization of a system is controlled by the source terms of the hydrodynamic fields. In addition, time dependent core-corona separation of hot regions is implemented. We show that our new model describes multiplicities and mean transverse mass in heavy-ion collisions within a beam energy region of $3<\sqrt{s_{NN}}<30$ GeV. Good agreement of the beam energy dependence of the $K^+/π^+$ ratio is obtained, which is explained by the fact that a part of the system is not thermalized in our core-corona approach. △ Less

Submitted 15 August, 2018; v1 submitted 23 May, 2018; originally announced May 2018.

Comments: 8 pages, 8 figures, new figures added; published in Physical Review C

Report number: J-PARC-TH-0131

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

Particle multiplicities in the central region of high-energy collisions from $k_T$-factorization with running coupling corrections

Authors: Adrian Dumitru, Andre V. Giannini, Matthew Luzum, Yasushi Nara

Abstract: Horowitz and Kovchegov have derived a $k_T$-factorization formula for particle production at small $x$ which includes running coupling corrections. We perform a first numerical analysis to confront the theory with data on the energy and centrality dependence of particle multiplicities at midrapidity in high-energy p+A (and A+A) collisions. Moreover, we point out a strikingly different dependence o… ▽ More Horowitz and Kovchegov have derived a $k_T$-factorization formula for particle production at small $x$ which includes running coupling corrections. We perform a first numerical analysis to confront the theory with data on the energy and centrality dependence of particle multiplicities at midrapidity in high-energy p+A (and A+A) collisions. Moreover, we point out a strikingly different dependence of the multiplicity per participant on $N_\text{part}$ in p+Pb vs.\ Pb+Pb collisions at LHC energies, and argue that the observed behavior follows rather naturally from the convolution of the gluon distributions of an asymmetric vs. symmetric projectile and target. △ Less

Submitted 27 August, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 8 pages, 8 figures; v3: 2 new figures, matches the version accepted by Physics Letters B

Beam energy dependence of squeeze-out effect on the directed and elliptic flow in Au+Au collisions at high baryon density region

Authors: Chao Zhang, Jiamin Chen, Xiaofeng Luo, Feng Liu, Yasushi Nara

Abstract: We present a detailed analysis of the beam energy dependence of the mechanisms for the generation of directed and elliptic flows in Au+Au collisions focusing on the role of hadronic rescattering and spectator shadowing within a microscopic transport model JAM with different equation of state. A systematic study of the beam energy dependence is performed for Au+Au collisions at… ▽ More We present a detailed analysis of the beam energy dependence of the mechanisms for the generation of directed and elliptic flows in Au+Au collisions focusing on the role of hadronic rescattering and spectator shadowing within a microscopic transport model JAM with different equation of state. A systematic study of the beam energy dependence is performed for Au+Au collisions at $\sqrt{s_{NN}} =2.3 - 62.4$ GeV. The transition of the dynamical origin of the directed flow is observed. We find that the initial Glauber type nucleon-nucleon collisions generate negative $v_1$ for nucleons at midrapidity due to the presence of spectator matter, and this negative nucleon $v_1$ is turned to be positive by the meson-baryon interactions at the beam energy region of $\sqrt{s_{NN}} < 30$ GeV. In contrast, above 30 GeV there is no spectator shadowing at midrapidity, and initial nucleon-nucleon collisions do not generate directed flow, but subsequent rescatterings among produced particles generate negative $v_1$ for nucleons. It is demonstrated that negative pion-directed flows are mostly generated by the interaction with the spectator matter. It is also shown that the squeeze-out effect is largely suppressed in the case of softening, which leads to the enhancement of elliptic flow around $\sqrt{s_{NN}}=5-7$ GeV. The elliptic flow at midrapidity above 10 GeV is not influenced by the squeeze-out due to spectator matter, while its effect is seen at the forward rapidity range of $y/y_\mathrm{c.m.}>0.5$, which decreases as beam energy increases. △ Less

Submitted 25 June, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: 14 pages, 17 figures

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

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)

Enhancement of elliptic flow can signal a first order phase transition in high energy heavy ion collisions

Authors: Yasushi Nara, Harri Niemi, Akira Ohnishi, Jan Steinheimer, Xiaofeng Luo, Horst Stoecker

Abstract: The beam energy dependence of the elliptic flow,$v_2$, is studied in mid-central Au+Au collisions in the energy range of $3\leq \sqrt{s_{NN}} \leq 30$ GeV within the microscopic transport model JAM. The results of three different modes of JAM are compared; cascade-,hadronic mean field-, and a new mode with modified equations of state, with a first order phase transition (1.O.P.T.) and with a cross… ▽ More The beam energy dependence of the elliptic flow,$v_2$, is studied in mid-central Au+Au collisions in the energy range of $3\leq \sqrt{s_{NN}} \leq 30$ GeV within the microscopic transport model JAM. The results of three different modes of JAM are compared; cascade-,hadronic mean field-, and a new mode with modified equations of state, with a first order phase transition (1.O.P.T.) and with a crossover transition. The standard hadronic mean field suppresses $v_2$, while the inclusion of the effects of a 1.O.P.T. (and also of a crossover transition) does enhance $v_2$ at $\sqrt{s_{NN}}<30$ GeV. The enhancement or suppression of the scaled energy flow, dubbed "elliptic flow"is understood as being due to out of plane- flow, i.e. $v_2<0$, dubbed out of plane - "squeeze-out", which occurs predominantly in the early, compression stage. Subsequently, the in-plane flow dominates, in the expansion stage, $v_2 > 0$. The directed flow, dubbed "bounce- off", is an independent measure of the pressure, which quickly builds up the transverse momentum transfer in the reaction plane. When the spectator matter leaves the participant fireball region, where the highest compression occurs, a hard expansion leads to larger $v_2$. A combined analysis of the three transverse flow coefficients, radial $v_0$-, directed $v_1$- and elliptic $v_2$- flow, in the beam energy range of $3\leq\sqrt{s_{NN}}\leq10$ GeV, distinguishes the different compression and expansion scenarios: a characteristic dependence on the early stage equation of state is observed. The enhancement of both the elliptic and the transverse radial flow and the simultaneous collapse of the directed flow of nucleons offers a clear signature if 1.O.P.T. is realized at the highest baryon densities created in high energy heavy-ion collisions. △ Less

Submitted 18 August, 2017; originally announced August 2017.

Comments: 12 pages,10 figures

Report number: YITP-17-80

Effects of mean-field and the softening of equation of state on elliptic flow in Au+Au collision at $\sqrt{s_{NN}}$ = 5 GeV from JAM model

Authors: Jiamin Chen, Xiaofeng Luo, Feng Liu, Yasushi Nara

Abstract: We perform a systematic study of elliptic flow ($v_2$) in Au+Au collision at $\sqrt{s_{NN}}=5$ GeV by using a microscopic transport model JAM. The centrality, pseudorapidity, transverse momentum and beam energy dependence of $v_2$ for charged as well as identified hadrons are studied. We investigate the effects of both the hadronic mean-field and the softening of equation of state (EoS) on ellipti… ▽ More We perform a systematic study of elliptic flow ($v_2$) in Au+Au collision at $\sqrt{s_{NN}}=5$ GeV by using a microscopic transport model JAM. The centrality, pseudorapidity, transverse momentum and beam energy dependence of $v_2$ for charged as well as identified hadrons are studied. We investigate the effects of both the hadronic mean-field and the softening of equation of state (EoS) on elliptic flow. The softening of EoS is realized by imposing attractive orbits in two body scattering, which can reduce the pressure of the system. We found that the softening of EoS leads to the enhancement of $v_2$, while the hadronic mean-field suppresses $v_2$ relative to the cascade mode. It indicates that elliptic flow at high baryon density regions is highly sensitive to the EoS and the enhancement of $v_2$ may probe the signature of a first-order phase transition in heavy-ion collisions at beam energies of a strong baryon stopping region. △ Less

Submitted 19 December, 2017; v1 submitted 2 January, 2017; originally announced January 2017.

Comments: 7 pages, 9 figures, to appear in Chinese Physics C

Journal ref: Chinese Physics C 42, 024001 (2018)

Equation of state dependence of directed flow in a microscopic transport model

Authors: Yasushi Nara, Harri Niemi, Jan Steinheimer, Horst Stoecker

Abstract: We study the sensitivities of the directed flow in Au+Au collisions on the equation of state (EoS), employing the transport theoretical model JAM. The EoS is modified by introducing a new collision term in order to control the pressure of a system by appropriately selecting an azimuthal angle in two-body collisions according to a given EoS. It is shown that this approach is an efficient method to… ▽ More We study the sensitivities of the directed flow in Au+Au collisions on the equation of state (EoS), employing the transport theoretical model JAM. The EoS is modified by introducing a new collision term in order to control the pressure of a system by appropriately selecting an azimuthal angle in two-body collisions according to a given EoS. It is shown that this approach is an efficient method to modify the EoS in a transport model. The beam energy dependence of the directed flow of protons is examined with two different EoS, a first-order phase transition and crossover. It is found that our approach yields quite similar results as hydrodynamical predictions on the beam energy dependence of the directed flow; Transport theory predicts a minimum in the excitation function of the slope of proton directed flow and does indeed yield negative directed flow, if the EoS with a first-order phase transition is employed. Our result strongly suggests that the highest sensitivity for the critical point can be seen in the beam energy range of $4.7\leq\srtNN\leq11.5$ GeV. △ Less

Submitted 23 November, 2016; originally announced November 2016.

Comments: 7 pages, 5 figures

Effects of Nuclear Potential on the Cumulants of Net-Proton and Net-Baryon Multiplicity Distributions in Au+Au Collisions at $\sqrt{s_{\text{NN}}} = 5\,\text{GeV}$

Authors: Shu He, Xiaofeng Luo, Yasushi Nara, ShinIchi Esumi, Nu Xu

Abstract: We analyze the rapidity and transverse momentum dependence for the cumulants of the net-proton and net-baryon distributions in Au+Au collisions at $\sqrt{s_{\text{NN}}} = 5\,\text{GeV}$ with a microscopic hadronic transport (JAM) model. To study the effects of mean field potential and softening of equation of state (EoS) on the fluctuations of net-proton (baryon) in heavy-ion collisions, the calcu… ▽ More We analyze the rapidity and transverse momentum dependence for the cumulants of the net-proton and net-baryon distributions in Au+Au collisions at $\sqrt{s_{\text{NN}}} = 5\,\text{GeV}$ with a microscopic hadronic transport (JAM) model. To study the effects of mean field potential and softening of equation of state (EoS) on the fluctuations of net-proton (baryon) in heavy-ion collisions, the calculations are performed with two different modes. The softening of EoS is realized in the model by implementing the attractive orbit in the two-body scattering to introduce a reduction pressure of the system. By comparing the results from the two modes with the results from default cascade, we find the mean field potential and softening of EoS have strong impacts on the rapidity distributions ($\text{d}N/\text{d}y$) and the shape of the net-proton (baryon) multiplicity distributions. The net-proton (baryon) cumulants and their ratios calculated from all of the three modes are with similar trends and show significant suppression with respect to unity, which can be explained by the presence of baryon number conservations. It indicates that the effects of mean field potential and softening of EoS might be not the ingredients that are responsible to the observed strong enhancement in the most central Au+Au collisions at 7.7 GeV measured by the STAR experiment at RHIC. △ Less

Submitted 15 November, 2016; v1 submitted 21 July, 2016; originally announced July 2016.

Comments: 6 pages, 5 figures

Journal ref: Physics Letter B 762, 296 (2016)

Predictions for $p+$Pb Collisions at $\sqrt{s_{NN}} = 5$ TeV: Comparison with Data

Authors: J. L. Albacete, F. Arleo, G. G. Barnafoldi, J. Barrette, W. -T. Deng, A. Dumitru, K. J. Eskola, E. G. Ferreiro, F. Fleuret, H. Fujii, M. Gyulassy, S. M. Harangozi, I. Helenius, Z. Kang, P. Kotko, K. Kutak, J. -P. Lansberg, P. Levai, Z. -W. Lin, Y. Nara, A. Rakotozafindrabe, G. Papp, H. Paukkunen, S. Peigne, M. Petrovici , et al. (13 additional authors not shown)

Abstract: Predictions made in Albacete {\it et al} prior to the LHC $p+$Pb run at $\sqrt{s_{NN}} = 5$ TeV are compared to currently available data. Some predictions shown here have been updated by including the same experimental cuts as the data. Some additional predictions are also presented, especially for quarkonia, that were provided to the experiments before the data were made public but were too late… ▽ More Predictions made in Albacete {\it et al} prior to the LHC $p+$Pb run at $\sqrt{s_{NN}} = 5$ TeV are compared to currently available data. Some predictions shown here have been updated by including the same experimental cuts as the data. Some additional predictions are also presented, especially for quarkonia, that were provided to the experiments before the data were made public but were too late for the original publication are also shown here. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: 55 pages 35 figures

Ultra-forward particle production from CGC+Lund fragmentation

Authors: Javier L. Albacete, Pablo Guerrero Rodríguez, Yasushi Nara

Abstract: We present an analysis of data on single inclusive pion production measured by the LHCf collaboration in high-energy proton-proton and proton-nucleus at ultra-forward rapidities, $8.8 \leq\! y\leq \!10.8$. We also analyse forward RHIC data for calibration purposes. Our analysis relies on the use of a Monte Carlo event generator that combines a perturbative description of the elementary scattering… ▽ More We present an analysis of data on single inclusive pion production measured by the LHCf collaboration in high-energy proton-proton and proton-nucleus at ultra-forward rapidities, $8.8 \leq\! y\leq \!10.8$. We also analyse forward RHIC data for calibration purposes. Our analysis relies on the use of a Monte Carlo event generator that combines a perturbative description of the elementary scattering process at partonic level based on the hybrid formalism of the Color Glass Condensate with an implementation of hadronization in the framework of the Lund string fragmentation model. This procedure allows us to reach values of the momenta of the produced particles as low as detected experimentally $p_t\sim0.1$ GeV. We achieve a good description of single inclusive spectra of charged particles and neutral pions at RHIC and the LHC respectively, and nuclear modification factors for proton-lead collisions at the LHC. Our results add evidence to the idea that particle production in the domain of very small Bjorken-$x$ is dominated by the saturation effects encoded in the unintegrated gluon distribution of the target. Being forward particle production of key importance in the development of air showers, we stress that this approach allows for a theoretically controlled extrapolation of our results to the scale of ultra-high energy cosmic rays, thus serving as starting point for future works on this topic. △ Less

Submitted 26 May, 2016; originally announced May 2016.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. D 94, 054004 (2016)

Examination of directed flow as a signature of the softest point of the equation of state in QCD matter

Authors: Yasushi Nara, Harri Niemi, Akira Ohnishi, Horst Stoecker

Abstract: We analyze the directed flow of protons and pions in high-energy heavy-ion collisions in the incident energy range from $\sqrt{s_{\scriptscriptstyle NN}}=7.7$ to 27 GeV within a microscopic transport model. Standard hadronic transport approaches do not describe the collapse of directed flow below $\sqrt{s_{\scriptscriptstyle NN}}\simeq 20$ GeV. By contrast, a model which simulates effects of a sof… ▽ More We analyze the directed flow of protons and pions in high-energy heavy-ion collisions in the incident energy range from $\sqrt{s_{\scriptscriptstyle NN}}=7.7$ to 27 GeV within a microscopic transport model. Standard hadronic transport approaches do not describe the collapse of directed flow below $\sqrt{s_{\scriptscriptstyle NN}}\simeq 20$ GeV. By contrast, a model which simulates effects of a softening of the equation of state, well describes the behavior of directed flow data recently obtained by the STAR Collaboration~\cite{STARv1}. We give a detailed analysis of how directed flow is generated. Particularly, we found that softening of effective equation of state at the overlapping region of two nuclei, i.e. the reaction stages where the system reaches high baryon density state, is needed to explain the observed collapse of proton directed flow within a hadronic transport approach. △ Less

Submitted 26 September, 2016; v1 submitted 28 January, 2016; originally announced January 2016.

Comments: 11 pages, 12 figures. Typos fixed, version accepted in PRC

Report number: YITP-15-120

Journal ref: Phys. Rev. C 94, 034906 (2016)

Probing neutron-proton dynamics by pions

Authors: Natsumi Ikeno, Akira Ono, Yasushi Nara, Akira Ohnishi

Abstract: In order to investigate the nuclear symmetry energy at high density, we study the pion production in central collisions of neutron-rich nuclei ${}^{132}\mathrm{Sn}+{}^{124}\mathrm{Sn}$ at 300 MeV/nucleon using a new approach by combining the antisymmetrized molecular dynamics (AMD) and a hadronic cascade model (JAM). The dynamics of neutrons and protons is solved by AMD, and then pions and $Δ$ res… ▽ More In order to investigate the nuclear symmetry energy at high density, we study the pion production in central collisions of neutron-rich nuclei ${}^{132}\mathrm{Sn}+{}^{124}\mathrm{Sn}$ at 300 MeV/nucleon using a new approach by combining the antisymmetrized molecular dynamics (AMD) and a hadronic cascade model (JAM). The dynamics of neutrons and protons is solved by AMD, and then pions and $Δ$ resonances in the reaction process are handled by JAM. We see the mechanism how the $Δ$ resonance and pions are produced reflecting the dynamics of neutrons and protons. We also investigate the impacts of cluster correlations as well as of the high-density symmetry energy on the nucleon dynamics and consequently on the pion ratio. We find that the $Δ^-/Δ^{++}$ production ratio agrees very well with the neutron-proton squared ratio $(N/Z)^2$ in the high-density and high-momentum region. We show quantitatively that $Δ$ production ratio, and therefore $(N/Z)^2$, are directly reflected in the $π^-/π^+$ ratio, with modification in the final stage of the reaction. △ Less

Submitted 5 April, 2018; v1 submitted 27 January, 2016; originally announced January 2016.

Comments: 14 pages, 10 figures; Figures 3-8 are updated with corrected numerical results. No change in the main conclusions

Report number: YITP-16-6

Violation of mass ordering for multi-strange hadrons at RHIC and LHC

Authors: Shiori Takeuchi, Koichi Murase, Tetsufumi Hirano, Pasi Huovinen, Yasushi Nara

Abstract: We study effects of the hadronic rescattering on final observables especially for multi-strange hadrons such as $φ$, $Ξ$ and $Ω$ in high-energy heavy-ion collisions within an integrated dynamical approach. In this approach, (3+1)-dimensional ideal hydrodynamics is combined with a microscopic transport model, JAM. We simulate the collisions with or without hadronic rescatterings and compare observa… ▽ More We study effects of the hadronic rescattering on final observables especially for multi-strange hadrons such as $φ$, $Ξ$ and $Ω$ in high-energy heavy-ion collisions within an integrated dynamical approach. In this approach, (3+1)-dimensional ideal hydrodynamics is combined with a microscopic transport model, JAM. We simulate the collisions with or without hadronic rescatterings and compare observables between these two options so that we quantify the effects of the hadronic rescattering. We find that the mean transverse momentum and the elliptic flow parameter of multi-strange hadrons are less affected by hadronic rescattering and, as a result, the mass ordering of the $p_T$-differential elliptic flow parameter $v_2(p_T)$ is violated: At the RHIC and the LHC energies the $v_2(p_T)$ for $φ$-mesons is larger than that for protons in the low-$p_T$ regions. △ Less

Submitted 18 January, 2016; originally announced January 2016.

Comments: 4 pages, 2 figures, proceedings of Quark Matter 2015

Journal ref: Nucl. Phys. A956 (2016) 457-460

Does non-monotonic behavior of directed flow signal the onset of deconfinement?

Authors: Yasushi Nara, Akira Ohnishi

Abstract: We investigate the effects of nuclear mean-field as well as the formation and decay of nuclear clusters on the directed flow $v_1$ in high energy nucleus-nucleus collisions from $\sqrt{s_{NN}}=7.7$ GeV to 27 GeV incident energies within a transport model. Specifically, we use the JAM transport model in which potentials are implemented based on the framework of the relativistic quantum molecular dy… ▽ More We investigate the effects of nuclear mean-field as well as the formation and decay of nuclear clusters on the directed flow $v_1$ in high energy nucleus-nucleus collisions from $\sqrt{s_{NN}}=7.7$ GeV to 27 GeV incident energies within a transport model. Specifically, we use the JAM transport model in which potentials are implemented based on the framework of the relativistic quantum molecular dynamics. Our approach reproduces the rapidity dependence of directed flow data up to $\sqrt{s_{NN}}\approx 8$ GeV showing the significant importance of mean-field. However, the slopes of $dv_1/dy$ at mid-rapidity are calculated to be positive at $\sqrt{s_{NN}}=11.7$ and 19.6 GeV, and becomes negative above 27 GeV. Thus the result from the JAM hadronic transport model with nuclear mean-field approach is incompatible with the data. Therefore within our approach, we conclude that the excitation function of the directed flow cannot be explained by the hadronic degree of freedom alone. △ Less

Submitted 19 December, 2015; originally announced December 2015.

Comments: Proceedings of the 25th International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2015), Kobe, Japan, 4 pages, 3figures

Report number: YITP-15-119

Effects of hadronic rescattering on multistrange hadrons in high-energy nuclear collisions

Authors: Shiori Takeuchi, Koichi Murase, Tetsufumi Hirano, Pasi Huovinen, Yasushi Nara

Abstract: We study the effects of hadronic rescattering on hadron distributions in high-energy nuclear collisions by using an integrated dynamical approach. This approach is based on a hybrid model combining (3+1)-dimensional ideal hydrodynamics for the quark gluon plasma (QGP), and a transport model for the hadron resonance gas. Since the hadron distributions are the result of the entire expansion history… ▽ More We study the effects of hadronic rescattering on hadron distributions in high-energy nuclear collisions by using an integrated dynamical approach. This approach is based on a hybrid model combining (3+1)-dimensional ideal hydrodynamics for the quark gluon plasma (QGP), and a transport model for the hadron resonance gas. Since the hadron distributions are the result of the entire expansion history of the system, understanding the QGP properties requires investigating how rescattering during the hadronic stage affects the final distributions of hadrons. We include multistrange hadrons in our study, and quantify the effects of hadronic rescattering on their mean transverse momenta and elliptic flow. We find that multistrange hadrons scatter less during the hadronic stage than nonstrange particles, and thus their distributions reflect the properties of the system in an earlier stage than the distributions of nonstrange particles. △ Less

Submitted 18 October, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 14 pages, 9 figures, v2: published version, title is changed

Journal ref: Phys. Rev. C 92, 044907 (2015)

Forward Hadron Productions in high energy pp collisions from a Monte-Carlo generator for Color Glass Condensate

Authors: Wei-Tian Deng, Hirotsugu Fujii, Kazunori Itakura, Yasushi Nara

Abstract: We develop a Monte-Carlo event generator based on combination of a parton production formula including the effects of parton saturation (called the DHJ formula) and hadronization process due to the Lund string fragmentation model. This event generator is designed for the description of hadron productions at forward rapidities and in a wide transverse momentum range in high-energy proton-proton col… ▽ More We develop a Monte-Carlo event generator based on combination of a parton production formula including the effects of parton saturation (called the DHJ formula) and hadronization process due to the Lund string fragmentation model. This event generator is designed for the description of hadron productions at forward rapidities and in a wide transverse momentum range in high-energy proton-proton collisions. We analyze transverse momentum spectra of charged hadrons as well as identified particles; pion, kaon, (anti-)proton at RHIC energy, and ultra-forward neutral pion spectra from LHCf experiment. We compare our results to those obtained in other models based on parton-hadron duality and fragmentation functions. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 10 pages, 9 figures

Structure of chromomagnetic fields in the glasma

Authors: T. Lappi, A. Dumitru, Y. Nara

Abstract: The initial stage of a heavy ion collision is dominated by nonperturbatively strong chromoelectric and -magnetic fields. The spatial Wilson loop provides a gauge invariant observable to probe the dynamics of the longitudinal chromomagnetic field. We discuss recent results from a real time lattice calculation of the area-dependence of the expectation value of the spatial Wilson loop. We show that a… ▽ More The initial stage of a heavy ion collision is dominated by nonperturbatively strong chromoelectric and -magnetic fields. The spatial Wilson loop provides a gauge invariant observable to probe the dynamics of the longitudinal chromomagnetic field. We discuss recent results from a real time lattice calculation of the area-dependence of the expectation value of the spatial Wilson loop. We show that at relatively early times after the collision, a universal scaling as a function of the area emerges at large distances for very different initial conditions, with a nontrivial critical exponent. A similar behavior has earlier been seen in calculations of the gluon transverse momentum spectrum, which becomes independent of the initial spectrum of gauge fields. We also show the distribution of eigenvalues of the spatial Wilson loop and the fluctuations of its real and imaginary parts. △ Less

Submitted 31 July, 2014; originally announced July 2014.

Comments: 4 pages, 3 figures, talk by TL at the Quark Matter 2014 conference, Darmstadt, Germany, May 19-24

Structure of longitudinal chromomagnetic fields in high energy collisions

Authors: A. Dumitru, T. Lappi, Y. Nara

Abstract: We compute expectation values of spatial Wilson loops in the forward light cone of high-energy collisions. We consider ensembles of gauge field configurations generated from a classical Gaussian effective action as well as solutions of high-energy renormalization group evolution with fixed and running coupling. The initial fields correspond to a color field condensate exhibiting domain-like struct… ▽ More We compute expectation values of spatial Wilson loops in the forward light cone of high-energy collisions. We consider ensembles of gauge field configurations generated from a classical Gaussian effective action as well as solutions of high-energy renormalization group evolution with fixed and running coupling. The initial fields correspond to a color field condensate exhibiting domain-like structure over distance scales of order the saturation scale. At later times universal scaling emerges at large distances for all ensembles, with a nontrivial critical exponent. Finally, we compare the results for the Wilson loop to the two-point correlator of magnetic fields. △ Less

Submitted 29 April, 2014; v1 submitted 16 January, 2014; originally announced January 2014.

Comments: RevTeX, 7 pages, 8 figures. V2: minor clarifications in text, Fig 6 revised

Journal ref: Phys. Lett. B 734 (2014) 7-12

Magnetic screening in high-energy heavy-ion collisions

Authors: Adrian Dumitru, Hirotsugu Fujii, Yasushi Nara

Abstract: We show that classical chromomagnetic fields produced coherently in the initial stage of a heavy-ion collision exhibit screening. From the two-point field strength correlator we determine the magnetic mass for SU(2) to be m ~ 5 times the saturation scale. Magnetic screening leads to an intuitive understanding of the area law scaling of spatial Wilson loops observed previously. The presence of scre… ▽ More We show that classical chromomagnetic fields produced coherently in the initial stage of a heavy-ion collision exhibit screening. From the two-point field strength correlator we determine the magnetic mass for SU(2) to be m ~ 5 times the saturation scale. Magnetic screening leads to an intuitive understanding of the area law scaling of spatial Wilson loops observed previously. The presence of screening effects in the initial state provides a basis for defining kinetic processes in the early stage of heavy-ion collisions, with electric and magnetic masses of the same order. △ Less

Submitted 13 May, 2013; originally announced May 2013.

Report number: RBRC 1020

Journal ref: Phys. Rev. D 88, 031503(R) (2013)

Magnetic flux loop in high-energy heavy-ion collisions

Authors: Adrian Dumitru, Yasushi Nara, Elena Petreska

Abstract: We consider the expectation value of a magnetic flux loop in the immediate forward light cone of collisions of heavy nuclei at high energies. Such collisions are characterized by a non-linear scale Q_s where color fields become strong. We find that loops of area greater than ~2/Q_s^2 exhibit area law behavior, which determines the scale of elementary flux excitations ("vortices"). We also estimate… ▽ More We consider the expectation value of a magnetic flux loop in the immediate forward light cone of collisions of heavy nuclei at high energies. Such collisions are characterized by a non-linear scale Q_s where color fields become strong. We find that loops of area greater than ~2/Q_s^2 exhibit area law behavior, which determines the scale of elementary flux excitations ("vortices"). We also estimate the magnetic string tension, sigma_M = 0.12 Q_s^2. By the time t ~ 1/Q_s even small loops satisfy an area law. We describe corrections to the propagator of semi-hard gluons at very early times in the background of fluctuating magnetic fields. △ Less

Submitted 16 September, 2013; v1 submitted 8 February, 2013; originally announced February 2013.

Comments: 4 pages, 5 figures; v2: added plot of Z(N) part of loop in fig.1; v3: added magnetic loop for asymmetric projectile/target saturation momenta, estimate of vortex density, and an appendix; v4: corrected outline of perturbative calculation of Wilson loop; to appear in PRD

Report number: RBRC 1008

Predictions for $p+$Pb Collisions at sqrt s_NN = 5 TeV

Authors: J. L. Albacete, N. Armesto, R. Baier, G. G. Barnafoldi, J. Barrette, S. De, W. -T. Deng, A. Dumitru, K. Dusling, K. J. Eskola, R. Fries, H. Fujii, F. Gelis, M. Gyulassy, I. Helenius, Z. -B. Kang, B. Z. Kopeliovich, K. Kutak, P. Levai, Z. Lin, A. H. Mueller, Y. Nara, J. Nemchik, G. Papp, M. Petrovici , et al. (16 additional authors not shown)

Abstract: Predictions for charged hadron, identified light hadron, quarkonium, photon, jet and gauge bosons in p+Pb collisions at sqrt s_NN = 5 TeV are compiled and compared. When test run data are available, they are compared to the model predictions. Predictions for charged hadron, identified light hadron, quarkonium, photon, jet and gauge bosons in p+Pb collisions at sqrt s_NN = 5 TeV are compiled and compared. When test run data are available, they are compared to the model predictions. △ Less

Submitted 22 January, 2013; v1 submitted 15 January, 2013; originally announced January 2013.

Comments: 84 pages

Report number: LLNL-JRNL-609944

Journal ref: Int. J. Mod. Phys. E Vol. 22 (2013) 1330007

CGC predictions for p+Pb collisions at the LHC

Authors: Javier L. Albacete, Adrian Dumitru, Hirotsugu Fujii, Yasushi Nara

Abstract: We present predictions for total multiplicities and single inclusive particle production in proton-lead collisions at the LHC. The main dynamical input in our calculations is the use of solutions of the running coupling BK equation tested in e+p data. We use a Monte-Carlo treatment of the nuclear geometry and either $k_t$-factorization or the hybrid formalisms to describe particle production in th… ▽ More We present predictions for total multiplicities and single inclusive particle production in proton-lead collisions at the LHC. The main dynamical input in our calculations is the use of solutions of the running coupling BK equation tested in e+p data. We use a Monte-Carlo treatment of the nuclear geometry and either $k_t$-factorization or the hybrid formalisms to describe particle production in the central and forward rapidity regions, respectively. △ Less

Submitted 25 September, 2012; v1 submitted 10 September, 2012; originally announced September 2012.

Comments: 17 pages, 13 figures; v1: draft of predictions submitted in time for the upcoming p+Pb test run at the LHC; a more polished version will follow soon. v2: text and discussion cleaned up. v3: added spectra at 5 TeV (pp and p+Pb mb), final version as submitted to Nucl Phys A

Report number: RBRC 969

Journal ref: Nucl. Phys. A 897 (2013) p. 1 - 27

Integrated Dynamical Approach to Relativistic Heavy Ion Collisions

Authors: Tetsufumi Hirano, Pasi Huovinen, Koichi Murase, Yasushi Nara

Abstract: We review integrated dynamical approaches to describe heavy ion reaction as a whole at ultrarelativistic energies. Since final observables result from all the history of the reaction, it is important to describe all the stages of the reaction to obtain the properties of the quark gluon plasma from experimental data. As an example of these approaches, we develop an integrated dynamical model, which… ▽ More We review integrated dynamical approaches to describe heavy ion reaction as a whole at ultrarelativistic energies. Since final observables result from all the history of the reaction, it is important to describe all the stages of the reaction to obtain the properties of the quark gluon plasma from experimental data. As an example of these approaches, we develop an integrated dynamical model, which is composed of a fully (3+1) dimensional ideal hydrodynamic model with the state-of-the-art equation of state based on lattice QCD, and subsequent hadronic cascade in the late stage. Initial conditions are obtained employing Monte Carlo versions of the Kharzeev-Levin-Nardi model (MC-KLN) or the Glauber model (MC-Glauber). Using this integrated model, we first simulate relativistic heavy ion collisions at the RHIC and LHC energies starting from conventional smooth initial conditions. We next utilise each Monte Carlo samples of initial conditions on an event-by-event basis and perform event-by-event dynamical simulations to accumulate a large number of minimum bias events. A special attention is paid to performing the flow analysis as in experiments toward consistent comparison of theoretical results with experimental data. △ Less

Submitted 11 February, 2013; v1 submitted 25 April, 2012; originally announced April 2012.

Comments: v2: 70 pages, 59 figures, two sections on overviews included, conclusion unchanged, version accepted for publication in Progress in Particle and Nuclear Physics

Dynamical modeling of high energy heavy ion collisions

Authors: Tetsufumi Hirano, Yasushi Nara

Abstract: We present theoretical approaches to high energy nuclear collisions in detail putting a special emphasis on technical aspects of numerical simulations. Models include relativistic hydrodynamics, Monte-Carlo implementation of k_T-factorization formula, jet quenching in expanding fluids, a hadronic transport model and the Vlasov equation for colored particles. We present theoretical approaches to high energy nuclear collisions in detail putting a special emphasis on technical aspects of numerical simulations. Models include relativistic hydrodynamics, Monte-Carlo implementation of k_T-factorization formula, jet quenching in expanding fluids, a hadronic transport model and the Vlasov equation for colored particles. △ Less

Submitted 20 March, 2012; originally announced March 2012.

Comments: 36 pages, 2 figures; Submitted to Progress of Theoretical and Experimental Physics

KNO scaling of fluctuations in pp and pA, and eccentricities in heavy-ion collisions

Authors: Adrian Dumitru, Yasushi Nara

Abstract: Multiplicity fluctuations at midrapidity in p+p collisions at high energies are described by a negative binomial distribution and exhibit approximate Koba-Nielsen-Olesen (KNO) scaling. We find that these KNO fluctuations are important also for reproducing the multiplicity distribution in d+Au collisions observed at RHIC, adding to the Glauber fluctuations of the number of binary collisions or part… ▽ More Multiplicity fluctuations at midrapidity in p+p collisions at high energies are described by a negative binomial distribution and exhibit approximate Koba-Nielsen-Olesen (KNO) scaling. We find that these KNO fluctuations are important also for reproducing the multiplicity distribution in d+Au collisions observed at RHIC, adding to the Glauber fluctuations of the number of binary collisions or participants. We predict that the multiplicity distribution in p+Pb collisions at the LHC also deviates little from the KNO scaling function. Finally, we analyze various moments of the eccentricity of the particle production zone in A+A collisions at RHIC and LHC and find that particle production fluctuations increase fluctuation dominated moments such as the triangularity epsilon_3 substantially. △ Less

Submitted 9 March, 2012; v1 submitted 30 January, 2012; originally announced January 2012.

Comments: 7 pages, 6 figures; v2: updated references and a few minor modifications, to appear in PRC

Report number: RBRC 941

Journal ref: Phys. Rev. C 85, 034907 (2012)

Gluon saturation in $pA$ collisions at the LHC: KLN model predictions for hadron multiplicities

Authors: Adrian Dumitru, Dmitri E. Kharzeev, Eugene M. Levin, Yasushi Nara

Abstract: The upcoming p+Pb run at the LHC will probe the nuclear gluon distribution at very small Bjorken x (from $x \sim 10^{-4}$ at mid-rapidity down to $x \sim 10^{-6}$ in the proton fragmentation region) and will allow to test approaches based on parton saturation. Here, we present the predictions of the KLN model for hadron multiplicities and multiplicity distributions in p+Pb collisions at a center-o… ▽ More The upcoming p+Pb run at the LHC will probe the nuclear gluon distribution at very small Bjorken x (from $x \sim 10^{-4}$ at mid-rapidity down to $x \sim 10^{-6}$ in the proton fragmentation region) and will allow to test approaches based on parton saturation. Here, we present the predictions of the KLN model for hadron multiplicities and multiplicity distributions in p+Pb collisions at a center-of-mass energy of 4.4 TeV. We also compare the model to the existing pp, dA and AA data from RHIC and LHC. △ Less

Submitted 19 October, 2012; v1 submitted 13 November, 2011; originally announced November 2011.

Comments: 7 pages, 7 figures; Added a note with comparison to the pPb LHC data from the ALICE Collaboration

Report number: RBRC 928

Journal ref: Phys. Rev. C 85, 044920 (2012)

Harmonic moments of the gluon density distribution in AA collisions

Authors: Yasushi Nara

Abstract: By using Monte-Carlo implementations of $k_T$-factorization formula with running-coupling BK unintegrated gluon distributions for nucleus-nucleus collisions, we compute higher order harmonic moments of the initial density distribution for both RHIC(Au+Au@200GeV) and LHC(Pb+Pb@2.76TeV) collisions. We study their sensitivity to the size of the valence parton distribution in the nucleon. By using Monte-Carlo implementations of $k_T$-factorization formula with running-coupling BK unintegrated gluon distributions for nucleus-nucleus collisions, we compute higher order harmonic moments of the initial density distribution for both RHIC(Au+Au@200GeV) and LHC(Pb+Pb@2.76TeV) collisions. We study their sensitivity to the size of the valence parton distribution in the nucleon. △ Less

Submitted 1 November, 2011; v1 submitted 13 October, 2011; originally announced October 2011.

Comments: 4 pages, 2 figures, references added, figure 2 is replaced by the results of central collisions

Forward particle productions at RHIC and the LHC from CGC within local rcBK evolution

Authors: H. Fujii, K. Itakura, Y. Kitadono, Y. Nara

Abstract: In order to describe forward hadron productions in high-energy nuclear collisions, we propose a Monte-Carlo implementation of Dumitru-Hayashigaki-Jalilian-Marian formula with the unintegrated gluon distribution obtained numerically from the running-coupling BK equation. We discuss influence of initial conditions for the BK equation by comparing a model constrained by global fit of small-x HERA dat… ▽ More In order to describe forward hadron productions in high-energy nuclear collisions, we propose a Monte-Carlo implementation of Dumitru-Hayashigaki-Jalilian-Marian formula with the unintegrated gluon distribution obtained numerically from the running-coupling BK equation. We discuss influence of initial conditions for the BK equation by comparing a model constrained by global fit of small-x HERA data and a newly proposed one from the running coupling MV model. △ Less

Submitted 7 July, 2011; originally announced July 2011.

Comments: Talk given at conference Quark Matter 2011, 4 pages

Report number: UT-Komaba-11-4, KEK-TH-1473

CGC initial conditions at RHIC and LHC

Authors: Javier L. Albacete, Adrian Dumitru, Yasushi Nara

Abstract: Monte-Carlo implementations of kT -factorization formula with both KLN and running-coupling BK unintegrated gluon distributions for nucleus-nucleus collisions are used to analyze recent experimental data on the particle multiplicities from RHIC(Au+Au@200GeV) and LHC(Pb+Pb@2.76TeV). We also compare the predicted transverse energy at midrapidity to new data from ALICE. Monte-Carlo implementations of kT -factorization formula with both KLN and running-coupling BK unintegrated gluon distributions for nucleus-nucleus collisions are used to analyze recent experimental data on the particle multiplicities from RHIC(Au+Au@200GeV) and LHC(Pb+Pb@2.76TeV). We also compare the predicted transverse energy at midrapidity to new data from ALICE. △ Less

Submitted 6 June, 2011; originally announced June 2011.

Comments: 6 pages and 7 figures. Contribution to the proceedings of the 27th Winter Workshop on Nuclear Dynamics held in Winter Park, Colorado, 6-13 February 2011