“…It is known that the elliptic partial differential operator A = d i,j=1 ∂ ∂x i a ij (x) ∂ ∂x j , under the Dirichlet boundary condition, generates a bounded analytic semigroup on the space L q := L q (Ω) for all 1 < q < ∞; see [33,Theorem 2.4]. When u 0 = 0, the solution of (1.1) possesses maximal L p -regularity on L q , namely, ∂ t u L p (R + ;L q ) + Au L p (R + ;L q ) ≤ C p,q f L p (R + ;L q ) ∀ 1 < p, q < ∞, (1.4) which is an important tool in studying well-posedness and regularity of solutions to nonlinear parabolic PDEs; see [3,7,31]. In the numerical solution of parabolic equations, it is also desirable to have a discrete analogue of this estimate, which has a number of applications in analysis of stability and convergence of numerical methods for nonlinear parabolic problems, including semilinear parabolic equations with strong nonlinearities [11], quasi-linear parabolic equations with nonsmooth coefficients [28], optimal control and inverse problems [23,25], and so on.…”