Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation

Basic Information

• Paper ID: 2501.00783
• Title: Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation
• Authors: Zhenghua Duan, Meng Li*, Chunjie Zhou (School of Mathematics and Statistics, Zhengzhou University)
• Classification: math.NA cs.NA
• Publication Date: January 3, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2501.00783

Abstract

This paper investigates the solid-state dewetting phenomenon of axisymmetric thin films on axisymmetric curved-surface substrates. Under the axisymmetric assumption, the surface evolution problem on curved substrates is reduced to a curve evolution problem on a static curved substrate. Based on thermodynamic variational principles with anisotropic surface energy, the authors rigorously derive a sharp-interface model governed by anisotropic surface diffusion along with corresponding boundary conditions. The continuous system satisfies energy decay and volume conservation laws, which motivates the design of structure-preserving numerical algorithms. By introducing a symmetrized surface energy matrix, a novel symmetrized variational formulation is derived. Through careful discretization of boundary terms in the variational formulation, an unconditionally energy-stable parametric finite element approximation for the axisymmetric system is established. Employing clever correction techniques, a structure-preserving method that simultaneously maintains energy stability and volume conservation is further developed.

Research Background and Motivation

Significance of the Problem

Solid-state dewetting (SSD) refers to the instability phenomenon of solid thin films at temperatures far below the melting point, where films undergo dewetting or agglomeration, forming complex morphological changes and small particles on the substrate. This phenomenon has widespread applications in optical and magnetic devices, thin films, sensors, and catalyst formation.

Limitations of Existing Research

1. Dimensional Constraints: Existing theoretical studies primarily focus on planar substrates with limited attention to topographically patterned substrates
2. Two-Dimensional Limitation: Current research mainly addresses two-dimensional scenarios with limited investigation of SSD on three-dimensional curved-surface substrates
3. Insufficient Numerical Methods: Lack of efficient numerical methods that simultaneously preserve energy stability and volume conservation

Research Motivation

• Three-dimensional thin films and curved-surface substrates typically possess rotational symmetry, allowing simplification of complex three-dimensional SSD to curve-based systems on curved substrates
• Development of numerical methods capable of handling strong anisotropic effects is needed
• The geometric characteristics of curved substrates significantly influence thin film evolution behavior

Core Contributions

1. Theoretical Modeling: Through thermodynamic variational principles, the authors derive a sharp-interface model for three-dimensional axisymmetric thin films on axisymmetric curved-surface substrates undergoing SSD, encompassing both weak and strong anisotropy cases
2. Numerical Method Innovation:
   • Introduction of a symmetrized surface energy matrix and derivation of a novel symmetrized variational formulation
   • Design of unconditionally energy-stable parametric finite element approximation through careful discretization of boundary terms
   • Development of structure-preserving methods that simultaneously maintain volume conservation and energy stability
3. Algorithm Verification: Extensive numerical experiments validate the convergence and structure-preserving properties of the proposed methods
4. Physical Phenomenon Exploration: Investigation of diverse interesting phenomena, including particle migration on positive/negative curvature substrates, pinch-off events, and edge retraction

Detailed Methodology

Problem Definition

Investigation of solid-state dewetting evolution of axisymmetric thin films on axisymmetric curved-surface substrates, including:

• Input: Initial film morphology, substrate geometry, material parameters
• Output: Film evolution trajectory over time
• Constraints: Volume conservation, energy decay, contact line boundary conditions

Mathematical Model

Total Free Energy

The system's total free energy is expressed as: $W = \int\int_S \gamma_{FV}(N) dS + (\gamma_{FS} - \gamma_{VS}) A(\Gamma_o/\Gamma_i)$

where $S$ is the film/vapor interface, and $\Gamma_i$ and $\Gamma_o$ are the inner and outer contact lines.

Axisymmetric Reduction

Under the axisymmetric assumption, the surface is parameterized as: $S(s,\phi) := (r(s)\cos\phi, r(s)\sin\phi, z(s))$

The total energy simplifies to: $W = 2\pi\int_0^L \gamma(\theta)r ds + 2\pi(\gamma_{FS} - \gamma_{VS})\int_{c_l}^{c_r} \hat{x}(c)dc$

Sharp-Interface Model

Through thermodynamic variational principles, the governing equations are derived: $X_t \cdot n = \frac{1}{r}(r\mu_s)_s, \quad 0 < s < L(t), t > 0$

The chemical potential is defined as: $\mu = (\gamma(\theta) + \gamma''(\theta))\kappa - \frac{\gamma(\theta)z_s + \gamma'(\theta)r_s}{r}$

Boundary conditions include:

1. Contact Line Condition: Ensures the contact line remains on the curved substrate
2. Relaxed Contact Angle Condition: Describes the kinetics of contact angle evolution
3. Zero Mass Flux Condition: Guarantees volume conservation

Symmetrized Variational Formulation

Introduction of a symmetric matrix: $B(\theta) = \begin{pmatrix} \gamma(\theta) & -\gamma'(\theta) \\ \gamma'(\theta) & \gamma(\theta) \end{pmatrix} \begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix} + S(\theta)\left[\frac{1}{2}I - \frac{1}{2}\begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix}\right]$

Through this matrix, the chemical potential equation can be written as: $r\mu n = \partial_s[rB(\theta)\partial_s X] - \gamma(\theta)e_1$

Numerical Methods

Parametric Finite Element Method

The parametric finite element method (PFEM) is employed to solve the variational formulation:

• Temporal discretization: Semi-implicit scheme
• Spatial discretization: Linear finite elements
• Boundary treatment: Carefully designed boundary term discretization

Energy Stable Method

The designed numerical scheme satisfies discrete energy stability: $W(X^{m+1}) - W(X^m) \leq 0$

