Source author record

Michael Strickland

Michael Strickland 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

94works
11topics
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

94 published item(s)

preprint2023arXiv

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

preprint2023arXiv

Bottomonium suppression in the quark-gluon plasma -- From effective field theories to non-unitary quantum evolution

In this proceedings contribution I review recent work which computes the suppression of bottomonium production in heavy-ion collisions using open quantum systems methods applied within the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory. I discuss how the computation of bottomonium suppression can be reduced to solving a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) quantum master equation for the evolution of the b-bbar reduced density matrix. The open quantum systems approach used allows one to take into account the non-equilibrium dynamics and decoherence of bottomonium in the quark-gluon plasma. Finally, I present comparisons of phenomenological predictions obtained using a recently obtained next-to-leading-order GKSL equation with ALICE, ATLAS, and CMS experimental data for bottomonium suppression and elliptic flow.

preprint2022arXiv

${\cal N}=4$ supersymmetric Yang-Mills thermodynamics from effective field theory

The free energy density of ${\cal N}=4$ supersymmetric Yang-Mills theory in four space-time dimensions is derived through second order in the 't Hooft coupling $λ$ at finite temperature using effective-field theory methods. The contributions to the free energy density at this order come from the hard scale $T$ and the soft scale $\sqrtλ T$. The effects of the scale $T$ are encoded in the coefficients of an effective three-dimensional field theory that is obtained by dimensional reduction at finite temperature. The effects of the electric scale $\sqrtλ T$ are taken into account by perturbative calculations in the effective theory.

preprint2022arXiv

Effective field theory treatment of ${\cal N}=4$ supersymmetric Yang-Mills thermodynamics

At finite temperature the free energy density of ${\cal N}=4$ supersymmetric Yang-Mills can be calculated using resummed perturbation theory through the order $λ^{5/2}$. Effective field theory methods provide a useful alternative approach to streamline these calculations. In this proceedings contribution, I review recent work with my collaborators where we used effective field theory methods to calculate the free energy density of ${\cal N}=4$ supersymmetric Yang-Mills in four spacetime dimensions through second order in the 't Hooft coupling $λ$. At this order the contributions to the free energy density come from the hard scale $T$ and the soft scale $\sqrtλT$. The contribution from the scale $T$ enters through the coefficients in the effective Lagrangian obtained by dimensional reduction and the effects of the scale $gT$ can be calculated using perturbative methods in the effective theory.

preprint2022arXiv

Long and short distance behavior of the imaginary part of the heavy-quark potential

The imaginary part of the effective heavy-quark potential is related to the total in-medium decay width of heavy quark-antiquark bound states. We extract the static limit of this quantity using classical-statistical simulations of real-time Yang-Mills dynamics by measuring the temporal decay of Wilson loops. By performing the simulations on finer and larger lattices, we are able to show that the nonperturbative results follow the same form as the perturbative ones. For large quark-antiquark separations, we quantify the magnitude of the non-perturbative long-range corrections to the imaginary part of the heavy-quark potential. We present our results for a wide range of temperatures, lattice spacings, and lattice volumes. We also extract approximations for the short-distance behavior of the classical potential.

preprint2022arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing 0+1d Bjorken expansion which goes beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. We then demonstrate that the dynamics can be re-scaled at intermediate times and represented by universal exponents. Finally, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions

preprint2022arXiv

Non-equilibrium evolution of quarkonium in medium in the open quantum system approach

In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the quark-gluon plasma using an open quantum system (OQS) approach that is applied in the framework of the potential non-relativistic QCD (pNRQCD) effective field theory and coupled to realistic hydrodynamical backgrounds that have been tuned to soft hadron observables. I review how the computation of bottomonium suppression can be reduced to solving a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation for the evolution of the b-bbar reduced density matrix, which includes both singlet and octet states plus medium-induced transitions between them at next-to-leading order (NLO) in the binding energy over temperature. Finally, I present comparisons of phenomenological predictions of the NLO OQS+pNRQCD approach and experimental data for bottomonium suppression and elliptic flow in LHC 5.02 TeV Pb-Pb collisions obtained using both smooth and fluctuating hydrodynamic initial conditions.

preprint2022arXiv

QTRAJ 1.0: A Lindblad equation solver for heavy-quarkonium dynamics

We introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.

preprint2022arXiv

QuantumFDTD -- A computational framework for the relativistic Schrödinger equation

We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schrödinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground, first, and second excited states. We (a) include the case of the relativistic Schrödinger equation and (b) add two optimized FFT-based kinetic energy terms for the non-relativistic case. All the three new kinetic terms are computed using Fast Fourier Transform (FFT). We release the resulting code as version 3 of quantumfdtd. Finally, the code now supports arbitrary external file-based potentials and the option to project out distinct parity eigenstates from the solutions. Our goal is quark models used for phenomenological descriptions of QCD bound states, described by the three-dimensional Schrödinger equation. However, we target any field where solving either the non-relativistic or the relativistic three-dimensional Schrödinger equation is required.

preprint2022arXiv

Scheme dependence of two-loop HTLpt-resummed $\text{SYM}_{4,4}$ thermodynamics

The resummed thermodynamics of ${\cal N}=4$ supersymmetric Yang-Mills theory in four space-time dimensions ($\text{SYM}_{4,4}$) has been calculated previously to two loop order within hard thermal loop perturbation theory (HTLpt) using the canonical dimensional regularization (DRG) scheme. Herein, we revisit this calculation using the regularization by dimensional reduction (RDR) scheme. Since the RDR scheme manifestly preserves supersymmetry it is the preferred scheme, however, it is important to assess if and by how much the resummed perturbative results depend on the regularization scheme used. Comparing predictions for the scaled entropy obtained using the DRG and RDR schemes we find that for $λ\lesssim 6$ they are numerically very similar. We then compare the results obtained in both schemes with the strict perturbative result, which is accurate up to order $λ^2$, and a generalized Padé approximant constructed from the known large-$N_c$ weak- and strong-coupling expansions. Comparing the strict perturbative expansion of the two-loop HTLpt result with the perturbative expansion to order $λ^2$, we find that both the DRG and RDR HTLpt calculations result in the same scheme-independent predictions for the coefficients at order $λ$, $λ^{3/2}$, and $λ^2 \logλ$, however, at order $λ^2$ there is a residual regularization scheme dependence.

preprint2022arXiv

Some recent advances in the understanding of ${\cal N}=4$ supersymmetric Yang-Mills thermodynamics

The interest in the thermodynamics of supersymmetric Yang-Mills started after Maldacena proposed the duality between string theory on AdS backgrounds and the large-N limit of SYM theories. One of the motivations to study the thermal properties of ${\cal N}=4$ supersymmetric Yang-Mills in four dimensions (SYM$_{4,4}$) is that at high temperatures, the weak-coupling limit of this theory has many similarities with high temperature quantum chromodynamics (QCD). In this proceedings contribution, we review recent calculations of the resummed perturbative free energy of ${\cal N}=4$ supersymmetric Yang-Mills in four spacetime dimensions through second order in the 't Hooft coupling $λ$ at finite temperature and zero chemical potential. We compare our final result with prior results obtained in the weak and strong-coupling limits and construct a generalized Padé approximant that interpolates between the weak-coupling result and the large-$N_c$ strong-coupling result.

preprint2022arXiv

The effective complex heavy-quark potential in an anisotropic quark-gluon plasma

We introduce a method for reducing anisotropic heavy-quark potentials to isotropic potentials by using an effective screening mass that depends on the quantum numbers $l$ and $m$ of a given state. We demonstrate that, using the resulting 1D effective potential model, one can solve a 1D Schrödinger equation and reproduce the full 3D results for the energies and binding energies of low-lying heavy-quarkonium bound states to relatively high accuracy. This includes the splitting of different p-wave polarizations. The resulting 1D effective model provides a way to include momentum anisotropy effects in open quantum system simulations of heavy-quarkonium dynamics in the quark-gluon plasma.

preprint2022arXiv

The non-equilibrium attractor for kinetic theory in relaxation time approximation

I demonstrate that the concept of a non-equilibrium attractor can be extended beyond the lowest-order moments typically considered in hydrodynamic treatments. Using a previously obtained exact solution to the relaxation-time approximation Boltzmann equation for a transversally homogeneous and boost-invariant system subject to Bjorken flow, I derive an equation obeyed by all moments of the one-particle distribution function. Using numerical solutions, I show that, similar to the pressure anisotropy, all moments of the distribution function exhibit attractor-like behavior wherein all initial conditions converge to a universal solution after a short time with the exception of moments which are sensitive to modes with zero longitudinal momentum and high transverse momentum. In addition, I compute the exact solution for the distribution function itself on very fine lattices in momentum space and demonstrate that (a) an attractor for the full distribution function exists and (b) solutions with generic initial conditions relax to this solution, first at low momentum and later at high momentum.

