Energetics and coarsening analysis of a simplified non-linear surface growth model

• Published in: Discrete and continuous dynamical systems. Series S Vol. 15; no. 1; pp. 161 - 177
• Main Authors: Khalfi, Hamza, Aarab, Amal, Alaa, Nour Eddine
• Format: Journal Article
• Language: English
• Published: 01.01.2022
• Subjects: coarsening; convergence to equilibrium; D. C. Algorithm; interface width; Epitaxial growth; energy minimization
• ISSN: 1937-1179; 1937-1632; 1937-1179
• DOI: 10.3934/dcdss.2021014

We study a simplified multidimensional version of the phenomenological surface growth continuum model derived in [6]. The considered model is a partial differential equation for the surface height profile \begin{document}$ u $\end{document} which possesses the following free energy functional: \begin{document}$ E(u) = \int_{\Omega} \left[ \frac{1}{2} \ln\left(1+\left|\nabla u \right|^2\right) - \left|\nabla u \right| \arctan\left(\left|\nabla u \right|\right) + \frac{1}{2} \left|\Delta u \right|^2 \right] {\rm d}x, $\end{document} where \begin{document}$ \Omega $\end{document} is the domain of a fixed support on which the growth is carried out. The term \begin{document}$ \left|\Delta u \right|^2 $\end{document} designates the standard surface diffusion in contrast to the second order term which phenomenologically describes the growth instability. The energy above is mainly carried out in regions of higher surface slope \begin{document}$ \left( \left|\nabla u \right| \right) $\end{document}. Hence minimizing such energy corresponds to reducing surface defects during the growth process from a given initial surface configuration. Our analysis is concerned with the energetic and coarsening behaviours of the equilibrium solution. The existence of global energy minimizers and a scaling argument are used to construct a sequence of equilibrium solutions with different wavelength. We apply our minimum energy estimates to derive bounds in terms of the linear system size \begin{document}$ \left| \Omega \right| $\end{document} for the characteristic interface width and average slope. We also derive a stable numerical scheme based on the convex-concave decomposition of the energy functional and study its properties while accommodating these results with 1d and 2d numerical simulations.