Microergodicity implies orthogonality of Matérn fields on bounded domains in $\mathbb{R}^4$

Resumo Técnico: Ortogonalidade de Campos de Matérn em Domínios Limitados em $\mathbb{R}^4$

1. O Problema e o Contexto

O artigo aborda um problema fundamental na estatística espacial: a identificabilidade de parâmetros de covariância para campos aleatórios gaussianos estacionários do tipo Matérn sob o regime de assintótica de domínio fixo (infill asymptotics). Neste regime, observa-se uma única realização de um campo aleatório em um domínio limitado $D \subset \mathbb{R}^d$, enquanto a densidade de amostragem aumenta indefinidamente.

O foco central é determinar se dois modelos de Matérn com parâmetros diferentes induzem medidas gaussianas equivalentes (onde os parâmetros não podem ser consistentemente estimados a partir de uma única trajetória) ou mutuamente singulares (onde os parâmetros são distinguíveis).

Para campos Matérn com suavidade $\nu$ conhecida, o parâmetro microergódico é definido como $m = \sigma^2 \alpha^{2\nu}$, onde $\sigma^2$ é a variância e $\alpha$ é o parâmetro de alcance.

• Resultados conhecidos:
  • Se $d \leq 3$: Medidas com o mesmo $m$ são equivalentes (os parâmetros $\sigma$ e $\alpha$ não podem ser separados).
  • Se $d > 4$: Medidas com o mesmo $m$ mas diferentes $\alpha$ são mutuamente singulares (os parâmetros podem ser separados).
• O Caso Crítico Aberto: O caso $d = 4$ permanecia sem solução. A questão era se a discrepância entre os modelos, que cresce muito lentamente em dimensão 4, seria suficiente para gerar singularidade.

2. Metodologia e Abordagem

O autor desenvolve uma abordagem baseada no domínio espectral, utilizando coeficientes de Fourier localizados para detectar discrepâncias sutis na estrutura de covariância.

• Coeficientes de Fourier Localizados:
  Em vez de trabalhar diretamente com incrementos no espaço físico (como feito em trabalhos anteriores de Anderes para $d>4$), o autor define coeficientes de Fourier ponderados por uma função de corte suave ($\chi$) no domínio $D$:
  $$X_k = \int_D \chi(t) e^{-ik \cdot t} Y(t) dt$$
  Isso permite isolar o comportamento de alta frequência do campo.

• Análise do Mismatch Espectral:
  Sob a condição de igualdade do parâmetro microergódico ($\sigma_1^2 \alpha_1^{2\nu} = \sigma_2^2 \alpha_2^{2\nu}$), as densidades espectrais $f_1(\xi)$ e $f_2(\xi)$ coincidem no termo dominante de alta frequência ($|\xi| \to \infty$). A primeira discrepância aparece no termo de segunda ordem.
  O autor demonstra que a razão de variâncias normalizada $\delta_k = \frac{v_2(k)}{v_1(k)} - 1$ decai como $|\delta_k| \asymp |k|^{-2}$ em $\mathbb{R}^4$.

• Escala de Informação Logarítmica:
  A soma das discrepâncias ao longo de uma casca de frequências $\Lambda_N$ (onde $K_0 \leq \|k\|_\infty \leq N$) define a escala de informação $L_N$:
  $$L_N = \sum_{k \in \Lambda_N} \delta_k^2 \asymp \sum_{k \in \Lambda_N} |k|^{-4} \asymp \log N$$
  Em dimensões $d < 4$, essa soma convergiria (informação finita $\to$ equivalência). Em $d > 4$, divergiria polinomialmente. Em $d=4$, a divergência é logarítmica, o que é o ponto crítico que permite a identificação.