preprint2021arXiv

Bottomonium suppression and elliptic flow in heavy-ion collisions

In this proceedings contribution I review recent progress concerning the suppression of bottomonium production in the quark-gluon plasma. Making use of open quantum system methods applied to potential non-relativistic quantum chromodynamics one can show that the dynamics of heavy-quarkonium bound states satisfying the scale hierarchy 1/a_0 >> pi T ~ m_D >> E obey a Lindblad equation whose solution provides the quantum evolution of the heavy-quarkonium reduced density matrix. To solve the resulting Lindblad equation we use a quantum trajectories algorithm which allows one to include all possible angular momentum states of the quark-antiquark probe in a scalable manner. We solve the Lindblad equation using a tuned 3+1D dissipative hydrodynamics code for the background temperature evolution. We then consider a large number of Monte-Carlo sampled bottomonium trajectories embedded in this background. This allows us to extract the centrality- and p_T-dependence of the nuclear suppression factor R_AA[Upsilon] and elliptic flow v_2[Upsilon]. We find good agreement between our model predictions and available sqrt(s_NN) = 5.02 TeV Pb-Pb collision experimental data from the ALICE, ATLAS, and CMS collaborations.

preprint2021arXiv

Bottomonium suppression and flow in heavy-ion collisions

The strong suppression of bottomonia production in ultra-relativistic heavy-ion collisions is a smoking gun for the creation of a deconfined quark-gluon plasma (QGP). In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the QGP through the use of pNRQCD and an open quantum systems approach. This approach allows one to evolve the heavy-quarkonium reduced density matrix, taking into account non-unitary effective Hamiltonian evolution of the wave-function and quantum jumps between different angular momentum and color states. In the case of a strong coupled QGP in which E << T,m_D << 1/a_0, the corresponding evolution equation is Markovian and can therefore be mapped to a Lindblad evolution equation. To solve the resulting Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm and couple the non-abelian quantum evolution to a realistic 3+1D viscous hydrodynamical background. Using a large number of Monte-Carlo sampled bottomonium trajectories, we make predictions for bottomonium R_AA and elliptic flow as a function of centrality and transverse momentum and compare to data collected by the ALICE, ATLAS, and CMS collaborations.

preprint2021arXiv

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

preprint2020arXiv

An improved anisotropic hydrodynamics ansatz

We introduce an improved form for the anisotropic hydrodynamics distribution function which explicitly takes into account the free-streaming and equilibrating contributions separately. We demonstrate that with this improvement one can better reproduce exact results available in the literature for the evolution of moments of the distribution function, in particular, for moments which contain no powers of the longitudinal momentum in their definition (m=0 moments). Using the resulting dynamical equations, we extract the non-equilibrium attractor associated with our improved aHydro ansatz and demonstrate that the improvement also allows one to better reproduce the exact dynamical attractor obtained using kinetic theory in the relaxation time approximation, particularly at early rescaled times and for m=0 moments.

preprint2020arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing one-dimensional Bjorken expansion. We demonstrate that the resulting far-from-equilibrium evolution is insensitive to certain features of the initial condition, including both the initial momentum-space anisotropy and initial occupancy. We find that this insensitivity extends beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. Based on our results, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions.

preprint2020arXiv

Pseudothermalization of the quark-gluon plasma

In this proceedings contribution I review recent work in kinetic theory which demonstrates that, for system undergoing Bjorken expansion, there exists an attractor in all moments of the one-particle distribution function. I discuss how this attractor emerges in both exact solutions obtained in relaxation time approximation (RTA) and the effective kinetic theory approach to high-temperature quantum chromodynamics (QCD). The QCD effective kinetic theory collisional kernel used includes both elastic (2 <-> 2) and LPM-resummed inelastic (2 <-> 1) contributions. The results obtained indicate that a pseudothermal attractor exists in both RTA and QCD kinetic theory and that their respective attractors can be extended to early times when the system is far from equilibrium. Finally, I discuss how knowledge of the QCD effective kinetic theory attractor can be used to assess different hydrodynamic freeze-out prescriptions used in heavy-ion phenomenology. The results obtained show that improved freeze-out prescriptions such as anisotropic hydrodynamics perform better in conditions corresponding to those generated in high-multiplicity pA and pp collisions, e.g. short lifetime and high inverse Reynolds number.

preprint2020arXiv

Two-loop HTL-resummed thermodynamics for N=4 supersymmetric Yang-Mills theory

We compute the two-loop hard-thermal-loop (HTL) resummed thermodynamic potential for N=4 supersymmetric Yang-Mills (SYM). Our final result is manifestly gauge-invariant and was renormalized using only simple vacuum energy, gluon mass, scalar mass, and quark mass counter terms. The HTL mass parameters m_D, M_D, and m_q are then determined self-consistently using a variational prescription which results in a set of coupled gap equations. Based on this, we obtain the two-loop HTL-resummed thermodynamic functions of N=4 SYM. We compare our final result with known results obtained in the weak- and strong-coupling limits. We also compare to previously obtained approximately self-consistent HTL resummations and Padé approximants. We find that the two-loop HTL resummed results for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 <= lambda <~ 2 and is in agreement with previous approximately self-consistent HTL resummation results for lambda <~ 6.

preprint2019arXiv

Exact solution for the non-equilibrium attractor in number-conserving relaxation time approximation

We extend previous studies of the conformal 0+1d kinetic non-equilibrium attractor in relaxation time approximation by enforcing number conservation through the introduction of a dynamical fugacity (chemical potential). We derive two coupled integral equations for the effective temperature and fugacity which are then solved numerically to obtain the exact solution. The resulting solutions exhibit convergence to a unique non-equilibrium attractor when the scaled moments of the distribution function are plotted as a function of the rescaled time w = tau/tau_eq. This occurs even though the system is out of chemical equilibrium at late times. In addition, compared to the case where number conservation was not imposed, we find that the moments converge to their respective attractors more quickly, particularly for moments with m=0. Finally, we compare the resulting attractor moments with predictions from different hydrodynamic frameworks.

preprint2019arXiv

Photon production and elliptic flow from momentum-anisotropic quark-gluon plasma

The emission of real photons from a momentum-anisotropic quark-gluon plasma (QGP) is affected by both the collective flow of the radiating medium and the modification of local rest frame emission rate due to the anisotropic momentum distribution of partonic degrees of freedom. In this paper, we first calculate the photon production rate from an ellipsoidally momentum-anisotropic QGP including hard contributions from Compton scattering and quark pair annihilation and soft contribution calculated using the hard thermal loop (HTL) approximation. We introduce a parametrization of the nonequilibrium rate in order to facilitate its further application in yield and flow calculations. We convolve the anisotropic photon rate with the space-time evolution of QGP provided by 3+1d anisotropic hydrodynamics (aHydro) to obtain the yield and the elliptic flow coefficient $v_2$ of photons from QGP generated at Pb-Pb collisions at LHC at 2.76 TeV and Au-Au collisions at RHIC at 200 GeV. We investigate the effects of various parameters on the results. In particular we analyze the sensitivity of results to initial momentum anisotropy.

preprint2016arXiv

Dilepton rate and quark number susceptibility with the Gribov action

We use a recently obtained resummed quark propagator at finite temperature which takes into account both the chromoelectric scale gT and the chromomagnetic scale g^2T through the Gribov action. The electric scale generates two massive modes whereas the magnetic scale produces a new massless spacelike mode in the medium. Moreover, the non-perturbative quark propagator is found to contain no discontinuity in contrast to the standard perturbative hard thermal loop approach. Using this non-perturbative quark propagator and self-consistent vertices, we compute the non-perturbative dilepton rate at vanishing three-momentum at one-loop order. The resulting rate has a rich structure at low energies due to the inclusion of the non-perturbative magnetic scale. We also calculate the quark number susceptibility, which is related to the conserved quark number density fluctuation in the deconfined state. Both the dilepton rate and quark number susceptibility are compared with results from lattice quantum chromodynamics and the standard hard thermal loop approach. Finally, we discuss how the absence of a discontinuity in the imaginary part of the non-perturbative quark propagator makes the results for both dilepton production and quark number susceptibility dramatically different from those in perturbative approaches and seemingly in conflict with known lattice data.

preprint2016arXiv

