Source author record

Carsten Greiner

Carsten Greiner 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

80works
10topics
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

80 published item(s)

preprint2022arXiv

Incomplete electromagnetic response of hot QCD matter

The electromagnetic response of hot QCD matter to decaying external magnetic fields is investigated. We examine the validity of Ohm's law and find that the induced electric current increases from zero and relaxes towards the value from Ohm's law. The relaxation time is larger than the lifetime of the external magnetic field for the QCD matter in relativistic heavy-ion collisions. The lower than expected electric current significantly suppresses the induced magnetic field and makes the electromagnetic response incomplete. We demonstrate the incomplete electromagnetic response of hot QCD matter by calculations employing the parton transport model combined with the solution of Maxwell's equations. Our results show a strong suppression by two orders of magnitude in the magnetic field, relatively to calculations assuming the validity of Ohm's law. This may undermine experimental efforts to measure magnetic-field-related effects in heavy-ion collisions.

preprint2022arXiv

Multicomponent relativistic dissipative fluid dynamics from the Boltzmann equation

We derive multicomponent relativistic second-order dissipative fluid dynamics from the Boltzmann equations for a reactive mixture of $N_{\text{spec}}$ particle species with $N_q$ intrinsic quantum numbers (e.g. electric charge, baryon number, and strangeness) using the method of moments. We obtain the continuity equations for multiple conserved charges as well as the conservation equations for the total energy and momentum in the single-fluid approximation. These $4+N_q$ conservation laws are closed by deriving the second-order equations of motion for the dissipative quantities in the $(10+4N_q)$-moment approximation. The resulting fluid-dynamical equations are formally similar to those of a single-component system, but feature different thermodynamic relations and transport coefficients. We derive general relations for all transport coefficients and compute them explicitly in the ultrarelativistic limit.

preprint2021arXiv

Early Universe Thermodynamics and Evolution in Nonviscous and Viscous Strong and Electroweak epochs: Possible Analytical Solutions

Based on recent perturbative and non-perturbative lattice calculations with almost quark flavors and the thermal contributions from photons, neutrinos, leptons, electroweak particles, and scalar Higgs bosons, various thermodynamic quantities, at vanishing net-baryon densities, such as pressure, energy density, bulk viscosity, relaxation time, and temperature have been calculated up to the TeV-scale, i.e. covering hadron, QGP and electroweak (EW) phases in the early Universe. This remarkable progress motivated the present study to determine the possible influence of the bulk viscosity in the early Universe and to understand how this would vary from epoch to epoch. We have taken into consideration first- (Eckart) and second-order (Israel-Stewart) theories for the relativistic cosmic fluid and integrated viscous equations of state in Friedmann equations. Nonlinear nonhomogeneous differential equations are obtained as analytical solutions. For Israel-Stewart, the differential equations are very sophisticated to be solved. They are outlined here as road-maps for future studies. For Eckart theory, the only possible solution is the functionality, $H(a(t))$, where $H(t)$ is the Hubble parameter and $a(t)$ is the scale factor, but none of them so far could to be directly expressed in terms of either proper or cosmic time $t$. For Eckart-type viscous background, especially at finite cosmological constant, non-singular $H(t)$ and $a(t)$ are obtained, where $H(t)$ diverges for QCD/EW and asymptotic EoS. For non-viscous background, the dependence of $H(a(t))$ is monotonic. The same conclusion can be drawn for an ideal EoS. We also conclude that the rate of decreasing $H(a(t))$ with increasing $a(t)$ varies from epoch to epoch, at vanishing and finite cosmological constant. These results obviously help in improving our understanding of the nucleosynthesis and the cosmological large-scale structure.

preprint2021arXiv

Solving the puzzle of high temperature light (anti)-nuclei production in ultra-relativistic heavy ion collisions

The creation of loosely bound objects in heavy ion collisions, e.g. light clusters, near the phase transition temperature $\left(T_{\rm ch} \approx 155 \, \rm{MeV} \right)$ has been a puzzling observation that seems to be at odds with Big Bang nucleosynthesis suggesting that deuterons and other clusters are formed only below a temperature $T\approx 0.1-1 \, \text{MeV}$. We solve this puzzle by showing that the light cluster abundancies in heavy ion reactions stay approximately constant from chemical freeze-out to kinetic freeze-out. To this aim we develop an extensive network of coupled reaction rate equations including stable hadrons and hadronic resonances to describe the temporal evolution of the abundancies of light (anti-)(hyper-)nuclei in the late hadronic environment of an ultrarelativistic heavy ion collision. It is demonstrated that the chemical equilibration of the light nuclei occurs on a very short timescale as a consequence of the strong production and dissociation processes. However, because of the partial chemical equilibrium of the stable hadrons, including the nucleon feeding from $Δ$ resonances, the abundancies of the light nuclei stay nearly constant during the evolution and cooling of the hadronic phase. This solves the longstanding contradiction between the thermal fits and the late stage coalescence (and the Big Bang nucleosynthesis) and explains why the observed light cluster yields are compatible with both a high chemical production temperature and a late state emission as modelled by coalescence. We also note in passing that the abundancies of the light clusters in the present approach are in excellent agreement with those measured by ALICE at LHC.

preprint2020arXiv

Cross-conductivity: novel transport coefficients to constrain the hadronic degrees of freedom of nuclear matter

In general, the constituents of the bulk matter produced in heavy-ion collisions carry, besides electric charge, multiple other conserved quantum numbers like baryon number and strangeness. Therefore, an electric field will not only generate an electric current but, at the same time, also currents in baryon number and strangeness. We propose that the impact of the electric field on these conserved currents should be characterized by additional transport coefficients, which we call cross-conductivities. In this paper, we introduce and present a calculation of these cross-conductivities from the Green-Kubo formalism within the transport code SMASH for different chemical compositions of hadron resonance gases. We find that the coefficients underlie an ordering in the active degrees of freedom and that thus the chemical composition of the system plays a crucial role. Further, we argue that in future comparisons of lattice QCD calculations with these findings, one could constrain which degrees of freedom and their corresponding charge properties are relevant for the QCD dynamics of the system.

preprint2020arXiv

Effects of fluctuations and color-neutrality in a finite volume

We investigate properties of strongly interacting matter in a schematic model, based on the combined degrees of freedom of a non-interacting hadronic phase and a non-interacting deconfined phase. It is found that in a finite system both phases contribute to the thermodynamic state due to fluctuations and that signatures of critical behviour like the divergence of statistical quantities are damped. The constraint of color-neutrality leads to a volume-dependent shift of the effective critical temperature, which follows a scaling law, independent of the baryochemical potential. According to the model, observable baryon-number susceptibilities at a given $T$ and $μ_B$ strongly depend on the system size. Finally, we compare hadronization conditions from the model with hadrochemical fits to experimental collider data, where a qualitatively similar system size dependence is extracted.

preprint2020arXiv

Kinetic freeze-out temperature from yields of short-lived resonances

A method to determine the kinetic freeze-out temperature in heavy-ion collisions from measured yields of short-lived resonances is presented. The resonance production is treated in the framework of thermal model with an evolution between chemical and kinetic freeze-outs. The yields of many short-lived resonances are suppressed at $T = T_{\rm kin} < T_{\rm ch}$. We determine the values of $T_{\rm kin}$ and $T_{\rm ch}$ for various centralities in Pb--Pb collisions at $\sqrt{s_{_{NN}}} = 2.76$ TeV by fitting the abundances of both the stable hadrons and the short-lived resonances such as $ρ^0$ and $ \text{K}^{*0}$, that were measured by the ALICE collaboration. This allows to extract the kinetic freeze-out temperature from the measured hadron and resonance yields alone, independent of assumptions about the flow velocity profile and the freeze-out hypersurface. The extracted $T_{\rm ch}$ values exhibit a moderate multiplicity dependence whereas $T_{\rm kin}$ drops, from $T_{\rm kin} \simeq T_{\rm ch} \simeq 155$ MeV in peripheral collisions to $T_{\rm kin} \simeq 110$ MeV in 0-20% central collisions. Predictions for other short-lived resonances are presented. A potential (non-)observation of a suppressed $f_0(980)$ meson yield will allow to constrain the lifetime of that meson.

preprint2020arXiv

Open charm and dileptons from relativistic heavy-ion collisions

