Energy quadratization Runge-Kutta method for the modified phase field crystal equation

In this paper, we propose high order and unconditionally energy stable methods for a modified phase field crystal equation by applying the strategy of the energy quadratization Runge-Kutta methods. We transform the original model into an equivalent system with auxiliary variables and quadratic free energy. The modified system preserves the laws of mass conservation and energy dissipation with the associated energy functional. We present rigorous proofs of the mass conservation and energy dissipation properties of the proposed numerical methods and present numerical experiments conducted to demonstrate their accuracy and energy stability. Finally, we compare long-term simulations using an indicator function to characterize the pattern formation.