Text this: First- and second-order accurate, unconditionally energy gradient stable, uniquely solvable, and mass-preserving linear numerical schemes for Cahn–Hilliard equation with source term.