• Estatística de Separação (Score Normalizado):
  O autor constrói uma estatística $T_N$ baseada em um escore de Whittle normalizado:
  $$T_N = \frac{1}{L_N} \sum_{k \in \Lambda_N} \delta_k \left( \frac{|X_k|^2}{v_1(k)} - 1 \right)$$
  Sob o modelo 1, $T_N \to 0$; sob o modelo 2, $T_N \to 1$. A prova envolve controlar as covariâncias fora da diagonal (devido à localização espacial) para mostrar que a variância da estatística decai suficientemente rápido ($Var(T_N) \lesssim L_N^{-1}$), permitindo a convergência quase certa ao longo de uma subsequência esparsa.

3. Principais Contribuições

1. Resolução do Caso Crítico $d=4$: O artigo prova que, para campos de Matérn estacionários em $\mathbb{R}^4$, se dois modelos possuem o mesmo parâmetro microergódico $m$ mas diferentes parâmetros de alcance $\alpha$, as medidas gaussianas induzidas em qualquer conjunto de observação denso e contável são mutuamente singulares.
2. Analogia Espectral de Alta Ordem: Introduz uma técnica que é o análogo espectral do método de incrementos de ordem superior de Anderes. Enquanto Anderes isola coeficientes de covariância no espaço físico, Pillai isola o primeiro mismatch espectral não nulo em altas frequências.
3. Tratamento de Dependência Fora da Diagonal: Uma contribuição técnica crucial é o manejo rigoroso das covariâncias entre coeficientes de Fourier localizados que não são independentes (devido à janela de corte $\chi$), provando que o decaimento rápido da transformada de Fourier de $\chi$ permite controlar o erro de aproximação diagonal.
4. Conexão com o Limite de Feldman-Hájek: O trabalho conecta a divergência logarítmica da soma de informações ($L_N$) com o critério de Feldman-Hájek para singularidade de medidas gaussianas, mostrando que $d=4$ é exatamente a fronteira onde a soma das raízes quadradas das discrepâncias diverge.

4. Resultados Principais

• Teorema 1: Seja $D \subset \mathbb{R}^4$ um domínio limitado e $S \subset D$ um conjunto contável e denso. Se dois campos de Matérn $Y_1$ e $Y_2$ satisfazem $\sigma_1^2 \alpha_1^{2\nu} = \sigma_2^2 \alpha_2^{2\nu}$ e $\alpha_1 \neq \alpha_2$, então as leis das observações em $S$ são mutuamente singulares ($P_1^S \perp P_2^S$).
• Estudo de Simulação: O autor realiza simulações numéricas em $\mathbb{R}^4$ que confirmam que a estatística $T_N$ separa efetivamente os modelos, concentrando-se em 0 para o modelo gerador e em 1 para o modelo alternativo, mesmo com a separação fraca (logarítmica).
• Estimativa de Parâmetros: O artigo sugere que a mesma estrutura espectral pode ser usada para estimar o parâmetro $\alpha$ via pseudo-verossimilhança de Whittle, demonstrando empiricamente a consistência do estimador em dimensão 4.

5. Significado e Impacto

• Completude Teórica: O trabalho fecha a lacuna teórica sobre a identificabilidade de parâmetros de Matérn em espaços euclidianos, estabelecendo que a dimensão 4 é a fronteira exata onde a separação de parâmetros se torna possível, embora a taxa de convergência seja lenta (logarítmica).
• Implicações Práticas: Para estatísticos espaciais, isso implica que, em dimensões 4 ou superiores, é teoricamente possível estimar tanto a variância quanto o alcance de um campo de Matérn a partir de uma única realização densa, ao contrário do que ocorre em dimensões 1, 2 e 3.
• Novas Ferramentas: A metodologia de "sondagem espectral localizada" (localized spectral probing) oferece uma nova ferramenta para estudar a equivalidade e ortogonalidade de medidas gaussianas em outros modelos espaciais complexos, especialmente aqueles onde a análise no domínio físico é difícil.

Em resumo, o artigo demonstra que, embora a informação adicional em dimensão 4 seja marginal (crescendo apenas logaritmicamente), ela é suficiente para quebrar a equivalência entre modelos de Matérn com o mesmo parâmetro microergódico, estabelecendo a singularidade mútua das medidas.