“…In the discrete case, it is hence necessary to develop numerical schemes that preserve the energy decay and maximum-principle for model (3). Historically, there has been enormous numerical exploration for tFAC equations or phase field models or more generally, the nonlinear subdiffustion models, where some energy stable schemes were developed (see, e.g., [16,15,14,25,24,19,5,21,29,44,33,23,9]). However, the studies of preservation on discrete energy dissipation, or energy decay and maximum-principle are still limited.…”