BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Swagato Mukherjee

Swagato Mukherjee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phnucl-thhep-thnucl-exhep-exDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

64works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

64 published item(s)

preprint2024arXiv

Analysis of Static Wilson Line Correlators from Lattice QCD at Finite Temperature with $T$-matrix Approach

We employ a thermodynamic $T$-matrix approach to study Wilson line correlators (WLCs) for a static quark-antiquark pair within the quark-gluon plasma (QGP). By refining earlier determined input parameters, self-consistent $T$-matrix calculations can approximately reproduce the WLCs from 2+1-flavor lattice-QCD (lQCD) computations, while also satisfying constraints on the bulk medium from the QGP equation of state. The main difference from earlier results is a reduced screening in the input potential toward higher temperatures. When elaborating the consequences of this scenario on transport properties we find a more pronounced temperature dependence of the low-momentum relaxation rate of heavy quarks in the QGP. This, in particular, leads to a rather weak increase of the (temperature-scaled) spatial diffusion coefficient, which turns out to be in fair agreement with the most recent pertinent lQCD data.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Generalized Parton Distributions from Lattice QCD with Asymmetric Momentum Transfer: Unpolarized Quarks

Traditionally, lattice QCD computations of generalized parton distributions (GPDs) have been carried out in a symmetric frame, where the transferred momentum is symmetrically distributed between the incoming and outgoing hadrons. However, such frames are inconvenient since they require a separate calculation for each value of the momentum transfer, increasing significantly the computational cost. In this work, by focusing on the quasi-distribution approach, we lay the foundation for faster and more effective lattice QCD calculations of GPDs exploiting asymmetric frames, with freedom in the transferred momentum distribution. An important ingredient of our approach is the Lorentz covariant parameterization of the matrix elements in terms of Lorentz-invariant amplitudes, which allows one to relate matrix elements in different frames. We also use this amplitude approach to propose a new definition of quasi-GPDs that is frame-independent and, more importantly, may lead to smaller power corrections in the matching relations to the light-cone GPDs. We demonstrate the efficacy of the formalism through numerical calculations using one ensemble of $N_f$=2+1+1 twisted mass fermions with a clover improvement. The value of the light-quark masses lead to a pion mass of about 260 MeV. Concentrating on the proton, and limiting ourselves to a vanishing longitudinal momentum transfer to the target, we extract the invariant amplitudes from matrix element calculations in both the symmetric and asymmetric frame, and obtain results for the twist-2 light-cone GPDs for unpolarized quarks, that is, $H$ and $E$.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

GPDs in asymmetric frames

It is often taken for granted that Generalized Parton Distributions (GPDs) are defined in the "symmetric" frame, where the transferred momentum is symmetrically distributed between the incoming/outgoing hadrons. However, such frames pose computational challenges for the lattice QCD practitioners. In these proceedings, we lay the foundation for lattice QCD calculations of GPDs in "asymmetric" frames, where the transferred momentum is not symmetrically distributed between the incoming/outgoing hadrons. The novelty of our work relies on the parameterization of the matrix elements in terms of Lorentz-invariant amplitudes, which not only helps in establishing relations between the said frames but also helps in isolating higher-twist contaminations. As an example, we focus on the unpolarized GPDs for spin-1/2 particles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Analytical structure of the equation of state at finite density: Resummation versus expansion in a low energy model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Heavy quark potential in quark-gluon Plasma: Deep neural network meets lattice quantum chromodynamics

Bottomonium states are key probes for experimental studies of the quark-gluon plasma (QGP) created in high-energy nuclear collisions. Theoretical models of bottomonium productions in high-energy nuclear collisions rely on the in-medium interactions between the bottom and antibottom quarks. The latter can be characterized by the temperature ($T$) dependent potential, with real ($V_R(T,r)$) and imaginary ($V_I(T,r)$) parts, as a function of the spatial separation ($r$). Recently, the masses and thermal widths of up to $3S$ and $2P$ bottomonium states in QGP were calculated using lattice quantum chromodynamics (LQCD). Starting from these LQCD results and through a novel application of deep neural network, here, we obtain $V_R(T,r)$ and $V_I(T,r)$ in a model-independent fashion. The temperature dependence of $V_R(T,r)$ was found to be very mild between $T\approx0-334$~MeV. For $T=151-334$~MeV, $V_I(T,r)$ shows a rapid increase with $T$ and $r$, which is much larger than the perturbation-theory-based expectations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lattice QCD Determination of the Bjorken-$x$ Dependence of Parton Distribution Functions at Next-to-next-to-leading Order

We report the first lattice QCD calculation of pion valence quark distribution with next-to-next-to-leading order perturbative matching correction, which is done using two fine lattices with spacings $a=0.04$ fm and $0.06$ fm and valence pion mass $m_π=300$ MeV, at boost momentum as large as $2.42$ GeV. As a crucial step to control the systematics, we renormalize the pion valence quasi distribution in the recently proposed hybrid scheme, which features a Wilson-line mass subtraction at large distances in coordinate space, and develop a procedure to match it to the $\overline{\rm MS}$ scheme. We demonstrate that the renormalization and the perturbative matching in Bjorken-$x$ space yield a reliable determination of the valence quark distribution for $0.03\lesssim x \lesssim 0.80$ with 5-20\% uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lattice QCD Equation of State for Nonvanishing Chemical Potential by Resumming Taylor Expansion

Taylor expansion in powers of baryon chemical potential ($μ_B$) is an oft-used method in lattice QCD to compute QCD thermodynamics for $μ_B>0$. Based only upon the few known lowest order Taylor coefficients, it is difficult to discern the range of $μ_B$ where such an expansion around $μ_B=0$ can be trusted. We introduce a resummation scheme for the Taylor expansion of the QCD equation of state in $μ_B$ that is based on the $n$-point correlation functions of the conserved current ($D_n$). The method resums the contributions of the first $N$ correlation function $D_1,\dots,D_N$ to the Taylor expansion of the QCD partition function to all orders in $μ_B$. We show that the resummed partition function is an approximation to the reweighted partition function at $μ_B\ne0$. We apply the proposed approach to high-statistics lattice QCD calculations using 2+1 flavors of Highly Improved Staggered Quarks with physical quark masses on $32^3\times8$ lattices and for temperatures $T\approx145$-176 MeV. We demonstrate that, as opposed to the Taylor expansion, the resummed version not only leads to improved convergence but also reflects the zeros of the resummed partition function and severity of the sign problem, leading to its eventual breakdown. We also provide a generalization of our scheme to include resummation of powers of temperature and quark masses in addition to $μ_B$, and show that the alternative expansion scheme of [S. Borsányi et al., Phys. Rev. Lett. 126, 232001 (2021).] is a special case of this generalized resummation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pion form factor and charge radius from Lattice QCD at physical point