We study the dynamics of open charm production and the dilepton radiation of the semi-leptonic decays of correlated $D\bar D$ pairs versus the quark-gluon plasma (QGP) radiation and hadronic sources in relativistic heavy-ion collisions. Our study is based on the Parton-Hadron-String Dynamics (PHSD) transport approach employing a non-perturbative QCD description of the strongly interacting quark-gluon plasma (sQGP) in terms of dynamical quasiparticles and the EoS based on lattice QCD. We compare the PHSD results for charm observables with the calculations from BAMPS (Boltzmann Approach to Multi-Parton Scatterings) which is based on perturbative QCD with massless partons and interaction cross sections calculated in leading order of the QCD coupling. We compare the $p_T$ dependence of the ratio $R_{AA}$ of $D$-mesons in $A+A$ over $p+p$ collisions scaled by the number of binary collisions $N_{bin}$ as well as the elliptic flow $v_2$ of $D$-mesons calculated within the PHSD and BAMPS at LHC energies. In other study, based on the PHSD calculations we find that the dileptons from correlated $D-$meson semi-leptonic decays dominate the 'thermal' radiation from the QGP in central Pb+Pb collisions at the intermediate masses ($1.2 < M < 3$ GeV) for higher invariant energies However, for invariant energies $\sqrt{s_{NN}} < 40$ GeV the QGP radiation overshines the contribution from $D,{\bar D}$ decays such that one should observe a rather clear signal from the partonic dilepton radiation. This finding provides promising perspectives to measure the QGP radiation in the dilepton experiments at RHIC BES and the future FAIR/NICA facilities.

preprint2019arXiv

Calculation of anisotropic transport coefficients for an ultrarelativistic Boltzmann gas in a magnetic field within a kinetic approach

According to the Kubo formulas we employ the (3+1)-d parton cascade, Boltzmann approach of multiparton scatterings (BAMPS), to calculate the anisotropic transport coefficients (shear viscosity and electric conductivity) for an ultrarelativistic Boltzmann gas in the presence of a magnetic field. The results are compared with those recently obtained by using the Grad's approximation. We find good agreements between both results, which confirms the general use of the derived Kubo formulas for calculating the anisotropic transport coefficients of quark-gluon plasma in a magnetic field.

preprint2019arXiv

Critical point signatures in the cluster expansion in fugacities

The QCD baryon number density can formally be expanded into a Laurent series in fugacity, which is a relativistic generalization of Mayer's cluster expansion. We determine properties of the cluster expansion in a model with a phase transition and a critical point at finite baryon density, in which the Fourier coefficients of the expansion can be determined explicitly and to arbitrary order. The asymptotic behavior of Fourier coefficients changes qualitatively as one traverses the critical temperature and it is connected to the branch points of a thermodynamic potential associated with the phase transition. The results are discussed in the context of lattice QCD simulations at imaginary chemical potential. We argue that the location of a branch point closest to the imaginary chemical potential axis can be extracted through an analysis of an exponential suppression of Fourier coefficients. This is illustrated using the four leading coefficients both in a toy model as well as by using recent lattice QCD data.

preprint2019arXiv

Diffusion processes involving multiple conserved charges: a first study from kinetic theory and implications to the fluid-dynamical modeling of heavy ion collisions

The bulk nuclear matter produced in heavy ion collisions carries a multitude of conserved quantum numbers: electric charge, baryon number, and strangeness. Therefore, the diffusion processes associated to these conserved charges cannot occur independently and must be described in terms of a set of coupled diffusion equations. This physics is implemented by replacing the traditional diffusion coefficients for each conserved charge by a diffusion coefficient matrix, which quantifies the coupling between the conserved quantum numbers. The diagonal coefficients of this matrix are the usual charge diffusion coefficients, while the off-diagonal entries describe the diffusive coupling of the charge currents. In this paper, we show how to calculate this diffusion coefficient matrix from kinetic theory and provide results for a hadron resonance gas and a gas of partons. We further find that the off-diagonal entries can reach similar magnitudes compared to the diagonal entries. In order to provide some insight on the influence that the coupling between the net charge diffusion currents can have on heavy ion observables, we present first results for the diffusive evolution of a hadronic system in a simple (1+1)D-fluid dynamics approach, and study different configurations of the diffusion matrix.

preprint2019arXiv

Polyakov linear-sigma model in mean-field approximation and optimized perturbation theory

We compare results from the Polyakov linear-sigma model (PLSM) in optimized perturbation theory (OPT) with the mean-field approximation (MFA). At finite temperatures and chemical potentials, the chiral condensates and the decofinement order parameters, the thermodynamic pressure, the pseudo-critical temperatures, the subtracted condensates, the second- and high-order moments of various conserved charges (cumulants) obtained in MFA are compared with OPT and also confronted to available lattice QCD simulations. We conclude that when moving from lower- to higher-order moments of various quantum charges, OPT becomes more closer to QCD.

preprint2019arXiv

The Shear Viscosity to Entropy Density Ratio of Hagedorn States

The fireball concept of Rolf Hagedorn, developed in the 1960's, is an alternative description of hadronic matter. Using a recently derived mass spectrum, we use the transport model GiBUU to calculate the shear viscosity of a gas of such Hagedorn states, applying the Green-Kubo method to Monte-Carlo calculations. Since the entropy density is rising ad infinitum near $T_H$, this leads to a very low shear viscosity to entropy density ratio near $T_H$. Further, by comparing our results with analytic expressions, we find a nice extrapolation behavior, indicating that a gas of Hagedorn states comes close or even below the boundary $1/4π$ from AdS-CFT.

preprint2016arXiv

Thermal Charm and Charmonium Production in Quark Gluon Plasma

We study the effect of thermal charm production on charmonium regeneration in high energy nuclear collisions. By solving the kinetic equations for charm quark and charmonium distributions in Pb+Pb collisions, we calculate the global and differential nuclear modification factors $R_{AA}(N_{part})$ and $R{AA}(p_t)$ for $J/Ψ$s. Due to the thermal charm production in hot medium, the charmonium production source changes from the initially created charm quarks at SPS, RHIC and LHC to the thermally produced charm quarks at Future Circular Collider (FCC), and the $J/Ψ$ suppression ($R_{AA}<1$) observed so far will be replaced by a strong enhancement ($R_{AA}>1$) at FCC at low transverse momentum.

preprint2015arXiv

A Monte Carlo framework for noncontinuous interactions between particles and classical fields

Particles and fields are standard components in numerical simulations like transport simulations in nuclear physics and have very well understood dynamics. Still, a common problem is the interaction between particles and fields due to their different formal description. Particle interactions are discrete, point-like events while fields have purely continuous equations of motion. A workaround is the use of effective theories like the Langevin equation with the drawback of energy conservation violation. We present a new method, which allows to model non-continuous interactions between particles and scalar fields, allowing us to simulate scattering-like interactions which exchange energy and momentum quanta between fields and particles obeying full energy and momentum conservation and control over interaction strengths and times. In this paper we apply this method to different model systems, starting with a simple scalar harmonic oscillator which is damped by losing discrete energy quanta. The second and third system is a scalar oscillator and a one dimensional field which are both damped by discrete energy loss and which are coupled to a stochastic force, leading to equilibrium states which correspond to statistical Langevin-like systems. The last example is a scalar field in 3D which is coupled to a microcanonical ensemble of particles by incorporating particle production and annihilation processes. Obeying the detailed-balance principle, the system equilibrates to thermal and chemical equilibrium with dynamical fluctuations on the fields, generated dynamically by the discrete interactions.

preprint2015arXiv

Constraints on the Jet-Medium Coupling from Measurements at RHIC and LHC

The measured data on the nuclear modification factor for pions and reconstructed jets as well as on the high-pT elliptic flow at RHIC and LHC energies are compared to results from a linear pQCD and a highly non-linear hybrid AdS holographic model of jet-energy loss. We find that the high-pT ellitic flow requires to include realistic medium transverse flow fields and a jet-medium coupling including the effects of the energy of the jet, the temperature of the bulk medium, and non-equilibrium effects close to the phase transition. We extend our jet-energy loss model that is coupled to state-of-the-art hydrodynamic prescriptions to backgrounds generated by the parton cascade BAMPS. We demonstrate that the results for the hydrodynamic and the parton-cascade backgrounds show a remarkable similarity. Unfortunately, the results for both the pion and a parton-jet nuclear modification factor are insensitive to the jet-path dependence of the models considered.

preprint2015arXiv

Glueballs amass at RHIC and LHC Colliders! - The early quarkless 1st order phase transition at $T=270$ MeV - from pure Yang-Mills glue plasma to GlueBall-Hagedorn states

