Hermite interpolation polynomials allow a function to be interpolated in terms of its values and derivatives at certain fixed points [1]. In the finite element setting this allows the number of degrees of freedom to be reduced by enforcing greater regularity at element interfaces [2]. Using these polynomials on a rectangular grid combined with the cutFEM framework [3, 4] for the wave equation and using the deal.II software package [5], a high order solution scheme is achieved with similar CFL values to an uncut Hermite scheme.