Source author record

Guy D. Moore

Guy D. Moore 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

71works
13topics
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

71 published item(s)

preprint2022arXiv

Lattice QCD noise reduction for bosonic correlators through blocking

We propose a method to substantially improve the signal-to-noise ratio of lattice correlation functions for bosonic operators or other operator combinations with disconnected contributions. The technique is applicable for correlations between operators on two planes (zero momentum correlators) when the dimension of the plane is larger than the separation between the two planes which are correlated. In this case, the correlation arises primarily from points whose in-plane coordinates are close, but noise arises from all pairs of points. By breaking each plane into bins and computing bin-bin correlations, it is possible to capture these short-distance correlators exactly while replacing (small) correlators at large spatial extent with a fit, with smaller uncertainty than the data. The cost is only marginally larger than averaging each plane before correlating, but the improvement in signal-to-noise can be substantial. We test the method on correlators of the gradient-flowed topological charge density and squared field strength, finding noise reductions by a factor of $\sim$ 3$-$7 compared to the conventional approach on the same ensemble of configurations.

preprint2022arXiv

Spacelike thermal correlators are almost time-independent

We show that, in relativistic field theories, the thermal correlation function of N bosonic operators <O_1(x_1,t_1) O_2(x_2,t_2) .. O_N(x_N,t_N)> at sufficiently spacelike-separated points shows exponentially weak dependence on the time variables t_1,...,t_N, when the space separations are held fixed. For classical thermal field theory, the time dependence vanishes when all points are spacelike separated

preprint2022arXiv

Spectral reconstruction details of a gradient-flowed color-electric correlator

In a recently published work we provide a proof-of-concept of a novel method to extract the heavy quark momentum diffusion coefficient from color-electric correlators on the lattice using gradient flow. The transport coefficient can be found in the infrared limit of the corresponding spectral function which is reconstructed through perturbative model fits of the correlator data. In this proceedings report we want to give more detailed insights into the systematic uncertainties of this procedure and compare our results with other studies.

preprint2021arXiv

Heavy quark momentum diffusion from the lattice using gradient flow

We apply the gradient flow on a color-electric two-point function that encodes the heavy quark momentum diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at $1.5\,T_c$. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations of this correlation function under Wilson flow are used to enhance the extrapolations of the non-perturbative lattice correlator. The final estimate for the continuum correlator at zero flow time largely agrees with one obtained from a previous study using the multi-level algorithm. We perform a spectral reconstruction based on perturbative model fits to estimate the heavy quark momentum diffusion coefficient. The approach we present here yields high-precision data for the correlator and is also applicable for actions with dynamical fermions.

preprint2020arXiv

Global cosmic string networks as a function of tension

We investigate the properties of global cosmic string networks as a function of the ratio of string tension to Goldstone-field coupling, and as a function of the Hubble damping strength. Our results show unambiguously that the string density is sensitive to this ratio. We also find that existing semi-analytical (one-scale) models must be missing some important aspect of the network dynamics. Our results point the way towards improving such models.

preprint2020arXiv

Improved Reweighting for QCD Topology at High Temperature

In a previous paper \cite{Jahn:2018dke} we presented a methodology for computing the topological susceptibility of QCD at temperatures where it is small and standard methods fail. Here we improve on this methodology by removing two barriers to the reweighting method's moving between topological sectors. We present high-statistics, continuum-extrapolated results for the susceptibility of pure-glue QCD up to $7 T_c$. We show that the susceptibility varies with temperature as $T^{-6.7\pm 0.3}$ between $T=2.5 T_c$ and $T=7 T_c$, in good agreement with expectations based on the dilute instanton gas approximation.

preprint2020arXiv

Thermal Heavy Quark Self-Energy from Euclidean Correlators

Brambilla, Escobedo, Soto, and Vairo have derived an effective description of quarkonium with two parameters; a momentum diffusion term and a real self-energy term. We point out that there is a similar real self-energy term for a single open heavy flavor and that it can be expressed directly in terms of Euclidean electric field correlators along a Polyakov line. This quantity can be directly studied on the lattice without the need for analytical continuation. We show that Minkowski-space calculations of this correlator correspond with the known NLO Euclidean value of the relevant electric field two-point function and that it differs from the real self-energy term for quarkonium.

preprint2016arXiv

Axion String Dynamics I: 2+1D

If the axion exists and if the initial axion field value is uncorrelated at causally disconnected points, then it should be possible to predict the efficiency of cosmological axion production, relating the axionic dark matter density to the axion mass. The main obstacle to making this prediction is correctly treating the axion string cores. We develop a new algorithm for treating the axionic string cores correctly in 2+1 dimensions. When the axionic string cores are given their full physical string tension, axion production is about twice as efficient as in previous simulations. We argue that the string network in 2+1 dimensions should behave very differently than in 3+1 dimensions, so this result cannot be simply carried over to the physical case. We outline how to extend our method to 3+1D axion string dynamics.

preprint2016arXiv

Challenges for models with composite states

Composite states of electrically charged and QCD-colored hyperquarks (HQs) in a confining SU(N_HC) hypercolor gauge sector are a plausible extension of the standard model at the TeV scale, and have been widely considered as an explanation for the tentative LHC diphoton excess. Additional new physics is required to avoid a stable charged hyperbaryon in such theories. We classify renormalizable models allowing the decay of this unwanted relic directly into standard model states, showing that they are significantly restricted if the new scalar states needed for UV completion are at the TeV scale. Alternatively, if hyperbaryon number is conserved, the charged relic can decay into a neutral hyperbaryon. Such theories are strongly constrained by direct detection, if the neutral constituent hyperquark carries color or weak isospin, and by LHC searches for leptoquarks if it is a color singlet. We show that the neutral hyperbaryon can have the observed relic abundance if the confinement scale and the hyperquark mass are above TeV scale, even in the absence of any hyperbaryon asymmetry.

