In this paper we consider the numerical analysis of a parabolic variational inequality arising from a deep quench limit of a model for phase separation in a binary mixture due to Cahn and Hilliard. Stability, convergence and error bounds for a finite element approximation are proven. Numerical simulations in one and two space dimensions are presented.