We present a C \(^0\) interior penalty finite element method for the sixth-order modified phase field crystal equation. We demonstrate that the numerical scheme is uniquely solvable, unconditionally energy stable, and convergent. Additionally, the error analysis presented develops a detailed methodology for analyzing time dependent problems utilizing the C \(^0\) interior penalty method. We furthermore support the theory with several numerical experiments.
