“…Thus, (46) holds with 𝜆 = 𝜇 and M = 𝛼 (−𝜇h 𝛼 ) . (b) Many authors use the following Mittag-Leffler stability (see Li et al[8] and their followers): There exists 𝜆 > 0, Because the Mittag-Leffler function admits the asymptotic behavior as in Lemma 1, it means that their Mittag-Leffler stability requires only a certain power decaying rate; meanwhile, our strict Mittag-Leffler stability needs the decay rate to be at least O(t −𝛼 ) where 𝛼 is the order of fractional derivative.Utilizing the argument in Ke and Thuy[28, Theorem 4.1], we gain the following strict Mittag-Leffler stability. Assume that the nonlinearity 𝑓 satisfies the Lipschitz condition||𝑓 (t, v 1 , w 1 ) − 𝑓 (t, v 2 , w 2 )|| ≤ 𝛽||v 1 − v 2 || + 𝜅||w 1 − w 2 ||, for all t ≥ 0, v i , w i ∈ L 2 (Ω), i ∈ {1, 2}, where 𝛽 ≥ 0, 𝜅 > 0 such that 𝛽 + 𝜅 < 𝜆 𝛾 1 ,where 𝜆 1 is the first eigenvalue of the Dirichlet Laplacian.…”