Topological susceptibility in finite temperature QCD with physical <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:mo stretchy="false">/</mml:mo><mml:mi>d</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math> domain-wall quarks
Yu‐Chih Chen, Ting-Wai Chiu, Tung‐Han Hsieh
Abstract
We perform hybrid Monte-Carlo simulation of lattice QCD with ${N}_{f}=2+1+1$ domain-wall quarks at the physical point, on the ${64}^{3}\ifmmode\times\else\texttimes\fi{}(64,20,16,12,10,8,6)$ lattices, each with three lattice spacings. The lattice spacings and the bare quark masses are determined on the ${64}^{4}$ lattices. The resulting gauge ensembles provide a basis for studying finite temperature QCD with ${N}_{f}=2+1+1$ domain-wall quarks at the physical point. In this paper, we determine the topological susceptibility of the QCD vacuum for $T>{T}_{c}\ensuremath{\sim}150\text{ }\text{ }\mathrm{MeV}$. The topological charge of each gauge configuration is measured by the clover charge in the Wilson flow at the same flow time in physical units, and the topological susceptibility ${\ensuremath{\chi}}_{t}(a,T)$ is determined for each ensemble with lattice spacing $a$ and temperature $T$. Using the topological susceptibility ${\ensuremath{\chi}}_{t}(a,T)$ of 15 gauge ensembles with three lattice spacings and different temperatures in the range $T\ensuremath{\sim}155--516\text{ }\text{ }\mathrm{MeV}$, we extract the topological susceptibility ${\ensuremath{\chi}}_{t}(T)$ in the continuum limit. To compare our results with others, we survey the continuum extrapolated ${\ensuremath{\chi}}_{t}(T)$ in lattice QCD with ${N}_{f}=2+1(+1)$ dynamical quarks at/near the physical point and discuss their discrepancies. Moreover, a detailed discussion on the reweighting method for the domain-wall fermion is presented.