preprint2016arXiv

Condensates in Relativistic Scalar Theories

Scalar field theory with large infrared initial occupancy develops very large deep-infrared occupancy, which locally resembles a Bose-Einstein condensate. We study the structure and spatial coherence of this condensate. The O(N) symmetric theory with N>1 is qualitatively different than N=1. We explain the thermodynamical reason why, for N>1, the condensate locally carries nearly maximal conserved charge density. We also show how this property impedes the condensate's decay, and we show that it prevents the condensate from ever becoming fully spatially homogeneous. For N <= 4 the condensate can carry topological defects, but these do not appear to control the large-k tail in its power spectrum, which is the same for N=8 where there are no topological defects.

preprint2016arXiv

Intercommutation of U(1) global cosmic strings

Global strings (those which couple to Goldstone modes) may play a role in cosmology. In particular, if the QCD axion exists, axionic strings may control the efficiency of axionic dark matter abundance. The string network dynamics depend on the string intercommutation efficiency (whether strings re-connect when they cross). We point out that the velocity and angle in a collision between global strings "renormalize" between the network scale and the microscopic scale, and that this plays a significant role in their intercommutation dynamics. We also point out a subtlety in treating intercommutation of very nearly antiparallel strings numerically. We find that the global strings of a O(2)-breaking scalar theory do intercommute for all physically relevant angles and velocities.

preprint2016arXiv

Jet-Medium Interactions at NLO in a Weakly-Coupled Quark-Gluon Plasma

We present an extension to next-to-leading order in the strong coupling constant $g$ of the AMY effective kinetic approach to the energy loss of high momentum particles in the quark-gluon plasma. At leading order, the transport of jet-like particles is determined by elastic scattering with the thermal constituents, and by inelastic collinear splittings induced by the medium. We reorganize this description into collinear splittings, high-momentum-transfer scatterings, drag and diffusion, and particle conversions (momentum-preserving identity-changing processes). We show that this reorganized description remains valid to NLO in $g$, and compute the appropriate modifications of the drag, diffusion, particle conversion, and inelastic splitting coefficients. In addition, a new kinematic regime opens at NLO for wider-angle collinear bremsstrahlung. These semi-collinear emissions smoothly interpolate between the leading order high-momentum-transfer scatterings and collinear splittings. To organize the calculation, we introduce a set of Wilson line operators on the light-cone which determine the diffusion and identity changing coefficients, and we show how to evaluate these operators at NLO.

preprint2015arXiv

Axion dark matter: strings and their cores

Axions constitute a well-motivated dark matter candidate, and if PQ symmetry breaking occurred after inflation, it should be possible to make a clean prediction for the relation between the axion mass and the axion dark matter density. We show that axion (or other global) string networks in 3D have a network density that depends logarithmically on the string separation-to-core ratio. This logarithm would be about 10 times larger in axion cosmology than what we can achieve in numerical simulations. We simulate axion production in the early Universe, finding that, for the separation-to-core ratios we can achieve, the changing density of the network has little impact on the axion production efficiency.

preprint2014arXiv

3.5 keV X-rays as the "21 cm line" of dark atoms, and a link to light sterile neutrinos

The recently discovered 3.5 keV X-ray line from extragalactic sources may be evidence of dark matter scatterings or decays. We show that dark atoms can be the source of the emission, through their hyperfine transitions, which would be the analog of 21 cm radiation from a dark sector. We identify two families of dark atom models that match the X-ray observations and are consistent with other constraints. In the first, the hyperfine excited state is long-lived compared to the age of the universe, and the dark atom mass is relatively unconstrained; dark atoms could be strongly self-interacting in this case. In the second, the excited state is short-lived and viable models are parameterized by the value of the dark proton-to-electron mass ratio $R$: for $R = 10^2-10^4$, the dark atom mass is predicted be in the range $350-1300$ GeV, with fine structure constant $α'\cong 0.1-0.6$. In either class of models, the dark photon must be massive with $m_{γ'} \gtrsim$ 1 MeV and decay into $e^+ e^-$. Evidence for the model could come from direct detection of the dark atoms. In a natural extension of this framework, the dark photon could decay predominantly into invisible particles, for example $\sim 0.5$ eV sterile neutrinos, explaining the extra radiation degree of freedom recently suggested by data from BICEP2, while remaining compatible with BBN.

preprint2014arXiv

Low Mass Thermal Dilepton Production at NLO in a Weakly Coupled Quark-Gluon Plasma

We present a computation, within weakly-coupled thermal QCD, of the production rate of low invariant mass ($M^2 \sim g^2 T^2$) dileptons, at next-to-leading order (NLO) in the coupling (which is $O(g^3 e^2 T^2)$). This involves extending the NLO calculation of the photon rate which we recently presented to the case of small nonzero photon invariant mass. Numerical results are discussed and tabulated forms and code are provided for inclusion in hydrodynamical models. We find that NLO corrections can increase the dilepton rate by up to 30-40% relative to leading order. We find that the electromagnetic response of the plasma for real photons and for small invariant mass but high energy dilepton pairs (e.g., $M^2 < (300\:\mathrm{MeV})^2$ but $p_T > 1 \: \mathrm{GeV}$) are close enough that dilepton pair measurements really can serve as Ersatz photon measurements. We also present a matching a la Ghisoiu and Laine between our results and results at larger invariant masses.

preprint2014arXiv

nPI Resummation in 3D SU(N) Higgs Theory

We test the utility of the nPI formalism for solving nonperturbative dynamics of gauge theories by applying it to study the phase diagram of SU(N) Higgs theory in 3 Euclidean spacetime dimensions. Solutions reveal standard signatures of a first order phase transition with a critical endpoint leading to a crossover regime, in qualitative agreement with lattice studies. The location of the critical endpoint, x sim 0.14 for SU(2) with a fundamental Higgs, is in rough but not tight quantitative agreement with the lattice. We end by commenting on the overall effectiveness and limitations of an nPI effective action based study. In particular, we have been unable to find an nPI gauge-fixing procedure which can simultaneously display the right phase structure and correctly handle the large-VEV Higgs region. We explain why doing so appears to be a serious challenge.

