The local Turnpike Property in Mean Field Control and Games with quadratic Hamiltonian

1. Formulação do Problema

O artigo investiga o comportamento de longo prazo das soluções dos sistemas de Jogos de Campo Médio (MFG) e dos problemas de Controle de Campo Médio (MFC) com um Hamiltoniano quadrático à medida que o horizonte temporal $T \to +\infty$.

O sistema específico estudado é:
$$\begin{cases} -\partial_t u - \Delta u + \frac{|Du|^2}{2} - f(x, m) + \delta u = 0 & \text{em } (0, T) \times \Omega \\ \partial_t m - \Delta m - \text{div}(m Du) = 0 & \text{em } (0, T) \times \Omega \\ m(0, x) = m_0(x), \quad u(T, x) = u_T(x) & \text{em } \Omega \end{cases}$$
onde $\Omega$ é ou o toro plano $\mathbb{T}^n$ ou $\mathbb{R}^n$, e $\delta \geq 0$ é um fator de desconto.

A Questão Central:
Dado um equilíbrio estacionário $(\bar{u}, \bar{m})$ que resolve o sistema estacionário associado, sob quais condições este equilíbrio atrai soluções dinâmicas $(u, m)$ quando $T \to \infty$? Especificamente, os autores visam estabelecer uma propriedade do túnel local: se os dados inicial e terminal estiverem suficientemente próximos do estado estacionário, a solução passará a maior parte do tempo exponencialmente próxima de $(\bar{u}, \bar{m})$?

O Desafio:
Os resultados padrão de túnel na literatura de MFG geralmente dependem da condição de monotonicidade de Lasry-Lions ($f_m \geq c > 0$), que garante unicidade global do estado estacionário e estabilidade global. Este artigo aborda o caso mais difícil onde a monotonicidade está ausente, permitindo múltiplos estados estacionários e potencial não-unicidade. O objetivo é provar estabilidade para soluções estacionárias específicas sem assumir unicidade global.

2. Metodologia

Os autores empregam uma abordagem perturbativa combinada com estimativas de energia e um argumento de ponto fixo.

A. Linearização e Simetria

A análise gira em torno das perturbações $\mu = m - \bar{m}$ e $v = u - \bar{u}$. Os autores linearizam o sistema em torno do estado estacionário. Uma observação crucial é que, para o Hamiltoniano quadrático, o sistema linearizado possui uma propriedade de simetria específica devido à identidade $D\bar{m} = -\bar{m} D\bar{u}$.

Eles reformulam o sistema linearizado como um sistema Hamiltoniano envolvendo operadores $A$ e $BB^*$:
$$\begin{pmatrix} \dot{\tilde{\mu}} \\ \dot{\tilde{v}} \end{pmatrix} = \begin{pmatrix} A - \frac{\delta}{2}I & -BB^* \\ -Q & -A^* + \frac{\delta}{2}I \end{pmatrix} \begin{pmatrix} \tilde{\mu} \\ \tilde{v} \end{pmatrix}$$
onde $Q$ representa o termo de acoplamento $f_m$. Nos casos monótonos padrão, $Q$ é definido positivo. Aqui, $Q$ pode não ser positivo.

B. A Hipótese Chave (S)

Para substituir a condição global de monotonicidade, os autores introduzem uma hipótese de estabilidade local (S) sobre a solução estacionária $(\bar{u}, \bar{m})$. Esta hipótese consiste em três partes:

1. Identidade: $D\bar{m} = -\bar{m} D\bar{u}$ (vale para soluções estacionárias de Fokker-Planck).
2. Desigualdade de Poincaré: $\bar{m}$ satisfaz uma desigualdade de Poincaré ponderada.
3. Positividade Estrita de Segunda Ordem: Existe $\eta \in (0, 1)$ tal que para todo $\mu$ com média zero:
   $$\int_\Omega \left( f_m(x, \bar{m})\mu^2 + \bar{m} \left| D\left(\frac{\mu}{\bar{m}}\right) \right|^2 + \delta \frac{\mu^2}{\bar{m}} \right) dx \geq \eta \int_\Omega \bar{m} \left| D\left(\frac{\mu}{\bar{m}}\right) \right|^2 dx$$
   Esta condição é interpretada como a positividade do autovalor principal do operador estacionário linearizado. É mais fraca que a monotonicidade global, mas forte o suficiente para garantir não-degenerescência local.

C. Construção da Função de Lyapunov

Em vez de depender da monotonicidade do acoplamento, os autores constroem um funcional de Lyapunov modificado (denotado $\Phi$ ou $\tilde{\Phi}$):
$$\tilde{\Phi}(t) = \int_\Omega e^{-\delta t} \left( \mu v + \frac{\mu^2}{\bar{m}} \right) dx$$
Ao diferenciar $\tilde{\Phi}$ e utilizar a simetria dos operadores e a hipótese (S), eles derivam uma desigualdade diferencial da forma:
$$-\dot{\tilde{\Phi}}(t) \geq c \tilde{\Phi}(t)$$
Esta desigualdade implica o decaimento exponencial do funcional, o que impulsiona o comportamento de túnel. A prova envolve o tratamento cuidadoso de termos não lineares de ordem inferior e do fator de desconto $\delta$.

