Before revision:  The numerical scheme is constructed based on a recently developed energetic variational formulation, in which the reaction part is reformulated in terms of reaction trajectories , and both the reaction and the diffusion parts dissipate the same free energy . This subtle fact opens a path of an  operator splitting scheme for these systems.
After revision:  The numerical scheme is constructed based on a recently developed energetic variational formulation, in which the reaction part is reformulated in terms of reaction trajectories , and both the reaction and the diffusion parts dissipate the same free energy opens a path of an operator splitting scheme for these systems.