We present our results on the electromagnetic form factor of pion over a wide range of $Q^2$ using lattice QCD simulations with Wilson-clover valence quarks and HISQ sea quarks. We study the form factor at the physical point with a lattice spacing $a=0.076$ fm. To study the lattice spacing and quark mass effects, we also present results for 300 MeV pion at two different lattice spacings $a=0.04$ and 0.06 fm. The lattice calculations at the physical quark mass appear to agree with the experimental results. Through fits to the form factor, we estimate the charge radius of pion for physical pion mass to be $\langle r_π^2 \rangle=0.42(2)~{\rm fm}^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searching for the QCD critical point along the pseudo-critical/freeze-out line using Padé-resummed Taylor expansions of cumulants of conserved charge fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, $T_{pc}\simeq 156.5$~MeV, we find the estimator of the radius of convergence to be $μ_B/T \gtrsim\ 3$ for strangeness-neutral matter. We also present results for the pole structure of the Padé approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Padé is consistent with not having a critical point at temperatures larger than $135~$MeV and a baryon chemical potential smaller than $μ_B/T \sim \ 2.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Static quark anti-quark interactions at non-zero temperature from lattice QCD

We study the interactions of a static quark antiquark pair at non-zero temperature using realistic 2+1 flavor lattice QCD calculations. The study consists of two parts: the first investigates the properties of Wilson line correlators in Coulomb gauge and compares to predictions of hard-thermal loop perturbation theory. As a second step we extract the spectral functions underlying the correlators using four conceptually different methods: spectral function fits, a HTL inspired fit for the correlation function, Padé rational approximation and the Bayesian BR spectral reconstruction. We find that our high statistics Euclidean lattice data are amenable to different hypotheses for the shapes of the spectral function and we compare the implications of each analysis method for the existence and properties of a well defined ground state spectral peak.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Static quark anti-quark interactions at non-zero temperature from lattice QCD

We present results on the in-medium interactions of static quark anti-quark pairs using realistic 2+1 HISQ flavor lattice QCD. Focus is put on the extraction of spectral information from Wilson line correlators in Coulomb gauge using four complementary methods. Our results indicate that on HISQ lattices, the position of the dominant spectral peak associated with the real-part of the interquark potential remains unaffected by temperature. This is in contrast to prior work in quenched QCD and we present follow up comparisons to newly generated quenched ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Taylor expansions and Padé approximants for cumulants of conserved charge fluctuations at non-vanishing chemical potentials

Using high statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we present results for low order cumulants of net baryon-number fluctuations at non-zero values of the baryon chemical potential. We calculate Taylor expansions for the pressure (zeroth order cumulant), net baryon-number density (first order cumulant) and the variance of the distribution on net-baryon number fluctuations (second order cumulant). We obtain series expansions from an eighth order expansion of the pressure and compare these to diagonal Padé approximants. This allows us to estimate the range of values for the baryon chemical potential in which these expansions are reliable. We find $μ_B/T\le 2.5$, $2.0$ and $1.5$ for the zeroth, first and second order cumulants, respectively. We furthermore, construct estimators for the radius of convergence of the Taylor series of the pressure. In the vicinity of the pseudo-critical temperature, $T_{pc}\simeq 156.5$ MeV, we find $μ_B/T \gtrsim\ 2.9$ at vanishing strangeness chemical potential and somewhat larger values for strangeness neutral matter. These estimates are temperature dependent and range from $μ_B/T \gtrsim\ 2.2$ at $T=135$ MeV to $μ_B/T\ \gtrsim\ 3.2$ at $T=165$ MeV. The estimated radius of convergences is the same for any higher order cumulant.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Conserved charge fluctuations at vanishing net-baryon density from Lattice QCD

We present here continuum extrapolated results for all 2nd order cumulants using the most recent results obtained by the HotQCD collaboration in (2+1)-flavor QCD. We constrain the applicability of various HRG models by presenting a detailed comparison of our results based on different sets of hadron spectra as well as with virial expansion based model calculations. A comparison with our lattice QCD results for 2nd order cumulants with models that parametrize repulsive interactions among baryons and anti-baryons in a hadron resonance gas through a single excluded volume parameter (EVHRG) is also shown.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlated Dirac Eigenvalues and Axial Anomaly in Chiral Symmetric QCD

We introduce novel relations between the derivatives [$\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$] of the Dirac eigenvalue spectrum [$ρ(λ, m_{l})$] with respect to the light sea quark mass ($m_{l}$) and the $(n+1)$-point correlations among the eigenvalues ($λ$) of the massless Dirac operator. Using these relations we present lattice QCD results for $\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$ ($n=1, 2, 3$) for $m_{l}$ corresponding to pion masses $m_π=160-55$ MeV, and at a temperature of about 1.6 times the chiral phase transition temperature. Calculations were carried out using (2+1) flavors of highly improved staggered quarks with the physical value of strange quark mass, three lattice spacings $a=0.12, 0.08, 0.06$ fm, and lattices having aspect ratios $4-9$. We find that $ρ(λ\to0, m_{l})$ develops a peaked structure. This peaked structure arises due to non-Poisson correlations within the infrared part of the Dirac eigenvalue spectrum, becomes sharper as $a\to0$, and its amplitude is proportional to $m_{l}^2$. We demonstrate that this $ρ(λ\to0,m_l)$ is responsible for the manifestations of axial anomaly in two-point correlation functions of light scalar and pseudoscalar mesons. After continuum and chiral extrapolations we find that axial anomaly remains manifested in two-point correlation functions of scalar and pseudoscalar mesons in the chiral limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bethe-Salpeter amplitudes of Upsilons

