“A Positivity-Preserving Second-Order BDF Scheme for the Cahn-Hilliard Equation With Variable Interfacial Parameters”. Communications in Computational Physics, vol. 28, no. 3, July 2020, pp. 967-98, https://doi.org/10.4208/cicp.OA-2019-0037.