skip to main content

Title: 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.
; ; ; ; ; ; ; ;
David, G.; Garg, P.; Kalweit, A.; Mukherjee, S.; Ullrich, T.; Xu, Z.; Yoo, I.-K.
Award ID(s):
Publication Date:
Journal Name:
EPJ Web of Conferences
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  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 »at energies much larger than the mass of the heavy quark, we solve the renormalization group equation perturbatively to two-loop order for the $$ {}^1{S}_0^{\left[8\right]} $$ 1 S 0 8 case.« less
  3. We study the phase transitions at finite temperature and density of the magnetic dual chiral density wave (MDCDW) phase. This spatially inhomogeneous phase emerges in cold, dense QCD in the presence of a strong magnetic field. Starting from the generalized Ginzburg-Landau (GL) expansion of the free energy, we derive several analytical formulas that enable fast numerical computation of the expansion coefficients to arbitrary order, allowing high levels of precision in the determination of the physical dynamical parameters, as well as in the transition curves in the temperature vs chemical potential plane at different magnetic fields. At magnetic fields and temperatures compatible with neutron star (NS) conditions, the MDCDW remains favored over the symmetric ground state at all densities. The phase’s “resilience” manifests in (1) a region of small but nonzero remnant mass and significant modulation at intermediate densities, originating in part from the nontrivial topology of the lowest Landau level, and (2) a region of increasing condensate parameters at high densities. Our analysis suggests the MDCDW condensate remains energetically favored at densities and temperatures much higher than previously considered, opening the possibility for this phase to be a viable candidate for the matter structure of even young neutron stars producedmore »by binary neutron star (BNS) mergers.« less
  4. null (Ed.)
    In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility H − 1 H^{-1} gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, n n and p p , is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of 0 0 prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough errormore »estimate (to establish the ℓ ∞ \ell ^\infty bound for n n and p p ), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme.« less
  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 »due to the strong nonlinearity and system size. On the other hand, explicit time-stepping schemes that are not also asymptotic-preserving in the highly collisional limit require resolving the mean-free path between collisions, making such schemes prohibitively expensive. In this work we develop a numerical method that is based on a nodal discontinuous Galerkin discretization in space, coupled with a semi-implicit discretization in time. In particular, we make use of a second order explicit Runge-Kutta scheme for the streaming term and an implicit Euler scheme for the material coupling term. Furthermore, in order to solve the material energy equation implicitly after each predictor and corrector step, we linearize the temperature term using a Taylor expansion; this avoids the need for an iterative procedure, and therefore improves efficiency. In order to reduce unphysical oscillation, we apply a slope limiter after each time step. Finally, we conduct several numerical experiments to verify the accuracy, efficiency, and robustness of the H$^T_N$ ansatz and the numerical discretizations.« less