Based on lattice non-relativistic QCD (NRQCD) studies we present results for Bethe-Salpeter amplitudes for $Υ(1S)$, $Υ(2S)$ and $Υ(3S)$ in vacuum as well as in quark-gluon plasma. Our study is based on 2+1 flavor $48^3 \times 12$ lattices generated using the Highly Improved Staggered Quark (HISQ) action and with a pion mass of $161$ MeV. At zero temperature the Bethe-Salpeter amplitudes follow the expectations based on non-relativistic potential models. At non-zero temperatures, the interpretation of Bethe-Salpeter amplitudes turns out to be more nuanced, but consistent with our previous lattice QCD study of excited Upsilons in quark-gluon plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dirac Eigenvalue spectrum of $N_f$=2+1 QCD towards the chiral limit using HISQ fermions

We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number to determine the eigenvalue spectrum. Simulations of (2 + 1)-flavor QCD are performed using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ$ = 8 lattices with aspect ratios $N_σ/N_τ$ ranging from 5 to 7. The strange quark mass is fixed to its physical value $m_{s}^{\rm phy}$, and the light quark masses $m_{l}$ are varied from $m_{s}^{\rm phy}/40$ to $m_{s}^{\rm phy}/160$ which correspond to pion mass $m_π$ ranging from 110 MeV to 55 MeV in the continuum limit. We compute the chiral condensate and $χ_π - χ_δ$ through the eigenvalue spectrum obtained from the the eigenvalue filtering method. We compare these results with those obtained from a direct calculation of the observables which involves inversions of the fermion matrix using the stochastic "noise vector" method. We find that these approaches yield consistent results. Furthermore, we also investigate the quark mass and temperature dependences of the Dirac eigenvalue density at zero eigenvalues to gain more insights about the $U_A(1)$ symmetry breaking in QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Excited bottomonia in quark-gluon plasma from lattice QCD

We present the first lattice QCD study of up to $3S$ and $2P$ bottomonia at non-zero temperatures. Correlation functions of bottomonia were computed using novel bottomonium operators and a variational technique, within the lattice non-relativistic QCD framework. We analyzed the bottomonium correlation functions based on simple physically-motivated spectral functions. We found evidence of sequential in-medium modifications, in accordance with the sizes of the bottomonium states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Higher order cumulants of net baryon-number distributions at non-zero $μ_B$

Using recent results on higher order cumulants of conserved charge fluctuations from lattice QCD, we construct mean, variance, skewness, kurtosis, hyper-skewness and hyper-kurtosis of net-baryon number distributions for small baryon chemical potentials $μ_B$. For the strangeness neutral case ($μ_S=0$) at fixed ratio of electric charge to baryon number density ($\frac{n_Q}{n_B}=0.4$), which is appropriate for a comparison with heavy ion collisions, we present results for $κ_B σ_B^2$, $S_B σ_B^3/M_B$, $κ^{H}_{B}σ_{B}^4$ and $S^{H}_{B}σ^5_{B}/M_{B}$ on the crossover line for the chiral transition, $T_{pc}(μ_B)$. Continuum extrapolations for this pseudo-critical transition line have recently been reported by HotQCD up to baryon chemical potentials $μ_B\simeq 300$ MeV [arXiv:1812.08235]. These cumulant ratios are of direct relevance for comparisons with corresponding ratios measured by STAR in the BES-I and II runs at beam energies $\sqrt{s_{NN}}\ge 20$ GeV. In particular, we point out that recent high statistics results on skewness and kurtosis of net-baryon number distributions obtained by STAR at $\sqrt{s_{NN}} = 54.4$ GeV put strong constraints on freeze-out parameters and are consistent with predictions from thermal QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Meson masses in external magnetic fields with HISQ fermions

We studied the temporal correlation function of mesons in the pseudo-scalar channel in (2+1)-flavor QCD in the presence of external magnetic fields at zero temperature. The simulations were performed on $32^3 \times 96$ lattices using the Highly Improved Staggered Quarks (HISQ) action with $m_π \approx $ 230 MeV. The strength of magnetic fields $|eB|$ ranges from 0 to around 3.3 GeV$^2$ ($\sim 60 m_π^2$). We found that the masses of neutral pseudo-scalar particles, e.g. neutral pion and kaon, monotonouslly decrease as the magnetic field grows and then saturate at a nonzero value. It is observed that heavier neutral pseudo-scalars are less affected by magnetic fields. Moreover, we found a non-monotonous behavior of charged pion and kaon mass in magnetic field for the first time. In the case of small magnetic field (0 $\leq~|eB| \lesssim$ 0.3 GeV$^2~\sim 6m_π^2$ ) the mass of charged pseudo-scalar grows with magnetic field and can be well described by the Lowest Landau Level approximation, while for $|eB|$ larger than 0.3 GeV$^2$ the mass starts to decrease. The possible connection between $|eB|$ dependences of neutral pion mass and the decreasing behavior of pseudo-critical temperature in magnetic field is discussed. Due to the nonzero value of neutral pion mass our simulation indicates that the superconducting phase of QCD does not exist in the current window of magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pion valence quark PDF from lattice QCD

We present lattice results on the valence-quark structure of the pion using a coordinate space method within the framework of Large Momentum Effective Theory (LaMET). In this method one relies on the matrix elements of a Euclidean correlator in boosted hadronic states, which have an operator product expansion at short distance that allows us to extract the moments of PDFs. We renormalize the Euclidean correlator by forming the reduced Ioffe-time distribution (rITD), and reconstruct the second and fourth moments of the pion PDF by taking into account of QCD evolution effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Skewness, kurtosis and the 5th and 6th order cumulants of net baryon-number distributions from lattice QCD confront high-statistics STAR data

