The formations of individual antiferromagnetic (AF) skyrmions and AF skyrmionic lattices on two-dimensional (2D) magnets with square crystal structure are debatable in recent years, for only an isolated skyrmion can be generated in such systems if classical Monte Carlo (CMC) method is employed. For the sake, we apply here an optimized quantum Monte Carlo approach to a 2D square magnet where the AF Heisenberg exchange (HE) and Dzyaloshinskii-Moriya (DM) interactions co-exist. Consequently, the computing program converges to the equilibrium states with appreciable computational speed, and the re...