🔢 mathematics

Existence of minimizers for the SDRI model in $\mathbb{R}^n$ : Wetting and dewetting regimes with mismatch strain

This paper extends the existence and regularity results for the stress-driven rearrangement instabilities (SDRI) model from two dimensions to arbitrary dimensions $n \geq 2$ , establishing minimizers for both wetting and dewetting regimes with mismatch strain by proving compactness and lower semicontinuity of the energy functional through a novel analysis of energy-bounded sequences.

Original authors: Shokhrukh Kholmatov, Paolo Piovano

Published 2026-01-29

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Shokhrukh Kholmatov, Paolo Piovano

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a crystal not as a perfect, rigid gem, but as a living, breathing material that is constantly trying to find its most comfortable shape. This paper is about understanding how and why these crystals change their shape when they are stressed, and proving mathematically that there is always a "best" shape they can settle into.

Here is the story of the paper, broken down into simple concepts and analogies.

The Big Picture: The Crystal's Dilemma

Think of a crystal sitting on a table (the "substrate"). The crystal wants to be happy, but it has two conflicting desires:

The Elastic Desire (The Stretch): The crystal's atoms are trying to line up perfectly with the table. But the table's atoms are spaced slightly differently than the crystal's atoms. This mismatch creates a "stretch" or tension, like a rubber band being pulled tight. The crystal wants to relieve this tension by deforming or moving.
The Surface Desire (The Smoothness): The crystal also hates having a rough, jagged surface. It wants to be smooth and compact to save energy, like a water droplet trying to become a perfect sphere.

The Conflict: If the crystal stretches to fix the mismatch, it might become jagged and rough. If it stays smooth, it stays stretched and tense. The "optimal shape" is the perfect compromise between these two forces.

The Problem: Too Many Shapes to Count

In the past, mathematicians could only solve this puzzle for flat, 2D shapes (like a drawing on a piece of paper) or for very simple scenarios where the crystal was glued tightly to the table.

But in the real world (3D), things get messy:

The crystal might peel off the table (dewetting).
It might form holes, bubbles, or sharp corners.
It might merge with other crystals or break apart.

These complex shapes cannot be described by simple graphs (like a hill or a valley). They are topological nightmares. The authors asked: "Does a perfect, stable shape actually exist for these messy, 3D, peeling, hole-filled crystals?"

The Solution: A New Mathematical Toolkit

The authors, Kholmatov and Piovano, developed a new way to look at this problem. They didn't just guess; they built a rigorous mathematical framework to prove that a solution always exists, no matter how weird the shape gets.

Here is how they did it, using everyday analogies:

1. The "Energy Scorecard"

Imagine every possible shape the crystal could take has a "score."

Elastic Score: How much tension is the crystal feeling? (High tension = bad score).
Surface Score: How rough is the surface? (Roughness = bad score).
Total Score: The sum of both. The crystal wants the lowest possible score.

The paper proves that if you keep lowering the score, you will eventually hit a bottom. You won't keep falling forever; there is a "floor" (a minimum) where the crystal settles.

2. The "Ghost" Problem (Handling the Mess)

In the real world, if you watch a crystal change shape, parts of it might suddenly disappear, reappear, or split. In math, this causes "ghosts" or gaps where the numbers break down.

The authors used a clever trick involving GSBD functions (a fancy type of mathematical function). Think of this as a way to track the crystal even when it tears or jumps.

The Analogy: Imagine trying to track a flock of birds. Sometimes they fly apart, sometimes they merge. Standard math loses track of them. The authors' method is like a special camera that can see the birds even when they split into two groups or fly through a wall, ensuring no bird is ever "lost" in the calculation.

3. The "Rigid Displacement" Fix

Sometimes, a chunk of the crystal just slides or rotates without changing its shape. In math, this can make the energy look infinite or undefined.

The Fix: The authors realized they could "subtract" these slides and rotations. It's like watching a car drive down the highway. If you only care about the shape of the car, you can ignore the fact that the car is moving forward. They mathematically "anchored" the moving parts so they could focus on the shape changes that actually matter.

