Search

These lecture notes, prepared for the 2024 XQCD PhD, provide an introduction to the analytic structure of an equation of state near a second-order phase transition and its most prominent landmark: the Yang-Lee edge singularity. In addition to discussing general properties, the notes review recent theoretical progress in locating the QCD critical point by tracking the trajectory of the Yang-Lee edge singularity., Comment: 37 pages, 19 figures

We present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of Quantum Chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar density, which is known as model H. We show that the observed dynamical scaling behavior depends on the correlation length and the shear viscosity of the fluid. As the correlation length is increased or the viscosity is decreased we observe a cross-over from the dynamical exponent of critical diffusion, $z\simeq 4$, to the expected scaling exponent of model H, $z\simeq 3$. We use our method to investigate time-dependent correlation function of non-Gaussian moments $M^n(t)$ of the order parameter. We find that the relaxation time depends in non-trivial manner on the power $n$., Comment: 5 pages, 4 figures

We discuss the scaling of the Yang-Lee singularity (YLs) and show how the universal scaling can be used to locate phase transitions in QCD. We describe two complementary methods to extract the location of the Yang-Lee singularity from lattice QCD data of the baryon-number density and higher order cumulants of the baryon number, obtained at imaginary chemical potential. The first method (multi-point Pad\'e resummation) is used to determine the Roberge-Weiss phase transition temperature. Our continuum extrapolated result is $T_{RW}=211.1\pm3.1$ MeV. The second method is based on the asymptotic behaviour of the Fourier coefficients of the baryon-number density. We discuss the derivation of a fitting function and demonstrate that the procedure can successfully locate the YLs in the Quark Meson model., Comment: 7 pages, 4figures, Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

In this paper we study the universal properties of the baryon chemical potential Fourier coefficients in Quantum Chromodynamics. We show that by following a well-defined strategy, the Fourier coefficients can be used to locate Yang-Lee edge singularities associated with chiral phase transition (and by extension with the Roberge-Weiss) in the complex chemical potential plane. We comment on the viability of performing this analysis using lattice QCD data., Comment: 18 pages, 7 figures

We investigate the effect of gluon Bose enhancement in the nuclear wave function on the dijet production in incoherent diffractive processes in DIS and ultraperipheral collisions. We demonstrate that Bose enhancement leads to an enhancement of diffractive dijet production cross section when the transverse momenta of the two jets are aligned at zero relative angle. This enhancement is maximal when the magnitude of the transverse momenta of the two jets are equal, and disappears rather quickly as a function of the ratio of the two momenta. We study both the dilute limit and fully nonlinear dense regime where the nuclear wave function is evolved with the leading order JIMWLK equation. In both cases we observe a visible effect, with it being enhanced by the evolution due to the dynamical generation of the color neutralization scale., Comment: 15 pages, 4 figures

This paper introduces a transverse-momentum dependent (TMD) factorization scheme designed to unify both large and small Bjorken-x regimes. We compute the next-to-leading order (NLO) quantum chromodynamics (QCD) corrections to the gluon TMD operator for an unpolarized hadron within this proposed scheme. This leads to the emergence of a new TMD evolution, incorporating those in transverse momentum, rapidity, and Bjorken-x. When matched to the collinear factorization scheme, our factorization scheme faithfully reproduces the well-established Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) and Collins-Soper-Sterman (CSS) evolutions. Conversely, matching with high-energy factorization not only yields the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution but also reveals distinctive signatures of CSS logarithms. The development of this novel TMD factorization scheme, capable of seamlessly reconciling disparate Bjorken-x regimes and faithfully reproducing established QCD evolution equations, has the potential to significantly advance our comprehension of high-energy processes and three-dimensional parton structures of hadrons., Comment: 38 pages, 7 figures; v2: Published version

Various conventional running coupling prescriptions reproducing $\beta_ 0$-dependent terms of NLO JIMWLK are reviewed and found to be theoretically inconsistent: the JIMWLK evolution Hamiltonian with running coupling violates the requirement of positive semidefiniteness. This requirement appears to be tightly related to the possibility of having a Langevin formulation for the evolution. We also review the scheme that attributes a part of $\beta_0$-dependent terms to the DGLAP evolution of the projectile. The remaining $\beta_0$-dependent contributions sum up into so-called the ``daughter dipole'' prescription, which leads to a manifestly positive semidefinite Hamiltonian., Comment: 28 pages; version accepted for publication