Structure-Preserving Method

Exact volume conservation is achieved through correction term $\delta f^{m+\frac{1}{2}}$: $V(X^{m+1}) - V(X^m) = 0$

Experimental Setup

Test Cases

The paper considers three typical initial configurations:

1. Spherical Film: Spherical film with radius 1.5 on a hemispherical substrate
2. Annular Film: Axisymmetric annular film with thickness 0.5 on a hemispherical substrate
3. Spherical Film: Spherical film with radius 1.8 on a spherical substrate

Material Parameters

• Contact line mobility: $\eta = 100$
• Surface energy parameter: $\sigma = -\frac{\sqrt{3}}{2}$
• Anisotropic surface energy: $\gamma(\theta) = 1 + \beta\cos(4\theta)$
  • $\beta = 0$: Isotropic
  • $0 \leq \beta \leq \frac{1}{15}$: Weak anisotropy
  • $\beta > \frac{1}{15}$: Strong anisotropy

Numerical Parameters

• Mesh size: $h = 2^{-7}$ to $2^{-8}$
• Time step: $\Delta t = 2^{-9}$ to $2^{-10}$
• Convergence tolerance: $10^{-8}$

Experimental Results

Convergence Verification

Numerical error is measured using manifold distance metric: $e^{h,\Delta t}(t) := |(\Omega^{h,\Delta t}(t)\backslash\Omega^r(t)) \cup (\Omega^r(t)\backslash\Omega^{h,\Delta t}(t))|$

Results demonstrate:

• Spatial Convergence Rate: Second-order convergence, consistent with theoretical predictions
• Temporal Convergence Rate: First-order convergence
• Maintains good convergence in isotropic, weakly anisotropic, and strongly anisotropic cases

Structure-Preserving Property Verification

Energy Stability

Numerical experiments confirm monotonic energy decay:

• Energy curves show decreasing trends under different time steps
• Anisotropy strength does not affect energy stability

Volume Conservation

• Energy Stable Method: Relative volume error on the order of $10^{-3}$
• Structure-Preserving Method: Achieves machine precision volume conservation (error $\sim 10^{-15}$)

Physical Phenomenon Investigation

Particle Migration Phenomenon

Observed on curved substrates:

• Positive Curvature Substrate: Small particles tend to migrate toward regions of smaller curvature
• Negative Curvature Substrate: Particle behavior influenced by anisotropy strength
• Anisotropic Effects: Strong anisotropy significantly affects migration rate

Pinch-Off Phenomenon

Long thin films undergo pinch-off during evolution:

• After pinch-off, reinitialization as two independent parts
• Eventually form smaller island and annular films

Edge Retraction

Edge retraction of semi-infinite step films:

• Films gradually climb over substrate corners
• Stronger anisotropy results in slower retraction rate

Related Work

Theoretical Modeling

• Srolovitz & Safran (1986): First proposed sharp-interface model
• Jiang et al. (2012): Phase-field method for handling topological changes
• Zhao et al. (2019): Dimension reduction model for axisymmetric geometry

Numerical Methods

• Barrett, Garcke & Nürnberg: Parametric finite element method (BGN method)
• Bao & Zhao (2021): Structure-preserving PFEM
• Li et al. (2023): Symmetrized methods

Curved Substrate Research

• Jiang et al. (2018): SSD on two-dimensional curved substrates
• Bao et al. (2024): Structure-preserving methods for two-dimensional curved substrates

This paper extends these methods to three-dimensional axisymmetric curved substrates, representing important progress in the field.

Conclusions and Discussion

Main Conclusions

1. Successfully established a complete mathematical model for SSD on axisymmetric curved substrates
2. Developed efficient numerical methods that simultaneously preserve energy stability and volume conservation
3. Revealed the significant influence of curved substrate geometry on thin film evolution behavior

Limitations

1. Axisymmetric Assumption: Restricts the geometric shapes that can be handled
2. Surface Diffusion Dominance: Neglects elastic effects and chemical reactions
3. Computational Complexity: Computational cost for three-dimensional problems remains high

Future Directions

1. Extension to non-axisymmetric geometries
2. Consideration of multi-physics coupling effects
3. Development of adaptive mesh methods
4. Application to practical industrial problems

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Derivation from first thermodynamic principles is complete and rigorous
2. Method Innovation: Symmetrized variational formulation and structure-preserving schemes are innovative
3. Numerical Stability: Unconditional energy stability avoids time step restrictions
4. Physical Significance: Numerical results reveal rich physical phenomena
5. Comprehensive Verification: Includes convergence, stability, and applicability validation

Weaknesses

1. Application Scope: Axisymmetric assumption limits method generalizability
2. Computational Efficiency: Computational cost may be high for large-scale three-dimensional problems
3. Parameter Sensitivity: Selection of certain material parameters lacks detailed discussion
4. Experimental Validation: Lacks comparison with actual experimental data

Impact

1. Academic Contribution: Provides important theoretical and numerical tools for SSD research on curved substrates
2. Practical Value: Applicable to nanodevice design and thin film engineering
3. Methodological Significance: Structure-preserving numerical methods can be generalized to other geometric evolution problems

Applicable Scenarios

• Thin film device design with axisymmetric characteristics
• Nanostructure fabrication on curved surfaces
• Morphology control of catalyst particles
• Surface engineering of optical and magnetic devices

References

The paper includes 55 high-quality references covering major works in SSD theoretical modeling, numerical methods, and experimental research, demonstrating comprehensive and authoritative literature review.

Overall Assessment: This is a high-quality numerical analysis paper with significant contributions in theoretical modeling, numerical methods, and application verification. The methods possess solid mathematical foundations and excellent numerical properties, providing powerful tools for solid-state dewetting research on three-dimensional curved substrates.