Non-boost-invariant dissipative hydrodynamics

The one-dimensional non-boost-invariant evolution of the quark-gluon plasma, presumably produced during the early stages of heavy-ion collisions, is analyzed within the frameworks of viscous and anisotropic hydrodynamics. We neglect transverse dynamics and assume homogeneous conditions in the transverse plane but, differently from Bjorken expansion, we relax longitudinal boost invariance in order to study the rapidity dependence of various hydrodynamical observables. We compare the results obtained using several formulations of second-order viscous hydrodynamics with a recent approach to anisotropic hydrodynamics, which treats the large initial pressure anisotropy in a non-perturbative fashion. The results obtained with second-order viscous hydrodynamics depend on the particular choice of the second-order terms included, which suggests that the latter should be included in the most complete way. The results of anisotropic hydrodynamics and viscous hydrodynamics agree for the central hot part of the system, however, they differ at the edges where the approach of anisotropic hydrodynamics helps to control the undesirable growth of viscous corrections observed in standard frameworks.

preprint2016arXiv

Photon production from a non-equilibrium quark-gluon plasma

We calculate leading-order medium photon yields from a quark-gluon plasma using (3+1)-dimensional anisotropic hydrodynamics. Non-equilibrium corrections to the photon rate are taken into account using a self-consistent modification of the particle distribution functions and the corresponding anisotropic hard-loop fermionic self-energies. We present predictions for the high-energy photon spectrum and photon elliptic flow as a function of transverse momentum, shear viscosity, and initial momentum-space anisotropy. Our findings indicate that high-energy photon production is sensitive to the assumed level of initial momentum-space anisotropy of the quark-gluon plasma. As a result, it may be possible to experimentally constrain the early-time momentum-space anisotropy of the quark-gluon plasma generated in relativistic heavy-ion collisions using high-energy photon yields.

preprint2016arXiv

Predictions for bottomonia suppression in 5.023 TeV Pb-Pb collisions

We compute the suppression of the bottomonia states Upsilon(1S), Upsilon(2S), Upsilon(3S), chi_b(1P), chi_b(2P), and chi_b(3P) states in Large Hadron Collider (LHC) sqrt(s_NN)) = 5.023 TeV Pb-Pb collisions. For the background evolution we use 3+1d anisotropic hydrodynamics with conditions extrapolated from sqrt(s_NN) = 2.76 TeV and we self-consistently compute bottomonia decay rates including non-equilibrium corrections to the interaction potential. For our final results, we make predictions for R_AA as function of centrality, rapidity, and p_T for the Upsilon(1S) and Upsilon(2S) states, including feed down effects. In order to assess the dependence on some of the model assumptions, we vary the shear viscosity-to-entropy density ratio, 4 pi eta/s = (1, 2, 3), and the initial momentum-space anisotropy parameter, xi_0 = (0, 10, 50), while holding the total light hadron multiplicity fixed.

preprint2016arXiv

Quark self-energy in an ellipsoidally anisotropic quark-gluon plasma

We calculate the quark self-energy in a quark-gluon plasma that possesses an ellipsoidal momentum-space anisotropy in the local rest frame. By introducing additional transverse momentum anisotropy parameters into the parton distribution functions, we generalize previous results which were obtained for the case of a spheroidal anisotropy. Our results demonstrate that the presence of anisotropies in the transverse directions affects the real and imaginary parts of quark self-energy and, consequently, the self-energy depends on both the polar and azimuthal angles in the local rest frame of the matter. Our results for the quark self-energy set the stage for the calculation of the effects of ellipsoidal momentum-space anisotropy on quark-gluon plasma photon spectra and collective flow.

preprint2016arXiv

Three-loop hard-thermal-loop perturbation theory thermodynamics at finite temperature and finite baryonic and isospin chemical potential

In a previous paper (JHEP {\bf 05} (2014) 27), we calculated the three-loop thermodynamic potential of QCD at finite temperature $T$ and quark chemical potentials $μ_q$ using the hard-thermal-loop perturbation theory (HTLpt) reorganization of finite temperature and density QCD. The result allows us to study the thermodynamics of QCD at finite temperature and isospin chemical potential $μ_I$. We calculate the pressure, energy density, and entropy density, the trace anomaly, and the speed of sound at zero and nonzero $μ_I$. The second, fourth, and sixth-order isospin susceptibilities are calculated at zero $μ_I$. Our results can be directly compared to lattice QCD without Taylor expansions around $μ_q=0$ since QCD has no sign problem at finite isospin chemical potential.

preprint2015arXiv

Anisotropic hydrodynamics for conformal Gubser flow

In this proceedings contribution, we review the exact solution of the anisotropic hydrodynamics equations for a system subject to Gubser flow. For this purpose, we use the leading-order anisotropic hydrodynamics equations which assume that the distribution function is ellipsoidally symmetric in local-rest-frame momentum. We then prove that the SO(3)_q symmetry in de Sitter space constrains the anisotropy tensor to be of spheroidal form with only one independent anisotropy parameter remaining. As a consequence, the exact solution reduces to the problem of solving two coupled non-linear differential equations. We show that, in the limit that the relaxation time goes to zero, one obtains Gubser's ideal hydrodynamic solution and, in the limit that the relaxation time goes to infinity, one obtains the exact free streaming solution obtained originally by Denicol et al. For finite relaxation time, we solve the equations numerically and compare to the exact solution of the relaxation-time-approximation Boltzmann equation subject to Gubser flow. Using this as our standard, we find that anisotropic hydrodynamics describes the spatio-temporal evolution of the system better than all currently known dissipative hydrodynamics approaches.

preprint2015arXiv

Bottomonia suppression in 2.76 TeV Pb-Pb collisions

We compute the QGP suppression of Upsilon(1s), Upsilon(2s), Upsilon(3s), chi_b1, and chi_b2 states in sqrt(s_NN)=2.76 TeV Pb-Pb collisions. Using the suppression of each of these states, we estimate the inclusive R_AA for the Upsilon(1s) and Upsilon(2s) states as a function of N_part, y, and p_T including the effect of excited state feed down. We find that our model provides a reasonable description of preliminary CMS results for the N_part-, y-, and p_T-dependence of R_AA for both the Upsilon(1s) and Upsilon(2s). Comparing to our previous model predictions, we find a flatter rapidity dependence, thereby reducing some of the tension between our model and ALICE forward-rapidity results for Upsilon(1s) suppression.

preprint2015arXiv

Dilepton production from the quark-gluon plasma using (3+1)-dimensional anisotropic dissipative hydrodynamics

We compute dilepton production from the deconfined phase of the quark-gluon plasma using leading-order (3+1)-dimensional anisotropic hydrodynamics. The anisotropic hydrodynamics equa- tions employed describe the full spatiotemporal evolution of the transverse temperature, spheroidal momentum-space anisotropy parameter, and the associated three-dimensional collective flow of the matter. The momentum-space anisotropy is also taken into account in the computation of the dilepton production rate, allowing for a self-consistent description of dilepton production from the quark-gluon plasma. For our final results, we present predictions for high-energy dilepton yields as a function of invariant mass, transverse momentum, and pair rapidity. We demonstrate that high- energy dilepton production is extremely sensitive to the assumed level of initial momentum-space anisotropy of the quark-gluon plasma. As a result, it may be possible to experimentally constrain the early-time momentum-space anisotropy of the quark-gluon plasma generated in relativistic heavy ion collisions using high-energy dilepton yields.

preprint2015arXiv

Dilepton production from the quark-gluon plasma using leading-order (3+1)D anisotropic hydrodynamics

Dilepton production from the quark-gluon plasma (QGP) phase of ultra-relativistic heavy-ion collisions is computed using the leading-order (3+1)-dimensional anisotropic hydrodynamics. It is shown that high-energy dilepton spectrum is sensitive to the initial local-rest-frame momentum-space anisotropy of the QGP. Our findings suggest that it may be possible to constrain the early-time momentum-space anisotropy in relativistic heavy-ion collisions using high-energy dilepton yields.

preprint2015arXiv

Leading-order anisotropic hydrodynamics for central collisions

