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

Resumo Técnico: Existência de Minimizadores para o Modelo SDRI em $\mathbb{R}^n$

Enunciado do Problema
Este artigo aborda a modelagem matemática de Instabilidades de Rearranjo Impulsionadas por Tensão (SDRI, do inglês Stress-Driven Rearrangement Instabilities) em materiais cristalinos. Os fenômenos de SDRI ocorrem quando efeitos elásticos desestabilizam o equilíbrio cristalino livre ideal, levando a mudanças morfológicas como rugosidade superficial, fissuração interna, nucleação de discordâncias e delaminação. Os autores investigam um modelo variacional onde a energia total é a soma das contribuições de energia elástica e de superfície. Diferentemente de estudos anteriores restritos a duas dimensões ($n=2$) ou regimes específicos (por exemplo, apenas molhagem, ou modelos que dependem de condições de contorno de Dirichlet), este trabalho visa estabelecer a existência e a regularidade parcial de minimizadores para o modelo SDRI em qualquer dimensão $n \geq 2$.

Uma característica crítica do modelo é a inclusão de uma deformação de desajuste ($M_0$) para representar a tensão motriz causada pelo desajuste de rede entre o cristal livre e os materiais adjacentes (o substrato). Esta abordagem evita prescrever condições de contorno de Dirichlet, que são frequentemente fisicamente irreais para cristais livres ou cenários específicos de molhagem/desmolhagem (wetting/dewetting). O funcional de energia $F(A, u)$ depende de pares configuracionais $(A, u)$, onde $A \subset \Omega$ é um conjunto de perímetro finito representando a região cristalina, e $u$ é um campo de deslocamento pertencente ao espaço de Funções Especiais de Deformação Limitada Generalizadas ($GSBD^2$).

Metodologia
Os autores empregam o método direto do cálculo de variações. Os principais desafios neste cenário são a falta de compacidade para sequências de configurações admissíveis e a semicontinuidade inferior do funcional de energia, particularmente devido ao acoplamento entre as energias elástica e de superfície e à presença da deformação de desajuste.

1. Topologia e Convergência: A análise é conduzida com relação a uma topologia específica $\tau_C$ no espaço de configurações admissíveis $\mathcal{C}$. A convergência $(A_k, u_k) \to (A, u)$ é definida pela convergência $L^1$ dos conjuntos $A_k$ e convergência quase em toda parte (q.t.p.) dos campos de deslocamento $u_k$ em $\Omega \cup S$.
2. Tratamento da Deformação de Desajuste: A deformação de desajuste $M_0$ é descontínua (não nula no recipiente $\Omega$, zero no substrato $S$). Isso impede a aplicação direta de técnicas utilizadas em modelos com condições de contorno de Dirichlet (onde deslocamentos rígidos podem ser controlados na fronteira).
3. Estratégia de Compacidade: Para provar a $\tau_C$-compacidade de sequências equibundadas em energia, os autores desenvolvem uma estratégia de partição inovadora. Eles particionam o substrato $S$ em componentes conexas e a região cristalina $A$ em "fases" $F^i$ que interagem com componentes específicas do substrato. Ao aplicar a desigualdade de Poincaré-Korn para cada fase, eles identificam sequências de deslocamentos rígidos $a^i_k$ tais que $u_k - a^i_k$ converge q.t.p. na fase correspondente.
   • Uma inovação técnica fundamental envolve a construção de "buracos" (conjuntos abertos $G^\delta_k$) ao redor das interfaces entre essas fases e o substrato. Isso permite que os autores redefinam o campo de deslocamento nas regiões de vácuo e próximo às fronteiras do substrato para lidar com a potencial divergência de $u_k$ (o "conjunto excepcional" onde $|u_k| \to \infty$) sem violar os limites de energia.
   • Esta abordagem contorna a necessidade de condições de Dirichlet, que eram anteriormente exigidas para controlar o conjunto excepcional em resultados de compacidade de $GSBD$.
