  The family of monolayer two-dimensional (2D) materials hosts a wide range of
interesting phenomena, including superconductivity, charge density waves,
topological states and ferromagnetism, but direct evidence for
antiferromagnetism in the monolayer has been lacking. Nevertheless,
antiferromagnets have attracted enormous interest recently in spintronics due
to the absence of stray fields and their terahertz resonant frequency. Despite
the great advantages of antiferromagnetic spintronics, controlling and
detecting N\'eel vectors have been limited in bulk materials. In this work, we
developed a sensitive second harmonic generation (SHG) microscope and detected
long-range N\'eel antiferromagnetic (AFM) order and N\'eel vector switching
down to the monolayer in MnPSe$_3$. Temperature-dependent SHG measurement in
repetitive thermal cooling surprisingly collapses into two curves, which
correspond to the switching of an Ising type N\'eel vector reversed by the
time-reversal operation, instead of a six-state clock ground state expected
from the threefold rotation symmetry in the structure. We imaged the spatial
distribution of the N\'eel vectors across samples and rotated them by an
arbitrary angle irrespective of the lattice in the sample plane by applying
strain. By studying both a Landau theory and a microscopic model that couples
strain to nearest-neighbor exchange, we conclude that the phase transition of
the XY model in the presence of strain falls into the Ising universality class
instead of the XY one, which could explain the extreme strain tunability.
Finally, we found that the 180{\deg} AFM domain walls are highly mobile down to
the monolayer after thermal cycles, paving the way for future control of the
antiferromagnetic domains by strain or external fields on demand for
ultra-compact 2D AFM terahertz spintronics.
