In this paper, we give a brief overview of recent analytical developments in the study of the deposition of a crystalline film onto a substrate, with the atoms of the film occupying the substrate’s natural lattice positions. This process is called epitaxial growth. Here we are interested in heteroepitaxy, that is, epitaxy when the film and the substrate have different crystalline structures. At the onset of the deposition, the film’s atoms tend to align themselves with those of the substrate because the energy gain associated with the chemical bonding effect is greater than the film’s strain due to the mismatch between the lattice parameters. As the film continues to grow, the stored strain energy per unit area of the interface increases with the film thickness, rendering the film’s flat layer morphologically unstable or metastable after the thickness reaches a critical value. As a result, the film’s free surface becomes corrugated, and the material agglomerates into clusters or isolated islands on the substrate. The formation of islands in systems such as In-GaAs/GaAs or SiGe/Si has essential high-end technology applications, such as modern semiconductor electronic and optoelectronic devices (quantum dots laser). The Stranski–Krastanow (SK) growth mode occurs when the islands are separated by a thin wetting layer, while the Volmer–Weber (VW) growth mode refers to the case when the substrate is exposed between islands.
In what follows, we adopt the variational model considered by Spencer in [41 B. J. Spencer, Asymptotic derivation of the glued-wetting-layer model and contact-angle condition for Stranski–Krastanow islands. Phys. Rev. B, 59, 2011 (1999) ] (see also [36 R. V. Kukta and L. B. Freund, Minimum energy configuration of epitaxial material clusters on a lattice-mismatched substrate. J. Mech. Phys. Solids45, 1835–1860 (1997) , 42 B. J. Spencer and J. Tersoff, Equilibrium shapes and properties of epitaxially strained islands. Phys. Rev. Lett.79, 4858 (1997) ], and the references contained therein). To be precise, the free energy functional associated with the physical system is given by
Here is the function whose graph describes the profile of the film, assumed to be -periodic, with , for some , is the region occupied by the film, i.e., writing ,
is displacement of the material, is the symmetric part of . Also, the elastic energy density is a positive definite quadratic form defined on the space of symmetric matrices
with a positive definite fourth-order tensor, so that for all , is an anisotropic surface energy density evaluated at the unit normal to , and denotes the two-dimensional Hausdorff measure. We suppose that is positively one-homogeneous and of class away from the origin, so that, in particular,
for some constant .
The substrate and the film admit different natural states corresponding to the mismatch between their respective crystalline structures. To be precise, a natural state for the substrate is given by , while a natural state for the film is given by for some nonzero matrix . Our models will reflect this mismatch, either by setting the elastic bulk energy as , where
or by imposing the Dirichlet boundary condition .
In the two-dimensional static case, existence of equilibrium solutions and their qualitative properties, including regularity, were studied in [3 M. Bonacini, Epitaxially strained elastic films: The case of anisotropic surface energies. ESAIM Control Optim. Calc. Var.19, 167–189 (2013) , 4 M. Bonacini, Stability of equilibrium configurations for elastic films in two and three dimensions. Adv. Calc. Var.8, 117–153 (2015) , 5 E. Bonnetier and A. Chambolle, Computing the equilibrium configuration of epitaxially strained crystalline films. SIAM J. Appl. Math.62, 1093–1121 (2002) , 15 E. Davoli and P. Piovano, Analytical validation of the Young–Dupré law for epitaxially-strained thin films. Math. Models Methods Appl. Sci.29, 2183–2223 (2019) , 16 E. Davoli and P. Piovano, Derivation of a heteroepitaxial thin-film model. Interfaces Free Bound.22, 1–26 (2020) , 17 B. De Maria and N. Fusco, Regularity properties of equilibrium configurations of epitaxially strained elastic films. In Topics in modern regularity theory, CRM Series 13, Ed. Norm., Pisa, 169–204 (2012) , 20 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Equilibrium configurations of epitaxially strained crystalline films: existence and regularity results. Arch. Ration. Mech. Anal.186, 477–537 (2007) , 24 I. Fonseca, G. Leoni and M. Morini, Equilibria and dislocations in epitaxial growth. Nonlinear Anal.154, 88–121 (2017) , 26 N. Fusco, Equilibrium configurations of epitaxially strained thin films. Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl.21, 341–348 (2010) , 29 N. Fusco and M. Morini, Equilibrium configurations of epitaxially strained elastic films: second order minimality conditions and qualitative properties of solutions. Arch. Ration. Mech. Anal.203, 247–327 (2012) , 33 S. Y. Kholmatov and P. Piovano, A unified model for stress-driven rearrangement instabilities. Arch. Ration. Mech. Anal.238, 415–488 (2020) ]. The variational techniques and analytical arguments developed in these papers have been used to treat other materials phenomena, such as voids and cavities in elastic solids [9 G. M. Capriani, V. Julin and G. Pisante, A quantitative second order minimality criterion for cavities in elastic bodies. SIAM J. Math. Anal.45, 1952–1991 (2013) , 19 I. Fonseca, N. Fusco, G. Leoni and V. Millot, Material voids in elastic solids with anisotropic surface energies. J. Math. Pures Appl. (9)96, 591–639 (2011) ].
The scaling regimes of the minimal energy in epitaxial growth were identified in [2 P. Bella, M. Goldman and B. Zwicknagl, Study of island formation in epitaxially strained films on unbounded domains. Arch. Ration. Mech. Anal.218, 163–217 (2015) , 30 M. Goldman and B. Zwicknagl, Scaling law and reduced models for epitaxially strained crystalline films. SIAM J. Math. Anal.46, 1–24 (2014) ] in terms of the parameters of the problem. The shape of the islands under the constraint of faceted profiles was addressed in [25 I. Fonseca, A. Pratelli and B. Zwicknagl, Shapes of epitaxially grown quantum dots. Arch. Ration. Mech. Anal.214, 359–401 (2014) ]. A variational model that takes into account the formation of misfit dislocations was introduced in [23 I. Fonseca, N. Fusco, G. Leoni and M. Morini, A model for dislocations in epitaxially strained elastic films. J. Math. Pures Appl. (9)111, 126–160 (2018) ].
The effect of atoms freely diffusing on the surface (called adatoms) was studied in [10 M. Caroccia, R. Cristoferi and L. Dietrich, Equilibria configurations for epitaxial crystal growth with adatoms. Arch. Ration. Mech. Anal.230, 785–838 (2018) ], where the model involves only surface energies.
A discrete-to-continuum analysis for free-boundary problems related to crystalline films deposited on substrates was undertaken in [35 L. C. Kreutz and P. Piovano, Microscopic validation of a variational model of epitaxially strained crystalline films. SIAM J. Math. Anal.53, 453–490 (2021) , 38 P. Piovano and I. Velčić, Microscopical justification of solid-state wetting and dewetting. arXiv:2010.08787 (2020) ].
The three-dimensional static case was studied in [6 A. Braides, A. Chambolle and M. Solci, A relaxation result for energies defined on pairs set-function and applications. ESAIM Control Optim. Calc. Var.13, 717–734 (2007) , 12 A. Chambolle and M. Solci, Interaction of a bulk and a surface energy with a geometrical constraint. SIAM J. Math. Anal.39, 77–102 (2007) ] in the case in which the symmetrized gradient is replaced by the gradient (see also [4 M. Bonacini, Stability of equilibrium configurations for elastic films in two and three dimensions. Adv. Calc. Var.8, 117–153 (2015) ]). More recently, new developments in the theory of , i.e., generalized special functions of bounded deformation (see [13 V. Crismale and M. Friedrich, Equilibrium configurations for epitaxially strained films and material voids in three-dimensional linear elasticity. Arch. Ration. Mech. Anal.237, 1041–1098 (2020) , 14 G. Dal Maso, Generalised functions of bounded deformation. J. Eur. Math. Soc. (JEMS)15, 1943–1997 (2013) ], and the references therein) have led to considerable progress on the relaxation of the functional (1) in the three dimensional case (see [13 V. Crismale and M. Friedrich, Equilibrium configurations for epitaxially strained films and material voids in three-dimensional linear elasticity. Arch. Ration. Mech. Anal.237, 1041–1098 (2020) ]). The regularity of equilibrium solutions remains an open problem. A local minimality sufficiency criterion, based on the strict positivity of the second variation, was established in [4 M. Bonacini, Stability of equilibrium configurations for elastic films in two and three dimensions. Adv. Calc. Var.8, 117–153 (2015) ], based on the work [29 N. Fusco and M. Morini, Equilibrium configurations of epitaxially strained elastic films: second order minimality conditions and qualitative properties of solutions. Arch. Ration. Mech. Anal.203, 247–327 (2012) ].
To study the morphological evolution of anisotropic epitaxially strained films, we assume that the surface evolves by surface diffusion under the influence of a chemical potential . To be precise, according to the Einstein–Nernst relation, the evolution is governed by the volume preserving equation
where , denotes the normal velocity of the evolving interface , stands for the tangential laplacian, and the chemical potential is given by the first variation of the underlying free-energy functional. In our context, this becomes (assuming )
where stands for the tangential divergence along , and is the elastic equilibrium in , i.e., the minimizer of the elastic energy under the prescribed periodicity and boundary conditions (see (7) below).
If the surface energy density is highly anisotropic, there may be directions for which
fails, see for instance [18 A. Di Carlo, M. E. Gurtin and P. Podio-Guidugli, A regularized equation for anisotropic motion-by-curvature. SIAM J. Appl. Math.52, 1111–1119 (1992) , 40 M. Siegel, M. J. Miksis and P. W. Voorhees, Evolution of material voids for highly anisotropic surface energy. J. Mech. Phys. Solids52, 1319–1353 (2004) ]. In this case, the evolution equation (4) is backward parabolic, and to overcome the ill-posedness of the problem we consider the following singular perturbation of the surface energy
where , stands for the sum of the principal curvatures of , and is a small positive constant (see [18 A. Di Carlo, M. E. Gurtin and P. Podio-Guidugli, A regularized equation for anisotropic motion-by-curvature. SIAM J. Appl. Math.52, 1111–1119 (1992) , 31 M. E. Gurtin and M. E. Jabbour, Interface evolution in three dimensions with curvature-dependent energy and surface diffusion: interface-controlled evolution, phase transitions, epitaxial growth of elastic films. Arch. Ration. Mech. Anal.163, 171–208 (2002) , 32 C. Herring, Some theorems on the free energies of crystal surfaces. Phys. Rev.82, 87 (1951) ]). The restriction in is motivated by the fact that the profile of the film will belong to , where , so that is continuously embedded into . This regularity is strongly used to prove existence of solutions. In contrast, in we can assume since is embedded in .
The regularized free-energy functional becomes
and (3) is replaced by
Coupling this evolution equation on the profile of the film with the elastic equilibrium elliptic system holding in the film, and parametrizing using , we obtain the following Cauchy system of equations with initial and natural boundary conditions:
where and is a -periodic function.
One can find in the literature sixth-order evolution equations of this type (see, e.g., [31 M. E. Gurtin and M. E. Jabbour, Interface evolution in three dimensions with curvature-dependent energy and surface diffusion: interface-controlled evolution, phase transitions, epitaxial growth of elastic films. Arch. Ration. Mech. Anal.163, 171–208 (2002) ] for the case without elasticity, see [40 M. Siegel, M. J. Miksis and P. W. Voorhees, Evolution of material voids for highly anisotropic surface energy. J. Mech. Phys. Solids52, 1319–1353 (2004) ] for the evolution of voids in elastically stressed materials, and [7 M. Burger, F. Haußer, C. Stöcker and A. Voigt, A level set approach to anisotropic flows with curvature regularization. J. Comput. Phys.225, 183–205 (2007) , 39 A. Rätz, A. Ribalta and A. Voigt, Surface evolution of elastically stressed films under deposition by a diffuse interface model. J. Comput. Phys.214, 187–208 (2006) ]).
We use the gradient flow structure of (7) with respect to a suitable -metric (see, e.g., [8 J. W. Cahn and J. E. Taylor, Overview no. 113, surface motion by surface diffusion. Acta Metall. Mater.42, 1045–1063 (1994) ]) to solve the equation via a minimizing movement scheme (see [1 L. Ambrosio, Minimizing movements. Rend. Accad. Naz. Sci. XL Mem. Mat. Appl. (5)19, 191–246 (1995) ]), i.e., we discretize the problem in time and solve suitable minimum incremental problems.
If instead of we used the gradient flow with respect to an -metric, we would obtain a fourth order evolution equation describing motion by evaporation-condensation (see [8 J. W. Cahn and J. E. Taylor, Overview no. 113, surface motion by surface diffusion. Acta Metall. Mater.42, 1045–1063 (1994) , 31 M. E. Gurtin and M. E. Jabbour, Interface evolution in three dimensions with curvature-dependent energy and surface diffusion: interface-controlled evolution, phase transitions, epitaxial growth of elastic films. Arch. Ration. Mech. Anal.163, 171–208 (2002) , 37 P. Piovano, Evolution of elastic thin films with curvature regularization via minimizing movements. Calc. Var. Partial Differ. Equ.49, 337–367 (2014) ]).
The short time existence of solutions to (7) established in [22 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Motion of three-dimensional elastic films by anisotropic surface diffusion with curvature regularization. Anal. PDE8, 373–423 (2015) ] is the first such result for geometric surface diffusion equations with elasticity in three-dimensions. In the recent paper [28 N. Fusco, V. Julin and M. Morini, The surface diffusion flow with elasticity in three dimensions. Arch. Ration. Mech. Anal.237, 1325–1382 (2020) ] (see also [27 N. Fusco, V. Julin and M. Morini, The surface diffusion flow with elasticity in the plane. Comm. Math. Phys.362, 571–607 (2018) ] for the two-dimensional case), the authors proved short-time existence of a smooth solution without the additional curvature regularization. They also showed asymptotic stability of strictly stable stationary sets.
The results summarized here can be found in the papers [20 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Equilibrium configurations of epitaxially strained crystalline films: existence and regularity results. Arch. Ration. Mech. Anal.186, 477–537 (2007) , 21 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Motion of elastic thin films by anisotropic surface diffusion with curvature regularization. Arch. Ration. Mech. Anal.205, 425–466 (2012) , 22 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Motion of three-dimensional elastic films by anisotropic surface diffusion with curvature regularization. Anal. PDE8, 373–423 (2015) ].
1 2D quasistatic equilibrium of epitaxially strained elastic films
In the following sections we assume self-similarity with respect to a planar axis and reduce the context to a two-dimensional framework. To be precise, we suppose that the material fills the infinite strip
where is a Lipschitz function representing the free profile of the film, which occupies the open set
The line corresponds to the film/substrate interface.
We assume that the mismatch strain corresponding to different natural states of the material in the substrate and in the film, respectively, is represented by
with . We will suppose that the film and the substrate share material properties, with homogeneous elasticity positive definite fourth-order tensor . Hence, bearing in mind the mismatch, the elastic energy per unit area is given by , where
for all symmetric matrices .
In turn, the interfacial energy density has a step discontinuity at , i.e.,
where the property
will favor the SK growth mode over the VW mode. For the case , and for different crystalline materials stress tensors for the substrate and for the film, we refer to [15 E. Davoli and P. Piovano, Analytical validation of the Young–Dupré law for epitaxially-strained thin films. Math. Models Methods Appl. Sci.29, 2183–2223 (2019) , 16 E. Davoli and P. Piovano, Derivation of a heteroepitaxial thin-film model. Interfaces Free Bound.22, 1–26 (2020) ].
The total energy of the system is given by
where represents the free surface of the film, that is,
Since the functional is not lower semicontinuous, and thus, in general, does not admit minimizers, we are led to study its relaxation. Let
and
where stands for the pointwise variation of the function . Note that coincides with the pointwise variation of the function , and so
For define
The following equalities hold:
We refer to [20 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Equilibrium configurations of epitaxially strained crystalline films: existence and regularity results. Arch. Ration. Mech. Anal.186, 477–537 (2007) ] for a proof.
Next we study regularity properties of minimizers of in . As customary in constrained variational problems, in order to have more flexibility in the choice of test functions, we prove that the volume constraint can be replaced by a volume penalization.
Let be a minimizer of the functional defined in (17) with . Then there exists such that for every integer , is a minimizer of the penalized functional
over all .
Proof. An argument similar to that of the proof of Theorem 1 guarantees that for every there exists a minimizer of . If for all sufficiently large, then
and so is a minimizer of .
Assume now that there is a subsequence, not relabeled, such that for all . If
for countably many , define
where has been chosen so that . Note that . Indeed, for every partition , we have that
for all . Hence,
which is a contradiction. Therefore, for all sufficiently large
Since
it follows from (18) and (20) that as and that . In turn, by (16), for some constant independent of .
Let be so large that for all . Then
and the function , , satisfies
Consider a partition . Then
where we used the fact that . Hence,
and so, by (20),
We deduce that
For define
By a change of variables and (10), we have
where is the matrix whose entries are
Observe that
Since is a positive definite quadratic form over the symmetric matrices (see (11)), we have that
for all symmetric matrices and . Hence by (1), (10) and (LABEL:601)
where depends only on the ellipticity constants of and . By (20), (21), and (24), we have that
Thus, if
we get a contradiction, and this completes the proof. ∎
To prove the regularity of the free boundary we use the following internal sphere condition.
Let be a minimizer of the functional defined in (17). Then there exists with the property that for every there exists an open ball , with , such that
This result was first proved in a slightly different context by Chambolle and Larsen [11 A. Chambolle and C. J. Larsen, C∞ regularity of the free boundary for a two-dimensional optimal compliance problem. Calc. Var. Partial Differential Equations18, 77–94 (2003) ] (see also [9 G. M. Capriani, V. Julin and G. Pisante, A quantitative second order minimality criterion for cavities in elastic bodies. SIAM J. Math. Anal.45, 1952–1991 (2013) , 20 I. Fonseca, N. Fusco, G. Leoni and M. Morini, Equilibrium configurations of epitaxially strained crystalline films: existence and regularity results. Arch. Ration. Mech. Anal.186, 477–537 (2007) ]). The argument is entirely two-dimensional and its extension to three dimensions is open.
Remark 4. Note that if is the outward unit normal to at , then . Thus, the set
is nonempty.
In the next theorem we prove that admits a left and right derivative at all but countably many points.
Let be a minimizer of the functional defined in (17