“A Positivity-Preserving Second-Order BDF Scheme for the Cahn-Hilliard Equation with Variable Interfacial Parameters” (2020) Communications in Computational Physics, 28(3), pp. 967–998. doi:10.4208/cicp.OA-2019-0037.