We reanalyze the origin of the large transverse logarithms associated with the QCD one loop beta-function coefficient in the NLO JIMWLK Hamiltonian. We show that some of these terms are not associated with the running of the QCD coupling constant but rather with the DGLAP evolution. The DGLAP-like resummation of these logarithms is mandatory within the JIMWLK Hamiltonian, as long as the color correlation length in the projectile is larger than that in the target. This regime in fact covers the whole range of rapidities at which JIMWLK evolution is supposed to be applicable. We derive the RG equation that resums these logarithms to all orders in alpha_s in the JIMWLK Hamiltonian. This is a nonlinear equation for the eikonal scattering matrix S(x). We solve this equation and perform the DGLAP resummation in two simple cases: the dilute limit, where both the projectile and the target are far from saturation, and the saturated regime, where the target correlation length also determines its saturation momentum., Comment: 43 pages

We extend the results of [Phys.Rev.Lett. 128 (2022) 18], where we argued that in the controlled environment of the Deep Inelastic Scattering experiments, Bose-Einstein correlation between gluons in a hadronic wave function can be accessed through the production of the diffractive dijet plus a third jet. In this observable, Bose-Einstein correlation causes the enhancement of the production cross sections at the zero relative angle between the transverse momentum imbalance of the photon-going dijet and the transverse momentum of the gluon jet, when the magnitude of the momentum imbalance is about the same as the magnitude of the produced gluon. In the present paper, we account for multiple scattering and non-linear effect in the target wave function. Although our equations can be applied to any high-energy DIS kinematics, to make them tractable numerically, we consider the high-momentum limit (momentum larger than $Q_s$) for the total momentum of the dijet, momentum imbalance, and the momentum of the produced gluon. By performing explicit numerical calculations, we confirm that the signal is present after accounting for multiple scattering., Comment: 31 pages, 7 figures

We calculate the entanglement entropy of a model proton wave function in coordinate space by integrating out degrees of freedom outside a small circular region $\bar A$ of radius $L$, where $L$ is much smaller than the size of the proton. The wave function provides a nonperturbative distribution of three valence quarks. In addition, we include the perturbative emission of a single gluon and calculate the entanglement entropy of gluons in $\bar A$. For both, quarks and gluons we obtain the same simple result: $S_E =-\int\frac{dx}{\Delta x}\, N_{L^2}(x)\log[N_{a^2}(x)]$, where $a$ is the UV cutoff in coordinate space and $\Delta x$ is the longitudinal resolution scale. Here $N_{S}(x)$ is the number of partons (of the appropriate species) with longitudinal momentum fraction $x$ inside an area $S$. It is related to the standard parton distribution function (PDF) by $N_S(x)=\frac{S}{A_p}\, \Delta x\, F(x)$, where $A_p$ denotes the transverse area of the proton., Comment: 24 pages, 1 figure

Motivated by recent discussions of entanglement in the context of high energy scattering, we consider the relation between the entanglement entropy of a highly excited state of a quantum system and the classical entanglement entropy of the corresponding classical system. We show on the example of two weakly coupled harmonic oscillators, that the two entropies are equal. Quantum mechanically, the reduced density matrix which yields this entropy is close to the maximally entangled state. We thus observe that the nature of entanglement in this type of state is purely classical., Comment: 19 pages

Employing the functional renormalization group approach at next-to-leading order of the derivative expansion, we refine our earlier findings for the location of the Yang-Lee edge singularity in classic O(N) universality classes. For the universality classes of interest to QCD, in three dimensions, we found $|z_c|/R_\chi^{1/\gamma} = 1.612(9),\ 1.597(3)$ for $N=2$, $4$ correspondingly. We also established $|z_c| = 2.04(8),\ 1.69(3)$ for $N=2$, $4$ albeit with greater systematic error., Comment: 36 pages; 6 figures; version accepted for publication