We use leading-order anisotropic hydrodynamics to study an azimuthally-symmetric boost-invariant quark-gluon plasma. We impose a realistic lattice-based equation of state and perform self-consistent anisotropic freeze-out to hadronic degrees of freedom. We then compare our results for the full spatiotemporal evolution of the quark-gluon plasma and its subsequent freeze-out to results obtained using 1+1d Israel-Stewart second-order viscous hydrodynamics. We find that for small shear viscosities, 4 pi eta/s ~ 1, the two methods agree well for nucleus-nucleus collisions, however, for large shear viscosity to entropy density ratios or proton-nucleus collisions we find important corrections to the Israel-Stewart results for the final particle spectra and the total number of charged particles. Finally, we demonstrate that the total number of charged particles produced is a monotonically increasing function of 4 pi eta/s in Israel-Stewart viscous hydrodynamics whereas in anisotropic hydrodynamics it has a maximum at 4 pi eta/s ~ 10. For all 4 pi eta/s > 0, we find that for Pb-Pb collisions Israel-Stewart viscous hydrodynamics predicts more dissipative particle production than anisotropic hydrodynamics.

preprint2015arXiv

Photon and dilepton production from a non-equilibrium quark-gluon plasma

We calculate leading-order medium photon and dilepton yields from a quark-gluon plasma using (3+1)-dimensional anisotropic hydrodynamics. In the case of dileptons, at leading-order it is sufficient to take into account non-equilibrium corrections to the rate through the use of anisotropic distribution functions. In the case of photons, non-equilibrium corrections to the rate are taken into account using a self-consistent modification of the distribution functions and the corresponding anisotropic hard-loop quark self-energies. We present predictions for the high-energy photon and dilepton spectrum and the photon elliptic flow as a function of invariant mass, transverse momentum, shear viscosity, and initial momentum-space anisotropy. Our findings indicate that both high-energy photon and dilepton production are sensitive to the assumed level of initial momentum-space anisotropy of the quark-gluon plasma. In addition, we find that the photon elliptic flow depends on the initial momentum-space anisotropy.

preprint2015arXiv

Quasiparticle equation of state for anisotropic hydrodynamics

We present a new method for imposing a realistic equation of state in anisotropic hydrodynamics. The method relies on the introduction of a single finite-temperature quasiparticle mass which is fit to lattice data. By taking moments of the Boltzmann equation, we obtain a set of coupled partial differential equations which can be used to describe the 3+1d spacetime evolution of an anisotropic relativistic system. We then specialize to the case of a 0+1d system undergoing boost-invariant Bjorken expansion and subject to the relaxation-time approximation collisional kernel. Using this setup, we compare results obtained using the new quasiparticle equation of state method with those obtained using the standard method for imposing the equation of state in anisotropic hydrodynamics. We demonstrate that the temperature evolution obtained using the two methods is nearly identical and that there are only small differences in the pressure anisotropy. However, we find that there are significant differences in the evolution of the bulk pressure correction.

preprint2015arXiv

Relativistic quantum transport coefficients for second-order viscous hydrodynamics

We express the transport coefficients appearing in the second-order evolution equations for bulk viscous pressure and shear stress tensor using Bose-Einstein, Boltzmann, and Fermi-Dirac statistics for the equilibrium distribution function and Grad's 14-moment approximation as well as the method of Chapman-Enskog expansion for the non-equilibrium part. Specializing to the case of transversally homogeneous and boost-invariant longitudinal expansion of the viscous medium, we compare the results obtained using the above methods with those obtained from the exact solution of the massive 0+1d relativistic Boltzmann equation in the relaxation-time approximation. We show that compared to the 14-moment approximation, the hydrodynamic transport coefficients obtained by employing the Chapman-Enskog method leads to better agreement with the exact solution of the relativistic Boltzmann equation.

preprint2015arXiv

Testing different formulations of leading-order anisotropic hydrodynamics

A recently obtained set of the equations for leading-order (3+1)D anisotropic hydrodynamics is tested against exact solutions of the Boltzmann equation with the collisional kernel treated in the relaxation time approximation. In order to perform the detailed comparisons, the new anisotropic hydrodynamics equations are reduced to the boost-invariant and transversally homogeneous case. The agreement with the exact solutions found using the new anisotropic hydrodynamics equations is similar to that found using previous, less general, formulations of anisotropic hydrodynamics. In addition, we find that, when compared to a state-of-the-art second-order viscous hydrodynamics framework, leading-order anisotropic hydrodynamics better reproduces the exact solution for the pressure anisotropy and gives comparable results for the bulk pressure evolution. Finally, we compare the transport coefficients obtained using linearized anisotropic hydrodynamics with results obtained using second-order viscous hydrodynamics.

preprint2015arXiv

Three-loop HTLpt thermodynamics at finite temperature and chemical potential

In this proceedings we present a state-of-the-art method of calculating thermodynamic potential at finite temperature and finite chemical potential, using Hard Thermal Loop perturbation theory (HTLpt) up to next-to-next-leading-order (NNLO). The resulting thermodynamic potential enables us to evaluate different thermodynamic quantities including pressure and various quark number susceptibilities (QNS). Comparison between our analytic results for those thermodynamic quantities with the available lattice data shows a good agreement.

preprint2014arXiv

A new exact solution of the relativistic Boltzmann equation and its hydrodynamic limit

We present an exact solution of the relativistic Boltzmann equation for a system undergoing boost-invariant longitudinal and azimuthally symmetric transverse flow ("Gubser flow"). The resulting exact non-equilibrium dynamics is compared to 1st- and 2nd-order relativistic hydrodynamic approximations for various shear viscosity to entropy density ratios. This novel solution can be used to test the validity and accuracy of different hydrodynamic approximations in conditions similar to those generated in relativistic heavy-ion collisions.

preprint2014arXiv

Anisotropic hydrodynamics for conformal Gubser flow

We derive the equations of motion for a system undergoing boost-invariant longitudinal and azimuthally-symmetric transverse "Gubser flow" using leading-order anisotropic hydrodynamics. This is accomplished by assuming that the one-particle distribution function is ellipsoidally-symmetric in the momenta conjugate to the de Sitter coordinates used to parameterize the Gubser flow. We then demonstrate that the SO(3)_q symmetry in de Sitter space further constrains the anisotropy tensor to be of spheroidal form. The resulting system of two coupled ordinary differential equations for the de Sitter-space momentum scale and anisotropy parameter are solved numerically and compared to a recently obtained exact solution of the relaxation-time-approximation Boltzmann equation subject to the same flow. We show that anisotropic hydrodynamics describes the spatio-temporal evolution of the system better than all currently known dissipative hydrodynamics approaches. In addition, we prove that anisotropic hydrodynamics gives the exact solution of the relaxation-time approximation Boltzmann equation in the ideal, eta/s -> 0, and free-streaming, eta/s -> infinity, limits.

preprint2014arXiv

Anisotropic Hydrodynamics: Motivation and Methodology

In this proceedings contribution I review recent progress in our understanding of the bulk dynamics of relativistic systems that possess potentially large local rest frame momentum-space anisotropies. In order to deal with these momentum-space anisotropies, a reorganization of relativistic viscous hydrodynamics can been made around an anisotropic background, and the resulting dynamical framework has been dubbed "anisotropic hydrodynamics." I also discuss expectations for the degree of momentum-space anisotropy of the quark gluon plasma generated in relativistic heavy ion collisions at RHIC and LHC from second-order viscous hydrodynamics, strong-coupling approaches, and weak-coupling approaches.

preprint2014arXiv

Anisotropic Hydrodynamics: Three lectures

Anisotropic hydrodynamics is a non-perturbative reorganization of relativistic hydrodynamics that takes into account the large momentum-space anisotropies generated in ultrarelativistic heavy-ion collisions. As a result, it allows one to extend the regime of applicability of hydrodynamic treatments to situations that can be quite far from isotropic thermal equilibrium. In this paper, I review the material presented in a series of three introductory lectures. I review the derivation of ideal and second-order viscous hydrodynamics from kinetic theory. I then show how to extend the methods used to a system that can be highly anisotropic in local-rest-frame momenta. I close by discussing recent work on this topic and then present an outlook to the future.

preprint2014arXiv

Bulk viscous evolution within anisotropic hydrodynamics

We derive a system of moment-based dynamical equations that describe the 1+1d space-time evolution of a cylindrically symmetric massive gas undergoing boost-invariant longitudinal expansion. Extending previous work, we introduce an explicit degree of freedom associated with the bulk pressure of the system. The resulting form generalizes the ellipsoidal one-particle distribution function appropriate for massless particles to massive particles. Using this generalized form, we obtain a system of partial differential equations that can be solved numerically. In order to assess the performance of this scheme, we compare the resulting anisotropic hydrodynamics solutions with the exact solution of the 0+1d Boltzmann equation in the relaxation time approximation. We find that the inclusion of the bulk degree of freedom improves agreement between anisotropic hydrodynamics and the exact solution for a massive gas.

preprint2014arXiv

Equation of State for QCD at finite temperature and density. Resummation versus lattice data

