We present an error analysis for a fully discrete finite difference scheme for the three-component Macromolecular Microsphere Composite (MMC) hydrogels system, a ternary Cahn–Hilliard system with a Flory–Huggins–deGennes free energy potential. The numerical scheme was recently proposed, and the positivity-preserving property and unconditional energy stability were theoretically established. In this paper, we rigorously prove first order convergence in time and second order convergence in space for the numerical scheme, in the LΔt∞(0,T;Hh−1)∩LΔt2(0,T;Hh1) norm. Many highly non-standard estimates have to be involved, due to the nonlinear and singular nature of the surface diffusion coefficients. The combination of (i) a higher order asymptotic expansion of the numerical solution (up to second order temporal accuracy); (ii) a rough error estimate (to establish the LΔt∞ bound for the phase variables); (iii) and a refined error estimate have to be carried out to conclude such a convergence result. To our knowledge, it will be the first work to provide an optimal rate convergence estimate for a ternary phase field system with singular energy coefficients.