Hostname: page-component-76fb5796d-5g6vh Total loading time: 0 Render date: 2024-04-25T08:45:02.101Z Has data issue: false hasContentIssue false
  • English
  • Français

Simulations of Nonlinear Strongly Anisotropic, Misfitting Crystals and Thin Films

Published online by Cambridge University Press:  01 February 2011

Solmaz Torabi ,
Steven Wise ,
Shuwang Li ,
Axel Voigt ,
John Lowengrub  and
Peng Zhou
Solmaz Torabi
Affiliation:
storabi@uci.edu, University of California, Irvine, Material Science and Engineering, 8826 Palo Verde Rd, Irvine, CA, 92617, United States, 9498365388
Steven Wise
Affiliation:
swise@math.utk.edu, University of Tennessee, Knoxville, Mathematics, Knoxville, TN, 37996, United States
Shuwang Li
Affiliation:
lis@math.uci.edu, University of California, Irvine, Mathematics, Irvine, CA, 92697, United States
Axel Voigt
Affiliation:
Axel.Voigt@mailbox.tu-dresden.de, Technische Universität Dresden, Institute für Wissenschaftliches Rechnen, Dresden, N/A, Germany
John Lowengrub
Affiliation:
lowengrb@math.uci.edu, University of California, Irvine, Mathematics, Irvine, CA, 92697, United States
Peng Zhou
Affiliation:
zhou@mrs.org, University of California, Irvine, Mathematics, Irvine, CA, 92697, United States
Get access

Abstract

We present a new approach for modeling strongly anisotropic crystal and epitaxial growth using regularized, anisotropic Cahn-Hilliard-type equations as a model for the growth and coarsening of thin films. When the surface anisotropy is sufficiently strong, sharp corners form and unregularized anisotropic Cahn-Hilliard equations become ill-posed. Our models contain a high order Willmore regularization to remove the ill posedness at the corners. A key feature of our approach is the development of a new formulation in which the interface thickness is independent of crystallographic orientation. In our previous work, we have provided matched asymptotic analysis to show the convergence of our diffuse interface model to the analytical sharp interface model. In previous models there was no such convergence to sharp interface model when the Willmore energy was considered. We present 2D numerical results using an adaptive, nonlinear multigrid finite-difference method. In particular, we find excellent agreement between the computed shapes using the Cahn-Hilliard approach, with a finite but small Willmore regularization, with dynamical numerical simulations of a sharp interface model. The equilibrium shapes from our diffuse model are compared with an analytical sharp-interface theory recently developed by Spencer [1] at the corners, and there is excellent match. Away from the corners there is an excellent agreement between the diffuse model and the classical Wulff shape. Finally, in order to model the misfit and displacement strains, we add the elastic energy and corresponding forces to our diffuse model. We analyze numerically the effect of elastic stress on the corner regularization in terms of two parameters: one parameter that describes the relative strength of the elastic energy to surface energy and the second that characteristics the strength of the surface energy anisotropy. Adding elastic energy modifies the equilibrium shape and in particular affects the shape of the corners. We can predict different Qdot shapes, such as pyramids and domes, based on the strength of the elastic interactions.

Keywords