We present new results on up to $6^{th}$ order cumulants of net baryon-number fluctuations at small values of the baryon chemical potential, $μ_B$, obtained in lattice QCD calculations with physical values of light and strange quark masses. Representation of the Taylor expansions of higher order cumulants in terms of the ratio of the two lowest order cumulants, $M_B/σ_B^2=χ_1^B(T,μ_B)/χ_2^B(T,μ_B)$, allows for a parameter free comparison with data on net proton-number cumulants obtained by the STAR Collaboration in the Beam Energy Scan at RHIC. We show that recent high statistics data on skewness and kurtosis ratios of net proton-number distributions, obtained at beam energy $\sqrt{s_{_{NN}}}=54.4$ GeV, agree well with lattice QCD results on cumulants of net baryon-number fluctuations close to the pseudo-critical temperature, $T_{pc}(μ_B)$, for the chiral transition in QCD. We also present first results from a next-to-leading order expansion of $5^{th}$ and $6^{th}$ order cumulants on the line of pseudo-critical temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal Broadening of Bottomonia: Lattice Non-Relativistic QCD with Extended Operators

We present lattice non-relativistic QCD calculations of bottomonium correlation functions at temperatures $T \simeq 150-350$ MeV. The correlation functions were computed using extended bottomonium operators, and on background gauge-field configurations for 2+1-flavor QCD having physical kaon and nearly-physical pion masses. We analyzed these correlation functions based on simple theoretically-motivated parameterizations of the corresponding spectral functions. The results of our analyses are compatible with significant in-medium thermal broadening of the ground state S- and P-wave bottomonia.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universality driven analytic structure of QCD crossover: radius of convergence and QCD critical point

Recent lattice QCD calculations show strong indications that the crossover of QCD at zero baryon chemical potential ($μ_B$) is a remnant of the second order chiral phase transition. The non-universal parameters needed to map temperature $T$ and $μ_B$ to the universal properties of the second order chiral phase transition were determined by lattice QCD calculations. Motivated by these advances, first, we discuss the analytic structure of the partition function -- the so-called Yang-Lee edge singularity -- in the QCD crossover regime, solely based on universal properties. Then, utilizing the lattice calculated non-universal parameters, we map this singularity to the real $T$ and complex $μ_B$ plane, in order to find the radius of convergence for a Taylor series expansion of QCD partition function around $μ_B=0$ in the QCD crossover regime. Our most important findings are: (i) An universality-based estimate of the radius of convergence around $μ_B=0$; (ii) Universality and lattice QCD based constraints on the location of the QCD critical point in the $T-μ_B$ plane.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Chiral crossover in QCD at zero and non-zero chemical potentials