The perturbative series for finite-temperature field theories has very poor convergence properties and one needs a way to reorganize it. In this talk, I review two ways of reorganizing the perturbative series for field theories at finite temperature and chemical potential, namely hard-thermal-loop perturbation theory (HTLpt) and dimensional reduction (DR). I will present results for the pressure, trace anomaly, speed of sound and the quark susceptibilities from a 3-loop HTLpt calculation and for the quark susceptibilities using DR at four loops. A careful comparison with available lattice data shows good agreement for a number of physical quantities.

preprint2014arXiv

Exact solution of the (0+1)-dimensional Boltzmann equation for a massive gas

We solve the one-dimensional boost-invariant kinetic equation for a relativistic massive system with the collision term treated in the relaxation time approximation. The result is an exact integral equation which can be solved numerically by the method of iteration to arbitrary precision. We compare predictions for the shear and bulk viscosities of a massive system with those obtained from the exact solution. Finally, we compare the time evolution of the bulk pressure obtained from our exact solution with results obtained from the dynamical equations of second-order viscous hydrodynamics.

preprint2014arXiv

Leading-order anisotropic hydrodynamics for systems with massive particles

The framework of anisotropic hydrodynamics is generalized to include finite particle masses. Two schemes are introduced and their predictions compared with exact solutions of the kinetic equation in the relaxation time approximation. The first formulation uses the zeroth and first moments of the kinetic equation, whereas the second formulation uses the first and second moments. For the case of one-dimensional boost-invariant expansion, our numerical results indicate that the second formulation yields much better agreement with the exact solutions.

preprint2014arXiv

Non-boost-invariant anisotropic dynamics

We study the non-boost-invariant evolution of a quark-gluon plasma subject to large early-time momentum-space anisotropies. Rather than using the canonical hydrodynamical expansion of the distribution function around an isotropic equilibrium state, we expand around a state which is anisotropic in momentum space and parameterize this state in terms of three proper-time and spatial-rapidity dependent parameters. Deviations from the Bjorken scaling solutions are naturally taken into account by the time evolution of the spatial-rapidity dependence of the anisotropic ansatz. As a result, we obtain three coupled partial differential equations for the momentum-space anisotropy, the typical momentum of the degrees of freedom, and the longitudinal flow. Within this framework (0+1)-dimensional Bjorken expansion is obtained as an asymptotic limit. Finally, we make quantitative comparisons of the temporal and spatial-rapidity evolution of the dynamical parameters and resulting pressure anisotropy in both the strong and weak coupling limits.

preprint2014arXiv

Second-order (2+1)-dimensional anisotropic hydrodynamics

We present a complete formulation of second-order (2+1)-dimensional anisotropic hydrodynamics. The resulting framework generalizes leading-order anisotropic hydrodynamics by allowing for deviations of the one-particle distribution function from the spheroidal form assumed at leading order. We derive complete second-order equations of motion for the additional terms in the macroscopic currents generated by these deviations from their kinetic definition using a Grad-Israel-Stewart 14-moment ansatz. The result is a set of coupled partial differential equations for the momentum-space anisotropy parameter, effective temperature, the transverse components of the fluid four-velocity, and the viscous tensor components generated by deviations of the distribution from spheroidal form. We then perform a quantitative test of our approach by applying it to the case of one-dimensional boost-invariant expansion in the relaxation time approximation (RTA) in which case it is possible to numerically solve the Boltzmann equation exactly. We demonstrate that the second-order anisotropic hydrodynamics approach provides an excellent approximation to the exact (0+1)-dimensional RTA solution for both small and large values of the shear viscosity.

preprint2014arXiv

Shear-bulk coupling in nonconformal hydrodynamics

We compute the temporal evolution of the pressure anisotropy and bulk pressure of a massive gas using second-order viscous hydrodynamics and anisotropic hydrodynamics. We then compare our results with an exact solution of the Boltzmann equation for a massive gas in the relaxation time approximation. We demonstrate that, within second-order viscous hydrodynamics, the inclusion of the full set of kinetic coefficients, particularly the shear-bulk couplings, is necessary to properly describe the time evolution of the bulk pressure. We also compare the results of second-order hydrodynamics with those obtained using the anisotropic hydrodynamics approach. We find that anisotropic hydrodynamics and second-order viscous hydrodynamics including the shear-bulk couplings are both able to reproduce the exact evolution with comparable accuracy.

preprint2014arXiv

Studying the validity of relativistic hydrodynamics with a new exact solution of the Boltzmann equation

We present an exact solution to the Boltzmann equation which describes a system undergoing boost-invariant longitudinal and azimuthally symmetric radial expansion for arbitrary shear viscosity to entropy density ratio. This new solution is constructed by considering the conformal map between Minkowski space and the direct product of three dimensional de Sitter space with a line. The resulting solution respects SO(3)_q x SO(1,1) x Z_2 symmetry. We compare the exact kinetic solution with exact solutions of the corresponding macroscopic equations that were obtained from the kinetic theory in ideal and second-order viscous hydrodynamic approximations. The macroscopic solutions are obtained in de Sitter space and are subject to the same symmetries used to obtain the exact kinetic solution.

preprint2014arXiv

Three loop HTL perturbation theory at finite temperature and chemical potential

In this proceedings contribution we present a recent three-loop hard-thermal-loop perturbation theory (HTLpt) calculation of the thermodynamic potential for a finite temperature and chemical potential system of quarks and gluons. We compare the resulting pressure, trace anomaly, and diagonal/off-diagonal quark susceptibilities with lattice data. We show that there is good agreement between the three-loop HTLpt analytic result and available lattice data.

preprint2014arXiv

Three-loop HTLpt Pressure and Susceptibilities at Finite Temperature and Density

We present results of a three-loop hard-thermal-loop perturbation theory calculation of the thermodynamical potential of a finite temperature and baryon chemical potential system of quarks and gluons. We compare the resulting pressure and diagonal quark susceptibilities with available lattice data. We find reasonable agreement between our analytic results and lattice data at both zero and finite chemical potential.

preprint2014arXiv

Three-loop HTLpt thermodynamics at finite temperature and chemical potential

We calculate the three-loop thermodynamic potential of QCD at finite temperature and chemical potential(s) using the hard-thermal-loop perturbation theory (HTLpt) reorganization of finite temperature and density QCD. The resulting analytic thermodynamic potential allows us to compute the pressure, energy density, and entropy density of the quark-gluon plasma. Using these we calculate the trace anomaly, speed of sound, and second-, fourth-, and sixth-order quark number susceptibilities. For all observables considered we find good agreement between our three-loop HTLpt calculations and available lattice data for temperatures above approximately 300 MeV.

preprint2014arXiv

Transport coefficients for bulk viscous evolution in the relaxation time approximation

We derive the form of the viscous corrections to the phase-space distribution function due to the bulk viscous pressure and shear stress tensor using the iterative Chapman-Enskog method. We then calculate the transport coefficients necessary for the second-order hydrodynamic evolution of the bulk viscous pressure and the shear stress tensor. We demonstrate that the transport coefficients obtained using the Chapman-Enskog method are different than those obtained previously using the 14-moment approximation for a finite particle mass. Specializing to the case of boost-invariant and transversally homogeneous longitudinal expansion, we show that the transport coefficients obtained using the Chapman-Enskog method result in better agreement with the exact solution of the Boltzmann equation in the relaxation-time approximation compared to results obtained in the 14-moment approximation. Finally, we explicitly confirm that the time evolution of the bulk viscous pressure is significantly affected by its coupling to the shear stress tensor.

preprint2014arXiv

Viscous hydrodynamics for strongly anisotropic expansion

A new formulation of second-order viscous hydrodynamics, based on an expansion around a locally anisotropic momentum distribution, is presented. It generalizes the previously developed formalism of anisotropic hydrodynamics (aHydro) to include a complete set of dissipative currents for which equations of motion are derived by solving the Boltzmann equation in the 14-moment approximation. By solving the vaHydro equations for a transversally homogeneous, longitudinally boost-invariant system ((0+1)-dimensional expansion) and comparing with the exact solution of the Boltzmann equation in relaxation-time approximation we show that vaHydro performs much better than all other known second-order viscous hydrodynamic approximations.

preprint2013arXiv

Anisotropic hydrodynamics -- basic concepts