The early stage of high multiplicity pp, pA and AA collider is represented by a nearly quarkless, hot, deconfined pure gluon plasma. According to pure Yang-Mills Lattice Gauge Theory, this hot pure glue matter undergoes, at a high temperature, $T_c = 270$ MeV, a first order phase transition into a confined Hagedorn-GlueBall fluid. These new scenario should be characterized by a suppression of high $p_T$ photons and dileptons, baryon suppression and enhanced strange meson production. We propose to observe this newly predicted class of events at LHC and RHIC.

preprint2015arXiv

Influence of multiple in-medium scattering processes on the momentum imbalance of reconstructed di-jets

Experimental data measured in sqrt(s)=2.76 TeV Pb+Pb collisions at the LHC show a significant enhancement of events with an unbalanced pair of reconstructed jet momenta in comparison with p+p collisions. This enhancement of momentum imbalance is supposed to be caused by the different momentum loss of the initial back-to-back di-partons by scatterings within the created dense medium. For investigating the underlying partonic momentum loss we employ the on-shell transport model BAMPS (Boltzmann Approach for Multi-Parton Scattering) for full heavy-ion collisions, which numerically solves the 3+1D Boltzmann equation based on 2->2 as well as inelastic 2<->3 scattering processes, together with PYTHIA initial state conditions for the parton showers. Due to the employed test-particle approach jet reconstruction within BAMPS events is not trivial. We introduce a method that nevertheless allows the microscopic simulation of the full evolution of the shower particles, recoiled medium particles, and the underlying bulk medium in one common microscopic framework. With this method it is possible to investigate the role of the medium recoil for the momentum imbalance A_J while using well-established background subtraction algorithms. Due to the available particle information in configuration as well as momentum space within BAMPS, it is additionally possible to reproduce the entire evolution of the reconstructed jets within the medium. With this information we investigate the sensitivity of the jet momentum loss from the difference in the partonic in-medium path lengths.

preprint2015arXiv

Kinetic approach to a relativistic Bose-Einstein condensate

We apply a Boltzmann approach to the kinetic regime of a relativistic Bose-Einstein condensate of scalar bosons by decomposing the one-particle distribution function in a condensate part and a non-zero momentum part of excited modes, leading to a coupled set of evolution equations which are then solved efficiently with an adaptive higher order Runge-Kutta scheme. We compare our results to the partonic cascade Monte-Carlo simulation BAMPS for a critical but far from equilibrium case of massless bosons. Motivated by the color glass condensate initial conditions in QCD with a strongly overpopulated initial glasma state, we also discuss the time evolution starting from an overpopulated initial distribution function of massive scalar bosons. In this system a self-similar evolution of the particle cascade with a non-relativistic turbulent scaling in the infrared sector is observed as well as a relativistic exponent for the direct energy cascade, confirming a weak wave turbulence in the ultraviolet region.

preprint2015arXiv

Monte-Carlo approach to particle-field interactions and the kinetics of the chiral phase transition

The kinetics of the chiral phase transition is studied within a linear quark-meson-$σ$ model, using a Monte-Carlo approach to semiclassical particle-field dynamics. The meson fields are described on the mean-field level and quarks and antiquarks as ensembles of test particles. Collisions between quarks and antiquarks as well as the $q\overline{q}$ annihilation to $σ$ mesons and the decay of $σ$ mesons is treated, using the corresponding transition-matrix elements from the underlying quantum field theory, obeying strictly the rule of detailed balance and energy-momentum conservation. The approach allows to study fluctuations without making ad hoc assumptions concerning the statistical nature of the random process as necessary in Langevin-Fokker-Planck frameworks.

preprint2015arXiv

Simulation of stationary Gaussian noise with regard to the Langevin equation with memory effect

We present an efficient method for simulating a stationary Gaussian noise with an arbitrary covariance function and then study numerically the impact of time-correlated noise on the time evolution of a 1 + 1 dimensional generalized Langevin equation by comparing also to analytical results. Finally, we apply our method to the generalized Langevin equation with an external harmonic and double-well potential.

preprint2015arXiv

Thermalization of gluons with Bose-Einstein condensation

We study the thermalization of gluons far from thermal equilibrium in relativistic kinetic theory. The initial distribution of gluons is assumed to resemble that in the early stage of ultrarelativistic heavy ion collisions. Only elastic scatterings in static, nonexpanding gluonic matter are considered. At first we show that the occurrence of condensation in the limit of vanishing particle mass requires a general constraint for the scattering matrix element. Then the thermalization of gluons with Bose-Einstein condensation is demonstrated in a transport calculation. We see a continuously increasing overpopulation of low energy gluons, followed by a decrease to the equilibrium distribution, when the condensation occurs. The times of the completion of the gluon condensation and of the entropy production are calculated. These times scale inversely with the energy density.

preprint2014arXiv

Elastic and radiative heavy quark interactions in ultra-relativistic heavy-ion collisions

Elastic and radiative heavy quark interactions with light partons are studied with the partonic transport description BAMPS (Boltzmann Approach to MultiParton Scatterings). After calculating the cross section of radiative processes for finite masses in the improved Gunion-Bertsch approximation and verifying this calculation by comparing to the exact result, we study elastic and radiative heavy quark energy loss in a static medium of quarks and gluons. Furthermore, the full 3+1D space-time evolution of gluons, light quarks, and heavy quarks in ultra-relativistic heavy-ion collisions at the BNL Relativistic Heavy-Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) are calculated with BAMPS including elastic and radiative heavy flavor interactions. Treating light and heavy particles on the same footing in the same framework, we find that the experimentally measured nuclear modification factor of charged hadrons and D mesons at the LHC can be simultaneously described. In addition, we calculate the heavy flavor evolution with an improved screening procedure from hard-thermal-loop calculations and confront the results with experimental data of the nuclear modification factor and the elliptic flow of heavy flavor particles at RHIC and LHC.

preprint2014arXiv

Electric conductivity of the quark-gluon plasma investigated using a perturbative QCD based parton cascade

Electric conductivity is sensitive to effective cross sections among the particles of the partonic medium. We investigate the electric conductivity of a hot plasma of quarks and gluons, solving the relativistic Boltzmann equation. In order to extract this transport coefficient, we employ the Green-Kubo formalism and, independently, a method motivated by the classical definition of electric conductivity. To this end we evaluate the static electric diffusion current upon the influence of an electric field. Both methods give identical results. For the first time, we obtain numerically the Drude electric conductivity formula for an ultrarelativistic gas of quarks and gluons employing constant isotropic binary cross sections. Furthermore, we extract the electric conductivity for a system of massless quarks and gluons including screened binary and inelastic, radiative $2\leftrightarrow 3$ perturbative QCD scattering. Comparing with recent lattice results, we find an agreement in the temperature dependence of the conductivity.

preprint2014arXiv

Elliptic flow and nuclear modification factor in ultra-relativistic heavy-ion collisions within a partonic transport model

The quark gluon plasma produced in ultra-relativistic heavy-ion collisions exhibits remarkable features. It behaves like a nearly perfect liquid with a small shear viscosity to entropy density ratio and leads to the quenching of highly energetic particles. We show that both effects can be understood for the first time within one common framework. Employing the parton cascade Boltzmann Approach to Multi-Parton Scatterings (BAMPS), the microscopic interactions and the space-time evolution of the quark gluon plasma are calculated by solving the relativistic Boltzmann equation. Based on cross sections obtained from perturbative QCD with explicitly taking the running coupling into account, we calculate the nuclear modification factor and elliptic flow in ultra-relativistic heavy-ion collisions. With only one single parameter associated with coherence effects of medium-induced gluon radiation, the experimental data of both observables can be understood on a microscopic level. Furthermore, we show that perturbative QCD interactions with a running coupling lead to a sufficiently small shear viscosity to entropy density ratio of the quark gluon plasma, which provides a microscopic explanation for the observations stated by hydrodynamic calculations.

preprint2014arXiv

Elliptic Flow Suppression due to Hadron Mass Spectrum

Hadron resonance gas models provide a good description of the equation of state of quantum chromodynamics determined by lattice QCD calculations at temperatures $T \sim 100-155$ MeV. In this paper we investigate the effects of an exponentially increasing hadron mass spectrum (Hagedorn spectrum) on the azimuthal anisotropy of the rapidly expanding matter formed in ultrarelativistic heavy ion collisions. If the temperature at which the conversion from fluid degrees of freedom to hadrons is sufficiently close to the Hagedorn temperature, the production of Hagedorn resonances suppresses the differential elliptic flow of all hadron species.

