Source author record

H. Niemi

H. Niemi 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

27works
6topics
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

27 published item(s)

preprint2016arXiv

Derivation of anisotropic dissipative fluid dynamics from the Boltzmann equation

Fluid-dynamical equations of motion can be derived from the Boltzmann equation in terms of an expansion around a single-particle distribution function which is in local thermodynamical equilibrium, i.e., isotropic in momentum space in the rest frame of a fluid element. However, in situations where the single-particle distribution function is highly anisotropic in momentum space, such as the initial stage of heavy-ion collisions at relativistic energies, such an expansion is bound to break down. Nevertheless, one can still derive a fluid-dynamical theory, called anisotropic dissipative fluid dynamics, in terms of an expansion around a single-particle distribution function, $\hat{f}_{0\bf k}$, which incorporates (at least parts of) the momentum anisotropy via a suitable parametrization. We construct such an expansion in terms of polynomials in energy and momentum in the direction of the anisotropy and of irreducible tensors in the two-dimensional momentum subspace orthogonal to both the fluid velocity and the direction of the anisotropy. From the Boltzmann equation we then derive the set of equations of motion for the irreducible moments of the deviation of the single-particle distribution function from $\hat{f}_{0\bf k}$. Truncating this set via the 14-moment approximation, we obtain the equations of motion of anisotropic dissipative fluid dynamics.

preprint2016arXiv

Hydrodynamic modeling of a pure-glue initial scenario in high-energy hadron and heavy-ion collisions

Partonic matter produced in the early stage of ultrarelativistic nucleus-nucleus collisions is assumed to be composed mainly of gluons, and quarks and antiquarks are produced at later times. The comparable hydrodynamic simulations of heavy-ion collisions for (2+1)-flavor and Yang-Mills equations of state performed by using three different hydrodynamic codes are presented. Assuming slow chemical equilibration of quarks, the spectra and elliptic flows of thermal dileptons and photons are calculated for central Pb+Pb collisions at the LHC energy of $\sqrt{s_{_{\rm NN}}} = 2.76$ TeV. It is shown that a suppression of quarks at early times leads to a significant reduction of the yield of the thermal dileptons, but only to a rather modest suppression of the $p_T$-distribution of direct photons. It is demonstrated that an enhancement of photon and dilepton elliptic flows might serve as a promising signature of the pure-glue initial state. Calculations based on Bjorken hydrodynamics suggest that collisions of small systems at intermediate energies available at RHIC or future FAIR facilities may show stronger effects associated with initial pure gluodynamic evolution.

preprint2015arXiv

Event-by-event fluctuations in perturbative QCD + saturation + hydro model: pinning down QCD matter shear viscosity in ultrarelativistic heavy-ion collisions

We introduce an event-by-event perturbative-QCD + saturation + hydro ("EKRT") framework for ultrarelativistic heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order perturbative QCD using a saturation conjecture to control soft particle production, and describe the space-time evolution of the QCD matter with dissipative fluid dynamics, event by event. We perform a simultaneous comparison of the centrality dependence of hadronic multiplicities, transverse momentum spectra, and flow coefficients of the azimuth-angle asymmetries, against the LHC and RHIC measurements. We compare also the computed event-by-event probability distributions of relative fluctuations of elliptic flow, and event-plane angle correlations, with the experimental data from Pb+Pb collisions at the LHC. We show how such a systematic multi-energy and multi-observable analysis tests the initial state calculation and the applicability region of hydrodynamics, and in particular how it constrains the temperature dependence of the shear viscosity-to-entropy ratio of QCD matter in its different phases in a remarkably consistent manner.

preprint2015arXiv

Pinning down QCD-matter shear viscosity in A+A collisions via EbyE fluctuations using pQCD + saturation + hydrodynamics

We compute the initial energy densities produced in ultrarelativistic heavy-ion collisions from NLO perturbative QCD using a saturation conjecture to control soft particle production, and describe the subsequent space-time evolution of the system with hydrodynamics, event by event. The resulting centrality dependence of the low-$p_T$ observables from this pQCD + saturation + hydro ("EKRT") framework are then compared simultaneously to the LHC and RHIC measurements. With such an analysis we can test the initial state calculation, and constrain the temperature dependence of the shear viscosity-to-entropy ratio $η/s$ of QCD matter. Using these constraints from the current RHIC and LHC measurements we then predict the charged hadron multiplicities and flow coefficients for the 5.023 TeV Pb+Pb collisions.

preprint2015arXiv

Pinning down QCD-matter shear viscosity in ultrarelativistic heavy-ion collisions via EbyE fluctuations using pQCD + saturation + hydrodynamics

