An energy stable, hexagonal finite difference scheme for the 2D phase field crystal amplitude equations

In this paper we construct an energy stable finite difference scheme for the amplitude expansion equations for the two-dimensional phase field crystal (PFC) model. The equations are formulated in a periodic hexagonal domain with respect to the reciprocal lattice vectors to achieve a provably unconditionally energy stable and solvable scheme. To our knowledge, this is the first such energy stable scheme for the PFC amplitude equations. The convexity of each part in the amplitude equations is analyzed, in both the semi-discrete and fully-discrete cases. Energy stability is based on a careful convexity analysis for the energy (in both the spatially continuous and discrete cases). As a result, unique solvability and unconditional energy stability are available for the resulting scheme. Moreover, we show that the scheme is point-wise stable for any time and space step sizes. An efficient multigrid solver is devised to solve the scheme, and a few numerical experiments are presented, including grain rotation and shrinkage and grain growth studies, as examples of the strength and robustness of the proposed scheme and solver.

[1]  Slobodan Nickovic,et al.  Geostrophic Adjustment on Hexagonal Grids , 2002 .

[2]  Cheng Wang,et al.  Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation , 2012, J. Comput. Phys..

[3]  Nikolas Provatas,et al.  Phase field crystal study of deformation and plasticity in nanocrystalline materials. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  M. Grant,et al.  Phase-field crystal modeling and classical density functional theory of freezing , 2007 .

[5]  Zhizhang Chen.  Unconditionally stable finite-difference time-domain methods and their applications , 2004, Proceedings. ICCEA 2004. 2004 3rd International Conference on Computational Electromagnetics and Its Applications, 2004..

[6]  Hirofumi Tomita,et al.  Shallow water model on a modified icosahedral geodesic grid by using spring dynamics , 2001 .

[7]  Tao Tang,et al.  An Adaptive Time-Stepping Strategy for the Molecular Beam Epitaxy Models , 2011, SIAM J. Sci. Comput..

[8]  Xesús Nogueira,et al.  An unconditionally energy-stable method for the phase field crystal equation , 2012 .

[9]  Badrinarayan P. Athreya,et al.  Renormalization Group Approach to Multiscale Modelling in Materials Science , 2005, cond-mat/0508671.

[10]  Badrinarayan P. Athreya,et al.  Renormalization-group theory for the phase-field crystal equation. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Yan Xu,et al.  Local Discontinuous Galerkin Method and High Order Semi-Implicit Scheme for the Phase Field Crystal Equation , 2016, SIAM J. Sci. Comput..

[12]  E Granato,et al.  Patterning of heteroepitaxial overlayers from nano to micron scales. , 2012, Physical review letters.

[13]  Steven M. Wise,et al.  Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation , 2012, SIAM J. Numer. Anal..

[14]  László Gránásy,et al.  Faceting and branching in 2D crystal growth. , 2011, Physical review letters.

[15]  Steven M. Wise,et al.  Unconditionally stable schemes for equations of thin film epitaxy , 2010 .

[16]  T Ala-Nissila,et al.  Phase-field-crystal models and mechanical equilibrium. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Tomohiro Takaki,et al.  Development of numerical scheme for phase field crystal deformation simulation , 2009 .

[18]  Yuan Ma,et al.  An adaptive time-stepping strategy for solving the phase field crystal model , 2013, J. Comput. Phys..

[19]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[20]  Peter W Voorhees,et al.  Stress-induced morphological instabilities at the nanoscale examined using the phase field crystal approach , 2009 .

[21]  Cheng Wang,et al.  A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation , 2014, Numerische Mathematik.

[22]  Ken Elder,et al.  Morphological instability, evolution, and scaling in strained epitaxial films: An amplitude-equation analysis of the phase-field-crystal model , 2010, 1003.2611.

[23]  Zhi-Feng Huang,et al.  Density-amplitude formulation of the phase-field crystal model for two-phase coexistence in two and three dimensions , 2010 .

[24]  Steven M. Wise,et al.  Solving the regularized, strongly anisotropic Cahn-Hilliard equation by an adaptive nonlinear multigrid method , 2007, J. Comput. Phys..

[25]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

[26]  Joel Berry,et al.  Modeling multiple time scales during glass formation with phase-field crystals. , 2011, Physical review letters.

[27]  N. Provatas,et al.  Phase-field crystals with elastic interactions. , 2006, Physical review letters.

[28]  E Granato,et al.  Phase diagram and commensurate-incommensurate transitions in the phase field crystal model with an external pinning potential. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  Nikolas Provatas,et al.  Phase-Field Methods in Materials Science and Engineering , 2010 .

[30]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[31]  Martin Grant,et al.  Modeling elasticity in crystal growth. , 2001, Physical review letters.

[32]  Zhonghua Qiao,et al.  An Adaptive Time-Stepping Strategy for the Cahn-Hilliard Equation , 2012 .

[33]  Daniel Lee.  Complete solution to seven-point schemes of discrete anisotropic Laplacian on regular hexagons , 2015 .

[34]  Alain Karma,et al.  Phase-field crystal study of grain-boundary premelting , 2008, 0807.5083.

[35]  Zhi-zhong Sun,et al.  Two finite difference schemes for the phase field crystal equation , 2015 .

[36]  T. Ala‐Nissila,et al.  Thermodynamics of bcc metals in phase-field-crystal models. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  Steven M. Wise,et al.  An Energy Stable and Convergent Finite-Difference Scheme for the Modified Phase Field Crystal Equation , 2011, SIAM J. Numer. Anal..

[38]  Joel Berry,et al.  Simulation of an atomistic dynamic field theory for monatomic liquids: freezing and glass formation. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[39]  S. M. Wise,et al.  Unconditionally Stable Finite Difference, Nonlinear Multigrid Simulation of the Cahn-Hilliard-Hele-Shaw System of Equations , 2010, J. Sci. Comput..

[40]  Badrinarayan P. Athreya,et al.  Renormalization group approach to multiscale simulation of polycrystalline materials using the phase field crystal model. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  E Granato,et al.  Thermal fluctuations and phase diagrams of the phase-field crystal model with pinning. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  Robert Spatschek,et al.  Amplitude equations for polycrystalline materials with interaction between composition and stress , 2010, 1002.1580.

[43]  Badrinarayan P. Athreya,et al.  Adaptive mesh computation of polycrystalline pattern formation using a renormalization-group reduction of the phase-field crystal model. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  Zhi-Feng Huang,et al.  Mesoscopic and microscopic modeling of island formation in strained film epitaxy. , 2008, Physical review letters.

[45]  Tapio Ala-Nissila,et al.  Modeling Self-Organization of Thin Strained Metallic Overlayers from Atomic to Micron Scales , 2013 .

[46]  Guohua Zhou,et al.  Fourier Analysis of Multigrid Methods on Hexagonal Grids , 2009, SIAM J. Sci. Comput..

[47]  E Granato,et al.  Dynamical transitions and sliding friction of the phase-field-crystal model with pinning. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  Cheng Wang,et al.  An Energy-Stable and Convergent Finite-Difference Scheme for the Phase Field Crystal Equation , 2009, SIAM J. Numer. Anal..

[49]  M. Grant,et al.  Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[50]  E Granato,et al.  Nonlinear driven response of a phase-field crystal in a periodic pinning potential. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[51]  Joel Berry,et al.  Melting at dislocations and grain boundaries: A phase field crystal study , 2008, 0801.0765.