thin filmsimulationnanostructure
Type
Research Article
Information
MRS Online Proceedings Library (OPL) , Volume 1087: Symposium V – Crystal-Shape Control and Shape-Dependent Properties – Methods, Mechanism, Theory, and Simulation , 2008 , 1087-V02-01
Copyright
Copyright © Materials Research Society 2008

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Spencer, B.J., Asymptotic solution for the equilibrium crystal shape with small corner energy regularization, Physical Review, 69, 2004 2557–2567.Google Scholar
[2]. Barth, J.V., Costantini, G., Kern, K., Engineering atomic and molecular nanostructures at surfaces, Nature, 437, 2005, p. 671–679.Google Scholar
[3] Costantini, G., Rastelli, A., Manzano, C., Songmuang, R., Schmidt, O.G., Kanel, H. von., Universal island shapes of self-organized quantum dots, Appl. Phys. Lett., 85, 2004, 5673–5675.Google Scholar
[4] Gray, J.L., Atha, S., Hull, R. and Florio, J.A., Hierarchical self-assembly of epitaxial semiconductor nanostructures, Nano Letters, 4, 2004, 2447–2450.Google Scholar
[5] , Ratz, Ribalta, A and Voigt, A, Surface evolution of elastically stressed films under deposition by a diffusive interface model, J. of Comp. Phys., 214, 2006, p 187–208.Google Scholar
[6] Wise, S.M., Kim, J.S., Lowengrub, J.S., Solving the regularized strongly anisotropic Cahn-Hilliard Equation by an adaptive multigrid method, J. of Comp. Phys., 2007 (in Press).Google Scholar
[7] Willmore, T.J., Riemannian Geometry, 1992, Oxford Science Publications, Oxford University Press.Google Scholar
[8] Burton, W., Cabrera, N., Frank, F., The growth of crystals and the equilibrium structure of their surfaces, Trans. R. Soc. Lond., A 243, 1951, p 299–358.Google Scholar
[9] Herring, C., Some theorems on the free energies of crystal surfaces, Phys. Rev., 82, 1951, p 8793 Google Scholar
[10] Gurtin, M., Jabbour, M., Interface evolution in three dimensions with curvature dependent energy and surface diffusion: Interface-controlled evolution, phase transitions, epitaxial growth of elastic films, Arch. Rational Mech. Anal., 163, 2002, p 171208.Google Scholar
[11] DiCarlo, , Gurtin, , Podio-Duidugli, , A regularized equation for anisotropic motion-by-curvature, SIAM J. Appl. Math., 52, 1992, p 1111–1118.Google Scholar
[12] Burger, M., Haußer, F., Stocker, C., Voigt, A., A level set approach to anisotropic flows with curvature regularization, J. Comput. Phys., 2007 (in Press).Google Scholar
[13] Kobayashi, R., Modeling and numerical simulations of dendritic crystal growth, Physica D, 63, 1993, p 410–423.Google Scholar
[14] Eggleston, J.J., McFadden, G.B. and Vorhees, P.W., A phase-field model for highly anisotropic interfacial energy, Physica D, 150, 2001, p 91–103.Google Scholar
[15] Haußer, F., Voigt, A., A discrete scheme for regularized anisotropic surface diffusion: A 6th order geometric evolution equation, Interfaces and Free Boundaries, 7, 2005, p 353–369.Google Scholar
[16] Wise, S., Lowengrub, J., Kim, J., Thornton, K., Voorhees, P., Johnson, W., Quantum dot formation on a strain-patterned epitaxial thin film, Appl. Phys. Lett. 87 (2005) 133102.Google Scholar
[17] Wheeler, A., Phase-field theory of edges in an anisotropic crystal, Proc. R. Soc. A, 462, 2006, p 3363–3384.Google Scholar
[18] Du, Q., Liu, C., Ryham, R. and Wang, X., A phase field formulation of the Willmore problem, Nonlinearity, 18, 2005, p 1249–1267.Google Scholar
[19] Torabi, S., Wise, S., Rätz, A., Lowengrub, J., Voigt, A., A New Method for Simulating Strongly Anisotropic Cahn-Hilliard Equations, in progress.Google Scholar
[20] Berger, M., Oliger, J., Adaptive mesh refinement for hyperbolic partial differential equations, J. Comput. Phys., 53, 1984, p 484–512.Google Scholar
[21] Berger, M., Colella, P., Local adaptive mesh refinement for shock hydrodynamics, J. Comput. Phys., 82, 1989, p 64–84.Google Scholar
[22] Colella, P., Graves, D. T., Ligocki, T. J., Martin, D. F., Modiano, D., Serafini, D. B., Straalen, B. V., CHOMBO software package for amr applications: Design document, Tech. rep., Lawrence Berkeley National Laboratory, Applied Numerical Algorithms Group; NERSC Division; Berkeley, CA, USA, 2003.Google Scholar
[23] Martin, D., Colella, P., A cell-centered adaptive projection method for the incompressible Euler equations, J. Comput. Phys., 163, 2000, p 271–312.Google Scholar
[24] Mitran, S., BEARCLAW a code for multiphysics applications with embedded boundaries: User's manual, Tech. rep., Department of Mathematics, University of North Carolina, http://www.amath.unc.edu/Faculty/mitran/bearclaw.html, 2006 Google Scholar
[25] Kim, J., Kang, K., Lowengrub, J., Conservative multigrid methods for Cahn-Hilliard fluids, J. Comput. Phys., 193, 2004, p 511–543.Google Scholar
[26] Brandt, A., Multi-level adaptive solutions to boundary-value problems, Math. Comput., 31, 1977, p 333–390.Google Scholar
[27] Trottenberg, U., Oosterlee, C., Schüller, A., Multigrid, Academic Press, London, UK, 2001.Google Scholar
[28] Garcke, H., Rumpf, M., Weikard, U., The Cahn-Hilliard equation with elasticity: Finite element approximation and qualitative studies, Interfaces and Free Boundarie, 3, 2001, p 101–118.Google Scholar
[29] Kay, D., Welford, R., A multigrid finite element solver for the Cahn-Hilliard equation, J. Comput. Phys., 212, 2006, p 288–304.Google Scholar
[30] Siegel, M, Miksis, MJ, Voorhees, PW, Evolution of material voids for highly anisotropic surface energy, Journal of the Mechanics and Physics of Solids, 52, 2004, p 1319–1353.Google Scholar