preprint2014arXiv

Renormalization of Null Wilson Lines in EQCD

Radiation and energy loss of a light, high-energy parton in a perturbative Quark-Gluon Plasma is controlled by transverse momentum exchange. The troublesome infrared contributions to transverse momentum exchange can be computed on the lattice using dimensional reduction to EQCD. However a novel extended operator, the Null Wilson Line of EQCD, is involved. We compute the renormalization properties of this object's lattice implementation to next-to-leading order, which should facilitate its efficient calculation on the lattice.

preprint2014arXiv

Scattering properties of dark atoms and molecules

There has been renewed interest in the possibility that dark matter exists in the form of atoms, analogous to those of the visible world. An important input for understanding the cosmological consequences of dark atoms is their self-scattering. Making use of results from atomic physics for the potentials between hydrogen atoms, we compute the low-energy elastic scattering cross sections for dark atoms. We find an intricate dependence upon the ratio of the dark proton to electron mass, allowing for the possibility to "design" low-energy features in the cross section. Dependences upon other parameters, namely the gauge coupling and reduced mass, scale out of the problem by using atomic units. We derive constraints on the parameter space of dark atoms by demanding that their scattering cross section does not exceed bounds from dark matter halo shapes. We discuss the formation of molecular dark hydrogen in the universe, and determine the analogous constraints on the model when the dark matter is predominantly in molecular form.

preprint2014arXiv

UV Cascade in Classical Yang-Mills via Kinetic Theory

We show that classical Yang-Mills theory with statistically homogeneous and isotropic initial conditions has a kinetic description and approaches a scaling solution at late times. We find the scaling solution by explicitly solving the Boltzmann equations, including all dominant processes (elastic and number-changing). Above a scale $p_{max} \propto t^{1/7}$ the occupancy falls exponentially in $p$. For asymptotically late times and sufficiently small momenta the occupancy scales as $f(p)\propto 1/p$, but this behavior sets in only at very late time scales. We find quantitative agreement of our results with lattice simulations, for times and momenta within the range of validity of kinetic theory.

preprint2013arXiv

3-loop 3PI effective action for 3D SU(3) QCD

The 3PI method is a technique to resum an infinite class of diagrams, which may be useful in studying nonperturbative thermodynamics and dynamics in quantum field theory. But it has never been successfully applied to gauge theories, where there are serious questions about gauge invariance breaking. We show how to perform the 3PI resummation of QCD in 3 Euclidean spacetime dimensions, a warmup problem to the 4 or 3+1 dimensional case. We present the complete details of the technical problems and how they are overcome. We postpone a comparison of gauge invariant correlation functions with their lattice-determined counterparts to a future publication.

preprint2013arXiv

Next-to-leading order thermal photon production in a weakly coupled quark-gluon plasma

We compute the next-to-leading order O(g) correction to the thermal photon production rate in a QCD plasma. The NLO contributions can be expressed in terms of gauge invariant condensates on the light cone, which are amenable to novel sum rules and Euclidean techniques. We expect these technologies to be generalizable to other NLO calculations. For the phenomenologically interesting value of alpha_s=0.3, the NLO correction represents a 20% increase and has a functional form similar to the LO result.

preprint2012arXiv

Composite magnetic dark matter and the 130 GeV line

We propose an economical model to explain the apparent 130 GeV gamma ray peak, found in the Fermi/LAT data, in terms of dark matter annihilation through a dipole moment interaction. The annihilating dark matter particles represent a subdominant component, with mass density 7-17% of the total DM density; and they only annihilate into gamma gamma, gamma Z, and ZZ, through a magnetic (or electric) dipole moment. Annihilation into other standard model particles is suppressed, due to a mass splitting in the magnetic dipole case, or to p-wave scattering in the electric dipole case. In either case, the observed signal requires a dipole moment of strength mu ~ 2/TeV. We argue that composite models are the preferred means of generating such a large dipole moment, and that the magnetic case is more natural than the electric one. We present a simple model involving a scalar and fermionic techniquark of a confining SU(2) gauge symmetry. We point out some generic challenges for getting such a model to work. The new physics leading to a sufficiently large dipole moment is below the TeV scale, indicating that the magnetic moment is not a valid effective operator for LHC physics, and that production of the strongly interacting constituents, followed by techni-hadronization, is a more likely signature than monophoton events. In particular, 4-photon events from the decays of bound state pairs are predicted.

preprint2012arXiv

UV Cascade in Classical Yang-Mills Theory

We study the real-time behavior of classical Yang-Mills theory under initial conditions with nonperturbatively large, infrared field amplitudes. Our lattice study confirms the cascade of energy towards higher momenta and lower occupancy, which occurs via a scaling solution $f[p,t_1] = (t_0/t_1)^{4/7}\, f[p (t_0/t_1)^{1/7},t_0]$. Above a characteristic scale p_{max}, f falls exponentially; below p_{max}, $f[p] \propto p^{-4/3}$. We find no evidence for different infrared exponents or for infrared occupancies in excess of those described by this scaling solution. We also investigate what the fate of large occupancies would be, both in the electric and the magnetic sector.

preprint2011arXiv

Bjorken Flow, Plasma Instabilities, and Thermalization

At asymptotically high energies, thermalization in heavy ion collisions can be described via weak-coupling QCD. We present a complete treatment of how thermalization proceeds, at the parametric weak-coupling level. We show that plasma instabilities dominate the dynamics, from immediately after the collision until well after the plasma becomes nearly in equilibrium. Initially they drive the system close to isotropy, but Bjorken expansion and increasing diluteness makes the system again become more anisotropic. At time τ~ α^(-12/5) Q^(-1) the dynamics become dominated by a nearly-thermal bath; and at time τ~ α^(-5/2) Q^(-1)$ the bath comes to dominate the energy density, completing thermalization. After this time there is a nearly isotropic and thermal Quark-Gluon Plasma.