We introduce an event-by-event pQCD + saturation + hydro ("EKRT") framework for high-energy heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order (NLO) perturbative QCD (pQCD) using saturation to control soft particle production, and describe the space-time evolution of the QCD matter with viscous hydrodynamics, event by event (EbyE). We compare the computed centrality dependence of hadronic multiplicities, p_T spectra and flow coefficients v_n against LHC and RHIC data. We compare also the computed EbyE probability distributions of relative fluctuations of v_n, as well as correlations of 2 and 3 event-plane angles, with LHC data. Our systematic multi-energy and -observable analysis not only tests the initial state calculation and applicability of hydrodynamics, but also makes it possible to constrain the temperature dependence of the shear viscosity-to-entropy ratio, eta/s(T), of QCD matter in its different phases. Remarkably, we can describe all these different flow observables and correlations consistently with eta/s(T) that is independent of the collision energy.

preprint2015arXiv

Predictions for 5.023 TeV Pb+Pb collisions at the LHC

We compute predictions for various low-transverse-momentum bulk observables in $\sqrt{s_{NN}} = 5.023$ TeV Pb+Pb collisions at the LHC from the event-by-event next-to-leading-order perturbative-QCD + saturation + viscous hydrodynamics ("EKRT") model. In particular, we consider the centrality dependence of charged hadron multiplicity, flow coefficients of the azimuth-angle asymmetries and correlations of event-plane angles. The centrality dependencies of the studied observables are predicted to be very similar to those at 2.76 TeV, and the magnitudes of the flow coefficients and event-plane angle correlations are predicted to be close to those at 2.76 TeV. The flow coefficients may, however, offer slightly more discriminating power on the temperature dependence of QCD matter viscosity than the 2.76 TeV measurements. Our prediction for the multiplicity in the 0-5\% centrality class, obtained using the two temperature-dependent shear-viscosity-to-entropy ratios that give the best overall fit to RHIC and LHC data is $dN_{\rm ch}/dη\big|_{|η|\le 0.5} =1876\dots2046$. We also predict a power-law increase from 200 GeV Au+Au collisions at RHIC to 2.76 and 5.023 TeV Pb+Pb collisions at the LHC, $dN_{\rm ch}/dη\big|_{|η|\le 0.5} \propto s^{0.164\dots0.174}$.

preprint2014arXiv

How large is the Knudsen number reached in fluid dynamical simulations of ultrarelativistic heavy ion collisions?

We investigate the applicability of fluid dynamics in ultrarelativistic heavy ion (AA) collisions and high multiplicity proton nucleus (pA) collisions. In order for fluid dynamics to be applicable the microscopic and macroscopic distance/time scales of the system have to be sufficiently separated. The degree of separation is quantified by the ratio between these scales, usually referred to as the Knudsen number. In this work, we calculate the Knudsen numbers reached in fluid dynamical simulations of AA and pA collisions at RHIC and LHC energies. For this purpose, we consider different choices of shear viscosity parametrizations, initial states and initialization times. We then estimate the values of shear viscosity for which the fluid dynamical description of ultrarelativistic AA and pA collisions breaks down. In particular, we study how such values depend on the centrality, in the case of AA collision, and multiplicity, in the case of pA collision. We found that the maximum viscosity in AA collisions is of the order $η/s \sim 0.1 \ldots 0.2$, which is similar in magnitude to the viscosities currently employed in simulations of heavy ion collisions. For pA collisions, we found that such limit is significantly lower, being less than $η/s=0.08$

preprint2014arXiv

Next-to-leading order improved perturbative QCD + saturation + hydrodynamics model for A+A collisions

We calculate initial conditions for the hydrodynamical evolution in ultrarelativistic heavy-ion collisions at the LHC and RHIC in an improved next-to-leading order perturbative QCD + saturation framework. Using viscous relativistic hydrodynamics, we show that we obtain a good simultaneous description of the centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow at the LHC and at RHIC. In particular, we discuss how the temperature dependence of the shear viscosity is constrained by these data.

preprint2014arXiv

On the relative importance of second-order terms in relativistic dissipative fluid dynamics

In Denicol et al., Phys. Rev. D 85, 114047 (2012), the equations of motion of relativistic dissipative fluid dynamics were derived from the relativistic Boltzmann equation. These equations contain a multitude of terms of second order in Knudsen number, in inverse Reynolds number, or their product. Terms of second order in Knudsen number give rise to non-hyperbolic (and thus acausal) behavior and must be neglected in (numerical) solutions of relativistic dissipative fluid dynamics. The coefficients of the terms which are of the order of the product of Knudsen and inverse Reynolds numbers have been explicitly computed in the above reference, in the limit of a massless Boltzmann gas. Terms of second order in inverse Reynolds number arise from the collision term in the Boltzmann equation, upon expansion to second order in deviations from the single-particle distribution function in local thermodynamical equilibrium. In this work, we compute these second-order terms for a massless Boltzmann gas with constant scattering cross section. Consequently, we assess their relative importance in comparison to the terms which are of the order of the product of Knudsen and inverse Reynolds numbers.