preprint2014arXiv

Hadronic vs. partonic J/psi production in the Statistical Hadronization Model

Assuming the Statistical Hadronization Model and charm flavour conservation, one can quantify the deviation of charm and anti-charm quarks from chemical equilibrium both in the hadronic phase and in the quark-gluon plasma. By linking the hadronic fireball volume and the corresponding plasma source volume via entropy conservation as proposed by Grandchamp et al., a direct comparison between the $J/ψ$ yields in both pictures can be used to obtain limits for the charm quark mass in the medium.

preprint2014arXiv

Momentum imbalance of D mesons in ultra-relativistic heavy-ion collisions at LHC

As a new observable for heavy flavor correlations the momentum imbalance $A_D$ of D mesons is proposed. It is defined analogously to the jet momentum imbalance $A_J$ of fully reconstructed jets. However, since D mesons are flavor tagged particles, no jet reconstruction is necessary. $A_D$ quantifies the influence of the medium created in heavy-ion collisions on correlated charm pairs. We present results with the partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS), which describes well the nuclear modification factor and elliptic flow of all heavy flavor particles at RHIC and LHC. The $A_D$ distribution in heavy-ion collisions at LHC is shifted to larger values of $A_D$ compared to proton-proton collisions. We argue that this shift is due to medium effects and can be explained partially by a path length imbalance of charm pairs and partially by momentum fluctuations in the initial charm pair distribution.

preprint2014arXiv

Shear viscosity of an ultrarelativistic Boltzmann gas with isotropic inelastic scattering processes

We derive an analytic expression for the shear viscosity of an ultra-relativistic gas in presence of both elastic $2\to 2$ and inelastic $2\leftrightarrow 3$ processes with isotropic differential cross sections. The derivation is based on the entropy principle and Grad's approximation for the off-equilibrium distribution function. The obtained formula relates the shear viscosity coefficient $η$ to the total cross sections $σ_{22}$ and $σ_{23}$ of the elastic resp. inelastic processes. The values of shear viscosity extracted using the Green-Kubo formula from kinetic transport calculations are shown to be in excellent agreement with the analytic results which demonstrates the validity of the derived formula.

preprint2014arXiv

Spin Asymmetry of J/psi in Peripheral Pb+Pb Collisions at LHC

By generalizing the statistical model for particle production to the spin degree of freedom of initially produced J/psi, we study the spin projection J_y of J/psi perpendicular to the reaction plane in peripheral heavy ion collisions at the LHC energy that leads to a strong, albeit of short duration, magnetic field. We find that for J/psis produced directly from charm and anticharm quarks in the color singlet state, like that in the Color-Singlet Model, their yield in the presence of the magnetic field is larger for J_y=0 than for J_y=1 or -1. This leads to a spin asymmetry of finally produced J/ψeven after including their final-state scattering in the produced quark-gluon plasma.

preprint2014arXiv

Suppression of the LHC $p/π$ ratio due to the QCD mass spectrum

Recent measurements of the proton to pion ratio at $\sqrt{s}_{NN}=2.7$ TeV by the ALICE collaboration at the LHC have found it to be lower than predictions from thermal fits. In this paper we investigate the role that the extended mass spectrum via Hagedorn states- massive resonances that follow an exponential mass spectrum and possess large decay widths- play in the determination of particle ratios at LHC through a scenario of multi-particle reactions and dynamical chemical equilibrium within the hadron gas phase. We show that it is possible to describe the lower $p/π$ ratio at LHC while still obtaining the experimental ratio of $K/π$ and $Λ/π^+$ in the Hagedorn state scenario if the protons are underpopulated at the switching temperature from hydrodynamics to the hadron gas phase.

preprint2014arXiv

Understanding the $p/π$ ratio at LHC due to QCD mass spectrum

Thermal fits have consistently reproduced the experimental particles yields of heavy ion collisions, however, the proton to pion ratio from ALICE Pb+Pb $\sqrt{s_{NN}}=2.76$ TeV is over-predicted by thermal models- known at the $p/π$ puzzle. Here we test the relevance of the extended mass spectrum, i.e., include Hagedorn states (resonances that follow an exponential mass spectrum and have very short life times) on the $p/π$ puzzle. We find that the extended mass spectrum is able to reproduce particle ratios at both RHIC and the LHC as well as being able to match the lower $p/π$ ratio at the LHC through dynamical chemical equilibration.

preprint2013arXiv

Asymptotic description of finite lifetime effects on the photon emission from a quark-gluon plasma

Direct photons play an important role as electromagnetic probes from the quark-gluon plasma (QGP) which occurs during ultrarelativistic heavy-ion collisions. In this context, it is of particular interest how the finite lifetime of the QGP affects the resulting photon production. Earlier investigations on this question were accompanied by a divergent contribution from the vacuum polarization and by the remaining contributions not being integrable in the ultraviolet (UV) domain. In this work, we provide a different approach in which we do not consider the photon number density at finite times, but for free asymptotic states obtained by switching the electromagnetic interaction according to the Gell-Mann and Low theorem. This procedure eliminates a possible unphysical contribution from the vacuum polarization and, moreover, renders the photon number density UV integrable. It is emphasized that the consideration of free asymptotic states is, indeed, crucial to obtain such physically reasonable results.

preprint2013arXiv

B_c meson enhancement and the momentum dependence in Pb+Pb collisions at LHC energy

B_c meson production in Pb+Pb collisions at sqrt{s}=2.76 A TeV is surveyed in both a statistical coalescence model and a transport model. The nuclear modification factor R_{AA} is predicted to be between 2 and 18 in the most central collisions, which can help to confirm the regeneration mechanism. In addition, the momentum dependence is also investigated as given by the transport model. A strong suppression of the transverse momentum is found in central collisions accompanying the enhancement in yield. The spectrum and elliptic flow of B_c are also discussed.

preprint2013arXiv

Collisional energy loss of heavy quarks

We develop a transport approach for heavy quarks in a quark-gluon plasma, which is based on improved binary collision rates taking into account quantum statistics, the running of the QCD coupling and an effective screening mass adjusted to hard-thermal loop calculations. We quantify the effects of in-medium collisions by calculating the heavy flavor nuclear modification factor and the elliptic flow for RHIC energies, which are comparable to radiative effects. We also derive an analytic formula for the mean collisional energy loss of an energetic heavy quark in a streaming quark gluon plasma.

preprint2013arXiv

Heavy vs. light flavor energy loss within a partonic transport model

The full space-time evolution of gluons, light and heavy quarks in ultra-relativistic heavy-ion collisions is studied within the partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS). We discuss for light and heavy quarks the elastic and radiative energy loss with a running coupling. Radiative processes, in particular, are implemented through an improved version of the Gunion-Bertsch matrix element, which is derived from comparisons to the exact result, explicitly taking finite heavy quark masses into account. Subsequently, we present results with and without radiative processes and compare them to experimental data at LHC. A focus is put on the nuclear modification factor of charged hadrons and D mesons.

preprint2013arXiv

Non-equilibrium photon production arising from the chiral mass shift

We investigate the photon emission arising from the shift of the quark/antiquark masses during the chiral phase transition in the early stage of ultrarelativistic heavy-ion collisions. As this mass shift leads to spontaneous creation of quark-antiquark pairs and thus contributes to the formation of the quark-gluon plasma, our investigations are relevant in the context of finite-lifetime effects on the photon emission from the latter. Earlier investigations on this topic by Boyanovsky et al. were plagued by a UV divergent contribution from the vacuum polarization and the non-integrability of the remaining contributions in the ultraviolet domain. In contrast to these investigations, we do not consider the photon number-density at finite times but for free asymptotic states obtained by an adiabatic switching of the electromagnetic interaction according to the Gell-Mann and Low theorem. This approach eliminates possible unphysical contributions from the vacuum polarization and renders the resulting photon spectra integrable in the ultraviolet domain. It is emphasized that the consideration of free asymptotic states is indeed crucial to obtain such physically reasonable results.

preprint2013arXiv

Nonequilibrium dynamics and transport near the chiral phase transition of a quark-meson model

Based on the 2PI quantum effective action of the linear sigma model with constituent quarks, we develop a transport approach to study systems out of equilibrium. In particular, we focus on the chiral phase transition as well as the critical point, where nonequilibrium effects near the phase transition give rise to critical behavior such as the fluctuation of the baryon number density. Predictions for long-range correlations and fluctuations of observables in our model could be used to study fundamental properties of the QCD phase transition. In the previous version of our transport model the chiral fields are implemented as mean fields, whereas quarks are treated as on-shell particles in the Vlasov equation with a dynamical force term. The current update includes also the distribution functions of sigma mesons and pions in a self-consistent way. On this basis a dissipation kernel between the mean fields and particle modes can be implemented.

