Maximal $L^p$-$L^q$ regularity of finite element methods for stochastic parabolic equations in polygonal domains

We prove that a spatial semidiscrete finite element method preserves the stochastic maximal regularity of the heat equation on bounded convex polygonal and polyhedral domains in dimensions two and three. A main difficulty is to establish suitable Green-function estimates for stochastic convolutions,