Topological susceptibility in finite temperature QCD with physical (u=d;s;c) domain-wall quarks

We perform hybrid Monte-Carlo simulation of lattice QCD with Nf = 2 + 1 + 1 domain-wall quarks at the physical point, on the 643 x (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 644 lattices. The resulting gauge ensembles provide a basis for studying finite temperature QCD with Nf = 2 + 1 + 1 domain-wall quarks at the physical point. In this paper, we determine the topological susceptibility of the QCD vacuum for T > Tc-, 150 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 chi t(a; T) is determined for each ensemble with lattice spacing a and temperature T. Using the topological susceptibility chi t(a; T) of 15 gauge ensembles with three lattice spacings and different temperatures in the range T-, 155-516 MeV, we extract the topological susceptibility chi t(T) in the continuum limit. To compare our results with others, we survey the continuum extrapolated chi t(T) in lattice QCD with Nf = 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.