preprint2013arXiv

Radiative parton processes in perturbative QCD - an improved version of the Gunion and Bertsch cross section from comparisons to the exact result

In this work we compare the Gunion-Bertsch approximation of the leading order perturbative QCD radiation matrix element to the exact result. To this end, we revisit the derivation of the Gunion-Bertsch approximation as well as perform extensive numerical comparisons of the Gunion-Bertsch and the exact result. We find that when employing the matrix elements to obtain rates or cross sections from phase space integration, the amplitude by Gunion and Bertsch deviates from the correct result in characteristic regions of the phase space. We propose an improved version of the Gunion-Bertsch matrix element which agrees very well with the exact result in all phase space regions.

preprint2012arXiv

Dissipative hydrodynamics for multi-component systems

Second-order dissipative hydrodynamic equations for each component of a multi-component system are derived using the entropy principle. Comparison of the solutions with kinetic transport results demonstrates validity of the obtained equations. We demonstrate how the shear viscosity of the total system can be calculated in terms of the involved cross sections and partial densities. Presence of the inter-species interactions leads to a characteristic time-dependence of the shear viscosity of the mixture, which also means that the shear viscosity of a mixture cannot be calculated using the Green-Kubo formalism the way it has been done recently. This finding is of interest for understanding of the shear viscosity of a quark-gluon-plasme extracted from comparisons of hydrodynamic simulations with experimental results from RHIC and LHC.

preprint2012arXiv

Hadron Mass Spectrum and the Shear Viscosity to Entropy Density Ratio of Hot Hadronic Matter

Lattice calculations of the QCD trace anomaly at temperatures $T<160$ MeV have been shown to match hadron resonance gas model calculations, which include an exponentially rising hadron mass spectrum. In this paper we perform a more detailed comparison of the model calculations to lattice data that confirms the need for an exponentially increasing density of hadronic states. Also, we find that the lattice data is compatible with a hadron density of states that goes as $ρ(m) \sim m^{-a}\exp(m/T_H) $ at large $m$ with $a> 5/2$ (where $T_H \sim 167$ MeV). With this specific subleading contribution to the density of states, heavy resonances are most likely undergo 2-body decay (instead of multi-particle decay), which facilitates their inclusion into hadron transport codes. Moreover, estimates for the shear viscosity and the shear relaxation time coefficient of the hadron resonance model computed within the excluded volume approximation suggest that these transport coefficients are sensitive to the parameters that define the hadron mass spectrum.

preprint2012arXiv

Off-equilibrium photon production during the chiral phase transition

In the early stage of ultrarelativistic heavy-ion collisions chiral symmetry is restored temporarily. During this so-called chiral phase transition, the quark masses change from their constituent to their bare values. This mass shift leads to the spontaneous non-perturbative creation of quark-antiquark pairs, which effectively contributes to the formation of the quark-gluon plasma. We investigate the photon production induced by this creation process. We provide an approach that eliminates possible unphysical contributions from the vacuum polarization and renders the resulting photon spectra integrable inthe ultraviolet domain. The off-equilibrium photon numbers are of quadratic order in the perturbative coupling constants while a thermal production is only of quartic order. Quantitatively, we find, however, that for the most physical mass-shift scenarios and for photon momenta larger than 1 GeV the off-equilibrium processes contribute less photons than the thermal processes.

preprint2012arXiv

Open heavy flavor and J/psi at RHIC and LHC

The suppression and collective flow of open and hidden heavy flavor in ultra-relativistic heavy ion collisions is studied with the partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS). Charm and bottom quarks interact with the rest of the medium via elastic scatterings with an improved Debye screening and a running coupling. Heavy flavor electron data from RHIC can only be described if the cross section is multiplied with K=3.5 to mimic the contribution from radiative processes. With the same value for K we compare to elliptic flow data of D mesons at the LHC and find a good agreement. However, the nuclear modification factor of D mesons, non-prompt J/psi, muons, and electrons at the LHC is underestimated in our calculations. In addition, we present first BAMPS results on the suppression and elliptic flow of J/psi mesons, which are treated in the same framework.

preprint2012arXiv

Open heavy flavor in Pb+Pb collisions at sqrt(s)=2.76 TeV within a transport model

The space-time evolution of open heavy flavor is studied in Pb+Pb collisions at sqrt(s)=2.76 TeV using the partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS). An updated version of BAMPS is presented which allows interactions among all partons: gluons, light quarks, and heavy quarks. Heavy quarks, in particular, interact with the rest of the medium via binary scatterings with a running coupling and a Debye screening which is matched by comparing to hard thermal loop calculations. The lack of radiative processes in the heavy flavor sector is accounted for by scaling the binary cross section with a phenomenological factor K=3.5, which describes well the elliptic flow v_2 and nuclear modification factor R_AA at RHIC. Within this framework we calculate in a comprehensive study the v_2 and R_AA of all interesting open heavy flavor particles at the LHC: electrons, muons, D mesons, and non-prompt J/psi from B mesons. We compare to experimental data, where it is already available, or make predictions. To do this accurately next-to-leading order initial heavy quark distributions are employed which agree well with proton-proton data of heavy flavor at sqrt(s)=7 TeV.

preprint2012arXiv

RHIC and LHC phenomena with an unified parton transport

We discuss recent applications of the partonic pQCD based cascade model BAMPS with focus on heavy-ion phenomeneology in hard and soft momentum range. The nuclear modification factor as well as elliptic flow are calculated in BAMPS for RHIC end LHC energies. These observables are also discussed within the same framework for charm and bottom quarks. Contributing to the recent jet-quenching investigations we present first preliminary results on application of jet reconstruction algorithms in BAMPS. Finally, collective effects induced by jets are investigated: we demonstrate the development of Mach cones in ideal matter as well in the highly viscous regime.

preprint2012arXiv

Soft gluon emission off a heavy quark revisited

An improved generalized suppression factor for gluon emission off a heavy quark is derived within perturbative QCD, which is valid for the full range of rapidity of the radiated gluon and also has no restriction on the scaled mass of the quark with its energy. In the appropriate limit it correctly reproduces the usual dead cone factor in the forward rapidity region. On the other hand, this improved suppression factor becomes close to unity in the backward direction. This indicates a small suppression of gluon emission in the backward region, which should have an impact on the phenomenology of heavy quark energy loss in the hot and dense matter produced in ultra-relativistic heavy-ion collisions.

preprint2012arXiv

Transition from ideal to viscous Mach cones in BAMPS

We investigate in a microscopical transport model the evolution of conical structures originating from the supersonic projectile moving through the 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. Furthermore, the two-particle correlations for different viscosities are extracted from the numerical calculations and we compare them to an analytical approximation. In addition, by adjusting he cross section we investigate the influence of the viscosity to the structure of Mach cones.

preprint2011arXiv

Dissipative hydrodynamics for relativistic multi-component systems

Second-order dissipative hydrodynamic equations for each component of a multi-component system are derived using the entropy principle. The shear viscosity of the whole system, appearing in the equation summed-up over all components, is related to the partial shear pressures and cannot be considered as an external parameter. We demonstrate that it is essential to solve hydrodynamic equations for each component, instead of treating a mixture as an effective one-component system with a free parameter $η/s$. Thus, extractions of the $η/s$ value of the QGP at RHIC and LHC have to be reexamined.

preprint2011arXiv

Elliptic and Triangular Flow and their Correlation in Ultrarelativistic High Multiplicity Proton Proton Collisions at 14 TeV

The spatial configuration of initial partons in high-multiplicity proton-proton scatterings at 14 TeV is assumed as three randomly positioned "hot spots". The parton momentum distribution in the hot spots is calculated by HIJING2.0 with some modifications. This initial condition causes not only large eccentricity $ε_2$ but also triangularity $ε_3$ and the correlation of $ε_2-ε_3$ event-plane angles. The final elliptic flow $v_2$, triangular flow $v_3$, and the correlation of $v_2-v_3$ event-plane angles are calculated by using the parton cascade model BAMPS to simulate the space-time parton evolution. Our results show that the $v_2-v_3$ correlation is different from that of $ε_2-ε_3$. This finding indicates that translations of different Fourier components of the initial spatial asymmetry to the final flow components are not independent. A dynamical correlation between the elliptic and triangular flow appears during the collective expansion.

