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 $.