“…Based on the idea of the Lagrange multiplier method, the original problem transforms to an equivalent system with quadratic energy form by introducing some arbitrary variables, and the IEQ approach aims to develop the linear, implicit and unconditionally energy stable scheme. Thanks to the advantages, the IEQ approach has been successfully applied to solve gradient flows, such as phase field elastic bending energy models [29,48], Allen-Cahn and Cahn-Hilliard equations [6,35], Peng-Robinson equation of state [23,36,51]. Shen et al [37] implemented the idea of IEQ approach using a scalar auxiliary variable and called it the SAV method to circumvent the disadvantages of the IEQ approach.…”