We perform hybrid Monte-Carlo simulation of $N_f=2+1+1$ lattice QCD with domain-wall quarks at the physical point.
The simulation is carried out on the $L^3 \times T = 64^3 \times 64$ lattice with lattice spacing $a \sim 0.064$ fm ($L > 4$ fm, and $M_\pi L > 3$), using the Nvidia DGX-1 (8 V100 GPUs interconnected by the NVLink).
To attain the maximal chiral symmetry for a finite extent ($N_s=16$) in the fifth dimension, we use the optimal domain-wall fermion for the quark action, together with the exact one-flavor action for domain-wall fermion.
We outline the salient features of our simulation (without topology freezing, and small residual masses), together with the preliminary result of the masses of $\pi^\pm$, $K^\pm$, and $D^\pm$.