Source author record

Yasushi Nara

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

33 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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}$

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Ultra-forward particle production from CGC+Lund fragmentation

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

Integrated Dynamical Approach to Relativistic Heavy Ion Collisions

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.

preprint2013arXiv

Magnetic flux loop in high-energy heavy-ion collisions

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.

preprint2013arXiv

Magnetic screening in high-energy heavy-ion collisions

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.

preprint2012arXiv

CGC predictions for p+Pb collisions at the LHC

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

Harmonic moments of the gluon density distribution in AA collisions

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.

preprint2010arXiv

Elliptic flow in Pb+Pb collisions at sqrt{s_{NN}} = 2.76 TeV: hybrid model assessment of the first data

We analyze the elliptic flow parameter v_2 in Pb+Pb collisions at sqrt{s_{NN}} = 2.76 TeV and in Au+Au collisions at sqrt{s_{NN}} =200 GeV using a hybrid model in which the evolution of the quark gluon plasma is described by ideal hydrodynamics with a state-of-the-art lattice QCD equation of state, and the subsequent hadronic stage by a hadron cascade model. For initial conditions, we employ Monte-Carlo versions of the Glauber and the Kharzeev-Levin-Nardi models and compare results with each other. We demonstrate that the differential elliptic flow v_2(p_T) hardly changes when the collision energy increases, whereas the integrated v_2 increases due to the enhancement of mean transverse momentum. The amount of increase of both v_2 and mean p_T depends significantly on the model of initialization.

preprint2010arXiv

Elliptic flow in U+U collisions at sqrt{s_{NN}} = 200 GeV and in Pb+Pb collisions at sqrt{s_{NN}} = 2.76 TeV: Prediction from a hybrid approach

We predict the elliptic flow parameter v_2 in U+U collisions at sqrt{s_{NN}}=200 GeV and in Pb+Pb collisions at sqrt{s_{NN}} = 2.76 TeV using a hybrid model in which the evolution of the quark gluon plasma is described by ideal hydrodynamics with a state-of-the-art lattice QCD equation of state, and the subsequent hadronic stage by a hadron cascade model.

preprint2009arXiv

Transverse momentum diffusion and jet energy loss in non-Abelian plasmas

We consider momentum broadening and energy loss of high momentum partons in a hot non-Abelian plasma due to collisions. We solve the coupled system of Wong-Yang-Mills equations on a lattice in real time, including binary hard elastic collisions among the partons. The collision kernel is constructed such that the total collisional energy loss and momentum broadening are lattice spacing independent. We find that the transport coefficient $\hat{q}$ corresponding to transverse momentum broadening receives sizable contributions from a power-law tail in the $p_\perp$-distribution of high-momentum partons. We establish the scaling of $\hat{q}$ and of $dE/dx$ with density, temperature and energy in the weak-coupling regime. We also estimate the nuclear modification factor $R_{AA}$ due to elastic energy loss of a jet in a classical Yang-Mills field.

preprint2003arXiv

Classical Gluodynamics of High Energy Nuclear Collisions: an Erratum and an Update

We comment on the relation of our previous work on the classical gluodynamics of high energy nuclear collisions to recent work by Lappi (hep-ph/0303076). While our results for the non-perturbative number liberation coefficient agree, those for the energy disagree by a factor of 2. This discrepancy can be traced to an overall normalization error in our non-perturbative formula for the energy. When corrected for, all previous results are in excellent agreement with those of Lappi. The implications of the results of these two independent computations for RHIC phenomenology are noted.

preprint2003arXiv

Hadrons as Signature of Black Hole Production at the LHC

In models with several large extra dimensions and fundamental Planck scale of the order of 1 TeV, black holes can be produced in large numbers at LHC energies. We compute the charged hadron spectra obtained from the decay of black holes created in p-p and Pb+Pb collisions at LHC. We show that hadrons from black hole decay dominate at transverse momenta p_T>30-100 GeV/c compared to usual QCD processes and black hole production signals are easy to identify in hadron transverse momentum spectra. Furthermore we show that a measurement of the charged hadron spectra probes Planck scales up to ~5 TeV for any number of extra dimensions.

preprint2003arXiv

The Cronin Effect, Quantum Evolution and the Color Glass Condensate

We show that the numerical solution of the classical SU(3) Yang-Mills equations of motion in the McLerran-Venugopalan model for gluon production in central heavy ion collisions leads to a suppresion at low $p_t$ and an enhancement at the intermediate $p_t$ region as compared to peripheral heavy ion and pp collisions at the same energy. Our results are compared to previous, Color Glass Condensate inspired calculations of gluon production in heavy ion collisions. We revisit the predictions of the Color Glass Condensate model for $pA$ ($dA$) collisions in Leading Order and show that quantum evolution--in particular the phenomenon of geometric scaling and change of anomalous dimensions--preserves the Cronin enhancement of $pA$ cross section (when normalized to the leading twist term) in the Leading Order approximation even though the $p_t$ spectrum can change. We comment on the case when gluon radiation is included.

preprint2002arXiv

Elliptic flow of colored glass in high energy heavy ion collisions

We compute the elliptic flow generated by classical gluon fields in a high energy nuclear collision. The classical gluon fields are described by a typical momentum scale, the saturation scale Lambda_s, which is, for RHIC energies, of the order of 1-2 GeV. A significant elliptic flow is generated only over time scales on the order of the system size R. The flow is dominated by soft modes $p_{T} \sim Lambda_s / 4$ which linearize at very late times $tau \sim R \gg 1/Lambda_s$. We discuss the implications of our result for the theoretical interpretation of the RHIC data.