“…Therefore, the MBP is an indispensable tool to study physical features of semilinear parabolic equations, including the aspects of mathematical analysis and numerical simulation. Up to now, great efforts have been made in developing MBP-preserving numerical methods for equations like (1.1), such as the stablized linear semi-implicit method [24,25], the nonlinear second-order method [9,10], the exponential time differencing method [7,8], the integrating factor method [13,16,17], the exponential cut-off method [15,29], and the exponential-SAV method [11,12]. As for the spatial discretizations, a partial list includes the works for finite element method [2,5,15,27,28,30], finite difference method [3,4,26], and finite volume method [21,22].…”