Summary: We present an unconditionally energy stable finite difference scheme for the modified phase field crystal equation, a generalized damped wave equation for which the usual phase field crystal equation is a special degenerate case. The method is based on a convex splitting of a discrete pseudoenergy and is semi-implicit. The equation at the implicit time level is nonlinear but represents the gradient of a strictly convex function and is thus uniquely solvable, regardless of time step-size. We present a local-in-time error estimate that ensures the pointwise convergence of the scheme.