Surfaces with parallel mean curvature in Sasakian space forms

We study the global geometry of surfaces in Sasakian space forms whose mean curvature vector is parallel in the normal bundle (these include the Riemannian Heisenberg space of dimension $2n+1$). We prove a codimension reduction theorem. We introduce two holomorphic quadratic differentials on anti-invariant such surfaces and use them to obtain classification theorems.