We provide the first calculation of two-gluon production at mid-rapidity in ultra-peripheral collisions in the Color Glass Condensate framework. To estimate systematic uncertainty associated with poor understanding of the wave function of the nearly real photon, we consider two diametrically different models: the dilute quark-antiquark dipole approximation and a vector meson, in which color charge density is approximated by McLerran-Venugopalan model. In the experimentally relevant range, the target nucleus can be faithfully approximated by a highly saturated state. This simplification enables us to perform efficient numerical simulations and extract the two-gluon correlation functions and the associated azimuthal harmonics., Comment: 35 pages,8 figures

We determine the universal location of the Yang-Lee edge singularity in the entire relevant domain of spatial dimensions $1\le d \le 4$ for the Ising universality class. To that end, we present analytical results for $d=1,2,4$ and near four dimensions. For $d=3$ and a set of fractional dimensions, we perform numerical calculations using a systematic Functional Renormalization Group approach., Comment: 42 pages, 4 figures, 2 tables; version accepted for publication

We give an elementary and pedagogical review of the phase diagrams which are possible in Quantum ChromoDynamics (QCD). Currently, the emphasis is upon the appearance of a critical endpoint, where disordered and ordered phases meet. In many models, though, a Lifshitz point also arises. At a Lifshitz point, three phases meet: disordered, ordered, and one where spatially inhomogeneous phases arise. At the level of mean field theory, the appearance of a Lifshitz point does not dramatically affect the phase diagram. We argue, however, that fluctuations about the Lifshitz point are very strong in the infrared, and significantly alter the phase diagram. We discuss at length the analogy to inhomogenous polymers, where the Lifshitz regime produces a bicontinuous microemulsion. We briefly mention the possible relevance to the phase diagram of QCD., Comment: Contribution to the International Conference on New Frontiers in Physics, 2018. v2: corrected acknowledgements for RDP

We revisit and extend the calculation of the density matrix and entanglement entropy of a Color Glass Condensate by including the leading saturation corrections in the calculation. We show that the density matrix is diagonal in the quasi particle basis, where it has the Boltzmann form. The quasi particles in a wide interval of momenta behave as massless two-dimensional bosons with the temperature proportional to the typical semi-hard scale $T=Q_s/\sqrt{\alpha_sN_c}$. Thus the semi-hard momentum region $Q_s

In high energy proton-nucleus collisions, the gluon saturation effects from the nucleus are fully incorporated into the light-like Wilson lines. The gluon saturation effects from the proton, which are anticipated to be important either in the extreme high energy limit or towards the dense-dense (nucleus-nucleus) collision regimes, have been studied perturbatively within the Color Glass Condensate effective theory in previous papers of this series. A configuration-by-configuration expression for the single inclusive semi-hard gluon production including the first saturation correction was obtained. In this paper, we perform ensemble averaging in the McLerran-Venugopalan model and the Dipole Approximation. We find that, in the saturation correction, the effects of the initial state interactions are negligible while the final state interactions play most important role and give a positive-valued contribution to the semi-hard gluon spectrum. Furthermore, we show that the single gluon spectrum scales approximately $1/k_{\perp}^{4}$ at small $k_{\perp}$, suggesting that a resummation of higher order saturation corrections is required to regulate the infrared region of the gluon spectrum., Comment: 37 pages, 12 figures

For theories plagued with a sign problem at finite density, a Taylor expansion in the chemical potential is frequently used for lattice gauge theory based computations of the equation of state. Recently, in arXiv:2106.03165, a new resummation scheme was proposed for such an expansion that resums contributions of correlation functions of conserved currents to all orders in the chemical potential. Here, we study the efficacy of this resummation scheme using a solvable low energy model, namely the mean-field quark-meson model. After adapting the scheme for a mean-field analysis, we confront the results of this scheme with the direct solution of the model at finite density as well as compare with results from Taylor expansions. We study to what extent the two methods capture the analytical properties of the equation of state in the complex chemical potential plane. As expected, the Taylor expansion breaks down as soon as the baryon chemical potential reaches the radius of convergence defined by the Yang-Lee edge singularity. Encouragingly, the resummation not only captures the location of the Yang-Lee edge singularity accurately, but is also able to describe the equation of state for larger chemical potentials beyond the location of the edge singularity for a wide range of temperatures., Comment: 9 pages, 9 figures; matches the published version