preprint2011arXiv

Chain Inflation Reconsidered

We investigate density perturbations in the chain inflation scenario, where the inflaton undergoes successive tunneling transitions along one field direction. First we show that when the bubble walls associated with such a phase transition meet, they induce the next phase transition and continue to propagate as phase interfaces. Then we present an analytical calculation of the density fluctuations and an estimate of non-Gaussianities such as f_{NL}, which we find to be small (of order 1). To get the right amplitude for the power spectrum, there have to be 10^8 phase transitions per Hubble time, a significant model building challenge. We find that working models of chain inflation must be rather strongly coupled, and thus have a very limited range of validity as effective field theories. We discuss generalizations to the multiple field case, the curvaton scenario, and noncanonical kinetic terms.

preprint2011arXiv

Exploring the Phase Diagram with Taylor Series: Epic Voyage or Just Another Bad Trip

It has been suggested in the literature that it may be possible to locate the QCD critical end point using the Taylor series of thermodynamic variables about the $μ=0$ axis. Since the phase transition at the critical end point is believed to be in the 3D Ising universality class, it would seem natural to test this method with the Ising Model, for which the answer is already known. The finding is that it is in fact possible to pinpoint the location of the Ising critical point using Taylor coefficients.

preprint2011arXiv

Kubo Formulae for Second-Order Hydrodynamic Coefficients

At second order in gradients, conformal relativistic hydrodynamics depends on the viscosity eta and on five additional "second-order" hydrodynamical coefficients tauPi, kappa, lambda1, lambda2, and lambda3. We derive Kubo relations for these coefficients, relating them to equilibrium, fully retarded 3-point correlation functions of the stress tensor. We show that the coefficient lambda3 can be evaluated directly by Euclidean means and does not in general vanish.

preprint2011arXiv

The Bulk Viscosity of a Pion Gas

We compute the bulk viscosity of a gas of pions at temperatures below the QCD crossover temperature, for the physical value of pion mass, to lowest order in chiral perturbation theory. Bulk viscosity is controlled by number-changing processes which become exponentially slow at low temperatures when the pions become exponentially dilute, leading to an exponentially large bulk viscosity zeta ~ (F_0^8/m_π^5) exp(2m_π/T), where F_0 = 93 MeV is the pion decay constant.

preprint2011arXiv

The stickiness of sound: An absolute lower limit on viscosity and the breakdown of second order relativistic hydrodynamics

Hydrodynamics predicts long-lived sound and shear waves. Thermal fluctuations in these waves can lead to the diffusion of momentum density, contributing to the shear viscosity and other transport coefficients. Within viscous hydrodynamics in 3+1 dimensions, this leads to a positive contribution to the shear viscosity, which is finite but inversely proportional to the microscopic shear viscosity. Therefore the effective infrared viscosity is bounded from below. The contribution to the second-order transport coefficient $τ_π$ is divergent, which means that second-order relativistic viscous hydrodynamics is inconsistent below some frequency scale. We estimate the importance of each effect for the Quark-Gluon Plasma, finding them to be minor if $η/s = 0.16$ but important if $η/s = 0.08$.

preprint2011arXiv

Thermalization in Weakly Coupled Nonabelian Plasmas

We investigate how relativistic, nonabelian plasmas approach equilibrium in a general context. Our treatment is entirely parametric and for small Yang-Mills coupling $α$. First we study isotropic systems with an initially nonequilibrium momentum distribution. We consider both the case of initially very high occupancy and initially very low occupancy. Then we consider systems which are anisotropic. We consider both weak anisotropy and large anisotropy, and allow the occupancy to be parametrically large or small. Writing the typical momentum of an initial excitation as Q and the final temperature as T, full equilibration occurs in a time t ~ α^{-2}/T for T > Q, and t ~ α^{-2} Q^{1/2} T^{-3/2} for T < Q, unless the initial system is sufficiently anisotropic and T > α^{2/3} Q, in which case equilibration occurs somewhat faster, t ~ α^{-13/7} Q^{5/7} T^{-12/7} (or α^{-2}/T if that is longer).

preprint2010arXiv

The Sphaleron Rate in SU(N) Gauge Theory

The sphaleron rate is defined as the diffusion constant for topological number NCS = int g^2 F Fdual/32 pi^2. It establishes the rate of equilibration of axial light quark number in QCD and is of interest both in electroweak baryogenesis and possibly in heavy ion collisions. We calculate the weak-coupling behavior of the SU(3) sphaleron rate, as well as making the most sensible extrapolation towards intermediate coupling which we can. We also study the behavior of the sphaleron rate at weak coupling at large Nc.

preprint2009arXiv

A way to estimate the heavy quark thermalization rate from the lattice

The thermalization rate of a heavy quark is related to its momentum diffusion coefficient. Starting from a Kubo relation and using the framework of the heavy quark effective theory, we argue that in the large-mass limit the momentum diffusion coefficient can be defined through a certain Euclidean correlation function, involving color-electric fields along a Polyakov loop. Furthermore, carrying out a perturbative computation, we demonstrate that the spectral function corresponding to this correlator is relatively flat at small frequencies. Therefore, unlike in the case of several other transport coefficients, for which the narrowness of the transport peak makes analytic continuation from Euclidean lattice data susceptible to severe systematic uncertainties, it appears that the determination of the heavy quark thermalization rate could be relatively well under control.

preprint2009arXiv

Can electroweak bubble walls run away?