We present results for pseudo-critical temperatures of QCD chiral crossovers at zero and non-zero values of baryon ($B$), strangeness ($S$), electric charge ($Q$), and isospin ($I$) chemical potentials $μ_{X=B,Q,S,I}$. The results were obtained using lattice QCD calculations carried out with two degenerate up and down dynamical quarks and a dynamical strange quark, with quark masses corresponding to physical values of pion and kaon masses in the continuum limit. By parameterizing pseudo-critical temperatures as $ T_c(μ_X) = T_c(0) \left[ 1 -κ_2^{X}(μ_{X}/T_c(0))^2 -κ_4^{X}(μ_{X}/T_c(0))^4 \right] $, we determined $κ_2^X$ and $κ_4^X$ from Taylor expansions of chiral observables in $μ_X$. We obtained a precise result for $T_c(0)=(156.5\pm1.5)\;\mathrm{MeV}$. For analogous thermal conditions at the chemical freeze-out of relativistic heavy-ion collisions, i.e., $μ_{S}(T,μ_{B})$ and $μ_{Q}(T,μ_{B})$ fixed from strangeness-neutrality and isospin-imbalance, we found $κ_2^B=0.012(4)$ and $κ_4^B=0.000(4)$. For $μ_{B}\lesssim300\;\mathrm{MeV}$, the chemical freeze-out takes place in the vicinity of the QCD phase boundary, which coincides with the lines of constant energy density of $0.42(6)\;\mathrm{GeV/fm}^3$ and constant entropy density of $3.7(5)\;\mathrm{fm}^{-3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Chiral phase transition in (2 + 1)-flavor QCD

The chiral phase transition temperature $T_{c}^{0}$ is a fundamental quantity of QCD. To determine this quantity we have performed simulations of (2 + 1)-flavor QCD using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ=6, 8$ and 12 lattices with aspect ratios $N_σ/N_τ$ ranging from 4 to 8. In our simulations the strange quark mass is fixed to its physical value $m_{s}^{\rm{phy}}$, and the values of two degenerate light quark masses $m_{l}$ are varied from $m_{s}^{\rm{phy}}/20$ to $m_{s}^{\rm{phy}}/160$ which correspond to a Goldstone pion mass $m_π$ ranging from 160 MeV to 55 MeV in the continuum limit. By investigating the light quark mass dependence and the volume dependence of various chiral observables, e.g. chiral susceptibilities and Binder cumulants, no evidence for a first order phase transition in our current quark mass window is found. Two estimators $T_{60}$ and $T_δ$ are proposed to extract the chiral phase transition temperature $T_{c}^{0}$ in the chiral and continuum limit and our current estimate for $T_{c}^{0}$ is $132_{-6}^{+3}$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Meson Screening Masses in (2+1)-Flavor QCD

We present lattice QCD results for mesonic screening masses in the temperature range 140 MeV $\lesssim T \lesssim$ 2500 MeV. Our calculations were carried out using (2+1)-flavors of the Highly Improved Staggered Quark (HISQ) action, with a physical value for the strange quark mass and two values of the light quark mass corresponding to pion masses of 160 MeV and 140 MeV. Continuum-extrapolated results were obtained using calculations with a variety of lattice spacings corresponding to temporal lattice extents $N_τ= 6 - 16$. We discuss the implications of these results for the effective restoration of various symmetries in the high temperature phase of QCD, as well as the approach toward the perturbative limit.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Thermal modifications of charmonia and bottomonia from spatial correlation functions

We present our study on the thermal modifications of charmonia and bottomonia from spatial correlation functions at zero and nonzero momenta in quenched QCD. To accommodate the heavy quarks on the lattice we performed simulations on very fine lattices at a fixed beta value corresponding to a lattice spacing $a^{-1}=22.8$ GeV on $192^3\times32$, $192^3\times48$, $192^3\times56$, $192^3\times64$ and $192^3\times96$ lattices using clover-improved Wilson fermions. These lattices correspond to temperatures of $2.25 T_c$, $1.50 T_c$, $1.25 T_c$, $1.10 T_c$ and $0.75 T_c$. To increase the signal to noise ratio in the axial-vector and scalar channels we used multi-sources for the measurement of spatial correlation functions. By investigating on the differences between spatial and temporal correlators as well as the temperature dependence of screening masses we will discuss the thermal effects in different channels of quarkonium states. Besides this the dispersion relation of the screening mass at different momenta is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Charm Degrees of Freedom in Quark Gluon Plasma

Lattice QCD studies on fluctuations and correlations of charm quantum number have established that deconfinement of charm degrees of freedom sets in around the chiral crossover temperature, $T_c$, i.e. charm degrees of freedom carrying fractional baryonic charge start to appear. By reexamining those same lattice QCD data we show that, in addition to the contributions from quark-like excitations, the partial pressure of charm degrees of freedom may still contain significant contributions from open-charm meson and baryon-like excitations associated with integral baryonic charges for temperatures up to $1.2~ T_c$. Charm quark-quasiparticles become the dominant degrees of freedom for temperatures $T>1.2~ T_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The curvature of the freeze-out line in heavy ion collisions

We calculate the mean and variance of net-baryon number and net-electric charge distributions from Quantum Chromodynamics (QCD) using a next-to-leading order Taylor expansion in terms of temperature and chemical potentials. We compare these expansions with experimental data from STAR and PHENIX, determine the freeze-out temperature in the limit of vanishing baryon chemical potential, and, for the first time, constrain the curvature of the freeze-out line through a direct comparison between experimental data on net-charge fluctuations and a QCD calculation. We obtain a bound on the curvature coefficient, kappa_2^f < 0.011, that is compatible with lattice QCD results on the curvature of the QCD transition line.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The topological structures in strongly coupled QGP with chiral fermions on the lattice

The nature of chiral phase transition for two flavor QCD is an interesting but unresolved problem. One of the most intriguing issues is whether or not the anomalous U(1) symmetry in the flavor sector is effectively restored along with the chiral symmetry. This may determine the universality class of the chiral phase transition. Since the physics near the chiral phase transition is essentially non-perturbative, we employ first principles lattice techniques to address this issue. We use overlap fermions, which have exact chiral symmetry on the lattice, to probe the anomalous U(1) symmetry violation of 2+1 flavor dynamical QCD configurations with domain wall fermions. The latter also optimally preserves chiral and flavor symmetries on the lattice, since it is known that the remnant chiral symmetry of the light quarks influences the scaling of the chiral condensate in the crossover transition region. We observe that the anomalous U(1) is not effectively restored in the chiral crossover region. We perform a systematic study of the finite size and cut-off effects since the signals of U(1) violation are sensitive to it. We also provide a glimpse of the microscopic topological structures of the QCD medium that are responsible for the strongly interacting nature of the quark gluon plasma phase. We study the effect of these microscopic constituents through our first calculations for the topological susceptibility of QCD at finite temperature, which could be a crucial input for the equation of state for anomalous hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal off-equilibrium scaling of critical cumulants in the QCD phase diagram

Exploiting the universality between the QCD critical point and the three dimensional Ising model, closed form expressions derived (arXiv:1506.00645 ) for non-equilibrium critical cumulants on the crossover side of the critical point reveal that they can differ both in magnitude and sign from equilibrium expectations. We demonstrate here that key elements of the Kibble-Zurek framework of non-equilibrium phase transitions can be employed to describe the dynamics of these critical cumulants. Our results suggest that observables sensitive to critical dynamics in heavy-ion collisions should be expressible as universal scaling functions, thereby providing powerful model independent guidance in searches for the QCD critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A stochastic approach to the reconstruction of spectral functions in lattice QCD

We present a Stochastic Optimization Method (SOM) for the reconstruction of the spectral functions (SPFs) from Euclidean correlation functions. In this approach the SPF is parameterized as a sum of randomly distributed boxes. By varying the width, location and height of the boxes stochastically an optimal SPF can be obtained. Using this approach we reproduce mock SPFs fairly well, which contain sharp resonance peaks, transport peaks and continuum spectra. We also analyzed the charmonium correlators obtained from $N_τ$=96, 48, 32 lattices using SOM and found similar conclusion on the dissociation temperatures of charmonium ground states as that obtained using the Maximum Entropy Method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Conserved Charge Fluctuations from Lattice QCD and the Beam Energy Scan

We discuss the next-to-leading order Taylor expansion of ratios of cumulants of net-baryon number fluctuations. We focus on the relation between the skewness ratio, $S_Bσ_B = χ_3^B/χ_1^B$, and the kurtosis ratio, $κ_Bσ_B^2 =χ_4^B/χ_2^B$. We show that differences in these two cumulant ratios are small for small values of the baryon chemical potential. The next-to-leading order correction to $κ_Bσ_B^2$ however is approximately three times larger than that for $S_Bσ_B$. The former thus drops much more rapidly with increasing beam energy, $\sqrt{s_{NN}}$. We argue that these generic patterns are consistent with current data on cumulants of net-proton number fluctuations measured by the STAR Collaboration at $\sqrt{s_{NN}}\ge 19.6$~GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diagonal and off-diagonal quark number susceptibilities at high temperatures

We present continuum extrapolated lattice QCD results for up to fourth-order diagonal and off-diagonal quark number susceptibilities in the high temperature region of 300-700 MeV. Lattice QCD calculations are performed using 2+1 flavors of highly improved staggered quarks with nearly physical quark masses and at four different lattice spacings. Comparisons of our results with recent weak coupling perturbative calculations yield reasonably good agreements for the entire temperature range.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

In-medium modifications of open and hidden strange-charm mesons from spatial correlation functions

We calculate spatial correlation functions of in-medium mesons consisting of strange--anti-strange, strange--anti-charm and charm--anti-charm quarks in (2+1)-flavor lattice QCD using the highly improved staggered quark action. A comparative study of the in-medium modifications of mesons with different flavor contents is performed. We observe significant in-medium modifications for the $ϕ$ and $D_s$ meson channels already at temperatures around the chiral crossover region. On the other hand, for the $J/ψ$ and $η_c$ meson channels in-medium modifications remain relatively small around the chiral crossover region and become significant only above 1.3 times the chiral crossover temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice QCD Thermodynamics with Physical Quark Masses

Over the past few years new physics methods and algorithms as well as the latest supercomputers have enabled the study of the QCD thermodynamic phase transition using lattice gauge theory numerical simulations with unprecedented control over systematic errors. This is largely a consequence of the ability to perform continuum extrapolations with physical quark masses. Here we review recent progress in lattice QCD thermodynamics, focussing mainly on results that benefit from the use of physical quark masses: the crossover temperature, the equation of state, and fluctuations of the quark number susceptibilities. In addition, we place a special emphasis on calculations that are directly relevant to the study of relativistic heavy ion collisions at RHIC and the LHC.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microscopic Origin of \boldmath{$U_A(1)$} Symmetry Violation in the High Temperature Phase of QCD

We investigate the low-lying eigenmodes of the Dirac matrix with the aim to gain more insight into the temperature dependence of the anomalous $U_A(1)$ symmetry. We use the overlap operator to probe dynamical QCD configurations generated with (2+1)-flavors of highly improved staggered quarks. We find no evidence of a gap opening up in the infrared region of the eigenvalue spectrum even at $1.5\,T_c$, $T_c$ being the chiral crossover temperature. Instead, we observe an accumulation of near-zero eigenmodes. We argue that these near-zero eigenmodes are primarily responsible for the anomalous breaking of the axial symmetry still being effective. At $1.5\,T_c$, these near-zero eigenmodes remain localized and their distribution is consistent with the dilute instanton gas picture. At this temperature, the average size of the instantons is $0.223(8)\,\text{fm}$ and their density is $0.147(7)\,\text{fm}^{-4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real time evolution of non-Gaussian cumulants in the QCD critical regime

We derive a coupled set of equations that describe the non-equilibrium evolution of cumulants of critical fluctuations for space-time trajectories on the cross-over side of the QCD phase diagram. In particular, novel expressions are obtained for the non-equilibrium evolution of non-Gaussian Skewness and Kurtosis cumulants. Utilizing a simple model of the space-time evolution of a heavy- ion collision, we demonstrate that, depending on the relaxation rate of critical fluctuations, Skewness and Kurtosis can differ significantly in magnitude as well as in sign from equilibrium expectations. Memory effects are important and shown to persist even for trajectories that skirt the edge of the critical regime. We use phenomenologically motivated parameterizations of freeze-out curves, and of the beam energy dependence of the net baryon chemical potential, to explore the implications of our model study for the critical point search in heavy-ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The U(1)A anomaly in high temperature QCD with chiral fermions on the lattice

The magnitude of the $U_A(1)$ symmetry breaking is expected to affect the nature of $N_f=2$ QCD chiral phase transition. The explicit breaking of chiral symmetry due to realistic light quark mass is small, so it is important to use chiral fermions on the lattice to understand the effect of $U_A(1)$ near the chiral crossover temperature, $T_c$. We report our latest results for the eigenvalue spectrum of 2+1 flavour QCD with dynamical Mobius domain wall fermions at finite temperature probed using the overlap operator on $32^3\times 8$ lattice. We check how sensitive the low-lying eigenvalues are to the sea-light quark mass. We also present a comparison with the earlier independent results with domain wall fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermodynamics of strong-interaction matter from Lattice QCD

We review results from lattice QCD calculations on the thermodynamics of strong-interaction matter with emphasis on input these calculations can provide to the exploration of the phase diagram and properties of hot and dense matter created in heavy ion experiments. This review is organized as follows: 1) Introduction, 2) QCD thermodynamics on the lattice, 3) QCD phase diagram at high temperature, 4) Bulk thermodynamics, 5) Fluctuations of conserved charges, 6) Transport properties, 7) Open heavy flavors and heavy quarkonia, 8) QCD in external magnetic fields, 9) Summary.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Additional Strange Hadrons from QCD Thermodynamics and Strangeness Freeze-out in Heavy Ion Collisions