4. Semicontinuidade Inferior: A prova de semicontinuidade inferior baseia-se na análise de medidas de Radon localizadas associadas à energia. Os autores estabelecem estimativas de densidade para essas medidas na fronteira reduzida do cristal, no conjunto de salto do deslocamento e na interface do substrato. Eles distinguem entre saltos "interiores" (onde o peso da tensão superficial é 2, correspondendo a fissuras) e saltos de fronteira, utilizando propriedades de fatiamento (slicing) de funções $GSBD$ e desigualdades isoperimétricas relativas.
5. Regularidade: A regularidade parcial dos minimizadores é estabelecida adaptando técnicas do modelo de fratura de Griffith. Os autores introduzem um funcional de energia localizado e provam estimativas de decaimento para os minimizadores, distinguindo entre "explosões" (blow-ups) dentro do cristal livre e aquelas nos vazios.

Principais Contribuições e Resultados

• Existência de Minimizadores (Teorema 1.1 / 2.4): O artigo prova que o problema de minimização para o funcional de energia $F$, sujeito a uma restrição de volume, admite uma solução em qualquer dimensão $n \geq 2$. Isso estende resultados anteriores de $n=2$ para dimensões superiores e inclui tanto regimes de molhagem quanto de desmolhagem.
• Generalização das Condições de Contorno: A metodologia lida com sucesso com a formulação de deformação de desajuste, que é fisicamente mais relevante para cristais livres do que as condições de Dirichlet. Os autores demonstram que seus argumentos também recuperam resultados de existência para modelos com condições de contorno de Dirichlet e vazios de material (Teorema 2.9, Observação 2.10).
• Compacidade em $GSBD$ sem Condições de Dirichlet: Uma contribuição técnica significativa é a prova de compacidade para sequências equibundadas em energia na ausência de condições de contorno fixas. Os autores mostram que mesmo quando o campo de deslocamento pode divergir para o infinito em um conjunto de medida positiva (a "fase pendente"), uma subsequência convergente pode ser extraída através da partição cuidadosa do domínio e da redefinição dos deslocamentos.
• Regularidade Parcial (Teorema 1.2 / 2.7): Os minimizadores $(A, u)$ possuem propriedades de regularidade parcial:
  • O conjunto de salto $J_u$ e a fronteira reduzida $\partial^* A$ são essencialmente fechados (salvo um conjunto $\mathcal{H}^{n-1}$-nulo).
  • Estimativas de densidade para o conjunto de salto são estabelecidas, fornecendo limites sobre a concentração local de energia.
  • Componentes "pendentes" do cristal (aqueles que não estão em contato com o substrato) mostram possuir um volume mínimo e um campo de deslocamento rígido em relação à deformação de desajuste.

Significância e Escopo
O artigo afirma estender o arcabouço matemático para modelos de SDRI para dimensões fisicamente relevantes ($n=3$) e regimes (desmolhagem) que eram anteriormente inacessíveis devido a limitações técnicas no tratamento de deformações de desajuste e à falta de condições de contorno. Ao evitar suposições baseadas em grafos (que não podem representar buracos, saliências ou eventos de estrangulamento/separação), o modelo captura mudanças topológicas complexas observadas em experimentos.

Os autores declaram explicitamente que seu trabalho estabelece existência e regularidade parcial para o cenário de duas fases (cristal e substrato). Eles reconhecem que diversas direções permanecem em aberto, incluindo:

• Extensão para modelos de múltiplas fases envolvendo múltiplas fases interagentes.
• Incorporação de condições de ângulo de contato na interface do substrato.
• Novos resultados de regularidade além da regularidade parcial (por exemplo, suavidade total).
• Análise de evoluções quase-estáticas e dinâmicas (fluxos de gradiente ou modelos dinâmicos com inércia).

O trabalho é apresentado como uma análise variacional rigorosa que faz a ponte entre as descrições físicas de instabilidades induzidas por tensão e a teoria de existência matemática, fornecendo uma base para investigações futuras sobre comportamentos de materiais mais complexos.