preprint2013arXiv

Fluid dynamics with saturated minijet initial conditions in ultrarelativistic heavy-ion collisions

Using next-to-leading order perturbative QCD and a conjecture of saturation to suppress the production of low-energy partons, we calculate the initial energy densities and formation times for the dissipative fluid dynamical evolution of the quark-gluon plasma produced in ultrarelativistic heavy-ion collisions. We identify the framework uncertainties and demonstrate the predictive power of the approach by a good global agreement with the measured centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow simultaneously for the Pb+Pb collisions at the LHC and Au+Au at RHIC. In particular, the shear viscosity in the different phases of QCD matter is constrained in this new framework simultaneously by all these data.

preprint2013arXiv

Transition From Ideal To Viscous Mach Cones In A Kinetic Transport Approach

Using a microscopic transport model we investigate the evolution of conical structures originating from the supersonic projectile moving through the hot matter of ultrarelativistic particles. Using different scenarios for the interaction between projectile and matter, and different transport properties of the matter, we study the formation and structure of Mach cones. Especially, a dependence of the Mach cone angle on the details and rate of the energy deposition from projectile to the matter is investigated. Furthermore, the two-particle correlations extracted from the numerical calculations are compared to an analytical approximation. We find that the propagation of a high energetic particle through the matter does not lead to the appearance of a double peak structure as observed in the ultrarelativistic heavy-ion collision experiments. The reason is the strongly forward-peaked energy and momentum deposition in the head shock region. In addition, by adjusting the cross section we investigate the influence of the viscosity to the structure of Mach cones. A clear and unavoidable smearing of the profile depending on a finite ratio of shear viscosity to entropy density is clearly visible.

preprint2012arXiv

Derivation of fluid dynamics from kinetic theory with the 14--moment approximation

We review the traditional derivation of the fluid-dynamical equations from kinetic theory according to Israel and Stewart. We show that their procedure to close the fluid-dynamical equations of motion is not unique. Their approach contains two approximations, the first being the so-called 14-moment approximation to truncate the single-particle distribution function. The second consists in the choice of equations of motion for the dissipative currents. Israel and Stewart used the second moment of the Boltzmann equation, but this is not the only possible choice. In fact, there are infinitely many moments of the Boltzmann equation which can serve as equations of motion for the dissipative currents. All resulting equations of motion have the same form, but the transport coefficients are different in each case.

preprint2012arXiv

Derivation of transient relativistic fluid dynamics from the Boltzmann equation

In this work we present a general derivation of relativistic fluid dynamics from the Boltzmann equation using the method of moments. The main difference between our approach and the traditional 14-moment approximation is that we will not close the fluid-dynamical equations of motion by truncating the expansion of the distribution function. Instead, we keep all terms in the moment expansion. The reduction of the degrees of freedom is done by identifying the microscopic time scales of the Boltzmann equation and considering only the slowest ones. In addition, the equations of motion for the dissipative quantities are truncated according to a systematic power-counting scheme in Knudsen and inverse Reynolds number. We conclude that the equations of motion can be closed in terms of only 14 dynamical variables, as long as we only keep terms of second order in Knudsen and/or inverse Reynolds number. We show that, even though the equations of motion are closed in terms of these 14 fields, the transport coefficients carry information about all the moments of the distribution function. In this way, we can show that the particle-diffusion and shear-viscosity coefficients agree with the values given by the Chapman-Enskog expansion.

preprint2012arXiv

Event-by-event distributions of azimuthal asymmetries in ultrarelativistic heavy-ion collisions

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.

preprint2012arXiv

Influence of a temperature-dependent shear viscosity on the azimuthal asymmetries of transverse momentum spectra in ultrarelativistic heavy-ion collisions

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

preprint2012arXiv

Solving the heat-flow problem with transient relativistic fluid dynamics

Israel-Stewart theory is a causal, stable formulation of relativistic dissipative fluid dynamics. This theory has been shown to give a decent description of the dynamical behavior of a relativistic fluid in cases where shear stress becomes important. In principle, it should also be applicable to situations where heat flow becomes important. However, it has been shown that there are cases where Israel-Stewart theory cannot reproduce phenomena associated with heat flow. In this paper, we derive a relativistic dissipative fluid-dynamical theory from kinetic theory which provides a good description of all dissipative phenomena, including heat flow. We explicitly demonstrate this by comparing this theory with numerical solutions of the relativistic Boltzmann equation.

preprint2011arXiv

Determination of the Shear Viscosity Relaxation Time at Weak and Strong Coupling