We compare lattice QCD results for appropriate combinations of net strangeness fluctuations and their correlations with net baryon number fluctuations with predictions from two hadron resonance gas (HRG) models having different strange hadron content. The conventionally used HRG model based on experimentally established strange hadrons fails to describe the lattice QCD results in the hadronic phase close to the QCD crossover. Supplementing the conventional HRG with additional, experimentally uncharted strange hadrons predicted by quark model calculations and observed in lattice QCD spectrum calculations leads to good descriptions of strange hadron thermodynamics below the QCD crossover. We show that the thermodynamic presence of these additional states gets imprinted in the yields of the ground-state strange hadrons leading to a systematic 5--8 MeV decrease of the chemical freeze-out temperatures of ground-state strange baryons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

HISQ inverter on Intel Xeon Phi and NVIDIA GPUs

The runtime of a Lattice QCD simulation is dominated by a small kernel, which calculates the product of a vector by a sparse matrix known as the "Dslash" operator. Therefore, this kernel is frequently optimized for various HPC architectures. In this contribution we compare the performance of the Intel Xeon Phi to current Kepler-based NVIDIA Tesla GPUs running a conjugate gradient solver. By exposing more parallelism to the accelerator through inverting multiple vectors at the same time we obtain a performance 250 GFlop/s on both architectures. This more than doubles the performance of the inversions. We give a short overview of both architectures, discuss some details of the implementation and the effort required to obtain the achieved performance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The equation of state in (2+1)-flavor QCD

We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent $N_τ=6,~8,~10$, and $12$. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range $(130-400)$ MeV. We compare our results with previous calculations, and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, $145~ {\rm MeV} \leq T \leq 163$ MeV, defined by the chiral transition, is $ε_c=(0.18-0.5)~{\rm GeV}/{\rm fm}^3$, $i.e.$, $(1.2-3.1)\ ε_{\rm nuclear}$. At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters $r_0$ from the static quark potential and $w_0$ from the gradient flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The melting and abundance of open charm hadrons

