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.