“…Hence there exists a constant C which depends only on θ and α such that, see Jin et al [21, (2.3)], (z α I + A) −1 ≤ C|z| −α , ∀ z ∈ Σ θ = {z = 0 : | arg z| < θ}. (1.3) Under the assumptions that the solutions of (1.1) are sufficiently smooth, for example u ∈ C 2 [0, T ], there are many time discretisation schemes in the literature, for example, Lubich's fractional multistep methods [7], [30], [55], [65], [61], [62], [64], the L1 scheme and its modification [28], [26], [18], [29], [47], the spectral method [58] [63], [4], [68], [57], nonpolynomial collocation method [14], discontinuous Galerkin method [39], etc.…”