BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
P. Huovinen contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We report on the application of a cascade + viscous hydro + cascade model for heavy ion collisions in the RHIC Beam Energy Scan range, $\sqrt{s_{\rm NN}}=6.3\dots200$ GeV. By constraining model parameters to reproduce the data we find that the effective(average) value of the shear viscosity over entropy density ratio $η/s$ decreases from 0.2 to 0.08 when collision energy grows from $\sqrt{s_{\rm NN}}\approx7$ to 39 GeV.
Hybrid approaches based on relativistic hydrodynamics and transport theory have been successfully applied for many years for the dynamical description of heavy ion collisions at ultrarelativistic energies. In this work a new viscous hybrid model employing the hadron transport approach UrQMD for the early and late non-equilibrium stages of the reaction, and 3+1 dimensional viscous hydrodynamics for the hot and dense quark-gluon plasma stage is introduced. This approach includes the equation of motion for finite baryon number, and employs an equation of state with finite net-baryon density to allow for calculations in a large range of beam energies. The parameter space of the model is explored, and constrained by comparison with the experimental data for bulk observables from SPS and the phase I beam energy scan at RHIC. The favored parameter values depend on energy, but allow to extract the effective value of the shear viscosity coefficient over entropy density ratio $η/s$ in the fluid phase for the whole energy region under investigation. The estimated value of $η/s$ increases with decreasing collision energy, which may indicate that $η/s$ of the quark-gluon plasma depends on baryochemical potential $μ_B$.
In most heavy ion collision simulations involving relativistic hydrodynamics, the Cooper-Frye formula is applied to transform the hydrodynamical fields to particles. In this article the so-called negative contributions in the Cooper-Frye formula are studied using a coarse-grained transport approach. The magnitude of negative contributions is investigated as a function of hadron mass, collision energy in the range of $E_{\rm lab} = 5$--$160A$ GeV, collision centrality and the energy density transition criterion defining the hypersurface. The microscopic results are compared to negative contributions expected from hydrodynamical treatment assuming local thermal equilibrium. The main conclusion is that the number of actual microscopic particles flying inward is smaller than the negative contribution one would expect in an equilibrated scenario. The largest impact of negative contributions is found to be on the pion rapidity distribution at midrapidity in central collisions. For this case negative contributions in equilibrium constitute 8--13\% of positive contributions depending on collision energy, but only 0.5--4\% in cascade calculation. The dependence on the collision energy itself is found to be non-monotonous with a maximum at 10-20$A$ GeV.
Many models of heavy ion collisions employ relativistic hydrodynamics to describe the system evolution at high densities. The Cooper-Frye formula is applied in most of these models to turn the hydrodynamical fields into particles. However, the number of particles obtained from the Cooper-Frye formula is not always positive-definite. Physically negative contributions of the Cooper-Frye formula are particles that stream backwards into the hydrodynamical region. We quantify the Cooper-Frye negative contributions in a coarse-grained transport approach, which allows to compare them to the actual number of underlying particles crossing the transition hypersurface. It is found that the number of underlying inward crossings is much smaller than the one the Cooper-Frye formula gives under the assumption of equilibrium distribution functions. The magnitude of Cooper-Frye negative contributions is also investigated as a function of hadron mass, collision energy in the range $E_{\rm lab} = 5-160A$ GeV, and collision centrality. The largest negative contributions we find are around 13% for the pion yield at midrapidity at $E_{\rm lab} = 20A$ GeV collisions.
We apply a 3+1D viscous hydrodynamic + cascade model to the heavy ion collision reactions with $\sqrt{s_{NN}}=6.3\dots39$ GeV. To accommodate the model for a given collision energy range, the initial conditions for hydrodynamic phase are taken from UrQMD, and the equation of state at finite baryon density is based on Chiral model coupled to the Polyakov loop. We study the collision energy dependence of pion and kaon rapidity distributions and $m_T$-spectra, as well as charged hadron elliptic flow and how shear viscosity affects them. The model calculations are compared to the data for Pb-Pb collisions at CERN SPS, as well as for Au-Au collisions in the Beam Energy Scan (BES) program energies at BNL RHIC. The data favours the value of shear viscosity $η/s\gtrsim0.2$ for this collision energy range.
We describe the details of 3+1 dimensional relativistic hydrodynamic code for the simulations of quark-gluon/hadron matter expansion in ultra-relativistic heavy ion collisions. The code solves the equations of relativistic viscous hydrodynamics in the Israel-Stewart framework. With the help of ideal-viscous splitting, we keep the ability to solve the equations of ideal hydrodynamics in the limit of zero viscosities using a Godunov-type algorithm. Milne coordinates are used to treat the predominant expansion in longitudinal (beam) direction effectively. The results are successfully tested against known analytical relativistic inviscid and viscous solutions, as well as against existing 2+1D relativistic viscous code.
Following the experimental program at BNL RHIC, we perform a similar "energy scan" using 3+1D viscous hydrodynamics coupled to the UrQMD hadron cascade, and study the collision energy dependence of pion and kaon rapidity distributions and $m_T$-spectra, as well as charged hadron elliptic flow. To this aim the equation of state for finite baryon density from a Chiral model coupled to the Polyakov loop is employed for hydrodynamic stage. 3D initial conditions from UrQMD are used to study gradual deviation from boost-invariant scaling flow. We find that the inclusion of shear viscosity in the hydrodynamic stage of evolution consistently improves the description of the data for Pb-Pb collisions at CERN SPS, as well as of the elliptic flow measurements for Au-Au collisions in the Beam Energy Scan (BES) program at BNL RHIC. The suggested value of shear viscosity is $η/s\ge0.2$ for $\sqrt{s_{NN}}=6.3\dots39$ GeV.
Relativistic dissipative fluid dynamics is a common tool to describe the space-time evolution of the strongly interacting matter created in ultrarelativistic heavy-ion collisions. For a proper comparison to experimental data, fluid-dynamical calculations have to be performed on an event-by-event basis. Therefore, fluid dynamics should be able to reproduce, not only the event-averaged momentum anisotropies, $<v_{n}>$, but also their distributions. In this paper, we investigate the event-by-event distributions of the initial-state and momentum anisotropies $ε_n$ and $v_n$, and their correlations. We demonstrate that the event-by-event distributions of relative $v_n$ fluctuations are almost equal to the event-by-event distributions of corresponding $ε_n$ fluctuations, allowing experimental determination of the relative anisotropy fluctuations of the initial state. Furthermore, the correlation $c(v_2,v_4)$ turns out to be sensitive to the viscosity of the fluid providing an additional constraint to the properties of the strongly interacting matter.
We study the influence of a temperature-dependent shear viscosity over entropy density ratio $η/s$, different shear relaxation times $τ_π$, as well as different initial conditions on the transverse momentum spectra of charged hadrons and identified particles. We investigate the azimuthal flow asymmetries as a function of both collision energy and centrality. The elliptic flow coefficient turns out to be dominated by the hadronic viscosity at RHIC energies. Only at higher collision energies the impact of the viscosity in the QGP phase is visible in the flow asymmetries. Nevertheless, the shear viscosity near the QCD transition region has the largest impact on the collective flow of the system. We also find that the centrality dependence of the elliptic flow is sensitive to the temperature dependence of $η/s$.
We compare the trace anomaly, strangeness and baryon number fluctuations calculated in lattice QCD with expectations based on hadron resonance gas model. We find that there is a significant discrepancy between the hadron resonance gas and the lattice data. This discrepancy is largely reduced if the hadron spectrum is modified to take into account the larger values of the quark mass used in lattice calculations as well as the finite lattice spacing errors. We also give a simple parametrization of QCD equation of state, which combines hadron resonance gas at low temperatures with lattice QCD at high temperatures. We compare this parametrization with other parametrizations of the equation of state used in hydrodynamical models and discuss differences in hydrodynamic flow for different equations of state.
We study the different nucleus-nucleus collisions, O+Au, S+S, S+Ag, S+Au and Pb+Pb, at the CERN-SPS energy in a one-fluid hydrodynamical approach using a parametrization based on baryon stopping in terms of the thickness of colliding nuclei. Good agreement with measured particle spectra is achieved. We deduce the mass number scaling behaviour of the initial energy density. We find that the equilibration time is nearly independent of the size of the colliding nuclei.