4. The "Hanging Phase"

One of the hardest parts was dealing with crystals that are not touching the table at all (dewetting).

The Analogy: Imagine a drop of water sitting on a table. If it peels off, it becomes a floating sphere. The math for a drop glued to a table is different from a floating drop.
The authors created a strategy to handle both cases simultaneously. They treated the "hanging" parts of the crystal as a separate zone where the rules change slightly, ensuring the math works even when the crystal is completely detached.

The Main Results (The "Takeaways")

Existence: They proved that for any amount of crystal and any level of stress, there is definitely a best shape. The crystal won't be stuck in an infinite loop of trying to find a shape; it will find one.
Regularity (The "Smoothness" Check): They also proved that while the shape might be complex, it isn't infinitely messy.
- The Analogy: The crystal might have sharp corners or cracks, but it won't have "fractal" edges that go on forever in a chaotic way. The boundaries are well-behaved.
- They showed that the "jumps" (cracks or tears) in the material are limited and predictable.

Why This Matters (According to the Paper)

The paper doesn't claim to cure diseases or build new phones. Instead, it solves a fundamental mathematical puzzle that physicists have been struggling with.

Before: We could only prove these shapes existed in simple, 2D, or "glued-down" scenarios.
Now: We know they exist in the full, messy, 3D world, whether the crystal is peeling off, forming holes, or interacting with different materials.

Summary

Think of this paper as the architect's blueprint for a crystal's shape. Before this, we knew the crystal wanted to find a perfect shape, but we couldn't prove it would actually find one in the complex, 3D world. The authors built a new mathematical bridge that proves the crystal always finds its "happy place," no matter how much it stretches, peels, or cracks. They did this by inventing a new way to count the energy of shapes that are constantly changing and breaking apart.

Technical Summary: Existence of Minimizers for the SDRI Model in $\mathbb{R}^n$

Problem Statement
This paper addresses the mathematical modeling of Stress-Driven Rearrangement Instabilities (SDRI) in crystalline materials. SDRI phenomena occur when elastic effects destabilize the optimal free-standing crystalline equilibrium, leading to morphological changes such as surface roughness, internal cracking, dislocation nucleation, and delamination. The authors investigate a variational model where the total energy is the sum of elastic and surface energy contributions. Unlike previous studies restricted to two dimensions ( $n=2$ ) or specific regimes (e.g., wetting only, or models relying on Dirichlet boundary conditions), this work aims to establish the existence and partial regularity of minimizers for the SDRI model in any dimension $n \geq 2$ .

A critical feature of the model is the inclusion of a mismatch strain ( $M_0$ ) to represent the driving stress caused by the lattice mismatch between the free crystal and adjacent materials (the substrate). This approach avoids prescribing Dirichlet boundary conditions, which are often physically unrealistic for free-standing crystals or specific wetting/dewetting scenarios. The energy functional $F(A, u)$ depends on configurational pairs $(A, u)$ , where $A \subset \Omega$ is a set of finite perimeter representing the crystal region, and $u$ is a displacement field belonging to the space of Generalized Special Functions of Bounded Deformation ( $GSBD^2$ ).

Methodology
The authors employ the direct method of the calculus of variations. The primary challenges in this setting are the lack of compactness for sequences of admissible configurations and the lower semicontinuity of the energy functional, particularly due to the coupling between elastic and surface energies and the presence of the mismatch strain.