Considering a dilute-dense situation suitable for pA collisions, we compute in the Color Glass Condensate the correlation between azimuthal asymmetries, specifically the squared second Fourier coefficient $v_2^2$, and the total multiplicity in the event. We also analyse the correlation between $v_2^2$ and the mean squared transverse momentum of particles in the event. In both cases, we find that the correlations are generally very small, consistent with the observations. We also note an interesting sharp change in the value of $v_2^2$ and its correlations as a function of the width of the transverse momentum bin, related with a change of the dominance of Bose and HBT quantum correlations., Comment: 6 pages, 5 figures; contribution to DIS2021; submission to SciPost

We study correlations originating from the quantum nature of gluons in a hadronic wave function. Bose-Einstein correlation between identical particles lead to the enhancement in the number of pairs of gluons with the same quantum numbers and small relative momentum. We show that these preexisting correlations can be probed in Deep Inelastic Scattering experiments at high energy. Specifically, we consider diffractive dijet plus a third jet production. The azimuthal dependence displays a peak at the zero relative angle between the transverse momentum imbalance of the photon-going dijet and the transverse momentum of the hadron-going jet. Our calculations explicitly show that the peak originates from Bose enhancement. Comparing electron-proton to electron-nucleus collisions, we demonstrate that the nuclear target enhances the relative strength of the peak. With the future high luminosity Electron-Ion Collider the proposed measurements of gluon Bose enhancement become experimentally feasible., Comment: 5 pages, 4 figures; updated version, corrected a mistake related to the Lipatov vertex, figures updated, qualitative predications and conclusions remain unchanged

Exploiting recently obtained analytic solutions of classical Yang-Mills equations for higher order perturbations in the field of the dilute object (proton), we derive the complete first saturation correction to the single inclusive semi-hard gluon production in high energy proton-nucleus collisions by applying the Lehmann-Symanzik-Zimmermann reduction formula. We thus finalize the program started by Balitsky (see Ref.~\cite{Balitsky:2004rr}) and independently by Chirilli, Kovchegov and Wertepny (see Ref.~\cite{Chirilli:2015tea}) albeit using a very different approach to carry out our calculations. We extracted the functional dependence of gluon spectrum on the color charge densities of the colliding objects; thus our results can be used to evaluate complete first saturation correction to the double/multiple inclusive gluon productions., Comment: 60 pages and 7 figures

