Abstract:In this paper we propose and analyze a second order accurate (in time) numerical scheme for the square phase field crystal equation, a gradient flow modeling crystal dynamics at the atomic scale in space but on diffusive scales in time. Its primary difference with the standard phase field crystal model is an introduction of the 4-Laplacian term in the free energy potential, which in turn leads to a much higher degree of nonlinearity. To make the numerical scheme linear while preserving the nonlinear energy stability, we make use of the scalar auxiliary variable (SAV) approach, in which a second order backward differentiation formula is applied in the temporal stencil. Meanwhile, a direct application of the SAV method faces certain difficulties, due to the involvement of the 4-Laplacian term, combined with a derivation of the lower bound of the nonlinear energy functional. In the proposed numerical method, an appropriate decomposition for the physical energy functional is formulated, so that the nonlinear energy part has a well-established global lower bound, and the rest terms lead to constant-coefficient diffusion terms with positive eigenvalues. In turn, the numerical scheme could be very efficiently implemented by constant-coefficient Poisson-like type solvers (via FFT), and energy stability is established by introducing an auxiliary variable, and an optimal rate convergence analysis is provided for the proposed SAV method. A few numerical experiments are also presented, which confirm the efficiency and accuracy of the proposed scheme.

Second-Order SAV Schemes for the Nonlinear Schrödinger Equation and Their Error Analysis

Error Estimates for the Scalar Auxiliary Variable (SAV) Schemes to the Modified Phase Field Crystal Equation

Split-step Spectral Galerkin Method for the Two-Dimensional Nonlinear Space-Fractional Schrödinger Equation

Numerical Methods for the Derivative Nonlinear Schrödinger Equation

Second-order linear adaptive time-stepping schemes for the fractional Allen–Cahn equation

Error Analysis of First- and Second-Order Linear, Unconditionally Energy-Stable Schemes for the Swift-Hohenberg Equation

SAV-based entropy-dissipative schemes for a class of kinetic equations

Two Second-Order and Linear Numerical Schemes for the Multi-Dimensional Nonlinear Time-Fractional Schrödinger Equation

Linear High-Order Energy-Preserving Schemes for the Nonlinear Schrödinger Equation with Wave Operator Using the Scalar Auxiliary Variable Approach

Unconditional Convergence and Optimal Error Estimates of the Euler Semi-Implicit Scheme for a Generalized Nonlinear Schrödinger Equation

Error Analysis with Polynomial Dependence on in SAV Methods for the Cahn-Hilliard Equation

A second-order L 2 - 1 σ difference scheme for the nonlinear time–space fractional Schrödinger equation

Error estimate of a fully decoupled numerical scheme based on the Scalar Auxiliary Variable (SAV) method for the Boussinesq system

Convergence and Error Analysis for the Scalar Auxiliary Variable (SAV) Schemes to Gradient Flows

Error analysis of the SAV-MAC scheme for the Navier-Stokes equations

Linearly Implicit Conservative Schemes for the Nonlocal Schrödinger Equation

Error Analysis of a First-Order IMEX Scheme for the Logarithmic Schrödinger Equation

Asymptotic error analysis for stochastic gradient optimization schemes with first and second order modified equations

Long-time stable SAV-BDF2 numerical schemes for the forced Navier-Stokes equations

A Second Order Accurate Scalar Auxiliary Variable (SAV) Numerical Method for the Square Phase Field Crystal Equation

Stability and error analysis of the SAV schemes for the inductionless MHD equations