Due to the rapid longitudinal expansion of the quark-gluon plasma created in relativistic heavy ion collisions, potentially large local rest frame momentum-space anisotropies are generated. The magnitude of these momentum-space anisotropies can be so large as to violate the central assumption of canonical viscous hydrodynamical treatments which linearize around an isotropic background. In order to better describe the early-time dynamics of the quark gluon plasma, one can consider instead expanding around a locally anisotropic background which results in a dynamical framework called anisotropic hydrodynamics. In this proceedings contribution we review the basic concepts of the anisotropic hydrodynamics framework presenting viewpoints from both the phenomenological and microscopic points of view.

preprint2013arXiv

Anisotropic Hydrodynamics for Rapidly Expanding Systems

We exactly solve the relaxation-time approximation Boltzmann equation for a system which is transversely homogeneous and undergoing boost-invariant longitudinal expansion. We compare the resulting exact numerical solution with approximate solutions available in the anisotropic hydrodynamics and second order viscous hydrodynamics frameworks. In all cases studied, we find that the anisotropic hydrodynamics framework is a better approximation to the exact solution than traditional viscous hydrodynamical approaches.

preprint2013arXiv

Bulk Properties of a Fermi Gas in a Magnetic Field

We calculate the number density, energy density, transverse pressure, longitudinal pressure, and magnetization of an ensemble of spin one-half particles in the presence of a homogenous background magnetic field. The magnetic field direction breaks spherical symmetry causing the pressure transverse to the magnetic field direction to be different than the pressure parallel to it. We present explicit formulae appropriate at zero and finite temperature for both charged and uncharged particles including the effect of the anomalous magnetic moment. We demonstrate that the resulting expressions satisfy the canonical relations, Omega = - P_parallel and P_perp = P_parallel - M B, with M = - d Omega/d B being the magnetization of the system. We numerically calculate the resulting pressure anisotropy for a gas of protons and a gas of neutrons and demonstrate that the inclusion of the anomalous magnetic increases the level of pressure anisotropy in both cases.

preprint2013arXiv

Charmonia and Bottomonia in a Magnetic Field

We study the effect of a static homogeneous external magnetic field on charmonium and bottomonium states. In an external magnetic field, quarkonium states do not have a conserved center-of-mass momentum. Instead there is a new conserved quantity called the pseudomomentum which takes into account the Lorentz force on the particles in the system. When written in terms of the pseudomomentum, the internal and center-of-mass motions do not decouple and, as a result, the properties of quarkonia depend on the states' center-of-mass momentum. We analyze the behavior of heavy particle-antiparticle pairs subject to an external magnetic field assuming a three dimensional harmonic potential and Cornell potential plus spin-spin interaction. In the case of the Cornell potential, we also take into account the mixing of the eta_c and J/psi states and eta_b and Upsilon states due to the background magnetic field. We then numerically calculate the dependence of the masses and mixing fractions on the magnitude of the background magnetic field and center-of-mass momentum of the state.

preprint2013arXiv

Equation of State of hot and dense QCD: Resummed perturbation theory confronts lattice data

We perform a detailed analysis of the predictions of resummed perturbation theory for the pressure and the second-, fourth-, and sixth-order diagonal quark number susceptibilities in a hot and dense quark-gluon plasma. First, we present an exact one-loop calculation of the equation of state within hard-thermal-loop perturbation theory (HTLpt) and compare it to a previous one-loop HTLpt calculation that employed an expansion in the ratios of thermal masses and the temperature. We find that this expansion converges reasonably fast. We then perform a resummation of the existing four-loop weak coupling expression for the pressure, motivated by dimensional reduction. Finally, we compare the exact one-loop HTLpt and resummed dimensional reduction results with state-of-the-art lattice calculations and a recent mass-expanded three-loop HTLpt calculation.

preprint2013arXiv

Hydrodynamics of anisotropic quark and gluon fluids

The recently developed framework of anisotropic hydrodynamics is generalized to describe the dynamics of coupled quark and gluon fluids. The quark and gluon components of the fluids are characterized by different dynamical anisotropy parameters. The dynamical equations describing such mixtures are derived from kinetic theory with the collisional kernel treated in the relaxation- time approximation allowing for different relaxation times for quarks and gluons. Baryon number conservation is enforced in the quark and anti-quark components of the fluid, but overall parton number non-conservation is allowed in the system. The resulting equations are solved numerically in the (0+1)-dimensional boost-invariant case at zero and finite baryon density.

preprint2013arXiv

Longitudinal thermalization via the chromo-Weibel instability

Non-Abelian plasma instabilities play an important role in the non-equilibrium dynamics of a weakly coupled quark-gluon plasma. Using the discretized hard loop framework we calculate the time evolution of soft gluonic fields in a longitudinally free streaming background. Extrapolating our results to energies probed in relativistic heavy-ion collisions we find a pressure anisotropy that persists for a few fm/c. However, the chromofields quickly develop a Boltzmann longitudi- nal energy spectrum, suggesting fast longitudinal thermalization of the quark gluon plasma even though it remains momentum-space anisotropic. In this proceedings contribution we review our recent numerical results, present new results for the scaling of the isotropization time with the initial current fluctuation amplitude, and present tests of the gauge invariance of the extracted longitudinal spectra.

preprint2013arXiv

Quark Number Susceptibilities from Two-Loop Hard Thermal Loop Perturbation Theory

We use the recently obtained two-loop hard thermal loop perturbation theory thermodynamics functions of a plasma of quarks and gluons to compute the diagonal second- and fourth-order quark number susceptibilities. The two-loop hard thermal loop perturbation theory thermodynamic functions used are reliable in the limit that the ratio of the quark chemical potential to temperature is small. Using this result, we are able to obtain (semi-)analytic expressions for the quark number susceptibilities at leading- and next-to-leading-order in hard thermal loop perturbation theory. We compare the hard thermal loop perturbation theory results with perturbative quantum chromodynamics calculations, a Polyakov-loop Nambu-Jona-Lasinio model calculation, and lattice quantum chromodynamics results.

preprint2013arXiv

Quarkonia in the Quark Gluon Plasma

In this paper we review recent progress towards understanding the nature of quarkonia in the quark-gluon plasma. We review the theory necessary to understand the melting of bound states due to color-screening, including lattice results for the heavy quark potential, lattice results on the correlation functions related to the relevant spectral functions, and the emergence of a complex-valued potential in high-temperature quantum chromodynamics. We close with a brief survey of phenomenological models of quarkonium suppression in relativistic heavy ion collisions.

preprint2013arXiv

Testing viscous and anisotropic hydrodynamics in an exactly solvable case

We exactly solve the one-dimensional boost-invariant Boltzmann equation in the relaxation time approximation for arbitrary shear viscosity. The results are compared with the predictions of viscous and anisotropic hydrodynamics. Studying different non-equilibrium cases and comparing the exact kinetic-theory results to the second-order viscous hydrodynamics results we find that recent formulations of second-order viscous hydrodynamics agree better with the exact solution than the standard Israel-Stewart approach. Additionally, we find that, given the appropriate connection between the kinetic and anisotropic hydrodynamics relaxation times, anisotropic hydrodynamics provides a very good approximation to the exact relaxation time approximation solution.

preprint2013arXiv

The chromo-Weibel instability in an expanding background

In this proceedings contribution we review recent calculations of the dynamics of the chromo-Weibel instability in the quark gluon plasma. This instability is present in gauge theories with a one-particle distribution function which is momentum-space anisotropic in the local rest frame. The conditions necessary for triggering this instability can be present already in the color-glass-condensate initial state or dynamically generated by the rapid longitudinal expansion of the matter created in a heavy-ion collision. Using the hard-loop framework we study the case that the one-particle distribution function possesses an arbitrary initial momentum anisotropy that increases in time due to longitudinal free streaming. The resulting three-dimensional dynamical equations for the chromofield evolution are solved numerically. We find that there is regeneration of the longitudinal pressure due to unstable plasma modes; nevertheless, the system seems to maintain a high-degree of momentum-space anisotropy. Despite this anisotropy, we find that there is rapid longitudinal thermalization of the plasma due to the non-linear mode couplings inherent in the unstable evolution.

preprint2013arXiv

Thermalization and isotropization in heavy-ion collisions

I review our current understanding of the processes driving the thermalization and isotropization of the quark gluon plasma (QGP) created in ultrarelativistic heavy ion collisions (URHICs). I begin by discussing the phenomenological evidence in favor of the creation of a thermal but momentum-space anisotropic QGP in URHICs. I then discuss the degree of isotropization using viscous (dissipative) hydrodynamics, weak-coupling approaches to QGP dynamics, and strong-coupling approaches to QGP dynamics. Finally, I report on recent progress in the area of real-time non-abelian gauge field simulations and non-abelian Boltzmann-Vlasov-based hard-loop simulations.

