“…Great attention has been devoted in the last decades to numerical approximations of evolutionary stochastic partial differential equations (SPDEs) (see, e.g. [2,4,9,11,12,21,23,25] and references therein). In the present work, we concentrate on a class of semi-linear SPDEs of second order with damping, described by du t = αLu t dt + Lu dt + F (u) dt + dW (t), in D × (0, T ], u(·, 0) = u 0 , u t (·, 0) = v 0 , in D, u = 0, on ∂D × (0, T ], The initial data u 0 , v 0 are assumed to be F 0 -measurable random variables.…”