Stable determination of the first order perturbation of the biharmonic operator from partial data

Resumo Técnico: Determinação Estável da Perturbação de Primeira Ordem do Operador Biharmônico a partir de Dados Parciais

Enunciado do Problema
Este artigo aborda um problema inverso de valor de contorno para o operador biharmônico perturbado em um domínio limitado $\Omega \subseteq \mathbb{R}^n$ ($n \ge 3$) com uma fronteira suave $\partial\Omega$. O operador é definido como $L_{A,q} = \Delta^2 + A(x) \cdot D + q(x)$, onde $D = i^{-1}\nabla$, $A$ é um campo vetorial em $W^{1,\infty}(\Omega, \mathbb{C}^n)$, e $q$ é um potencial em $L^\infty(\Omega, \mathbb{C})$. O estudo foca no cenário de dados parciais: recuperar os coeficientes $A$ e $q$ a partir de medições tomadas apenas em subconjuntos abertos arbitrariamente pequenos $\gamma_1, \gamma_2 \subset \partial\Omega$.

O problema assume que os coeficientes são a priori conhecidos em uma vizinhança $\omega_0$ da fronteira. Os dados são fornecidos por um mapa Dirichlet-para-Neumann (DtN) parcial $\Lambda_{\gamma_1, \gamma_2}^{A,q}$, que mapeia dados de Dirichlet $(f, g)$ suportados em $\gamma_1$ para as derivadas normais $(\partial_\nu u, \partial_\nu(\Delta u))$ medidas em $\gamma_2$. O objetivo principal é estabelecer estimativas de estabilidade (quantificando como erros nas medições afetam a reconstrução dos coeficientes) tanto para a perturbação de primeira ordem $A$ quanto para a perturbação de ordem zero $q$.

Metodologia
Os autores empregam o método de soluções de Óptica Geométrica Complexa (CGO), uma técnica padrão para problemas inversos de operadores elípticos. A estratégia de prova envolve vários componentes técnicos fundamentais:

1. Construção de Soluções CGO: Os autores constroem soluções $u_1$ e $u_2$ para as equações adjunta e direta, respectivamente, da forma $u(x) = e^{ix\cdot\zeta/h}(a(x) + r(x))$, onde $\zeta$ é um vetor complexo satisfazendo $\zeta \cdot \zeta = 0$, $h$ é um parâmetro semiclássico, e $r(x)$ é um termo de resto que desaparece quando $h \to 0$.
2. Identidade Integral e Fórmula de Green: Ao aplicar a fórmula de Green à diferença de dois operadores e utilizar o mapa DtN parcial, os autores derivam uma identidade integral que relaciona a diferença dos coeficientes ($A_1 - A_2$ e $q_1 - q_2$) com a diferença dos mapas DtN.
3. Continuidade Única Quantitativa Local: Uma inovação crucial neste trabalho é a derivação de uma estimativa de continuidade única quantitativa local para soluções da equação biharmônica não homogênea com condições de contorno homogêneas. Isso é alcançado usando estimativas de Carleman para o operador biharmônico semiclássico $h^4\Delta^2$. Este passo permite que os autores controlem as normas das soluções especificamente sobre o subconjunto de medição $\gamma_2$, em vez de toda a fronteira, o que é essencial para o cenário de dados parciais.
4. Decomposição de Campos Vetoriais: Para recuperar o campo vetorial $A$ a partir de sua derivada exterior $dA$ (que é naturalmente obtida via transformada de Fourier da identidade integral), os autores utilizam uma decomposição do tipo Helmholtz $A = A_{sol} + \nabla \phi$. Eles estimam a parte solenoide $A_{sol}$ usando o limite de $dA$ e a parte de potencial $\nabla \phi$ resolvendo um problema de valor de contorno para $\phi$ e aplicando desigualdades de imersão de Sobolev e interpolação.
5. Otimização de Parâmetros: As estimativas de estabilidade são derivadas otimizando o parâmetro semiclássico $h$ e o corte de frequência $\rho$ no domínio de Fourier, equilibrando os termos de erro provenientes dos restos CGO e os pesos exponenciais nas estimativas de Carleman.

Principais Contribuições e Resultados
O artigo estabelece estimativas de estabilidade do tipo logarítmico e do tipo duplo-logarítmico, que são características de problemas inversos mal-postos.

• Estabilidade para a Perturbação de Primeira Ordem ($A$):
  O Teorema 1.1 estabelece uma estimativa de estabilidade do tipo log para o campo vetorial $A$. Especificamente, a norma $L^\infty$ da diferença $A_1 - A_2$ é limitada por uma soma de uma potência da diferença do DtN e uma potência negativa do logaritmo da diferença do DtN:
  $$\|A_1 - A_2\|_{L^\infty} \lesssim \|\Lambda_1 - \Lambda_2\|^{\mu_1} + |\log \|\Lambda_1 - \Lambda_2\||^{-\mu_2}$$
  Este resultado complementa resultados de unicidade anteriores ao fornecer um limite quantitativo de estabilidade.

• Estabilidade para a Perturbação de Ordem Zero ($q$) com Termo de Primeira Ordem:
  O Teorema 1.2 estabelece uma estimativa de estabilidade do tipo log-log para o potencial $q$ quando o termo de primeira ordem $A$ está presente:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim |\log |\log \|\Lambda_1 - \Lambda_2\|||^{-\mu'}$$
  O Corolário 1.3 estende isso para uma estimativa $L^\infty$ sob suposições adicionais de regularidade em $q$.

• Estabilidade para a Perturbação de Ordem Zero ($q$) sem Termo de Primeira Ordem:
  O Teorema 1.4 mostra que, na ausência da perturbação de primeira ordem ($A=0$), a estabilidade para $q$ melhora de tipo log-log para tipo log:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim \|\Lambda_1 - \Lambda_2\|^{2/3} + |\log \|\Lambda_1 - \Lambda_2\||^{-2/(n+2)}$$
  Isso destaca a dificuldade específica introduzida pelo termo de primeira ordem no operador biharmônico.

Significância e Alegações
Os autores afirmam que este trabalho fornece os primeiros resultados de estabilidade para a perturbação de primeira ordem do operador biharmônico utilizando dados parciais. Embora a unicidade para este problema tenha sido estabelecida anteriormente (citando [42]), a análise de estabilidade é significativamente mais desafiadora.

O artigo enfatiza duas principais novidades técnicas:

1. Continuidade Única Local: A derivação de uma estimativa de continuidade única quantitativa local para operadores elípticos de ordem superior (especificamente o operador biharmônico) usando estimativas de Carleman. Os autores observam que este é o primeiro tal resultado para operadores de ordem superior neste contexto específico de dados parciais.
2. Tratamento de Dados Parciais: A capacidade de derivar estimativas de estabilidade usando medições em subconjuntos abertos arbitrariamente pequenos da fronteira, em vez de exigir dados em toda a fronteira ou configurações geométricas específicas (como partes inacessíveis planas).

Os autores posicionam seus resultados como um complemento aos teoremas de unicidade existentes e um passo à frente na análise de estabilidade de problemas inversos de ordem superior, traçando paralelos com métodos usados para os operadores de Schrödinger e de Schrödinger magnético, enquanto abordam os desafios únicos impostos pela estrutura do operador biharmônico e pela falta de invariância de calibre para sua perturbação de primeira ordem.