In extensions of the Standard Model with SU(2) singlet scalar fields, there can be regions of parameter space for which the electroweak phase transition is first order already at the mean-field level of analysis. We show that in this case the phase interface (bubble wall) can become ultra-relativistic, with the relativistic gamma factor gamma = (1-v_{wall}^2)^{-1/2} growing linearly with the wall's propagation distance. We provide a simple criterion for determining whether the bubble wall "runs away" in this way or if gamma approaches a terminal value.

preprint2009arXiv

Heavy Quark Thermalization in Classical Lattice Gauge Theory: Lessons for Strongly-Coupled QCD

Thermalization of a heavy quark near rest is controlled by the correlator of two electric fields along a temporal Wilson line. We address this correlator within real-time, classical lattice Yang-Mills theory, and elaborate on the analogies that exist with the dynamics of hot QCD. In the weak-coupling limit, it can be shown analytically that the dynamics on the two sides are closely related to each other. For intermediate couplings, we carry out non-perturbative simulations within the classical theory, showing that the leading term in the weak-coupling expansion significantly underestimates the heavy quark thermalization rate. Our analytic and numerical results also yield a general understanding concerning the overall shape of the spectral function corresponding to the electric field correlator, which may be helpful in subsequent efforts to reconstruct it from Euclidean lattice Monte Carlo simulations.

preprint2009arXiv

Jet energy loss, photon production, and photon-hadron correlations at RHIC

Jet energy loss, photon production and photon-hadron correlations are studied together at high transverse momentum in relativistic heavy-ion collisions at RHIC energies. The modification of hard jets traversing a hot and dense nuclear medium is evaluated by consistently taking into account induced gluon radiation and elastic collisions. The production of high transverse momentum photons in Au+Au collisions at RHIC is calculated by incorporating a complete set of photon-production channels. Comparison with experimental photon production and photon-hadron correlation data is performed, using a (3+1)-dimensional relativistic hydrodynamic description of the thermalized medium created in these collisions. Our results demonstrate that the interaction between the hard jets and the soft medium is important for the study of photon production and of photon-hadron correlation at RHIC.

preprint2008arXiv

Bulk viscosity and spectral functions in QCD

We examine the behavior of the spectral function for the trace of the stress tensor in QCD in the two regimes where it is possible to make analytical progress; weak coupling, and close to a second order QCD phase transition. We determine the behavior of the bulk viscosity in each regime. We discuss the problem of analytic continuation of the (lattice) Euclidean correlation function to determine the spectral function. In each case the spectral function has a narrow peak at small frequency; its shape would be challenging to extract accurately from lattice data with error bars.

preprint2008arXiv

Heavy quark diffusion in QCD and N=4 SYM at next-to-leading order