Topology and Convergence: The analysis is conducted with respect to a specific topology $\tau_C$ on the space of admissible configurations $\mathcal{C}$ . Convergence $(A_k, u_k) \to (A, u)$ is defined by $L^1$ -convergence of the sets $A_k$ and almost everywhere (a.e.) convergence of the displacement fields $u_k$ in $\Omega \cup S$ .
Handling the Mismatch Strain: The mismatch strain $M_0$ is discontinuous (non-zero in the container $\Omega$ , zero in the substrate $S$ ). This prevents the direct application of techniques used in models with Dirichlet boundary conditions (where rigid displacements can be controlled at the boundary).
Compactness Strategy: To prove the $\tau_C$ $τ_{C}$ -compactness of energy-equibounded sequences, the authors develop a novel partitioning strategy. They partition the substrate $S$ $S$ into connected components and the crystal region $A$ $A$ into "phases" $F^i$ $F^{i}$ interacting with specific substrate components. By applying the Poincaré-Korn inequality to each phase, they identify sequences of rigid displacements $a^i_k$ $a_{k}^{i}$ such that $u_k - a^i_k$ $u_{k} - a_{k}^{i}$ converges a.e. in the corresponding phase.
- A key technical innovation involves constructing "holes" (open sets $G^\delta_k$ ) around the interfaces between these phases and the substrate. This allows the authors to redefine the displacement field in the void regions and near the substrate boundaries to handle the potential divergence of $u_k$ (the "exceptional set" where $|u_k| \to \infty$ ) without violating the energy bounds.
- This approach circumvents the need for Dirichlet conditions, which were previously required to control the exceptional set in $GSBD$ compactness results.
Lower Semicontinuity: The proof of lower semicontinuity relies on analyzing localized Radon measures associated with the energy. The authors establish density estimates for these measures on the reduced boundary of the crystal, the jump set of the displacement, and the substrate interface. They distinguish between "interior" jumps (where the surface tension weight is 2, corresponding to cracks) and boundary jumps, utilizing slicing properties of $GSBD$ functions and relative isoperimetric inequalities.
Regularity: The partial regularity of minimizers is established by adapting techniques from the Griffith fracture model. The authors introduce a localized energy functional and prove decay estimates for minimizers, distinguishing between "blow-ups" inside the free crystal and those in the voids.

Key Contributions and Results

Existence of Minimizers (Theorem 1.1 / 2.4): The paper proves that the minimum problem for the energy functional $F$ subject to a volume constraint admits a solution in any dimension $n \geq 2$ . This extends previous results from $n=2$ to higher dimensions and includes both wetting and dewetting regimes.
Generalization of Boundary Conditions: The methodology successfully handles the mismatch strain formulation, which is physically more relevant for free crystals than Dirichlet conditions. The authors demonstrate that their arguments also recover existence results for models with Dirichlet boundary conditions and material voids (Theorem 2.9, Remark 2.10).
Compactness in $GSBD$ without Dirichlet Conditions: A significant technical contribution is the proof of compactness for energy-equibounded sequences in the absence of fixed boundary conditions. The authors show that even when the displacement field may diverge to infinity on a set of positive measure (the "hanging phase"), a convergent subsequence can be extracted by carefully partitioning the domain and redefining displacements.
Partial Regularity (Theorem 1.2 / 2.7): The minimizers $(A, u)$ $(A, u)$ possess partial regularity properties:
- The jump set $J_u$ and the reduced boundary $\partial^* A$ are essentially closed (up to an $\mathcal{H}^{n-1}$ -negligible set).
- Density estimates for the jump set are established, providing bounds on the local energy concentration.
- "Hanging" components of the crystal (those not in contact with the substrate) are shown to have a minimum volume and a rigid displacement field relative to the mismatch strain.

Significance and Scope
The paper claims to extend the mathematical framework for SDRI models to physically relevant dimensions ( $n=3$ ) and regimes (dewetting) that were previously inaccessible due to technical limitations in handling mismatch strains and the lack of boundary conditions. By avoiding graph-based assumptions (which cannot represent holes, overhangs, or pinch-off events), the model captures complex topological changes observed in experiments.

The authors explicitly state that their work establishes existence and partial regularity for the two-phase setting (crystal and substrate). They acknowledge that several directions remain open, including:

Extension to multiphase models involving multiple interacting phases.
Incorporation of contact-angle conditions at the substrate interface.
Further regularity results beyond partial regularity (e.g., full smoothness).
Analysis of quasistatic and dynamic evolutions (gradient flows or dynamic models with inertia).

The work is presented as a rigorous variational analysis that bridges the gap between physical descriptions of stress-driven instabilities and mathematical existence theory, providing a foundation for future investigations into more complex material behaviors.