“A Positivity-Preserving Second-Order BDF Scheme for the Cahn-Hilliard Equation with Variable Interfacial Parameters”, Commun. Comput. Phys., vol. 28, no. 3, pp. 967–998, Jul. 2020, doi: 10.4208/cicp.OA-2019-0037.