Ratios of cumulants of conserved net charge fluctuations are sensitive to the degrees of freedom that are carriers of the corresponding quantum numbers in different phases of strong interaction matter. Using lattice QCD with 2+1 dynamical flavors and quenched charm quarks we calculate second and fourth order cumulants of net charm fluctuations and their correlations with other conserved charges such as net baryon number, electric charge and strangeness. Analyzing appropriate ratios of these cumulants we probe the nature of charmed degrees of freedom in the vicinity of the QCD chiral crossover region. We show that for temperatures above the chiral crossover transition temperature, charmed degrees of freedom can no longer be described by an uncorrelated gas of hadrons. This suggests that the dissociation of open charm hadrons and the emergence of deconfined charm states sets in just near the chiral crossover transition. Till the crossover region we compare these lattice QCD results with two hadron resonance gas models --including only the experimentally established charmed resonances and also including additional states predicted by quark model and lattice QCD calculations. This comparison provides evidence for so far unobserved charmed hadrons that contribute to the thermodynamics in the crossover region.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD chiral transition, $\ua$ symmetry and the Dirac spectrum using domain wall fermions

We report on a study of the finite-temperature QCD transition region for temperatures between 139 and 196 MeV, with a pion mass of 200 MeV and two space-time volumes: $24^3\times8$ and $32^3\times8$, where the larger volume varies in linear size between 5.6 fm (at T=139 MeV) and 4.0 fm (at T=195 MeV). These results are compared with the results of an earlier calculation using the same action and quark masses but a smaller, $16^3\times8$ volume. The chiral domain wall fermion formulation with a combined Iwasaki and dislocation suppressing determinant ratio gauge action are used. This lattice action accurately reproduces the $\sua$ and $\ua$ symmetries of the continuum. Results are reported for the chiral condensates, connected and disconnected susceptibilities and the Dirac eigenvalue spectrum. We find a pseudo-critical temperature, $T_c$, of approximately 165 MeV consistent with previous results and strong finite volume dependence below $T_c$. Clear evidence is seen for $\ua$ symmetry breaking above $T_c$ which is quantitatively explained by the measured density of near-zero modes in accordance with the dilute instanton gas approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD phase transition with physical-mass, chiral quarks

We report on the first lattice calculation of the QCD phase transition using chiral fermions at physical values of the quark masses. This calculation uses 2+1 quark flavors, spatial volumes between (4 fm$)^3$ and (11 fm$)^3$ and temperatures between 139 and 196 MeV . Each temperature was calculated using a single lattice spacing corresponding to a temporal Euclidean extent of $N_t=8$. The disconnected chiral susceptibility, $χ_{\rm disc}$ shows a pronounced peak whose position and height depend sensitively on the quark mass. We find no metastability in the region of the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm. Each result is strong evidence that the QCD ``phase transition'' is not first order but a continuous cross-over for $m_π=135$ MeV. The peak location determines a pseudo-critical temperature $T_c = 155(1)(8)$ MeV. Chiral $SU(2)_L\times SU(2)_R$ symmetry is fully restored above 164 MeV, but anomalous $U(1)_A$ symmetry breaking is non-zero above $T_c$ and vanishes as $T$ is increased to 196 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chiral phase transition of $N_f$=2+1 QCD with the HISQ action

We present studies of universal properties of the chiral phase transition in $N_f$=2+1 QCD based on the simulations using Highly Improved Staggered fermions on lattices with temporal extent $N_τ$=6. We analyze the quark mass and volume dependence of the chiral condensates and chiral susceptibilities in QCD with two degenerate light quarks and a strange quark. The strange quark mass is chosen to be fixed to its physical value ($m^{phy}_s$) and five values of light quark masses ($m_l$) that are varied in the interval 1/20$\gtrsim m_l/m^{phy}_s \gtrsim$1/80. Here various quark masses correspond to pseudo Goldstone pion masses ranging from about 160 MeV to about 80 MeV. The O(N) scaling of chiral observables and the influence of universal scaling on physical observables in the region of physical quark mass values are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Deconfinement of strangeness and freeze-out from charge fluctuations

We use Lattice QCD calculations of fluctuations and correlations of various conserved charges to show that the deconfinement of strangeness takes place in the chiral crossover region of QCD; however, inside the quark-gluon plasma strange quarks remain strongly interacting at least up to temperatures twice the QCD crossover temperature. Further, we discuss how the freeze-out parameters of heavy-ion collisions can be determined in a model-independent way through direct comparisons between experimentally measured higher order cumulants of conserved charges and corresponding Lattice QCD calculations. Utilizing the preliminary data from the STAR and PHENIX experiments we illustrate this method. Although, the Lattice QCD based determinations of the freeze-out parameters utilizing data sets of different experiments and different observables are currently not consistent with each other, it is tantalizing to see that all the observed freeze-out parameters lie very close to the chiral/deconfinement crossover region of QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Investigation of the $U_A(1)$ in high temperature QCD on the lattice

In this project we study the effect of the $U_A(1)$ anomaly for (2+1)-flavour QCD at high temperature. We apply the overlap operator as a tool to probe the topological properties of gauge field configurations which have been generated within the Highly Improved Staggered Quark (HISQ) discretization scheme on lattices of size $32^3\times 8$ with $m_l/m_s=1/20$, commonly used for the study of QCD thermodynamics. Although we have at present, only results for one value of the quark masses and thus cannot monitor the change of the eigenvalue distributions with the light quark mass, the distribution of the low-lying eigenvalues of the overlap operator suggests that the $U_A(1)$ is not restored effectively even at 1.5 times the pseudo critical temperature. The corresponding low-lying eigenmodes show localization properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quark number susceptibilities at high temperatures

We calculate second and fourth order quark number susceptibilities for 2+1 flavor QCD in the high temperature region. In our study we use two improved staggered fermion formulations, namely the highly improved staggered quark formulation, and the so-called p4 formulation, as well as several lattice spacings. Second order quark number susceptibilities are calculated using both improved staggered fermion formulations, and we show that in the continuum limit the two formulations give consistent results. The fourth order quark number susceptibilities are studied only using the p4 formulation and at non-zero lattice spacings. We compare our results on quark number susceptibilities with recent weak coupling calculations, and find that these agree reasonably well with the lattice calculations within the estimated uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strangeness at high temperatures: from hadrons to quarks