preprint2013arXiv

Two-loop HTL pressure at finite temperature and chemical potential

We calculate the two-loop pressure of a plasma of quarks and gluons at finite temperature and chemical potential using the hard thermal loop perturbation theory (HTLpt) reorganization of finite temperature/density quantum chromodynamics. The computation utilizes a high temperature expansion through fourth order in the ratio of the chemical potential to temperature. This allows us to reliably access the region of high temperature and small chemical potential. We compare our final result for the two-loop HTLpt pressure at finite temperature and chemical potential with perturbative quantum chromodynamics (QCD) calculations and available lattice QCD results.

preprint2012arXiv

Boost-Invariant (2+1)-dimensional Anisotropic Hydrodynamics

We present results of the application of the anisotropic hydrodynamics (aHydro) framework to (2+1)-dimensional boost invariant systems. The necessary aHydro dynamical equations are derived by taking moments of the Boltzmann equation using a momentum-space anisotropic one-particle distribution function. We present a derivation of the necessary equations and then proceed to numerical solutions of the resulting partial differential equations using both realistic smooth Glauber initial conditions and fluctuating Monte-Carlo Glauber initial conditions. For this purpose we have developed two numerical implementations: one which is based on straightforward integration of the resulting partial differential equations supplemented by a two-dimensional weighted Lax-Friedrichs smoothing in the case of fluctuating initial conditions; and another that is based on the application of the Kurganov-Tadmor central scheme. For our final results we compute the collective flow of the matter via the lab-frame energy-momentum tensor eccentricity as a function of the assumed shear viscosity to entropy ratio, proper time, and impact parameter.

preprint2012arXiv

Bottomonia in the Quark Gluon Plasma

I review recent calculations of the suppression of bottomonium states in heavy ion collisions. A non-relativistic potential is used which is complex valued. This allows one to extract the binding energies and decay widths of the ground and excited states of bottomonium as a function of the typical plasma particle momentum and momentum-space anisotropy. The decay widths determined are used as input and integrated over space-time taking into account the dynamical evolution of the typical particle momentum and momentum-space anisotropy. The suppression of Upsilon(1s), Upsilon(2s), Upsilon(3s), chi_b1, and chi_b2 is obtained as a function of centrality, rapidity, and transverse momentum. The obtained results are compared with data from the STAR and CMS collaborations.

preprint2012arXiv

Chromoelectric oscillations in a dynamically evolving anisotropic background

We study the oscillations of a uniform longitudinal chromoelectric field in a dynamically-evolving momentum-space anisotropic background in the weak field limit. Evolution equations for the background are derived by taking moments of the Boltzmann equation in two cases: (i) a fixed relaxation time and (ii) a relaxation time that is proportional to the local inverse transverse momentum scale of the plasma. The second case allows us to reproduce 2nd-order viscous hydrodynamical dynamics in the limit of small shear viscosity to entropy ratio. We then linearize the Boltzmann-Vlasov equation in a dynamically-evolving background and obtain an integro-differential evolution equation for the chromoelectric field. We present numerical solutions to this integro-differential equation for a variety of different initial conditions and shear viscosity to entropy density ratios. The dynamical equations obtained are novel in that they include a non-trivial time-dependent momentum-space anisotropic background and the effect of collisional damping for the first time.

preprint2012arXiv

Highly anisotropic dissipative hydrodynamics

The quark gluon plasma generated in ultrarelativistic heavy ion collisions may possess sizable momentum-space anisotropies that cause the longitudinal and transverse pressures in the local rest frame to be significantly different. We review recent attempts to derive a dynamical framework that can reliably describe systems that possess a high degree of momentum-space anisotropy. The dynamical framework that has been developed can describe the evolution of the quark gluon plasma ranging from the longitudinal free-streaming limit to the ideal hydrodynamical limit.

preprint2012arXiv

Instabilities of an anisotropically expanding non-Abelian plasma: 3D+3V discretized hard-loop simulations

We study the (3+1)-dimensional evolution of non-Abelian plasma instabilities in the presence of a longitudinally expanding background of hard particles using the discretized hard loop framework. The free streaming background dynamically generates a momentum-space anisotropic distribution which is unstable to the rapid growth of chromomagnetic and chromoelectric fields. These fields produce longitudinal pressure that works to isotropize the system. Extrapolating our results to energies probed in ultrarelativistic heavy-ion collisions we find, however, that a pressure anisotropy persists for a few fm/c. In addition, on time scales relevant to heavy-ion collisions we observe continued growth of plasma instabilities in the strongly non-Abelian regime. Finally, we find that the longitudinal energy spectrum is well-described by a Boltzmann distribution with increasing temperature at intermediate time scales.

preprint2012arXiv

Thermal Bottomonium Suppression at RHIC and LHC

In this paper we consider the suppression of bottomonium states in ultrarelativistic heavy ion collisions. We compute the suppression as a function of centrality, rapidity, and transverse momentum for the states Upsilon(1s), Upsilon(2s), Upsilon(3s), chi_b1, and chi_b2. Using this information, we then compute the inclusive Upsilon(1s) suppression as a function of centrality, rapidity, and transverse momentum including feed down effects. Calculations are performed for both RHIC sqrt(s_NN)=200 GeV Au-Au collisions and LHC sqrt(s_NN)=2.76 TeV Pb-Pb collisions. From the comparison of our theoretical results with data available from the STAR and CMS Collaborations we are able to constrain the shear viscosity to entropy ratio to be in the range 0.08 < eta/S < 0.24. Our results are consistent with the creation of a high temperature quark-gluon plasma at both RHIC and LHC collision energies.

preprint2011arXiv

Non-Abelian plasma instabilities: SU(3) vs. SU(2)

We present the first 3+1 dimensional simulations of non-Abelian plasma instabilities in gauge-covariant Boltzmann-Vlasov equations for the QCD gauge group SU(3) as well as for SU(4) and SU(5). The real-time evolution of instabilities for a plasma with stationary momentum-space anisotropy is studied using a hard-loop effective theory that is discretized in the velocities of hard particles. We find that the numerically less expensive calculations using the group SU(2) essentially reproduce the nonperturbative dynamics of non-Abelian plasma instabilities with higher rank gauge groups provided the mass parameters of the corresponding hard-loop effective theories are the same. In particular we find very similar spectra for the turbulent cascade that forms in the strong-field regime, which is associated with an approximately linear growth of energy in collective fields. The magnitude of the linear growth however turns out to increase with the number of colors.

preprint2011arXiv

Quarkonium states in a complex-valued potential

We calculate quarkonium binding energies using a realistic complex-valued potential for both an isotropic and anisotropic quark-gluon plasma. We determine the disassociation temperatures of the ground and first excited states considering both the real and imaginary parts of the binding energy. We show that the effect of momentum-space anisotropy is smaller on the imaginary part of the binding energy than on the real part of the binding energy. In the case that one assumes an isotropic plasma, we find disassociation temperatures for the J/psi, Upsilon and chi_b of 1.6 T_c, 2.8 T_c, and 1.5 T_c, respectively. We find that a finite oblate momentum-space anisotropy increases the disassociation temperature for all states considered and results in a splitting of the p-wave states associated with the chi_b first excited state of bottomonium.

preprint2011arXiv

The QCD trace anomaly

In this brief report we compare the predictions of a recent next-to-next-to-leading order hard-thermal-loop perturbation theory (HTLpt) calculation of the QCD trace anomaly to available lattice data. We focus on the trace anomaly scaled by T^2 in two cases: N_f=0 and N_f=3. When using the canonical value of mu = 2 pi T for the renormalization scale, we find that for Yang-Mills theory (N_f=0) agreement between HTLpt and lattice data for the T^2-scaled trace anomaly begins at temperatures on the order of 8 T_c while when including quarks (N_f=3) agreement begins already at temperatures above 2 T_c. In both cases we find that at very high temperatures the T^2-scaled trace anomaly increases with temperature in accordance with the predictions of HTLpt.

preprint2011arXiv

Thermal Upsilon(1s) and chi_b1 suppression in sqrt(s_NN)=2.76 TeV Pb-Pb collisions at the LHC

I compute the thermal suppression of the Upsilon(1s) and chi_b1 states in sqrt(s_NN)=2.76 TeV Pb-Pb collisions. Using the suppression of each of these states I estimate the total R_AA for the Upsilon(1s) state as a function of centrality, rapidity, and transverse momentum. I find less suppression of the chi_b1 state than would be traditionally assumed; however, my final results for the total Upsilon(1s) suppression are in good agreement with recent preliminary CMS data.