preprint2011arXiv

Elliptic flow and energy loss of heavy quarks in ultra-relativistic heavy ion collisions

The space-time propagation of heavy quarks in ultra-relativistic heavy ion collisions is studied within the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). In this model heavy quarks interact with the partonic medium via binary scatterings. The cross sections for these interactions are calculated with leading order perturbative QCD, but feature a more precise Debye screening derived within the hard thermal loop approximation and obey the running of the coupling. Within this framework the elliptic flow and the nuclear modification factor of heavy quarks are computed for RHIC and LHC energies and compared to available experimental data. It is found that binary scatterings alone cannot reproduce the data and, therefore, radiative corrections have to be taken into account.

preprint2011arXiv

Extraction of shear viscosity in stationary states of relativistic particle systems

Starting from a classical picture of shear viscosity we construct a stationary velocity gradient in a microscopic parton cascade. Employing the Navier-Stokes ansatz we extract the shear viscosity coefficient $η$. For elastic isotropic scatterings we find an excellent agreement with the analytic values. This confirms the applicability of this method. Furthermore for both elastic and inelastic scatterings with pQCD based cross sections we extract the shear viscosity coefficient $η$ for a pure gluonic system and find a good agreement with already published calculations.

preprint2011arXiv

Hagedorn States and Thermalization in Heavy Ion Collisions

In recent years, Hagedorn states have been used to explain the equilibrium and transport properties of a hadron gas close to the QCD critical temperature. These massive resonances are shown to lower $η/s$ to near the AdS/CFT limit close to the phase transition. A comparison of the Hagedorn model to recent lattice results is made and it is found that the hadrons can reach chemical equilibrium almost immediately, well before the chemical freeze-out temperatures found in thermal fits for a hadron gas without Hagedorn states.

preprint2011arXiv

Heavy quarks and charmonium at RHIC and LHC within a partonic transport model

Heavy quark and charmonium production as well as their space-time evolution are studied in transport simulations of heavy-ion collisions at RHIC and LHC. In the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS) heavy quarks can be produced in initial hard parton scatterings or during the evolution of the quark-gluon plasma. Subsequently, they interact with the medium via binary scatterings with a running coupling and a more precise Debye screening which is derived from hard thermal loop calculations, participate in the flow and lose energy. We present results of the elliptic flow and nuclear modification factor of heavy quarks and compare them to available data. Furthermore, preliminary results on J/psi suppression at forward and mid-rapidity are reported for central and non-central collisions at RHIC. For this, we study cold nuclear matter effects and the dissociation as well as regeneration of J/psi in the quark-gluon plasma.

preprint2011arXiv

Jet quenching and elliptic flow at RHIC and LHC within a pQCD-based partonic transport model

Fully dynamic simulations of heavy ion collisions at RHIC and at LHC energies within the perturbative QCD-based partonic transport model BAMPS (Boltzmann Approach to Multi-Parton Scatterings) are presented, focusing on the simultaneous investigation of jet quenching and elliptic flow. The model features inelastic 2 <-> 3 processes based on the Gunion-Bertsch matrix element and has recently been extended to include light quark degrees of freedom, allowing for direct comparison to hadronic data on the nuclear modification factor via a fragmentation scheme for high-pT partons. The nuclear modification factor of neutral pions in central Au+Au collisions at RHIC energy is compared to experimental data. Furthermore first results on the nuclear modification factor and the integrated elliptic flow of charged hadrons in Pb+Pb collisions at LHC are presented and compared to recent ALICE data. These investigations are complemented by a study on the suppression of D-mesons at LHC based on elastic interactions with the medium.

preprint2011arXiv

Longitudinal correlation of the triangular flow event plane in a hybrid approach with hadron and parton cascade initial conditions

The longitudinal long-range correlations of the triangular flow event plane angles are calculated in a Boltzmann + hydrodynamics hybrid approach. The potential to disentangle different energy deposition scenarios is explored by utilizing two different transport approaches for the early non-equilibrium evolution. In the hadronic transport approach the particle production in high energy heavy ion reactions is mainly governed by string excitation and fragmentation processes which are absent in the parton cascade approach. We find that in both approaches the initial state shows a strong longitudinal correlation of the event plane angles which is diluted but still persists in the final state momentum space distributions of the produced particles. A ridge-like structure can also be caused by near-collinear gluon radiation in a parton cascade approach and does not necessarily prove longitudinal flux tubes in the initial state.

preprint2011arXiv

Open heavy flavor at RHIC and LHC in a partonic transport model

Heavy quarks are a unique probe to study the medium produced in ultra-relativistic heavy ion collisions. Within the partonic transport model Boltzmann approach to multi-parton scatterings (BAMPS) the production, energy loss and elliptic flow of heavy quarks are investigated in 3+1 dimensional simulations of the quark gluon plasma evolution. With only binary interactions between heavy quarks and particles from the medium calculated within perturbative QCD, the results on elliptic flow and the nuclear modification factor are not compatible with experimental data from RHIC and LHC. However, if the binary cross section is multiplied with K=4 both the elliptic flow and the nuclear modification factor are simultaneously described at RHIC and also LHC. Furthermore, preliminary results are presented that the implementation of radiative processes leads to a stronger suppression which agrees well with the measured nuclear modification factor at RHIC without the need of any K factor.

preprint2010arXiv

A relativistic dissipative hydrodynamic description for systems including particle number changing processes

Relativistic dissipative hydrodynamic equations are extended by taking into account particle number changing processes in a gluon system, which expands in one dimension boost-invariantly. Chemical equilibration is treated by a rate equation for the particle number density based on Boltzmann equation and Grad's ansatz for the off-equilibrium particle phase space distribution. We find that not only the particle production, but also the temperature and the momentum spectra of the gluon system, obtained from the hydrodynamic calculations, are sensitive to the rates of particle number changing processes. Comparisons of the hydrodynamic calculations with the transport ones employing the parton cascade BAMPS show the inaccuracy of the rate equation at large shear viscosity to entropy density ratio. To improve the rate equation, the Grad's ansatz has to be modified beyond the second moments in momentum.

preprint2010arXiv

Bulk matter evolution and extraction of jet transport parameter in heavy-ion collisions at RHIC

Within the picture of jet quenching induced by multiple parton scattering and gluon bremsstrahlung, medium modification of parton fragmentation functions and therefore the suppression of large transverse momentum hadron spectra are controlled by both the value and the space-time profile of the jet transport parameter along the jet propagation path. Experimental data on single hadron suppression in high-energy heavy-ion collisions at the RHIC energy are analyzed within the higher-twist (HT) approach to the medium modified fragmentation functions and the next-to-leading order (NLO) perturbative QCD (pQCD) parton model. Assuming that the jet transport parameter $\hat q$ is proportional to the particle number density in both QGP and hadronic phase, experimental data on jet quenching in deeply inelastic scattering (DIS) off nuclear targets can provide guidance on $\hat q_{h}$ in the hot hadronic matter. One can then study the dependence of extracted initial value of jet quenching parameter $\hat q_{0}$ at initial time $τ_{0}$ on the bulk medium evolution. Effects of transverse expansion, radial flow, phase transition and non-equilibrium evolution are examined. The extracted values are found to vary from $\hat q_{0}τ_{0}= 0.54$ GeV$^{2}$ in the (1+3)d ideal hydrodynamic model to 0.96 GeV$^{2}$ in a cascade model, with the main differences coming from the initial non-equilibrium evolution and the later hadronic evolution. The overall contribution to jet quenching from the hadronic phase, about 22-44%, is found to be significant. Therefore, realistic description of the early non-equilibrium parton evolution and later hadronic interaction will be critical for accurate extraction of the jet transport parameter in the strongly interacting QGP phase in high-energy heavy-ion collisions.

preprint2010arXiv

Dependence of elliptic flow on number of parton degrees of freedom

We calculate the elliptic flow parameter $v_2$ for Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV employing the parton cascade BAMPS (Boltzmann Approach of Multiparton Scatterings). Besides gluon interactions including the bremsstrahlung process, interactions with quarks are considered in an effective, but approximate way to investigate the dependence of the collective flow on the number of parton degrees of freedom. We find that $v_2$ as a function of the transverse momentum $p_T$ is sensitive to the number of parton degrees of freedom, whereas the $p_T$ averaged $v_2$ does not. When including quarks, $v_2(p_T)$ shifts to lower $p_T$, the parton transverse momentum spectra become softer and the mean parton transverse momenta decrease.

