Finite difference methods for a continuous-time heterogeneous agent model with recursive utility

Resumo Técnico: Métodos de Diferenças Finitas para um Modelo de Agentes Heterogêneos em Tempo Contínuo com Utilidade Recursiva

Enunciado do Problema
Este artigo aborda o desafio computacional de resolver modelos de tempo contínuo com agentes heterogêneos (HA), especificamente modelos do tipo Aiyagari-Bewley-Huggett, onde os agentes possuem utilidade recursiva de Epstein-Zin. Ao contrário da utilidade padrão de Aversão ao Risco Relativa Constante (CRRA), a utilidade de Epstein-Zin permite a separação da aversão ao risco ($\gamma$) e da elasticidade da substituição intertemporal (EIS, $\psi$). O modelo é formulado como um sistema de Equações Diferenciais Parciais (EDPs) acopladas: equações de Hamilton-Jacobi-Bellman (HJB) para as funções de valor e equações de Fokker-Planck-Kolmogorov (FPK) para a distribuição de riqueza.

As principais dificuldades na resolução das equações HJB neste contexto são duas:

1. O Hamiltoniano assume valores infinitos para momento negativo (correspondente às restrições de empréstimo).
2. O agregador de utilidade depende não linearmente da própria função de valor $v$, uma característica específica da utilidade recursiva. Essa dependência complica os métodos de solução iterativos padrão, como o algoritmo de Howard, que são tipicamente usados para a utilidade CRRA.

O comportamento do sistema muda fundamentalmente com base no parâmetro $\theta = \frac{1-\psi^{-1}}{1-\gamma}$. O caso $\theta \ge 1$ corresponde à "resolução tardia" da incerteza, enquanto $0 < \theta < 1$ corresponde à "resolução precoce". O artigo observa que, embora o caso $\theta \ge 1$ tenha sido estudado (ex: em [4]), o caso $0 < \theta < 1$ apresenta obstáculos teóricos e numéricos significativos, particularmente no que diz respeito à existência e unicidade das soluções.

Metodologia
Os autores propõem um método de diferenças finitas utilizando esquemas monótonos (upwind) para discretizar as equações HJB com restrição de estado em uma grade limitada. O núcleo da metodologia envolve o design de algoritmos iterativos específicos adaptados aos dois regimes de preferência distintos:

1. Caso de Resolução Tardia ($\theta \ge 1$):

   • Neste regime, o sistema satisfaz um princípio de comparação, garantindo a unicidade da solução de viscosidade.
   • Os autores propõem um algoritmo Howard-Newton. Este método envolve um laço externo que atualiza a política e um laço interno que resolve o sistema linearizado resultante usando o método de Newton.
   • O algoritmo baseia-se na concavidade estrita do Hamiltoniano em relação à função de valor neste regime.

2. Caso de Resolução Precoce ($0 < \theta < 1$):

   • Aqui, o princípio de comparação padrão pode não se manter, e o Hamiltoniano é convexo em relação à função de valor.
   • Os autores definem um mapa $\Gamma$ tal que uma solução para a equação HJB é um ponto fixo de $\Gamma$. Eles provam que este mapa é monótono e invariante dentro de um intervalo específico definido por sub e supersoluções discretas (barreiras).
   • Aproveitando o teorema do ponto fixo de Tarski-Kantorovich, eles propõem um algoritmo Howard-Tarski-Kantorovich. Este método utiliza um laço externo para atualizar a função de valor com base na monotonicidade de $\Gamma$ e um laço Howard interno para resolver a etapa de avaliação da política.
   • Duas variantes são apresentadas: uma variante "ascendente" (Algoritmo 5.1) partindo de uma subsolução para convergir para a solução mínima e uma variante "descendente" (Algoritmo 5.2) partindo de uma supersolução para convergir para a solução máxima.

Principais Contribuições

• Design Algorítmico: O artigo introduz e analisa rigorosamente algoritmos iterativos convergentes para equações HJB discretizadas com utilidade recursiva. Ele distingue entre a abordagem Howard-Newton para $\theta \ge 1$ e a abordagem Howard-Tarski-Kantorovich para $0 < \theta < 1$.
• Provas de Convergência: Os autores provam a convergência destes algoritmos iterativos. Para o caso $\theta \ge 1$, isso implica a existência e unicidade da solução das equações HJB discretizadas. Para o caso $0 < \theta < 1$, as garantias de convergência asseguram a existência de pelo menos uma solução (especificamente, soluções mínima e máxima).
• Estimativas de Erro: No caso de resolução tardia ($\theta \ge 1$), o artigo deriva estimativas de erro a priori entre a solução única da equação HJB contínua e sua contraparte discreta. Os autores estabelecem uma taxa de convergência de $O(\Delta x)$, utilizando a forte regularidade das soluções contínuas para evitar a técnica de "duplicação de variáveis" frequentemente exigida na teoria de soluções de viscosidade.
• Resultados de Existência: A análise fornece provas de existência para as soluções das equações discretizadas em ambos os regimes de preferência, um resultado não trivial dada a dependência não padrão da função de utilidade em relação à função de valor.

Resultos

• Validação Numérica: Os autores apresentam experimentos numéricos para três casos de teste com parâmetros variados ($\gamma, \psi$) cobrindo tanto $\theta > 1$ quanto $0 < \theta < 1$.
• Desempenho do Algoritmo:
  • Para $\theta \ge 1$, o algoritmo Howard-Newton produz uma sequência não decrescente de funções de grade, confirmando a monotonicidade teórica.
  • Para $0 < \theta < 1$, os algoritmos Howard-Tarski-Kantorovich (tanto a variante ascendente quanto a descendente) produzem sequências monótonas convergindo para as soluções mínima e máxima, respectivamente.
• Taxa de Convergência: Testes numéricos em um cenário de taxa de juros fixa (Teste 1 HJB) confirmam a taxa de convergência teórica de $O(\Delta x)$. As curvas de erro em uma escala log-log correm paralelas a uma linha de inclinação 1.
• Comportamento da Política: As simulações ilustram políticas de poupança e distribuições de ativos. Notavelmente, em um caso de teste específico (Teste 4 HJB) com baixo risco de transição, o modelo captura um cenário onde agentes de alta renda não poupam (a poupança é zero no limite de empréstimo), um comportamento consistente com previsões teóricas sobre motivos precaucionais.

Significância e Alegações
O artigo afirma fornecer um arcabouço computacional rigoroso para modelos de HA com utilidade de Epstein-Zin, uma classe de modelos que se tornou indispensável para quantificar os efeitos distributivos da política fiscal e monetária. Os autores enfatizam que seu trabalho preenche uma lacuna na literatura ao oferecer um método matematicamente sólido para o caso $0 < \theta < 1$, onde os métodos padrão falham devido à ausência de um princípio de comparação.

Embora o artigo não prove a unicidade da solução para o caso $0 < \theta < 1$ (notando que isso permanece uma questão em aberto), os autores observam em seus experimentos que tanto o algoritmo ascendente quanto o descendente convergem para a mesma solução. Eles sugerem modestamente que isso implica unicidade na prática, mesmo que não esteja estabelecido teoricamente. O trabalho estende a aplicabilidade dos métodos de diferenças finitas, anteriormente populares para a utilidade CRRA, para o cenário mais complexo da utilidade recursiva, fornecendo uma base para a futura análise macroeconômica quantitativa.