In high energy proton-nucleus collisions, the single- and double-inclusive soft gluon productions at the leading order have been calculated and phenomenologically studied in various approaches for many years. These studies do not take into account the saturation and multiple rescatterings in the field of the proton. The first saturation correction to these leading order results (the terms that are enhanced by the combination $\alpha_s^2 \mu^2$, where $\mu^2$ is the proton's color charge squared per unit transverse area) has not been completely derived despite recent attempts using a diagrammatic approach. This paper is the first in a series of papers towards analytically completing the first saturation correction to physical observables in high energy proton-nucleus collisions. Our approach is to analytically solve the classical Yang-Mills equations in the dilute-dense regime using the Color Glass Condensate effective theory and compute physical observables constructed from classical gluon fields. In the current paper, the Yang-Mills equations are solved perturbatively in the field of the dilute object (the proton). Next-to-leading order and next-to-next-to-leading order analytic solutions are explicitly constructed. A systematic way to obtain all higher order analytic solutions is outlined, Comment: 53 pages, 9 figures

Within the dense-dilute Color Glass Condensate approach, and using the Golec-Biernat-Wuesthoff model for the dipole scattering amplitude, we calculate $v_2^2$ as well as the correlations between $v_2^2$ and both the total multiplicity and the mean transverse momentum of produced particles. We find that the correlations are generally very small consistent with the observations. We note an interesting sharp change in the value of $v^2_2$ as well as of its correlations as a function of the width of the transverse momentum bin. This crossover is associated with the change from Bose enhancement dominance of the correlation for narrow bin to HBT dominated correlations for larger bin width., Comment: 28 pages, 4 figures

At a nonzero temperature T, a constant field $\overline{A}_0 \sim T/g$ generates nontrivial eigenvalues of the thermal Wilson line. We discuss contributions to the free energy of such a holonomous plasma when the coupling constant, $g$, is weak. We review the computation to $\sim g^2$ by several alternate methods, and show that gauge invariant sources, which are nonlinear in the gauge potential $A_0$, generate novel contributions to the gluon self energy at $\sim g^2$. These ensure the gluon self energy remains transverse to $\sim g^2$, and are essential in computing contributions to the free energy at $\sim g^3$ for small holonomy, $\overline{A}_0 \sim T$. We show that the contribution $\sim g^3$ from off-diagonal gluons is discontinuous as the holonomy vanishes. The contribution from diagonal gluons is continuous as the holonomy vanishes, but sharply constrains the possible sources which generate nonzero holonomy, and must involve an infinite number of Polyakov loops., Comment: 32 pages, 2 figures

We define the "entropy of ignorance" which quantifies the entropy associated with ability to perform only a partial set of measurement on a quantum system. For a parton model the entropy of ignorance is equal to a Boltzmann entropy of a classical system of partons. We analyze a calculable model used for describing low x gluons in Color Glass Condensate approach, which has similarities with the parton model of QCD. In this model we calculate the entropy of ignorance in the particle number basis as well as the entanglement entropy of the observable degrees of freedom. We find that the two are similar at high momenta, but differ by a factor of order unity at low momenta. This holds for the Renyi as well as von Neumann entropies. We conclude that the entanglement does not seem to play an important role in the context of the parton model., Comment: 11 pages, 1 figure; version accepted to Phys. Rev. D

We compute the term $\sim g^3$ in the free energy for a $SU(N)$ gauge theory with nonzero holonomy at nonzero temperature. If the holonomy is generated kinematically by the introduction of gauge invariant sources coupled to Polyakov loops, the contribution of charged (off-diagonal) gluons to the free energy at order $g^3$, ${\cal F}^{\left( 3\rm{:c.g.} \right)}$, is singular: ${\cal F}^{\left( 3\rm{:c.g.} \right)} \neq 0 $ without holonomy, but ${\cal F}^{\left( 3 \rm{:c.g.} \right)}= 0$ when the holonomy is nonzero, even infinitesimally. We show that the absence of the charged gluon contribution is required by gauge invariance alone and is therefore a universal feature., Comment: 14 pages

It was argued in arXiv:1805.09342 and arXiv:1807.00825 that the systematics of the azimuthal anisotropy coefficients $v_{2,3}$ measured in ultrarelativistic light-heavy ion collisions at RHIC and the LHC can be described in an initial state dilute-dense Color Glass Condensate (CGC) framework. We elaborate here on the discussion in these papers and provide further novel results that strengthen their conclusions. The underlying mathematical framework and numerical techniques employed are very similar to those in the CGC based IP-Glasma model used previously as initial conditions for heavy-ion collisions. The uncertainties in theory/data comparisons for small systems are discussed, with unknowns that are specific to the model distinguished from those that are generic to all models. We present analytical arguments that demonstrate that quantum interference effects such as Bose enhancement and Hanbury-Brown-Twiss correlations of gluons, as well as coherent multiple scattering of gluons in the projectile off color domains in the target, are enhanced in rare events. The quantum origins of the large anisotropies in small systems are corroborated by numerical results for deuteron-gold collisions that show that large anisotropies in rare configurations can occur when the nucleons in the projectile overlap significantly. This is at variance with the classical intuition of hydrodynamical models. We also comment on the consequences of ignoring the many-body color charge correlations of gluons in models that only consider geometrical fluctuations in the energy density., Comment: 25 pages, 13 figures, 1 appendix

We report on recent progress in understanding multiparticle correlations in small systems from the initial state. First, we consider a proof-of-principle parton model, which we use to demonstrate that many of the multiparticle correlations observed in light-heavy ion collisions, often ascribed to hydrodynamic collectivity, can be qualitatively reproduced in an initial state model. Then, we study the two-particle harmonics $v_2$ and $v_3$ for p/d/$^3$He+Au collisions at RHIC using the dilute-dense Color Glass Condensate Effective Field Theory framework. We show that this provides a viable alternative explanation to hydrodynamics, and elaborate on how such modeling can be improved., Comment: 4 pages, 2 figures. Proceedings for Quark Matter 2018 based on talk by M. Mace. v2: Minor edits, identical to published version

Simple power counting arguments in the dilute-dense framework of the Color Glass Condensate (CGC) Effective Field Theory predict that even and odd azimuthal anisotropy harmonics of two-particle correlations in proton-nucleus collisions at the LHC will respectively satisfy v^2_{2n} ~ N_{ch}^0 and v^2_{2n+1} ~ N_{ch}, where N_{ch} denotes the number of charged particles. We show that these expectations are borne out qualitatively, and even quantitatively, within systematic uncertainties, for v_2 and v_4 in comparisons with data from the ATLAS collaboration. We also observe that ATLAS data for the v_3 azimuthal harmonic are in excellent agreement with our qualitative expectation; quantitative comparisons are numerically challenging at present. The lessons from this study fully complement those gained by the recent comparison of the CGC dilute-dense framework [arXiv:1805.09342] to data from the PHENIX collaboration on small system collisions at RHIC., Comment: 6 pages, 4 figures

In the framework of dense-dilute CGC approach we study fluctuations in the multiplicity of produced particles in p-A collisions. We show that the leading effect that drives the fluctuations is the Bose enhancement of gluons in the proton wave function. We explicitly calculate the moment generating function that resums the effects of Bose enhancement. We show that it can be understood in terms of the Liouville effective action for the composite field which is identified with the fluctuating density, or saturation momentum of the proton. The resulting probability distribution turns out to be very close to the gamma-distribution. We also calculate the first correction to this distribution which is due to pairwise Hanbury Brown-Twiss correlations of produced gluons., Comment: 15 pages, 3 figures; version accepted for publication

We study analytically and numerically effects of spatial eccentricity of the projectile shape on the second flow harmonic in inclusive gluon production in p-A collisions in the CGC framework. Keeping the collision area fixed, we find that the two quantities are anti-correlated., Comment: 10 pages, 3 figures

We demonstrate that the striking systematics of two-particle azimuthal Fourier harmonics $v_2$ and $v_3$ in ultrarelativistic collisions of protons, deuterons and helium-3 ions off gold nuclei measured by the PHENIX Collaboration [arXiv:1805.02973] at the Relativistic Heavy Ion Collider (RHIC) is reproduced in the Color Glass Condensate (CGC) effective field theory. This contradicts the claim in [arXiv:1805.02973] that their data rules out initial state based explanations. The underlying systematics of the effect, as discussed previously in [arXiv:1705.00745,arXiv:1706.06260,arXiv:1801.09704], arise from the differing structure of strong color correlations between gluon domains of size $1/Q_S$ at fine ($p_\perp \gtrapprox Q_S$) or coarser ($p_\perp \lessapprox Q_S$) transverse momentum resolution. Further tests of the limits of validity of this framework can be carried out in light-heavy ion collisions at both RHIC and the Large Hadron Collider. Such measurements also offer novel opportunities for further exploration of the role of the surprisingly large short-range nuclear correlations measured at Jefferson Lab., Comment: 6 pages, 4 figures. v2: minor changes, version published in PRL

We show that, in the saturation/Color Glass Condensate framework, odd azimuthal harmonics of the two-gluon correlation function with a long-range separation in rapidity are generated by the higher-order saturation corrections in the interactions with the projectile and the target. At the very least, the odd harmonics require three scatterings in the projectile and three scatterings in the target. We derive the leading-order expression for the two-gluon production cross section which generates odd harmonics: the expression includes all-order interactions with the target and three interactions with the projectile. We evaluate the obtained expression both analytically and numerically, confirming that the odd-harmonics contribution to the two-gluon production in the saturation framework is non-zero., Comment: 31 pages, 8 figures

We consider the phase diagram of hadronic matter as a function of temperature, T , and baryon chemical potential, mu. Currently the dominant paradigm is a line of first order transitions which ends at a critical endpoint. In this work we suggest that spatially inhomogenous phases are a generic feature of the hadronic phase diagram at nonzero mu and low T . Familiar examples are pion and kaon condensates. At higher densities, we argue that these condensates connect onto chiral spirals in a quarkyonic regime. Both of these phases exhibit the spontaneous breaking of a global U(1) symmetry and quasi-long range order, analogous to smectic liquid crystals. We argue that there is a continuous line of first order transitions which separate spatially inhomogenous from homogenous phases, where the latter can be either a hadronic phase or a quark-gluon plasma. While mean field theory predicts that there is a Lifshitz point along this line of first order transitions, in three spatial dimensions strong infrared fluctuations wash out any Lifshitz point. Using known results from inhomogenous polymers, we suggest that instead there is a Lifshitz regime. Non-perturbative effects are large in this regime, where the momentum dependent terms for the propagators of pions and associated modes are dominated not by terms quadratic in momenta, but quartic. Fluctuations in a Lifshitz regime may be directly relevant to the collisions of heavy ions at (relatively) low energies, sqrt(s)/A : 1 to 20 GeV., Comment: 43 pages, 3 figures

We analyze the phase structure of $SU(\infty)$ gauge theory at finite temperature using matrix models. Our basic assumption is that the effective potential is dominated by double-trace terms for the Polyakov loops. As a function of the temperature, a background field for the Polyakov loop, and a quartic coupling, it exhibits a universal structure: in the large portion of the parameter space, there is a continuous phase transition analogous to the third-order phase transition of Gross, Witten and Wadia, but the order of phase transition can be higher than third. We show that different confining potentials give rise to drastically different behavior of the eigenvalue density and the free energy. Therefore lattice simulations at large $N$ could probe the order of phase transition and test our results., Comment: 7 pages, 2 figures, conference proceeding for Critical Point and Onset of Deconfinement - CPOD2017

We study phase transitions in $SU(\infty)$ gauge theories at nonzero temperature using matrix models. Our basic assumption is that the effective potential is dominated by double trace terms for the Polyakov loops. As a function of the various parameters, related to terms linear, quadratic, and quartic in the Polyakov loop, the phase diagram exhibits a universal structure. In a large region of this parameter space, there is a continuous phase transition whose order is larger than second. This is a generalization of the phase transition of Gross, Witten, and Wadia (GWW). Depending upon the detailed form of the matrix model, the eigenvalue density and the behavior of the specific heat near the transition differ drastically. We speculate that in the pure gauge theory, that although the deconfining transition is thermodynamically of first order, it can be nevertheless conformally symmetric at infinite $N$., Comment: 37 pages, 9 figures

We explore the influence of finite volume effects on baryon number fluctuations in a non-perturbative chiral model. In order to account for soft modes, we use the functional renormalization group in a finite volume, using a smooth regulator function in momentum space. We compare the results for a smooth regulator with those for a sharp (or Litim) regulator, and show that in a finite volume, the latter produces spurious artifacts. In a finite volume there are only apparent critical points, about which we compute the ratio of the fourth to the second order cumulant of quark number fluctuations. When the volume is sufficiently small the system has two apparent critical points; as the system size decreases, the location of the apparent critical point can move to higher temperature and lower chemical potential., Comment: minor text corrections, 24 pages, 10 figures

We consider how tetraquarks can affect the chiral phase transition in theories like QCD, with light quarks coupled to three colors. For two flavors the tetraquark field is an isosinglet, and its effect is minimal. For three flavors, however, the tetraquark field transforms in the same representation of the chiral symmetry group as the usual chiral order parameter, and so for very light quarks there may be two chiral phase transitions, which are both of first order. In QCD, results from the lattice indicate that any transition from the tetraquark condensate is a smooth crossover. In the plane of temperature and quark chemical potential, though, a crossover line for the tetraquark condensate is naturally related to the transition line for color superconductivity. For four flavors we suggest that a triquark field, antisymmetric in both flavor and color, combine to form hexaquarks., Comment: 34 pages, 4 figures

We modify the McLerran-Venugopalan model to include only a finite number of sources of color charge. We argue that Coulombic interactions between these color charges generates a source-source correlation function that properly includes the effects of color charge screening, a generalization of Debye screening for the Color Glass Condensate. Such a model may be useful for computing angular harmonics of flow measured in high energy hadron collisions for small systems. In this paper we provide a basic formulation of the problem on a lattice., Comment: 13 pages, 3 figures

A chiral matrix model applicable to QCD with 2+1 flavors is developed. This requires adding a SU(3)_L x SU(3)_R x Z(3)_A nonet of scalar fields, with both parities, and coupling these to quarks through a Yukawa coupling, y. Treating the scalar fields in mean field approximation, the effective Lagrangian is computed by integrating out quarks to one loop order. In addition to the usual symmetry breaking term, linear in the current quark mass m_qk, at a nonzero temperature T it is necessary to add a new term, ~ m_qk T^2. The parameters of the gluon part of the matrix model are identical to that for the pure glue theory without quarks. The parameters in the chiral matrix model are fixed by the values, at zero temperature, of the pion decay constant and the masses of the pions, kaons, eta, and eta'. The temperature for the chiral crossover at T_chi = 155 MeV is determined by adjusting the Yukawa coupling y. We find reasonable agreement with the results of numerical simulations on the lattice for the pressure and related quantities. In the chiral limit, besides the divergence in the chiral susceptibility there is also a milder divergence in the susceptibility between the Polyakov loop and the chiral order parameter, with critical exponent beta - 1. We compute derivatives with respect to a quark chemical potential to determine the susceptibilities for baryon number, the chi_2n. Especially sensitive tests are provided by chi_4 - chi_2 and by chi_6, which changes in sign about T_chi. The behavior of the susceptibilities in the chiral matrix model strongly suggests that as the temperature increases from T_chi, that the transition to deconfinement is significantly quicker than indicated by the measurements of the (renormalized) Polyakov loop on the lattice., Comment: 69 pages, 22 figures; version published in PRD

We compute the correlation functions of Polyakov loops in $SU(N_c)$ gauge theories by explicitly summing all diagrams at tree level in two special cases, for $N_c = 2$ and $N_c = \infty$. When $N_c =2$ we find the expected we find Coulomb-like behavior at short distances, $\sim 1/x$ as the distance $x \rightarrow 0$. In the planar limit at $N_c = \infty$ we find a weaker singularity, $\sim 1/\sqrt{x}$ as $x \rightarrow 0$. In each case, at short distances the behavior of the correlation functions between two Polyakov loops, and the corresponding Wilson loop, are the same. We suggest that such non-Coulombic behavior is an artifact of the planar limit., Comment: 23 pages, 1 figure

We apply the flow analysis for multi-particle correlations used in heavy ion collisions to multi-particle production from a Pomeron. We show that the n'th order angular harmonic arising from an m particle correlation v_n[m] satisfies v_n[m] ~ v_n[p] for n $\ge$ 1. We discuss some implications of this for the Color Glass Condensate description of high energy hadronic collisions., Comment: 11 pages, 3 figures

We investigate the process of Abelian pair production in the presence of strong inhomogeneous and time-dependent external electric fields. The spatial dependence of the external field is motivated by a non-Abelian color flux tube in heavy-ion collisions. We show that the inhomogeneity significantly increase the particle yield compared to that in the commonly used models with a constant and homogeneous field. Moreover our results indicate that in contrast to the latter, most of the particles are produced at the interface of the field profile in accordance with Heisenberg's prediction., Comment: 6 pages, 6 figures

We compute the collisional energy loss for a heavy quark above the critical temperature in Quantum ChromoDynamics (QCD). We work in the semi Quark-Gluon Plasma, which assumes that this region is dominated by the non-trivial holonomy of the thermal Wilson line. Relative to the result to leading order in perturbation theory, at a fixed value of the coupling constant we generically we find that collisional energy loss is suppressed by powers of the Polyakov loop, l < 1. For small values of the loop, this suppression is linear for the scattering off of light quarks, and quadratic for the scattering off of gluons, or for Compton scattering., Comment: 19 pages, 4 figures