Researcher profile

Swagato Mukherjee

Swagato Mukherjee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 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.

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$.

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.

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.

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.

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.

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.

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$.

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$.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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}$.

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.

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.

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.