D. Estimativas A Priori e Ponto Fixo

1. Estimativas $L^2$: Eles primeiro estabelecem estimativas de decaimento exponencial ponderado em $L^2$ para $\mu$ e $Dv$ usando o funcional de Lyapunov.
2. Estimativas $L^\infty$: Usando a teoria de regularidade parabólica (estimativas de Schauder, transformada de Cole-Hopf para a equação de Hamilton-Jacobi), eles elevam o decaimento em $L^2$ para decaimento em $L^\infty$.
3. Argumento de Ponto Fixo: Eles definem um mapeamento $F$ em um espaço de funções próximas à solução estacionária. Usando o teorema do ponto fixo de Leray-Schauder, eles provam a existência de uma solução para o sistema não linear que satisfaz as estimativas a priori derivadas, desde que os dados iniciais/terminais sejam suficientemente pequenos.

3. Contribuições Principais

1. Relaxação da Monotonicidade: O artigo estabelece com sucesso uma propriedade de túnel local sem a condição global de monotonicidade de Lasry-Lions. Isso permite o estudo de sistemas com múltiplos equilíbrios estacionários.
2. Estabilidade Local via Condições Espectrais: Os autores identificam que a positividade estrita do autovalor principal do sistema estacionário linearizado (Hipótese S) é suficiente para garantir estabilidade exponencial local. Isso conecta a teoria espectral do problema estacionário à dinâmica de longo prazo do problema dependente do tempo.
3. Exploração da Simetria: O trabalho destaca como a estrutura específica do Hamiltoniano quadrático permite uma simetria nos operadores linearizados que é essencial para construir o funcional de Lyapunov em configurações não monótonas.
4. Existência de Trajetórias Estáveis: O artigo prova que, para qualquer solução estacionária que satisfaça (S), existe pelo menos uma trajetória dinâmica (para horizontes finito e infinito) que permanece exponencialmente próxima dela, desde que as condições de contorno sejam suficientemente próximas.

4. Resultados Principais

• Teorema 1.1 (Estimativas A Priori): Sob a hipótese (S), se uma solução $(u, m)$ estiver inicialmente próxima de $(\bar{u}, \bar{m})$, ela satisfaz uma propriedade de contratividade: a distância ao estado estacionário decai exponencialmente no tempo (com taxas dependendo de $\delta$ e da dimensão).
• Teorema 1.2 (Existência em Horizonte Finito): No toro $\mathbb{T}^n$, se a distribuição inicial $m_0$ e o custo terminal $u_T$ estiverem suficientemente próximos de $\bar{m}$ e $\bar{u}$, existe uma solução clássica $(u_T, m_T)$ para o sistema MFG que satisfaz a propriedade de túnel exponencial local:
  $$\|m_T(t, \cdot) - \bar{m}(\cdot)\|_{L^\infty} \leq C (e^{-\sigma_1 t} + e^{-\sigma_2 (T-t)})$$
• Teorema 4.3 (Horizonte Infinito): Ao passar ao limite $T \to \infty$, os autores estabelecem a existência de um equilíbrio dinâmico em $[0, \infty)$ que converge exponencialmente para o estado estacionário $(\bar{u}, \bar{m})$.

5. Significado e Implicações

• Avanço Teórico: Este trabalho preenche a lacuna entre a teoria de controle ótimo (onde as propriedades de túnel são bem estudadas) e os Jogos de Campo Médio, abordando especificamente o regime não monótono, que é prevalente em muitos modelos físicos e econômicos (por exemplo, modelos com congestionamento ou atração-repulsão).
• Tratamento da Não-Unicidade: Fornece um quadro rigoroso para analisar a estabilidade em sistemas onde múltiplos equilíbrios coexistem. Esclarece que, mesmo que um sistema tenha múltiplos estados estacionários, específicos podem ser localmente estáveis e atrativos, desde que satisfaçam uma condição de não-degenerescência (S).
• Aplicações: Os resultados aplicam-se a modelos onde o acoplamento $f(x, m)$ não é monótono (por exemplo, $f_m$ pode ser negativo em algumas regiões), o que é comum em modelos de dinâmica de multidões com atração ou enxameamento biológico.
• Direções Futuras: Os autores observam que, embora o Hamiltoniano quadrático seja crucial para a simetria usada aqui, a abordagem pode se estender a Hamiltonianos mais gerais. Eles também destacam a questão em aberto de se uma propriedade de túnel global pode existir na presença de múltiplos estados estacionários.

Em resumo, Cirant e De Bernardi fornecem um quadro analítico robusto para entender a estabilidade de longo prazo dos Jogos de Campo Médio perto de estados estacionários, indo além da suposição restritiva de monotonicidade global ao aproveitar propriedades espectrais locais e a geometria específica dos Hamiltonianos quadráticos.