First study of $N_f=2+1+1$ lattice QCD with physical domain-wall quarks

Using 10-16 units of Nvidia DGX-1, we have generated the first gauge ensemble for $N_f=2+1+1$ lattice QCD with physical $(u/d, s, c)$ domain-wall quarks, on the $ 64^4 $ lattice with lattice spacing $a \sim 0.064 $~fm ($ L > 4 $ fm, and $ M_\pi L > 3 $). The salient feature of this gauge ensemble is that the chiral symmetry is preserved to a high precision and all topological sectors are sampled ergodically. In this paper, we present the first results of the topological susceptibility and the ground-state hadron mass spectra.