arXiv:2609.39996v2 Announce Type: replace Abstract: We develop a pressure-free virtual element method for the surface Stokes problem on polygonal approximations of closed surfaces of arbitrary genus. The method extends the planar VEM Stokes complex to discrete surfaces, yielding a nonconforming complex with an exactly tangential, $H(\operatorname{div})$-conforming velocity space, commuting interpolation, and the correct discrete cohomology. To obtain the weak consistency needed for optimal convergence, we introduce an edge correction that preserves the complex without adding degrees of freedom; the resulting construction reduces to the original complex in the planar setting. The virtual formulation requires no penalties, and a local divergence-preserving reconstruction provides a computable pressure-robust load. We also construct explicit macroelement representatives whose energy and load induce an alternative formulation on the VEM space, algebraically equivalent to a direct macroelement Galerkin method. The original virtual formulation and the macroelement method are stable and attain optimal first-order convergence in the broken $H^1$ norm and second-order convergence in the $L^2$ norm for the velocity. Numerical experiments illustrate these convergence rates, the role of the edge correction, and pressure robustness.