preprint2010arXiv

Energy loss in a partonic transport model including bremsstrahlung processes

A detailed investigation of the energy loss of gluons that traverse a thermal gluonic medium simulated within the perturbative QCD--based transport model BAMPS (a Boltzmann approach to multiparton scatterings) is presented in the first part of this work. For simplicity the medium response is neglected in these calculations. The energy loss from purely elastic interactions is compared to the case where radiative processes are consistently included based on the matrix element by Gunion and Bertsch. From this comparison gluon multiplication processes gg -> ggg are found to be the dominant source of energy loss within the approach employed here. The consequences for the quenching of gluons with high transverse momentum in fully dynamic simulations of Au+Au collisions at the RHIC energy of sqrt(s) = 200 AGeV are discussed in the second major part of this work. The results for central collisions as discussed in a previous publication are revisited and first results on the nuclear modification factor R_AA for non-central Au+Au collisions are presented. They show a decreased quenching compared to central collisions while retaining the same shape. The investigation of the elliptic flow v2 is extended up to non-thermal transverse momenta of 10 GeV, exhibiting a maximum v2 at roughly 4 to 5 GeV and a subsequent decrease. Finally the sensitivity of the aforementioned results on the specific implementation of the effective modeling of the Landau-Pomeranchuk-Migdal (LPM) effect via a formation time based cut-off is explored.

preprint2010arXiv

Hagedorn States and Thermalization

In recent years Hagedorn states have been used to explain the physics close to the critical temperature within a hadron gas. Because of their large decay widths these massive resonances lower $η/s$ to near the AdS/CFT limit within the hadron gas phase. A comparison of the Hagedorn model to recent lattice results is made and it is found that for both Tc =176 MeV and Tc=196 MeV, the hadrons can reach chemical equilibrium almost immediately, well before the chemical freeze-out temperatures found in thermal fits for a hadron gas without Hagedorn states. In this paper we also observe the effects of Hagedorn States on the $K^+/π^+$ horn seen at AGS, SPS, and RHIC.

preprint2010arXiv

Heavy quark production at RHIC and LHC within a partonic transport model

The production and space-time evolution of charm and bottom quarks in nucleus-nucleus collisions at RHIC and LHC are investigated with the partonic transport model BAMPS (Boltzmann Approach of MultiParton Scatterings). Heavy quarks, produced in primary hard parton scatterings during nucleon-nucleon collisions, are sampled using the Monte Carlo event generator PYTHIA or the leading order mini-jet model in conjunction with the Glauber model, revealing a strong sensitivity on the parton distribution functions, scales, and heavy quark mass. In a comprehensive study exploring different charm masses, K factors, and possible initial gluon conditions, secondary production and the evolution of heavy quarks are examined within a fully dynamic BAMPS simulation for central heavy ion collisions at RHIC and LHC. Although charm production in the quark-gluon plasma can be neglected at RHIC, it is significant at LHC but very sensitive to the initial conditions and the charm mass. Bottom production in the quark-gluon plasma, however, is negligible both at RHIC and LHC.

preprint2010arXiv

Heavy quarks at RHIC and LHC within a partonic transport model

Production and space-time evolution of heavy quarks in central and non-central heavy-ion collisions at RHIC and LHC are studied with the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). In addition to the initially created heavy quarks in hard parton scatterings during nucleon-nucleon collisions, secondary heavy quark production in the quark-gluon plasma is investigated and the sensitivity on various parameters is estimated. In BAMPS heavy quarks scatter with particles of the medium via elastic collisions, whose cross section is calculated with the running coupling and a more precise implementation of Debye screening. In this framework, we compute the elliptic flow and nuclear modification factor of heavy quarks and compare it to the experimental data.

preprint2010arXiv

Investigation of jet quenching and elliptic flow within a pQCD-based partonic transport model

The partonic transport model BAMPS (a Boltzmann approach to multiparton scatterings) is employed to investigate different aspects of heavy ion collisions within a common framework based on perturbative QCD. This report focuses on the joint investigation of the collective behavior of the created medium and the energy loss of high-pT gluons traversing this medium. To this end the elliptic flow and the nuclear modification factor of gluons in heavy ion collisions at 200 AGeV are simulated with BAMPS. The mechanism for the energy loss of high energy gluons within BAMPS is studied in detail. For this, purely elastic interactions are compared to radiative processes, gg -> ggg, that are implemented based on the matrix element by Gunion and Bertsch. The latter are found to be the dominant source of energy loss within the framework employed in this work.

preprint2010arXiv

Jets and flow within a pQCD-based partonic transport model

We present fully dynamic simulations of heavy ion collisions at RHIC energies within the perturbative QCD-based partonic transport model BAMPS, focusing on the simultaneous investigation of jet-quenching and elliptic flow. The model consistently features elastic and inelastic 2 -> 3 processes, the latter being based on the Gunion-Bertsch matrix element. We discuss first attempts to extend the model to include light quark degrees of freedom and study the energy loss of high energy gluons and quarks in a static partonic medium. The difference between gluons and quarks in inelastic processes is found to be weaker than expected from color factors, due to a self-quenching effect associated with a cut-off modeling the LPM effect.

preprint2010arXiv

Low Mass Dilepton Rate from the Deconfined Phase

We discuss low mass dilepton rates ($\le 1$ GeV) from the deconfined phase of QCD using both perturbative and non-perturbative models and compare with those from lattice gauge theory and in-medium hadron gas. Our analysis suggests that the rate at very low invariant mass ($ M\le 200$ MeV) using the nonperturbative gluon condensate in a semiempirical way within the Green function dominates over the Born-rate and independent of any uncertainty associated with the choice of the strong coupling in perturbation theory. On the other hand the rate from $ρ-q$ interaction in the deconfined phase is important between 200 MeV $\le M \le 1 GeV as it is almost of same order of the Born-rate as well as in-medium hadron gas rate. Also the higher order perturbative rate, leaving aside its various uncertainties, from HTL approximation becomes reliable at $M\ge 200$ MeV and also becomes comparable with the Born-rate and the lattice-rate for $M\ge 500$ MeV, constraining on the broad resonance structures in the dilepton rate at large invariant mass. We also discuss the lattice constraints on the low mass dilepton rate. Furthermore, we discuss a more realistic way to advocate the quark-hadron duality hypothesis based on the dilepton rates from QGP and hadron gas than it is done in the literature.

preprint2010arXiv

Production and elliptic flow of heavy quarks at RHIC and LHC within a partonic transport model

Production and elliptic flow of heavy quarks are investigated in nucleus-nucleus collisions at RHIC and LHC within the Boltzmann Approach of MultiParton Scatterings (BAMPS). The initial heavy quark yield is estimated with the leading order mini-jet model and compared to PYTHIA for several parton distribution functions. Secondary production of heavy quarks in the quark-gluon plasma is examined within the full 3+1 dimensional BAMPS simulation of heavy ion collisions. At RHIC this yield is negligible, but for LHC charm production in the QGP plays a significant role. In addition, we study the elliptic flow of charm quarks at RHIC with the result that leading order processes are not sufficient to describe the experimental data.

preprint2010arXiv

Production, elliptic flow and energy loss of heavy quarks in the quark-gluon plasma

Production, elliptic flow and the nuclear modification factor of charm and bottom quarks are studied in central and non-central heavy-ion collisions at RHIC and LHC using the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). Employing an initial heavy quark yield obtained with PYTHIA the full space-time evolution of charm and bottom quarks in the quark-gluon plasma (QGP) is carried out with BAMPS, taking also secondary production in the QGP into account. Only elastic collisions of heavy quarks with particles from the medium cannot describe the experimentally observed elliptic flow and nuclear modification factor. However, using an improved Debye screening and the running coupling yields a result which is much closer to data.

preprint2010arXiv

Relativistic Shock Waves and Mach Cones in Viscous Gluon Matter

To investigate the formation and the propagation of relativistic shock waves in viscous gluon matter we solve the relativistic Riemann problem using a microscopic parton cascade. We demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to entropy density ratio $η/s$. Furthermore we compare our results with those obtained by solving the relativistic causal dissipative fluid equations of Israel and Stewart (IS), in order to show the validity of the IS hydrodynamics. Employing the parton cascade we also investigate the formation of Mach shocks induced by a high-energy gluon traversing viscous gluon matter. For $η/s = 0.08$ a Mach cone structure is observed, whereas the signal smears out for $η/s \geq 0.32$.

preprint2010arXiv