preprint2011arXiv

Three-loop HTL QCD thermodynamics

The hard-thermal-loop perturbation theory (HTLpt) framework is used to calculate the thermodynamic functions of a quark-gluon plasma to three-loop order. This is the highest order accessible by finite temperature perturbation theory applied to a non-Abelian gauge theory before the high-temperature infrared catastrophe. All ultraviolet divergences are eliminated by renormalization of the vacuum, the HTL mass parameters, and the strong coupling constant. After choosing a prescription for the mass parameters, the three-loop results for the pressure and trace anomaly are found to be in very good agreement with recent lattice data down to $T \sim 2-3\,T_c$, which are temperatures accessible by current and forthcoming heavy-ion collision experiments.

preprint2010arXiv

A Parallel Algorithm for Solving the 3d Schrodinger Equation

We describe a parallel algorithm for solving the time-independent 3d Schrodinger equation using the finite difference time domain (FDTD) method. We introduce an optimized parallelization scheme that reduces communication overhead between computational nodes. We demonstrate that the compute time, t, scales inversely with the number of computational nodes as t ~ N_nodes^(-0.95 +/- 0.04). This makes it possible to solve the 3d Schrodinger equation on extremely large spatial lattices using a small computing cluster. In addition, we present a new method for precisely determining the energy eigenvalues and wavefunctions of quantum states based on a symmetry constraint on the FDTD initial condition. Finally, we discuss the usage of multi-resolution techniques in order to speed up convergence on extremely large lattices.

preprint2010arXiv

Dissipative Dynamics of Highly Anisotropic Systems

In this paper we present a method to improve the description of 0+1 dimensional boost invariant dissipative dynamics in the presence of large momentum-space anisotropies. We do this by reorganizing the canonical hydrodynamic expansion of the distribution function around a momentum-space anisotropic ansatz rather than an isotropic equilibrium one. At leading order the result obtained is two coupled ordinary differential equations for the momentum-space anisotropy and typical momentum of the degrees of freedom. We show that this framework can reproduce both the ideal hydrodynamic and free streaming limits. Additionally, we demonstrate that when linearized the differential equations reduce to 2nd order Israel-Stewart viscous hydrodynamics. Finally, we make quantitative comparisons of the evolution of the pressure anisotropy within our approach and 2nd order viscous hydrodynamics in both the strong and weak coupling limits.

preprint2010arXiv

Hard-thermal-loop QCD Thermodynamics

Naively resummed perturbative approximations to the thermodynamic functions of QCD do not converge at phenomenologically relevant temperatures. Here we review recent results of a three-loop hard-thermal-loop perturbation theory calculation of the thermodynamic functions of a quark-gluon plasma for general N_c and N_f. We show comparisons of our recent results with lattice data from both the hotQCD and Wuppertal-Budapest groups. We demonstrate that the three-loop hard-thermal-loop perturbation result for QCD thermodynamics agrees with lattice data down to temperatures T ~ 2 T_c.

preprint2010arXiv

Matching pre-equilibrium dynamics and viscous hydrodynamics

We demonstrate how to match pre-equilibrium dynamics of a 0+1 dimensional quark gluon plasma to 2nd-order viscous hydrodynamical evolution. The matching allows us to specify the initial values of the energy density and shear tensor at the initial time of hydrodynamical evolution as a function of the lifetime of the pre-equilibrium period. We compare two models for the pre-equilibrium quark-gluon plasma, longitudinal free streaming and collisionally-broadened longitudinal expansion, and present analytic formulas which can be used to fix the necessary components of the energy-momentum tensor. The resulting dynamical models can be used to assess the effect of pre-equilibrium dynamics on quark-gluon plasma observables. Additionally, we investigate the dependence of entropy production on pre-equilibrium dynamics and discuss the limitations of the standard definitions of the non-equilibrium entropy.

preprint2010arXiv

NNLO hard-thermal-loop thermodynamics for QCD

We calculate the thermodynamic functions of a quark-gluon plasma for general N_c and N_f to three-loop order using hard-thermal-loop perturbation theory. At this order, all the ultraviolet divergences can be absorbed into renormalizations of the vacuum, the HTL mass parameters, and the strong coupling constant.We show that at three loops, the results for the pressure and trace anomaly are in very good agreement with recent lattice data down to temperatures T~2T_c.

preprint2010arXiv

Three-loop HTL gluon thermodynamics at intermediate coupling

We calculate the thermodynamic functions of pure-glue QCD to three-loop order using the hard-thermal-loop perturbation theory (HTLpt) reorganization of finite temperature quantum field theory. We show that at three-loop order hard-thermal-loop perturbation theory is compatible with lattice results for the pressure, energy density, and entropy down to temperatures $T\simeq3\;T_c$. Our results suggest that HTLpt provides a systematic framework that can used to calculate static and dynamic quantities for temperatures relevant at LHC.

preprint2009arXiv

Hard-thermal-loop QED thermodynamics

The weak-coupling expansion for thermodynamic quantities in thermal field theories is poorly convergent unless the coupling constant is tiny. We discuss the calculation of the free energy for a hot gas of electrons and photons to three-loop order using hard-thermal-loop perturbation theory (HTLpt). We show that the hard-thermal-loop perturbation reorganization improves the convergence of the successive approximations to the QED free energy at large coupling, e ~ 2. The reorganization is gauge invariant by construction, and due to the cancellations among various contributions, we obtain a completely analytic result for the resummed thermodynamic potential at three loops.

preprint2009arXiv

Jet energy loss in the quark-gluon plasma by stream instabilities

We study the evolution of the plasma instabilities induced by two jets of particles propagating in opposite directions and crossing a thermally equilibrated non-Abelian plasma. In order to simplify the analysis we assume that the two jets of partons can be described with uniform distribution functions in coordinate space and by Gaussian distribution functions in momentum space. We find that while crossing the quark-gluon plasma, the jets of particles excite unstable chromomagnetic and chromoelectric modes. These fields interact with the particles (or hard modes) of the plasma inducing the production of currents; thus, the energy lost by the jets is absorbed by both the gauge fields and the hard modes of the plasma. We compare the outcome of the numerical simulations with the analytical calculation performed assuming that the jets of particles can be described by a tsunami-like distribution function. We find qualitative and semi-quantitative agreement between the results obtained with the two methods.

preprint2009arXiv

QED Thermodynamics at Intermediate Coupling

We discuss reorganizing finite temperature perturbation theory using hard-thermal-loop (HTL) perturbation theory in order to improve the convergence of successive perturbative approximations to the free energy of a gauge theory. We briefly review the history of the technique and present new results for the three-loop HTL-improved approximation for the free energy of QED. We show that the hard-thermal-loop perturbation reorganization improves the convergence of the successive approximations to the QED free energy at intermediate coupling, e ~ 2. The reorganization is gauge invariant by construction, and due to cancellation among various contributions, one can obtain a completely analytic result for the resummed thermodynamic potential at three loops.

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.

preprint2006arXiv

The chromo-Weibel instability

I discuss the physics of non-Abelian plasmas which are locally anisotropic in momentum space. Such momentum-space anisotropies are generated by the rapid longitudinal expansion of the matter created in the first 1 fm/c of an ultrarelativistic heavy ion collision. In contrast to locally isotropic plasmas anisotropic plasmas have a spectrum of soft unstable modes which are characterized by exponential growth of transverse chromo-magnetic/-electric fields at short times. This instability is the QCD analogue of the Weibel instability of QED. Parametrically the chromo-Weibel instability provides the fastest method for generation of soft background fields and dominates the short-time dynamics of the system. The existence of the chromo-Weibel instability has been proven using diagrammatic methods, transport theory, and numerical solution of classical Yang-Mills fields. I review the results obtained from each of these methods and discuss the numerical techniques which are being used to determine the late-time behavior of plasmas subject to a chromo-Weibel instability.

preprint1999arXiv

Hard-thermal-loop Resummation of the Thermodynamics of a Hot Gluon Plasma

We calculate the thermodynamic functions of a hot gluon plasma to leading order in hard-thermal-loop (HTL) perturbation theory. Effects associated with screening, gluon quasiparticles, and Landau damping are resummed to all orders. The ultraviolet divergences generated by the HTL propagator corrections can be cancelled by a counterterm that depends on the thermal gluon mass parameter. The HTL thermodynamic functions are compared to those from lattice gauge theory calculations and from quasiparticle models. For reasonable values of the HTL parameters, the deviations from lattice results for T > 2 T_c have the correct sign and roughly the correct magnitude to be accounted for by next-to-leading order corrections in HTL perturbation theory.