Exponential convergence of the hp virtual element method in presence of corner singularities

In the present work, we analyze the hp version of virtual element methods for the 2D Poisson problem. We prove exponential convergence of the energy error employing sequences of polygonal meshes geometrically refined, thus extending the classical choices for the decomposition in the hp finite element framework to very general decomposition of the domain. A new stabilization for the discrete bilinear form with explicit bounds in h and p is introduced. Numerical experiments validate the theoretical results. We also exhibit a numerical comparison between hp virtual elements and hp finite elements.