The quark-gluon-plasma phase transition diagram, Hagedorn matter and quark-gluon liquid

In order to study the nuclear matter in the relativistic heavy ion collisions and the compact stars, we need the hadronic density of states for the entire ($μ_B-T$) phase transition diagram. We present a model for the continuous high-lying mass (and volume) spectrum density of states that fits the Hagedorn mass spectrum. This model explains the origin of the tri-critical point besides various phenomena such as the quarkyonic matter and the quark-gluon liquid. The Hagedorn mass spectrum is derived for the color-singlet quark-gluon bag with various internal structures such as the unimodular unitary, orthogonal and color-flavor locked symplectic symmetry groups. The continuous high-lying hadronic mass spectrum is populated at first by the unitary Hagedorn states. Then the spectrum turns to be dominated by the colorless orthogonal states as the dilute system is heated up. Subsequently, the liquid/gas of orthogonal Hagedorn states undergoes higher order deconfinement phase transition to quark-gluon plasma. Under the deconfinement phase transition process, the color-singlet states is broken badly to form the colored $SU(N_c)$ symmetry group. On the other hand, when the hadronic matter is compressed to larger $μ_{B}$ and heated up, the colorless unitary states undergoes first order phase transition to explosive quark-gluon plasma. The tri-critical point emerges as a change in the characteristic behaviour of the matter and as an intersection among various phases with different internal symmetries. When the saturated hadronic matter is cooled down and compressed to higher density, it turns to be dominated by the colorless symplectic states. This matter exhibits the first order phase transition to quark-gluon plasma when it is heated up to higher temperature. The role of chiral phase transition is also discussed.

preprint2010arXiv

Thermalization through Hagedorn states - the importance of multiparticle collisions

Quick chemical equilibration times of hadrons within a hadron gas are explained dynamically using Hagedorn states, which drive particles into equilibrium close to the critical temperature. Within this scheme master equations are employed for the chemical equilibration of various hadronic particles like (strange) baryon and antibaryons. A comparison of the Hagedorn model to recent lattice results is made and it is found that for both Tc =176 MeV and Tc=196 MeV, the hadrons can reach chemical equilibrium almost immediately, well before the chemical freeze-out temperatures found in thermal fits for a hadron gas without Hagedorn states.

preprint2010arXiv

Third-order dissipative hydrodynamics from the entropy principle

We review the entropy based derivation of third-order hydrodynamic equations and compare their solutions in one-dimensional boost-invariant geometry with calculations by the partonic cascade BAMPS. We demonstrate that Grad's approximation, which underlies the derivation of both Israel-Stewart and third-order equations, describes the transverse spectra from BAMPS with high accuracy. At the same time solutions of third-order equations are much closer to BAMPS results than solutions of Israel-Stewart equations. Introducing a resummation scheme for all higher-oder corrections to one-dimensional hydrodynamic equation we demonstrate the importance of higher-order terms if the Knudsen number is large.

preprint2009arXiv

Shear viscosity and out of equilibrium dynamics

Using Grad's method, we calculate the entropy production and derive a formula for the second-order shear viscosity coefficient in a one-dimensionally expanding particle system, which can also be considered out of chemical equilibrium. For a one-dimensional expansion of gluon matter with Bjorken boost invariance, the shear tensor and the shear viscosity to entropy density ratio $η/s$ are numerically calculated by an iterative and self-consistent prescription within the second-order Israel-Stewart hydrodynamics and by a microscopic parton cascade transport theory. Compared with $η/s$ obtained using the Navier-Stokes approximation, the present result is about 20% larger at a QCD coupling $α_s \sim 0.3$(with $η/s\approx 0.18$) and is a factor of 2-3 larger at a small coupling $α_s \sim 0.01$. We demonstrate an agreement between the viscous hydrodynamic calculations and the microscopic transport results on $η/s$, except when employing a small $α_s$. On the other hand, we demonstrate that for such small $α_s$, the gluon system is far from kinetic and chemical equilibrium, which indicates the break down of second-order hydrodynamics because of the strong noneqilibrium evolution. In addition, for large $α_s$ ($0.3-0.6$), the Israel-Stewart hydrodynamics formally breaks down at large momentum $p_T\gtrsim 3$ GeV but is still a reasonably good approximation.

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.

preprint2008arXiv

QCD Matter Thermalization at RHIC and LHC

Employing the perturbative QCD inspired parton cascade, we investigate kinetic and chemical equilibration of the partonic matter created in central heavy ion collisions at RHIC and LHC energies. Two types of initial conditions are chosen. One is generated by the model of wounded nucleons using the PYTHIA event generator and Glauber geometry. Another is considered as a color glass condensate. We show that kinetic equilibration is almost independent on the chosen initial conditions, whereas there is a sensitive dependence for chemical equilibration. The time scale of thermalization lies between 1 and 1.5 fm/c. The final parton transverse energy obtained from BAMPS calculations is compared with the RHIC data and is estimated for the LHC energy.

preprint2008arXiv

Review of the "Bottom-Up" scenario

Thermalization of a longitudinally expanding color glass condensate with Bjorken boost invariant geometry is investigated within parton cascade BAMPS. Our main focus lies on the detailed comparison of thermalization, observed in BAMPS with that suggested in the Bottom-Up scenario. We demonstrate that the tremendous production of soft gluons via $gg \to ggg$, which is shown in the Bottom-Up picture as the dominant process during the early preequilibration, will not occur in heavy ion collisions at RHIC and LHC energies, because the back reaction $ggg\to gg$ hinders the absolute particle multiplication. Moreover, contrary to the Bottom-Up scenario, soft and hard gluons thermalize at the same time. The time scale of thermal equilibration in BAMPS calculations is of order $\as^{-2} (\ln \as)^{-2} Q_s^{-1}$. After this time the gluon system exhibits nearly hydrodynamic behavior. The shear viscosity to entropy density ratio has a weak dependence on $Q_s$ and lies close to the lower bound of the AdS/CFT conjecture.

preprint1998arXiv

Interpretation and resolution of pinch singularities in non-equilibrium quantum field theory

Ill-defined pinch singularities arising in a perturbative expansion in out of equilibrium quantum field theory have a natural analogue to standard scattering theory. We explicitly demonstrate that the occurrence of such terms is directly related to Fermi's golden rule known from elementary scattering theory and is thus of no mystery. We further argue that within the process of thermalization of a plasma one has to resum such contributions to all orders as the process itself is of non-perturbative nature. In this way the resummed propagators obtain a finite width. Within the Markov approximation of kinetic theory the actual phase space distribution at a given time of the evolution enters explicitly.

preprint1998arXiv

Stochastic interpretation of Kadanoff-Baym equations and their relation to Langevin processes

In this more pedagogical study we want to elucidate on stochastic aspects inherent to the (non-)equilibrium real time Green's function description (or `closed time path Green's function' -- CTPGF) of transport equations, the so called `Kadanoff-Baym equations'. As a toy model we couple a free scalar boson quantum field to an exemplaric heat bath with some given temperature T. It will be shown in detail that the emerging transport equations have to be understood as the ensemble average over stochastic equations of Langevin type. This corresponds to the equivalence of the influence functional approach by Feynman and Vernon and the CTP technique. The former, however, gives a more intuitive physical picture. In particular the physical role of (quantum) noise and the connection of its correlation kernel to the Kadanoff-Baym equations will be discussed. The inherent presence of noise and dissipation related by the fluctuation-dissipation theorem guarantees that the modes or particles become thermally populated on average in the long-time limit. For long wavelength modes with momenta much less than the temperature the emerging wave equation do behave nearly as classical. On the other hand, a kinetic transport description can be obtained in the semi-classical particle regime. Including fluctuations, its form resembles that of a phenomenological Boltzmann-Langevin description. However, we will point out some severe discrepancies in comparison to the Boltzmann- Langevin scheme. As a further byproduct we also note how the occurrence of so called pinch singularities is circumvented by a clear physical necessity of damping within the one-particle propagator.

preprint1996arXiv

Classical Fields Near Thermal Equilibrium

We discuss the classical limit for the long-distance (``soft'') modes of a quantum field when the hard modes of the field are in thermal equilibrium. We address the question of the correct semiclassical dynamics when a momentum cut-off is introduced. Higher order contributions leads to a stochastic interpretation for the effective action in analogy to Quantum Brownian Motion, resulting in dissipation and decoherence for the evolution of the soft modes. Particular emphasis is put on the understanding of dissipation. Our discussion focuses mostly on scalar fields, but we make some remarks on the extension to gauge theories.