We investigate the microscopic origin of the relaxation time coefficient in relativistic fluid dynamics. We show that the extraction of the shear viscosity relaxation time via the gradient expansion is ambiguous and in general fails to give the correct result. The correct value for the shear viscosity relaxation time is extracted from the slowest non-hydrodynamic pole of the corresponding retarded Green's function, if such a pole is purely imaginary. According to the AdS/CFT correspondence, in strongly-coupled $\mathcal{N}=4$ SYM the non-hydrodynamic poles of the shear stress tensor nearest to the origin have a nonzero real part, which implies that the transient fluid-dynamical equations for this gauge theory are not equivalent to the well-known Israel-Stewart equations.

preprint2011arXiv

Origin of the Relaxation Time in Dissipative Fluid Dynamics

We show how the linearized equations of motion of any dissipative current are determined by the analytical structure of the associated retarded Green's function. If the singularity of the Green's function, which is nearest to the origin in the complex-frequency plane, is a simple pole on the imaginary frequency axis, the linearized equations of motion can be reduced to relaxation-type equations for the dissipative currents. The value of the relaxation time is given by the inverse of this pole. We prove that, if the relaxation time is sent to zero, or equivalently, the pole to infinity, the dissipative currents approach the values given by the standard gradient expansion.

preprint2011arXiv

Second order dissipative fluid dynamics from kinetic theory

We derive the equations of second order dissipative fluid dynamics from the relativistic Boltzmann equation following the method of W. Israel and J. M. Stewart. We present a frame independent calculation of all first- and second-order terms and their coefficients using a linearised collision integral. Therefore, we restore all terms that were previously neglected in the original papers of W. Israel and J. M. Stewart.

preprint2010arXiv

Investigation of shock waves in the relativistic Riemann problem: A comparison of viscous fluid dynamics to kinetic theory

We solve the relativistic Riemann problem in viscous matter using the relativistic Boltzmann equation and the relativistic causal dissipative fluid-dynamical approach of Israel and Stewart. Comparisons between these two approaches clarify and point out the regime of validity of second-order fluid dynamics in relativistic shock phenomena. The transition from ideal to viscous shocks is demonstrated by varying the shear viscosity to entropy density ratio $η/s$. We also find that a good agreement between these two approaches requires a Knudsen number $Kn < 1/2$.

preprint2010arXiv

Mach Cones in Viscous Matter

Employing a microscopic transport model we investigate the evolution of high energetic jets moving through a viscous medium. For the scenario of an unstoppable jet we observe a clearly strong collective behavior for a low dissipative system $η/s \approx 0.005$, leading to the observation of cone-like structures. Increasing the dissipation of the system to $η/s \approx 0.32$ the Mach Cone structure vanishes. Furthermore, we investigate jet-associated particle correlations. A double-peak structure, as observed in experimental data, is even for low-dissipative systems not supported, because of the large influence of the head shock.

preprint2009arXiv

Numerical tests of causal relativistic dissipative fluid dynamics

We present numerical methods to solve the Israel-Stewart (IS) equations of causal relativistic dissipative fluid dynamics with bulk and shear viscosities. We then test these methods studying the Riemann problem in (1+1)-- and (2+1)-dimensional geometry. The numerical schemes investigated here are applicable to realistic (3+1)--dimensional modeling of a relativistic dissipative fluid.

preprint2009arXiv

Relativistic shock waves in viscous gluon matter

We solve the relativistic Riemann problem in viscous gluon matter employing a microscopic parton cascade. We demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to entropy density ratio $η/s$ from zero to infinity. We show that an $η/s$ ratio larger than 0.2 prevents the development of well-defined shock waves on timescales typical for ultrarelativistic heavy-ion collisions. Comparisons with viscous hydrodynamic calculations confirm our findings.

preprint2005arXiv

Transverse Spectra of Hadrons in Central $AA$ Collisions at RHIC and LHC from pQCD+Saturation+Hydrodynamics and from pQCD+Energy Losses

We study the transverse spectra of hadrons in nearly central $AA$ collisions at RHIC and LHC in a broad transverse momentum range Low-$p_T$ spectra are calculated by using boost-invariant hydrodynamics with initial energy and net-baryon densities from the EKRT pQCD+saturation model. High-$p_T$ spectra are obtained from pQCD jet calculation including the energy loss of the parton in the matter prior to its fragmentation to final hadrons.

preprint2004arXiv

Photon production from non-equilibrium QGP in heavy ion collisions

We present a calculation of thermal photon production i.e. photons from secondary interactions among particles produced in heavy ion collisions at collider energies. This is done within the framework of hydrodynamics. We take into account the lack of chemical equilibrium in QGP. It turns out that main effects from chemical non-equilibrium composition of QGP, reduction of particle number and increase in temperature, nearly cancel in photon spectrum.