Closely related to the HHO method are the virtual element methods (VEM) [1]}, [2]}, [3]}, which are also in principle arbitrary-order polytopal methods. A VEM Stokes complex is developed in [4]}, and used in the recent work [5]} to design a VEM scheme for MHD equations. However, only the lowest order VEM is considered and analysed in [5]}.
