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.