Stranski-Krastanov Growth Mode

• Stranski-Krastanov growth mode is a heteroepitaxial process characterized by an initial 2D wetting layer that transitions to 3D island nucleation due to strain and surface energy balance.
• It is governed by a thermodynamic framework where the wetting parameter and lattice mismatch determine the critical thickness for the 2D–3D transformation.
• The model enables quantitative predictions of critical misfit and island size, which are essential for engineering quantum dots and nanoscale surface structures.

The Stranski-Krastanov (SK) growth mode is a foundational mechanism in heteroepitaxial thin film synthesis, especially for systems with significant lattice mismatch. It describes a morphological transition where epitaxial growth begins as layer-by-layer (2D) then switches, above a critical thickness, to the nucleation and enlargement of three-dimensional (3D) islands atop a residual wetting layer. This article provides a rigorous definition, complete thermodynamic framework, quantitative criteria, and key experimental consequences, with explicit technical detail and reference to recent canonical works (Prieto et al., 2017).

1. Thermodynamic Foundation and Wetting Criterion

SK growth is governed by competition between surface/interfacial energies and lattice-mismatch-induced elastic strain. The equilibrium-shape analysis of Bauer specifies the excess surface energy for a coherently strained film A on a substrate S: $$\Delta \sigma = \sigma(A) + \sigma(A/S) - \sigma(S)$$ The wetting parameter $\Phi$ is central to the mode selection: $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$ where $\sigma$ is the film surface energy, $\sigma_s$ the substrate surface energy, $\sigma_i$ the film-substrate interfacial energy, and $\beta$ the adhesion energy per unit area.

The sign of $\Phi$ and the magnitude of the lattice misfit $f = (b-a)/a$ define the growth regime:

• $\Phi > 0$: Volmer-Weber (VW), direct 3D island formation.
• $\Phi$0: Frank–van der Merwe (FM), layer-by-layer (2D) growth.
• $\Phi$1: Stranski–Krastanov (SK), 2D wetting layer succeeded by 3D island nucleation.

The presence of both FM and VW sequences is fully explained by the wettability criterion, which is robust across metals, semiconductors, and advanced alloy systems.

2. Energetics and Critical Thickness

The driving force for 2D–3D transformation is the difference in chemical potential between an $\Phi$2-monolayer pseudomorphic film and bulk 3D material: $\Phi$3 with $\Phi$4, for shear modulus $\Phi$5, Poisson ratio $\Phi$6, monolayer thickness $\Phi$7, and misfit $\Phi$8.

Islands form when the 2D–3D driving force $\Phi$9 turns positive: $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$0

The transition thickness (critical thickness $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$1) satisfies $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$2: $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$3 This explicit dependence on both wetting parameter and misfit solidly quantifies the SK regime boundary.

3. Monolayer–Multilayer Transformations and Island Nucleation

SK island nucleation incorporates both nucleation-like and non-nucleation transformations depending on the nature of the overlayer. For mono–bilayer transitions, the critical island base area is

$$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$4

with nucleation barriers diverging as $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$5: $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$6 Critical nucleus shapes are universally compact rectangles or squares with one additional atom (Prieto et al., 2011). The monolayer is the necessary precursor for multilayer islands, with nucleation preferentially occurring at island corners or edges, not centers. In "stiff" materials and compressive overlayers, transformation proceeds via sequential nucleation; in "soft" materials (Pb, In) and tensile overlayers, direct multilayer formation dominates, producing "magic" thickness islands.

4. Phase Diagram and Material Dependence

Growth mode boundaries are sharply defined in the $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$7 plane. FM mode exists for $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$8 and $$\Phi = \frac{\sigma + \sigma_i - \sigma_s}{2\sigma} = 1 - \frac{\beta}{2\sigma}$$9, SK for $\sigma$0 and $\sigma$1, and VW for all $\sigma$2 if $\sigma$3.

Quantitative predictions are validated experimentally:

• Critical misfits $\sigma$4 observed: $\sigma$51.4% for Si$\sigma$6Ge$\sigma$7/Si, $\sigma$81.8% for InGaAs/GaAs.
• Monolayer–bilayer critical sizes diverge as $\sigma$9; $\sigma_s$0 atoms for $\sigma_s$1 when $\sigma_s$2 (1+1D model).
• Preferred island heights in soft, tensile metals: e.g., Pb/Si $\sigma_s$3 7 ML, In/Si $\sigma_s$4 4 ML, arising from minima in multilayer nucleation barriers.

5. Model Approximations and Implications

SK modeling typically employs anharmonic Morse-type interatomic potentials truncated to first neighbors in 2+1D or 1+1D, with rigid substrate and pseudomorphic wetting-layer assumptions—coherent, dislocation-free regime only. Step energies and step–step repulsion are treated atomistically.

"Stiff" overlayers with large force constants support classical nucleation-driven, layer-by-layer 2D–3D transformation (chiefly compressive cases). "Soft" overlayers favor direct multilayer 2D transformation and the emergence of stable, quantized ("magic-height") islands.

The derived relationships enable practical estimation of critical misfit, onset of 3D clustering, and quantum-dot size distributions as functions of key material parameters $\sigma_s$5 and lattice misfit $\sigma_s$6, facilitating rational design across semiconductor and metallic systems (Prieto et al., 2017).

6. Impact on Quantum Dot and Surface Nanostructure Engineering

The SK mechanism underpins self-assembled quantum dot synthesis, notably in III–V semiconductors (InAs/GaAs, InAs/InP), with defect-free, coherently strained islands nucleated above a thin wetting layer (Zhang et al., 2022, Berdnikov et al., 2023). The wetting layer, while necessary for SK growth, can introduce undesirable continuum electronic states; advanced engineering (e.g., AlAs monolayer insertion) eliminates these states for improved photonic performance (Löbl et al., 2018).

Surfactant-mediated SK transitions are now realized with atomic precision, enabling "on-demand" 3D islanding and QD morphology control (e.g., use of Bi or graphene) (Lewis et al., 2017, Lewis et al., 2019, Rougemaille et al., 2019). Model systems span metals, semiconductors, and complex alloys, with tunable density, size, and strain relaxation—a paradigm crucial for single-photon sources, quantum computation, and device nanostructure control.

7. Summary Table: Core SK Parameters (Selected Systems)

System	Typical $\sigma_s$7 (%)	h_c (nm/ML)	Overlayer Type
Si$\sigma_s$8Ge$\sigma_s$9/Si	$\sigma_i$01.4	$\sigma_i$11–2 ML	Stiff/semiconductor
InGaAs/GaAs	$\sigma_i$21.8	$\sigma_i$31–1.6 ML	Stiff/semiconductor
Pb/Si	Variable	7 ML ("magic height")	Soft/metal
Re/Nb	$\sigma_i$413	$\sigma_i$515 nm	Stiff/metal

For detailed definitions and formulae, see (Prieto et al., 2017).

---

Stranski-Krastanov growth remains a robust, quantitatively predictive framework for analyzing, controlling, and engineering strain-induced morphological transitions, enabling coherent self-assembly across a breadth of material systems. The doctrine of critical thickness, wetting parameter, and nucleation pathways grounded in explicit thermodynamic and atomistic modeling continues to inform experimental strategies and computational models in surface and interface science.