Appropriate combinations of up to fourth order cumulants of net strangeness fluctuations and their correlations with net baryon number and electric charge fluctuations, obtained from lattice QCD calculations, have been used to probe the strangeness carrying degrees of freedom at high temperatures. For temperatures up to the chiral crossover separate contributions of strange mesons and baryons can be well described by an uncorrelated gas of hadrons. Such a description breaks down in the chiral crossover region, suggesting that the deconfinement of strangeness takes place at the chiral crossover. On the other hand, the strangeness carrying degrees of freedom inside the quark gluon plasma can be described by a weakly interacting gas of quarks only for temperatures larger than twice the chiral crossover temperature. In the intermediate temperature window these observables show considerably richer structures, indicative of the strongly interacting nature of the quark gluon plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model

We calculate the quadratic fluctuations of net baryon number, electric charge and strangeness as well as correlations among these conserved charges in (2+1)-flavor lattice QCD at zero chemical potential. Results are obtained using calculations with tree level improved gauge and the highly improved staggered quark (HISQ) actions with almost physical light and strange quark masses at three different values of the lattice cut-off. Our choice of parameters corresponds to a value of 160 MeV for the lightest pseudo scalar Goldstone mass and a physical value of the kaon mass. The three diagonal charge susceptibilities and the correlations among conserved charges have been extrapolated to the continuum limit in the temperature interval 150 MeV <T < 250 MeV. We compare our results with the hadron resonance gas (HRG) model calculations and find agreement with HRG model results only for temperatures T<= 150 MeV. We observe significant deviations in the temperature range 160 MeV < T < 170 MeV and qualitative differences in the behavior of the three conserved charge sectors. At T < 160 MeV quadratic net baryon number fluctuations in QCD agree with HRG model calculations while, the net electric charge fluctuations in QCD are about 10% smaller and net strangeness fluctuations are about 20% larger. These findings are relevant to the discussion of freeze-out conditions in relativistic heavy ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Freeze-out Conditions in Heavy Ion Collisions from QCD Thermodynamics

We present a determination of chemical freeze-out conditions in heavy ion collisions based on ratios of cumulants of net electric charge fluctuations. These ratios can reliably be calculated in lattice QCD for a wide range of chemical potential values by using a next-to-leading order Taylor series expansion around the limit of vanishing baryon, electric charge and strangeness chemical potentials. From a computation of up to fourth order cumulants and charge correlations we first determine the strangeness and electric charge chemical potentials that characterize freeze-out conditions in a heavy ion collision and confirm that in the temperature range 150 MeV < T < 170 MeV the hadron resonance gas model provides good approximations for these parameters that agree with QCD calculations on the (5-15)% level. We then show that a comparison of lattice QCD results for ratios of up to third order cumulants of electric charge fluctuations with experimental results allows to extract the freeze-out baryon chemical potential and the freeze-out temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Signatures of charmonium modification in spatial correlation functions

We study spatial correlation functions of charmonium in 2+1 flavor QCD using an improved staggered formulation. Contrary to the temporal correlation functions the spatial correlation functions exhibit a strong temperature dependence above the QCD transition temperature. Above this temperature they are sensitive to temporal boundary conditions. Both features become significant at a temperature close to 1.5 Tc and suggest corresponding modifications of charmonium spectral functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The chiral transition and U(1)_A symmetry restoration from lattice QCD using Domain Wall Fermions

We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_τ= 8 and spatial extent N_σ= 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_π= 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

QCD Phase Transition in a Strong Magnetic Background

We investigate the properties of the deconfining/chiral restoring transition for two flavor QCD in presence of a uniform background magnetic field. We adopt a standard staggered discretization of the fermion action and a lattice spacing of the order of 0.3 fm. We explore different values of the bare quark mass, corresponding to pion masses in the range 200 - 480 MeV, and magnetic fields up to |e|B ~ 0.75 GeV^2. The deconfinement and chiral symmetry restoration temperatures remain compatible with each other and rise very slightly (< 2 % for our largest magnetic field) as a function of the magnetic field. On the other hand, the transition seems to become sharper as the magnetic field increases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Universal O(N) scaling and the chiral critical line in (2+1)-flavor QCD with small chemical potentials

We show that for small values of the chemical potential the curvature of the phase transition line can be deduced from an analysis of scaling properties of the chiral condensate and its susceptibilities. We make use of a recent analysis of the magnetic equation of state in (2+1)-flavor QCD where a connection between the QCD parameters and the universal scaling fields could be established. The remaining dependence of the reduced temperature on the chemical potential can be fixed by an analysis of a mixed susceptibility, obtained from a derivative with respect to quark mass and chemical potential. We extract this dependence which describes the curvature of the phase transition line, at two values of the cut-off, $aT=1/4$ and $1/8$. We find that cut-off effects are small for the curvature parameter and determine the transition line in the chiral limit to leading order in the light quark chemical potential. We obtain $T_c(μ_B)/T_c(0) = 1 - 0.00656(66) (μ_B/T)^2 +{\cal O}(μ_B^4)$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Lattice QCD equation of state : improving the differential method

We propose an improvement of the differential method for the computation of the equation of state of QCD from lattice simulations. In contrast to the earlier differential method our technique yields positive pressure for all temperatures including in the transition region. Employing it on temporal lattices of 8, 10 and 12 sites and by extrapolating to zero lattice spacing we obtained the pressure, energy density, entropy density, specific heat and speed of sound in quenched QCD for 0.9 < T/Tc < 3. A comparison of our results is made with those from the dimensional reduction approach and a conformal symmetric theory at high-temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Screening masses of mesons in 2+1 flavour QCD

We present results for screening masses of light and strange mesons in 2+1 flavour QCD using improved (p4fat3) staggered fermions on 6x24^3 lattices. We have studied the screening masses of scalar, pseudo-scalar, vector and axial-vector mesons along the line of constant physics, determined by a pion mass ~220 MeV and a kaon mass ~500 MeV. In order to investigate the cut-off and volume dependencies we have also performed studies of the meson screening correlators in the non-interacting theory using the p4 and the standard staggered discretizations.

0 citations0 reviewsNo code linkedOpen access