Finite chemical potential equation of state for QCD from an alternative expansion scheme
The Taylor expansion approach to the equation of state of QCD at finite chemical potential struggles to reach large chemical potential μ B . This is primarily due to the intrinsic diffculty in precisely determining higher order Taylor coefficients, as well as the structure of the temperature dependence of such observables. In these proceedings, we illustrate a novel scheme [1] that allows us to extrapolate the equation of state of QCD without suffering from the poor convergence typical of the Taylor expansion approach. We continuum extrapolate the coefficients of our new expansion scheme and show the thermodynamic observables up to μ B / T ≤ 3.5.
Authors:
; ; ; ; ; ; ; ;
Editors:
; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10336539
Journal Name:
EPJ Web of Conferences
Volume:
259
Page Range or eLocation-ID:
10015
ISSN:
2100-014X
1. Abstract Consider $$(X_{i}(t))$$ ( X i ( t ) ) solving a system of N stochastic differential equations interacting through a random matrix $${\mathbf {J}} = (J_{ij})$$ J = ( J ij ) with independent (not necessarily identically distributed) random coefficients. We show that the trajectories of averaged observables of $$(X_i(t))$$ ( X i ( t ) ) , initialized from some $$\mu$$ μ independent of  $${\mathbf {J}}$$ J , are universal, i.e., only depend on the choice of the distribution $$\mathbf {J}$$ J through its first and second moments (assuming e.g., sub-exponential tails). We take a general combinatorial approach to proving universality for dynamical systems with random coefficients, combining a stochastic Taylor expansion with a moment matching-type argument. Concrete settings for which our results imply universality include aging in the spherical SK spin glass, and Langevin dynamics and gradient flows for symmetric and asymmetric Hopfield networks.
2. A bstract We study heavy quarkonium production associated with gluons in e + e − annihilation as an illustration of the perturbative QCD (pQCD) factorization approach, which incorporates the first nonleading power in the energy of the produced heavy quark pair. We show how the renormalization of the four-quark operators that define the heavy quark pair fragmentation functions using dimensional regularization induces “evanescent” operators that are absent in four dimensions. We derive closed forms for short-distance coefficients for quark pair production to next-to-leading order ( $${\alpha}_s^2$$ α s 2 ) in the relevant color singlet and octet channels. Using non-relativistic QCD (NRQCD) to calculate the heavy quark pair fragmentation functions up to v 4 in the velocity expansion, we derive analytical results for the differential energy fraction distribution of the heavy quarkonium. Calculations for $${}^3{S}_1^{\left[1\right]}$$ 3 S 1 1 and $${}^1{S}_0^{\left[8\right]}$$ 1 S 0 8 channels agree with analogous NRQCD analytical results available in the literature, while several color-octet calculations of energy fraction distributions are new. We show that the remaining corrections due to the heavy quark mass fall off rapidly in the energy of the produced state. To explore the importance of evolutionmore »
5. The thermal radiative transfer (TRT) equations form an integro-differential system that describes the propagation and collisional interactions of photons. Computing accurate and efficient numerical solutions TRT are challenging for several reasons, the first of which is that TRT is defined on a high-dimensional phase space that includes the independent variables of time, space, and velocity. In order to reduce the dimensionality of the phase space, classical approaches such as the P$_N$ (spherical harmonics) or the S$_N$ (discrete ordinates) ansatz are often used in the literature. In this work, we introduce a novel approach: the hybrid discrete (H$^T_N$) approximation to the radiative thermal transfer equations. This approach acquires desirable properties of both P$_N$ and S$_N$, and indeed reduces to each of these approximations in various limits: H$^1_N$ $\equiv$ P$_N$ and H$^T_0$ $\equiv$ S$_T$. We prove that H$^T_N$ results in a system of hyperbolic partial differential equations for all $T\ge 1$ and $N\ge 0$. Another challenge in solving the TRT system is the inherent stiffness due to the large timescale separation between propagation and collisions, especially in the diffusive (i.e., highly collisional) regime. This stiffness challenge can be partially overcome via implicit time integration, although fully implicit methods may become computationally expensivemore »