We present the full details of a calculation at next-to-leading order of the momentum diffusion coefficient of a heavy quark in a hot, weakly coupled, QCD plasma. Corrections arise at O(g_s); physically they represent interference between overlapping scatterings, as well as soft, electric scale (p ~ gT) gauge field physics, which we treat using the hard thermal loop (HTL) effective theory. In 3-color, 3-flavor QCD, the momentum diffusion constant of a fundamental representation heavy quark at NLO is kappa = (16π/3) alpha_s^2 T^3 (log(1/g) + 0.07428 + 1.9026 g). We extend the computation to a heavy fundamental representation ``probe'' quark in large N_c, N=4 Super Yang-Mills theory, where the result is kappa^{SYM}= (lambda^2 T^3)(6π) (log(1/\sqrtλ) + 0.4304 + 0.8010 \sqrt{lambda}) (where lambda=g_s^2 N_c is the t'Hooft coupling). In the absence of some resummation technique, the convergence of perturbation theory is poor.

preprint2006arXiv

The Bulk Viscosity of High-Temperature QCD

We compute the bulk viscosity zeta of high-temperature QCD to leading order in powers of the running coupling alpha_s(T). We find that it is negligible compared to shear viscosity eta for any alpha_s that might reasonably be considered small. The physics of bulk viscosity in QCD is a bit different than in scalar phi^4 theory. In particular, unlike in scalar theory, we find that an old, crude estimate of zeta as 15 ((1/3)-v_s^2)^2 eta gives the correct order of magnitude, where v_s is the speed of sound. We also find that leading-log expansions of our result for zeta are not accurate except at very small coupling.

preprint2005arXiv

Apparent Thermalization due to Plasma Instabilities in Quark-Gluon Plasma

Hydrodynamical modeling of heavy ion collisions at RHIC suggests that the quark-gluon plasma (QGP) "thermalizes" in a remarkably short time scale, about 0.6 fm/c. We argue that this should be viewed as indicating fast isotropization, but not necessarily complete thermalization, of the non-equilibrium QGP. Non-Abelian plasma instabilities can drive local isotropization of an anisotropic QGP on a time scale which is faster than ordinary perturbative scattering processes. As a result, we argue that theoretical expectations based on weak coupling analysis are not necessarily in conflict with hydrodynamic modeling of the early part of RHIC collisions, provided one recognizes the key role of non-Abelian plasma instabilities.

preprint2005arXiv

Constraining the New Aether: Gravitational Cherenkov Radiation

We study the simplest concrete theory for spontaneous Lorentz violation, the ``New Aether Theory'' of Jacobson and Mattingly, which is a vector-tensor gravitational theory with a fixed-modulus condition on the vector field. We show that the observation of ultra-high energy cosmic rays (which implies the absence of energy loss via various Cherenkov type processes) places constraints on the parameters of this theory, which are much stronger than those previously found in the literature and are also stronger than the constraints generically arising when gravity displays sub-luminal propagation.

preprint2005arXiv

Energy Loss of Leading Hadrons and Direct Photon production in Evolving Quark-Gluon Plasma

We calculate the nuclear modification factor of neutral pions and the photon yield at high p_T in central Au-Au collisions at RHIC (\sqrt{s}=200 GeV) and Pb-Pb collisions at the LHC (\sqrt{s}=5500 GeV). A leading-order accurate treatment of jet energy loss in the medium has been convolved with a physical description of the initial spatial distribution of jets and a (1+1) dimensional expansion. We reproduce the nuclear modification factor of pion R_{AA} at RHIC, assuming an initial temperature T_i=370 MeV and a formation time τ_i=0.26 fm/c, corresponding to dN/dy=1260. The resulting suppression depends on the particle rapidity density dN/dy but weakly on the initial temperature. The jet energy loss treatment is also included in the calculation of high p_T photons. Photons coming from primordial hard N-N scattering are the dominant contribution at RHIC for p_T > 5 GeV, while at the LHC, the range 8<p_T<14 GeV is dominated by jet-photon conversion in the plasma.

preprint2005arXiv

QCD Plasma Instabilities: The Nonabelian Cascade

Magnetic plasma instabilities appear to play an important role in the early stages of quark-gluon plasma equilibration in the high energy (weak coupling) limit. Numerical studies of the growth of such instabilities from small seed fluctuations have found initial exponential growth in their energy, followed by linear growth once the associated color magnetic fields become so large that their non-abelian interaction are non-perturbative. In this paper, we use simulations to determine the nature of this linear energy growth. We find that the long-wavelength modes associated with the instability have ceased to grow, but that they cascade energy towards the ultraviolet in the form of plasmon excitations of ever increasing energy. We find a quasi-steady-state power-law distribution $f_k \propto k^{-ν}$ for this cascade, with spectral index $ν\simeq 2$.

preprint2005arXiv

The turbulent spectrum created by non-Abelian plasma instabilities

Recent numerical work on the fate of plasma instabilities in weakly-coupled non-Abelian gauge theory has shown the development of a cascade of energy from long to short wavelengths. This cascade has a steady-state spectrum, analogous to the Kolmogorov spectrum for turbulence in hydrodynamics or for energy cascades in other systems. In this paper, we theoretically analyze processes responsible for this cascade and find a steady-state spectrum f_k ~ k^-2, where f_k is the phase-space density of particles with momentum k. The exponent -2 is consistent with results from numerical simulations. We also discuss implications of the emerging picture of instability development on the "bottom-up" thermalization scenario for (extremely high energy) heavy ion collisions, emphasizing fundamental questions that remain to be answered.

preprint2004arXiv

QCD Plasma Instabilities and Bottom-Up Thermalization

We study the role of QCD plasma instabilities in non-equilibrium quark-gluon plasmas. First, we argue that such instabilities must drastically modify the "bottom-up" thermalization scenario for heavy-ion collisions. Second, we discuss conditions for the existence of instabilities in a more general context than previously treated in the QCD literature. We also give a thorough qualitative review of the origin of instabilities. We discuss some mechanisms whereby the growth of plasma instabilities saturates. Finally, we solve explicitly for instabilities and their growth rates for two extreme cases of anisotropic non-equilibrium plasmas that can be treated relatively simply and analytically: $f(\p) = F(p_\perp) δ(p_z)$ and $f(\p) = F(p_z) δ^{(2)}(p_\perp)$, where $f(\p)$ is the distribution of particles in momentum space.

preprint2003arXiv

Effective Kinetic Theory for High Temperature Gauge Theories

Quasiparticle dynamics in relativistic plasmas associated with hot, weakly-coupled gauge theories (such as QCD at asymptotically high temperature $T$) can be described by an effective kinetic theory, valid on sufficiently large time and distance scales. The appropriate Boltzmann equations depend on effective scattering rates for various types of collisions that can occur in the plasma. The resulting effective kinetic theory may be used to evaluate observables which are dominantly sensitive to the dynamics of typical ultrarelativistic excitations. This includes transport coefficients (viscosities and diffusion constants) and energy loss rates. We show how to formulate effective Boltzmann equations which will be adequate to compute such observables to leading order in the running coupling $g(T)$ of high-temperature gauge theories [and all orders in $1/\log g(T)^{-1}$]. As previously proposed in the literature, a leading-order treatment requires including both $2<->2$ particle scattering processes as well as effective ``$1<->2$'' collinear splitting processes in the Boltzmann equations. The latter account for nearly collinear bremsstrahlung and pair production/annihilation processes which take place in the presence of fluctuations in the background gauge field. Our effective kinetic theory is applicable not only to near-equilibrium systems (relevant for the calculation of transport coefficients), but also to highly non-equilibrium situations, provided some simple conditions on distribution functions are satisfied.

preprint2003arXiv

Monte Carlo simulation of O(2) phi^4 field theory in three dimensions

Using standard numerical Monte Carlo lattice methods, we study non-universal properties of the phase transition of three-dimensional phi^4 theory of a 2-component real field phi = (phi_1,phi_2) with O(2) symmetry. Specifically, we extract the renormalized values of <phi^2>/u and r/u^2 at the phase transition, where the continuum action of the theory is \int d^3x [ (1/2) |\gradϕ|^2 + \half r ϕ^2 + {u\over4!} ϕ^4 ]. These values have applications to calculating the phase transition temperature of dilute or weakly-interacting Bose gases (both relativistic and non-relativistic). In passing, we also provide perturbative calculations of various O(a) lattice-spacing errors in three-dimensional O(N) scalar field theory, where (a) is the lattice spacing.

preprint2003arXiv

Transport coefficients in high temperature gauge theories: (II) Beyond leading log

Results are presented of a full leading-order evaluation of the shear viscosity, flavor diffusion constants, and electrical conductivity in high temperature QCD and QED. The presence of Coulomb logarithms associated with gauge interactions imply that the leading-order results for transport coefficients may themselves be expanded in an infinite series in powers of 1/log(1/g); the utility of this expansion is also examined. A next-to-leading-log approximation is found to approximate the full leading-order result quite well as long as the Debye mass is less than the temperature.

preprint2002arXiv

Photon and Gluon Emission in Relativistic Plasmas

We recently derived, using diagrammatic methods, the leading-order hard photon emission rate in ultra-relativistic plasmas. This requires a correct treatment of multiple scattering effects which limit the coherence length of emitted radiation (the Landau-Pomeranchuk-Migdal effect). In this paper, we provide a more physical derivation of this result, and extend the treatment to the case of gluon radiation.

preprint2002arXiv

Transport Coefficients in Large $N_f$ Gauge Theory: Testing Hard Thermal Loops

We compute shear viscosity and flavor diffusion coefficients for ultra-relativistic gauge theory with many fermionic species, Nf >> 1, to leading order in 1/Nf. The calculation is performed both at leading order in the effective coupling strength g^2 Nf, using the Hard Thermal Loop (HTL) approximation, and completely to all orders in g^2 Nf. This constitutes a nontrivial test of how well the HTL approximation works. We find that in this context, the HTL approximation works well wherever the renormalization point sensitivity of the leading order HTL result is small.

preprint2001arXiv

Lower Bound on the Propagation Speed of Gravity from Gravitational Cherenkov Radiation

Recently, interesting 4-D Lorentz violating models have been proposed, in which all particles have a common maximum velocity $c$, but gravity propagates (in the preferred frame) with a different maximum velocity $c_g \neq c$. We show that the case $c_g < c$ is very tightly constrained by the observation of the highest energy cosmic rays. Assuming a galactic origin for the cosmic rays gives a conservative bound of $c-c_g < 2 \times 10^{-15} c$; if the cosmic rays have an extragalactic origin the bound is orders of magnitude tighter, of order $c-c_g < 2 \times 10^{-19} c$.

preprint2001arXiv

Non-perturbative computation of the bubble nucleation rate in the cubic anisotropy model

At first order phase transitions the transition proceeds through droplet nucleation and growth. We discuss a lattice method for calculating the droplet nucleation rate, including the complete dynamical factors. The method is especially suitable for very strongly suppressed droplet nucleation, which is often the case in physically interesting transitions. We apply the method to the 3-dimensional cubic anisotropy model in a parameter range where the model has a radiatively induced strong first order phase transition, and compare the results with analytical approaches.

preprint2001arXiv

Photon Emission from Quark-Gluon Plasma: Complete Leading Order Results

We compute the photon emission rate of an equilibrated, hot QCD plasma at zero chemical potential, to leading order in both alpha_{EM} and the QCD coupling g_s(T). This requires inclusion of near-collinear bremsstrahlung and inelastic pair annihilation contributions, and correct incorporation of Landau-Pomeranchuk-Migdal suppression effects for these processes. Analogous results for a QED plasma are also included.

preprint2001arXiv

Photon Emission from Ultrarelativistic Plasmas

The emission rate of photons from a hot, weakly coupled ultrarelativistic plasma is analyzed. Leading-log results, reflecting the sensitivity of the emission rate to scattering events with momentum transfers from $gT$ to $T$, have previously been obtained. But a complete leading-order treatment requires including collinearly enhanced, inelastic processes such as bremsstrahlung. These inelastic processes receive O(1) modifications from multiple scattering during the photon emission process, which limits the coherence length of the emitted radiation (the Landau-Pomeranchuk-Migdal effect). We perform a diagrammatic analysis to identify, and sum, all leading-order contributions. We find that the leading-order photon emission rate is not sensitive to non-perturbative $g^2 T$ scale dynamics. We derive an integral equation for the photon emission rate which is very similar to the result of Migdal in his original discussion of the LPM effect. The accurate solution of this integral equation for specific theories of interest will be reported in a companion paper.

preprint2001arXiv

T_c for homogeneous dilute Bose gases: a second-order result

The transition temperature for a dilute, homogeneous, three-dimensional Bose gas has the expansion T_c = T_0 {1 + c_1 a n^(1/3) + [c_2' ln(a n^(1/3)) + c_2''] a^2 n^(2/3) + O(a^3 n)}, where a is the scattering length, n the number density, and T_0 the ideal gas result. The first-order coefficient c_1 depends on non-perturbative physics. In this paper, we show that the coefficient c_2' can be computed perturbatively. We also show that the remaining second-order coefficient c_2'' depends on non-perturbative physics but can be related, by a perturbative calculation, to quantities that have previously been measured using lattice simulations of three-dimensional O(2) scalar field theory. Making use of those simulation results, we find T_c = T_0 {1 + (1.32+-0.02) a n^(1/3) + [19.7518 ln(a n^(1/3)) + (75.7+-0.4)] a^2 n^(2/3) + O(a^3 n)}.

preprint1999arXiv

Chern-Simons Number Diffusion and Hard Thermal Loops on the Lattice

We develop a discrete lattice implementation of the hard thermal loop effective action by the method of added auxiliary fields. We use the resulting model to measure the sphaleron rate (topological susceptibility) of Yang-Mills theory at weak coupling. Our results give parametric behavior in accord with the arguments of Arnold, Son, and Yaffe, and are in quantitative agreement with the results of Moore, Hu, and Muller.

preprint1999arXiv

Real Time Simulations in Lattice Gauge Theory

I review the study of real (Minkowski) time correlators in hot, weakly coupled Yang-Mills theory via lattice methods. I concentrate on the Minkowski time topological susceptibility, which is related to the efficiency of baryon number violation at high temperature. It can be computed by approximating the IR fields as classical and solving their dynamics nonperturbatively on the lattice. However it is essential to include the UV degrees of freedom. Their influence can be computed perturbatively and included in the lattice model by the addition of auxiliary fields.

preprint1999arXiv

Was the Electroweak Phase Transition Preceded by a Color-Broken Phase?

It has been suggested, in connection with electroweak baryogenesis in the Minimal Supersymmetric Standard Model (MSSM), that the right-handed top squark has a negative mass squared parameter, such that its field could condense prior to the electroweak phase transition (EWPT). Thus color and electric charge could have been broken just before the EWPT. Here we investigate whether the tunneling rate from the color-broken vacuum can ever be large enough for the EWPT to occur in this case. We find that, even when all parameters are adjusted to their most favorable values, the nucleation rate is many orders of magnitude too small. We conclude that, without additional physics beyond the MSSM, the answer to our title question is ``no.'' This gives constraints in the plane of the light stop mass versus parameters related to stop mixing. However it may be possible to get color breaking in extended models, such as those with R-parity violation.

preprint1998arXiv

Chern-Simons number diffusion with hard thermal loops

We construct an extension of the standard Kogut-Susskind lattice model for classical 3+1 dimensional Yang-Mills theory, in which ``classical particle'' degrees of freedom are added. We argue that this will correctly reproduce the ``hard thermal loop'' effects of hard degrees of freedom, while giving a local implementation which is numerically tractable. We prove that the extended system is Hamiltonian and has the same thermodynamics as dimensionally reduced hot Yang-Mills theory put on a lattice. We present a numerical update algorithm and study the abelian theory to verify that the classical gauge theory self-energy is correctly modified. Then we use the extended system to study the diffusion constant for Chern-Simons number. We verify the Arnold-Son-Yaffe picture that the diffusion constant is inversely proportional to hard thermal loop strength. Our numbers correspond to a diffusion constant of Gamma = 29 +- 6 alpha_w^5 T^4 for m_D^2 = 11 g^2 T^2/6.

preprint1998arXiv

Measuring the Broken Phase Sphaleron Rate Nonperturbatively

We present details for a method to compute the broken phase sphaleron rate (rate of hot baryon number violation below the electroweak phase transition) nonperturbatively, using a combination of multicanonical and real time lattice techniques. The calculation includes the ``dynamical prefactor,'' which accounts for prompt recrossings of the sphaleron barrier. The prefactor depends on the hard thermal loops, getting smaller with increasing Debye mass; but for realistic Debye masses the effect is not large. The baryon number erasure rate in the broken phase is slower than a perturbative estimate by about exp(-3.6). Assuming the electroweak phase transition has enough latent heat to reheat the universe to the equilibrium temperature, baryon number is preserved after the phase transition if the ratio of (``dimensionally reduced'' thermal) scalar to gauge couplings (lambda / g^2) is less than .037.

preprint1997arXiv

Classical Field Dynamics of the Electroweak Phase Transition

We investigate the thermodynamics and dynamics of the electroweak phase transition by modelling the infrared physics with classical Yang-Mills Higgs theory. We discuss the accuracy of this approach and conclude that, for quantities whose determination is dominated by the infrared, the classical method should be correct up to parametrically suppressed (ie O(alpha)) corrections. For a Higgs self-coupling which at tree level corresponds to m_H ~ 50 GeV, we determine the jump in the order parameter to be delta phi = 1.5 gT, the surface tension to be sigma = 0.07 g^4 T^3, and the friction coefficient on the moving bubble wall due to infrared bosons to be η= P/v_w = 0.03 \pm .004 g^6 T^4. We also investigate the response of Chern-Simons number to a spatially uniform chemical potential and find that it falls off a short distance inside the bubble wall, both in equilibrium and below the equilibrium temperature.

preprint1997arXiv

Lattice Chern-Simons Number Without Ultraviolet Problems

We develop a topological method of measuring Chern-Simons number change in the real time evolution of classical lattice SU(2) and SU(2) Higgs theory. We find that the Chern-Simons number diffusion rate per physical 4-volume is very heavily suppressed in the broken phase, and that it decreases with lattice spacing in pure Yang-Mills theory, although not as quickly as predicted by Arnold, Son, and Yaffe.

preprint1996arXiv

Fermion Fluctuation Determinant and Sphaleron Bound

We investigate analytically the fermionic fluctuation determinant at finite temperatures in the minimal standard model, including all operators up to dimension 6 and all contributions to the effective potential to all orders in the high $T$ expansion, to 1 loop. We apply the results to find corrections to the Sphaleron erasure rate in the broken phase. We conclude that the analytic treatment of fermions is very reliable and that there is a great deal of baryon erasure after the phase transition for any physical Higgs mass.

preprint1995arXiv

Bubble Wall Velocity at the Electroweak Phase Transition

We calculate the velocity and thickness of a bubble wall at the electroweak phase transition in the Minimal Standard Model. We model the wall with semiclassical equations of motion and show that friction arises from the deviation of massive particle populations from thermal equilibrium. We treat these with Boltzmann equations in a fluid approximation in the background of the wall. Our analysis improves on the previous work by using the two loop effective potential, accounting for particle transport, and determining the wall thickness dynamically. We find that the wall is significantly thicker than at phase equilibrium, and that the velocity is fairly high, $v_w \simeq 0.7c$, and quite weakly dependent on the Higgs mass.

preprint1995arXiv

How fast can the wall move? A study of the electroweak phase transition dynamics

We consider the dynamics of bubble growth in the Minimal Standard Model at the electroweak phase transition and determine the shape and the velocity of the phase boundary, or bubble wall. We show that in the semi-classical approximation the friction on the wall arises from the deviation of massive particle populations from thermal equilibrium. We treat these with Boltzmann equations in a fluid approximation. This approximation is reasonable for the top quarks and the light species while it underestimates the friction from the infrared $W$ bosons and Higgs particles. We use the two-loop finite temperature effective potential and find a subsonic bubble wall for the whole range of Higgs masses $0<m_H<90$GeV. The result is weakly dependent on $m_H$: the wall velocity $v_w$ falls in the range $0.36<v_w<0.44$, while the wall thickness is in the range $29> L T > 23 $. The wall is thicker than the phase equilibrium value because out of equilibrium particles exert more friction on the back than on the base of a moving wall. We also consider the effect of an infrared gauge condensate which may exist in the symmetric phase; modelling it simplemindedly, we find that the wall may become supersonic, but not ultrarelativistic.