“…Such an energy dissipation low plays an important role in developing stable numerical methods for dissipation systems due to its importance to the long time numerical simulation, see e.g., [12,13,14,44,36,49,16,46,42] and references therein. In recent years, fractional-type phase-field models have attracted more and more attentions [1,2,3,7,30,38,23]. For instance, consider a fractional type free energy [38] E α (φ) := Ω ε 2 2 |∇ α φ| 2 + F (φ) dx, (1.5) where ∇ α is the fractional gradient ∇ α = ( ∂ α ∂x1 , ..., ∂ α ∂x d ) and { ∂ α ∂x k } k are fractional derivatives.…”