Abstract:Computing equilibrium shapes of crystals (ESC) is a challenging problem in materials science that involves minimizing an orientation-dependent (i.e., anisotropic) surface energy functional subject to a prescribed mass constraint. The highly nonlinear and singular anisotropic terms in the problem make it very challenging from both the analytical and numerical aspects. Especially, when the strength of anisotropy is very strong (i.e., strongly anisotropic cases), the ESC will form some singular, sharp corners even if the surface energy function is smooth. Traditional numerical approaches, such as the $H^{-1}$ gradient flow, are unable to produce true sharp corners due to the necessary addition of a high-order regularization term that penalizes sharp corners and rounds them off. In this paper, we propose a new numerical method based on the Davis-Yin splitting (DYS) optimization algorithm to predict the ESC instead of using gradient flow approaches. We discretize the infinite-dimensional phase-field energy functional in the absence of regularization terms and transform it into a finite-dimensional constraint minimization problem. The resulting optimization problem is solved using the DYS method which automatically guarantees the mass-conservation and bound-preserving properties. We also prove the global convergence of the proposed algorithm. These desired properties are numerically observed. In particular, the proposed method can produce real sharp corners with satisfactory accuracy. Finally, we present numerous numerical results to demonstrate that the ESC can be well simulated under different types of anisotropic surface energies, which also confirms the effectiveness and efficiency of the proposed method.

Unconditionally energy-stable linear convex splitting algorithm for the L2 quasicrystals

Phase-field Modeling and Linearly Energy-Stable Runge–Kutta Algorithm of Colloidal Crystals on Curved Surfaces

New third-order convex splitting methods and analysis for the phase field crystal equation

A Second-Order, Uniquely Solvable, Energy Stable BDF Numerical Scheme for the Phase Field Crystal Model

Efficient numerical methods for computing the stationary states of phase field crystal models

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

A Linearly Second-Order, Unconditionally Energy Stable Scheme and Its Error Estimates for the Modified Phase Field Crystal Equation

Efficient Numerical Schemes with Unconditional Energy Stabilities for the Modified Phase Field Crystal Equation

Efficient Second-Order Unconditionally Stable Numerical Schemes for the Modified Phase Field Crystal Model with Long-Range Interaction

An Efficient Linear Second Order Unconditionally Stable Direct Discretization Method for the Phase-Field Crystal Equation on Surfaces

A New Lagrange Multiplier Method for the Mass-Conserved Allen–Cahn Type Square Phase-Field Crystal Model

Parallel Energy-Stable Phase Field Crystal Simulations Based on Domain Decomposition Methods.

An Unconditionally Energy Stable Method for the Anisotropic Phase-Field Crystal Model in Two Dimension

High-order energy stable schemes of incommensurate phase-field crystal model

Splitting Methods for Semi-Classical Hamiltonian Dynamics of Charge Transfer in Nonlinear Lattices

An operator-splitting optimization approach for phase-field simulation of equilibrium shapes of crystals

Stability and Error Estimate of the Operator Splitting Method for the Phase Field Crystal Equation

Quasicrystals as an intermediate form of matter between crystalline and amorphous solids

A new Allen–Cahn type two-model phase-field crystal model for fcc ordering and its numerical approximation

A New Space-Fractional Modified Phase Field Crystal Equation and Its Numerical Algorithm

A Second Order Unconditionally Stable Scheme for the Modified Phase Field Crystal Model with Elastic Interaction and Stochastic Noise Effect