The island dynamics model on parallel quadtree grids

Abstract We introduce an approach for simulating epitaxial growth by use of an island dynamics model on a forest of quadtree grids, and in a parallel environment. To this end, we use a parallel framework introduced in the context of the level-set method. This framework utilizes: discretizations that achieve a second-order accurate level-set method on non-graded adaptive Cartesian grids for solving the associated free boundary value problem for surface diffusion; and an established library for the partitioning of the grid. We consider the cases with: irreversible aggregation, which amounts to applying Dirichlet boundary conditions at the island boundary; and an asymmetric (Ehrlich–Schwoebel) energy barrier for attachment/detachment of atoms at the island boundary, which entails the use of a Robin boundary condition. We provide the scaling analyses performed on the Stampede supercomputer and numerical examples that illustrate the capability of our methodology to efficiently simulate different aspects of epitaxial growth. The combination of adaptivity and parallelism in our approach enables simulations that are several orders of magnitude faster than those reported in the recent literature and, thus, provides a viable framework for the systematic study of mound formation on crystal surfaces.

[1]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[2]  Hanan Samet,et al.  Applications of spatial data structures - computer graphics, image processing, and GIS , 1990 .

[3]  J. Wollschläger,et al.  Diffraction characterization of rough films formed under stable and unstable growth conditions , 1998 .

[4]  F. Hudda,et al.  Atomic View of Surface Self‐Diffusion: Tungsten on Tungsten , 1966 .

[5]  Frédéric Gibou,et al.  Geometric integration over irregular domains with application to level-set methods , 2007, J. Comput. Phys..

[6]  B. Merriman,et al.  Kinetic model for a step edge in epitaxial growth. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[7]  Frédéric Gibou,et al.  A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids , 2006, J. Comput. Phys..

[8]  Joachim Krug,et al.  Island nucleation in the presence of step-edge barriers: Theory and applications , 1999, cond-mat/9912410.

[9]  J. Krug,et al.  Islands, mounds and atoms : patterns and processes in crystal growth far from equilibrium , 2004 .

[10]  A. Chernov,et al.  THE SPIRAL GROWTH OF CRYSTALS , 1961 .

[11]  J. Langer.  Instabilities and pattern formation in crystal growth , 1980 .

[12]  Jianfeng Lu,et al.  Emergence of step flow from an atomistic scheme of epitaxial growth in 1+1 dimensions. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  R. Schwoebel.  Step motion on crystal surfaces , 1968 .

[14]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[15]  Frédéric Gibou,et al.  Rate equations and capture numbers with implicit islands correlations , 2001 .

[16]  Ernst,et al.  Observation of a growth instability during low temperature molecular beam epitaxy. , 1994, Physical review letters.

[17]  J. Wendelken,et al.  EVOLUTION OF MOUND MORPHOLOGY IN REVERSIBLE HOMOEPITAXY ON CU(100) , 1997 .

[18]  T. Michely,et al.  The homoepitaxial growth of Pt on Pt(111) studied with STM , 1992 .

[19]  J. Villain,et al.  Physics of crystal growth , 1998 .

[20]  T. S. P. S.,et al.  GROWTH , 1924, Nature.

[21]  W. K. Burton,et al.  The growth of crystals and the equilibrium structure of their surfaces , 1951, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[22]  G. Shortley,et al.  The Numerical Solution of Laplace's Equation , 1938 .

[23]  A. D. Yoffe,et al.  Low-dimensional systems: Quantum size effects and electronic properties of semiconductor microcrystallites (zero-dimensional systems) and some quasi-two-dimensional systems , 1993 .

[24]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[25]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[26]  Carsten Burstedde,et al.  p4est: Scalable Algorithms for Parallel Adaptive Mesh Refinement on Forests of Octrees , 2011, SIAM J. Sci. Comput..

[27]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[28]  Richard L. Schwoebel,et al.  Step Motion on Crystal Surfaces. II , 1966 .

[29]  Ellen D. Williams,et al.  Steps on surfaces: experiment and theory , 1999 .

[30]  Evans,et al.  Transition to Multilayer Kinetic Roughening for Metal (100) Homoepitaxy. , 1995, Physical review letters.

[31]  Frédéric Gibou,et al.  A Supra-Convergent Finite Difference Scheme for the Poisson and Heat Equations on Irregular Domains and Non-Graded Adaptive Cartesian Grids , 2007, J. Sci. Comput..

[32]  A Second Order Accurate Finite Difference Scheme for the Heat Equation on Irregular Domains and Adaptive Grids , 2006 .

[33]  Frédéric Gibou,et al.  A numerical scheme for the Stefan problem on adaptive Cartesian grids with supralinear convergence rate , 2009, J. Comput. Phys..

[34]  S. Osher,et al.  Island dynamics and the level set method for epitaxial growth , 1999 .

[35]  Frédéric Gibou,et al.  A Poisson-Boltzmann solver on irregular domains with Neumann or Robin boundary conditions on non-graded adaptive grid , 2011, J. Comput. Phys..

[36]  Per Brinch Hansen.  Numerical Solution of Laplace's Equation , 1992 .

[37]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[38]  Myung-joo Kang,et al.  Fluctuations and scaling in aggregation phenomena , 2000 .

[39]  J. Villain.  Continuum models of crystal growth from atomic beams with and without desorption , 1991 .

[40]  R. Caflisch,et al.  Level-set method for island dynamics in epitaxial growth , 2002 .

[41]  Frédéric Gibou,et al.  A level set approach for diffusion and Stefan-type problems with Robin boundary conditions on quadtree/octree adaptive Cartesian grids , 2013, J. Comput. Phys..

[42]  Cedric Boeckx,et al.  Islands , 2008, Lang. Linguistics Compass.

[43]  Frédéric Gibou,et al.  Capture numbers in rate equations and scaling laws for epitaxial growth , 2003 .

[44]  P. Smereka,et al.  A Remark on Computing Distance Functions , 2000 .

[45]  C. Thompson,et al.  Evolution of thin-film and surface microstructure , 1991 .

[46]  Bales,et al.  Dynamics of irreversible island growth during submonolayer epitaxy. , 1994, Physical review. B, Condensed matter.

[47]  Smilauer,et al.  Coarsening and slope evolution during unstable spitaxial growth. , 1995, Physical review. B, Condensed matter.

[48]  Sander,et al.  Coarsening of Unstable Surface Features during Fe(001) Homoepitaxy. , 1995, Physical review letters.

[49]  F Gibou,et al.  Singularities and spatial fluctuations in submonolayer epitaxy. , 2003, Physical review letters.

[50]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[51]  Hanan Samet,et al.  The Design and Analysis of Spatial Data Structures , 1989 .

[52]  Russel E. Caflisch,et al.  Growth, Structure and Pattern Formation for Thin Films , 2008, J. Sci. Comput..

[53]  Russel E. Caflisch,et al.  Level set simulation of directed self-assembly during epitaxial growth , 2006 .

[54]  Mohammad Mirzadeh,et al.  Parallel level-set methods on adaptive tree-based grids , 2016, J. Comput. Phys..

[55]  Frédéric Gibou,et al.  Efficient symmetric discretization for the Poisson, heat and Stefan-type problems with Robin boundary conditions , 2010, J. Comput. Phys..

[56]  Ronald Fedkiw,et al.  A level set method for thin film epitaxial growth , 2001 .