Numerical analysis for a unified 2 factor model of structural and reduced form types for corporate bonds with fixed discrete coupon

Resumo Técnico: Análise Numérica para um Modelo Unificado de Dois Fatores de Títulos Corporativos com Cupons Discretos Fixos

Definição do Problema
O artigo aborda a precificação de títulos corporativos sujeitos a default que pagam cupons discretos fixos. Embora exista um trabalho teórico significativo para títulos de cupom zero (Merton, 1977) e modelos estruturais para títulos com cupom (Geske, 1977; Agliardi, 2011), há uma lacuna em soluções analíticas para um modelo unificado que combine abordagens estruturais e de forma reduzida especificamente para títulos com cupons discretos fixos.

Os autores observam que, embora um modelo unificado de um fator para cupons discretos fixos possua uma fórmula analítica de precificação, e um modelo unificado de dois fatores para cupons discretos descontados tenha sido resolvido, o modelo unificado de dois fatores para cupons discretos fixos parece carecer de uma solução analítica. Isso se deve à complexidade das Equações Diferenciais Parciais (EDPs) regentes, que são equações parabólicas envolvendo três variáveis (valor da firma, taxa de juros de curto prazo e tempo) e derivadas mistas. Além disso, o modelo deve considerar dois tipos de default:

1. Default Esperado: Ocorre apenas em datas de cupons discretos, quando o patrimônio líquido da firma é insuficiente para cobrir a dívida e os cupons.
2. Default Inesperado: Pode ocorrer a qualquer momento, governado por uma intensidade de default (taxa de risco) constante.

Metodologia
Para resolver este problema, os autores propõem um modelo unificado de dois fatores sob a medida de martingale livre de risco, incorporando:

• Modelo Vasicek: Para a taxa de juros de curto prazo ($r$).
• Movimento Browniano Geométrico: Para o valor da firma ($V$), com um rendimento de dividendos constante.
• Correlação: Um coeficiente de correlação ($\rho$) entre os processos de Wiener que impulsionam a taxa de juros e o valor da firma.

O problema de precificação é formulado como um conjunto de problemas de valor inicial para EDPs (Equação 2.4). A solução envolve uma sequência de EDPs definidas sobre subintervalos entre as datas de cupom, com condições de contorno que alternam entre o preço do título mais o próximo cupom e o valor de recuperação, dependendo se o valor da firma cobre a dívida.

Como uma fórmula analítica é considerada intratável para o caso geral ($i < N-1$), os autores empregam um esquema de diferença finita explícito para análise numérica. As principais etapas metodológicas incluem:

1. Transformação de Variáveis: Transformar a variável do valor da firma $V$ para $x = \ln V$ para lidar com o domínio $V > 0$ como $x \in \mathbb{R}$.
2. Discretização: Utilizar diferenças progressivas para derivadas temporais e diferenças centrais para derivadas espaciais. Para derivadas mistas ($\partial^2 B / \partial x \partial r$), os autores exploram esquemas de diferença central, progressiva e mista progressiva-regressiva, dependendo do sinal da correlação $\rho$.
3. Análise de Estabilidade: Uma derivação rigorosa das condições de estabilidade para o esquema explícito. Os autores demonstram que o esquema explícito padrão é geralmente instável devido às derivadas mistas. Eles derivam restrições específicas para o passo de tempo ($\Delta t$) e passos espaciais ($\Delta x, \Delta r$) em relação aos coefos de difusão e à correlação $\rho$ para garantir a estabilidade.

Principais Contribuições

1. Formulação do Modelo Unificado: O artigo formaliza um modelo unificado de dois fatores que integra o default estrutural (travessia de barreira em tempos discretos) e o default de forma reduzida (impulsionado por intensidade em qualquer momento) para títulos de cupons discretos fixos.
2. Condições de Estabilidade: A principal contribuição teórica é a derivação das condições de estabilidade necessárias para tornar viável o esquema de diferença finita explícita aplicado a este sistema de EDP de 3 variáveis com derivadas mistas.
   • Para $\rho = 0$, a estabilidade requer $\mu_x + \mu_r < 1$ e limites específicos nos termos de drift.
   • Para $0 < \rho < 1/2$, os autores mostram que a estabilidade requer $0 < \rho < 1/2$ e razões específicas entre os tamanhos dos passos espaciais ($\sqrt{\mu_r/\mu_x}$).
   • Para $-1 < \rho < 0$, restrições semelhantes são derivadas usando uma diferença mista progressiva-regressiva para a derivada mista.
3. Implementação Numérica: Os autores fornecem um arcabouço computacional para calcular preços de títulos, spreads de crédito e durações onde fórmulas analíticas não estão disponíveis.

Resultados
Utilizando o esquema explícito derivado com $\rho = 0$, os autores apresentam resultados numéricos baseados em um conjunto específico de parâmetros (ex: $N=2$ cupons, parâmetros Vasicek, intensidades de default específicas $\lambda_0, \lambda_1$). A análise produz:

• Preços de Títulos:
  • Os preços diminuem conforme a taxa de juros de curto prazo ($r$) aumenta.
  • Os preços aumentam conforme o valor da firma ($V$) aumenta, exibindo convexidade em relação a $V$.
  • Os preços exibem saltos (descontinuidades) nas datas de cupom, refletindo o pagamento de cupons e o reset das condições de risco de default.
• Spread de Crédito (CS):
  • O CS diminui conforme o valor da firma aumenta.
  • O CS aumenta rapidamente conforme o vencimento se aproxima, se o valor da firma for baixo.
  • O CS diminui conforme a taxa de juros de curto prazo aumenta (devido ao efeito de desconto sobre o componente livre de risco).
  • O CS aumenta com intensidades de default ($\lambda$) mais altas.
  • O CS exibe saltos nas datas de cupom; a magnitude do salto aumenta com o tamanho do cupom.
• Duração:
  • A duração aumenta com o valor da firma e a taxa de juros.
  • A duração exibe saltos após os pagamentos de cupons, consistente com a definição de duração como uma média ponderada do tempo até os fluxos de caixa.

Significância e Alegações
O artigo alega modestamente fornecer uma ferramenta numérica necessária para um modelo que generaliza o modelo estrutural de Agliardi em um arcabouço unificado para títulos de cupons discretos fixos. Os autores afirmam que obter uma fórmula analítica para este modelo unificado de dois fatores específico é "difícil" e provavelmente impossível no caso geral. Portanto, a significância do trabalho reside em:

1. Estabelecer as condições de estabilidade necessárias para tornar viável o esquema de diferença finita explícita para este sistema de EDP complexo.
2. Demonstrar que os resultados numéricos são compatíveis com situações financeiras reais, como evidenciado pelo comportamento intuitivo dos preços de títulos, spreads de crédito e durações sob diversas condições de mercado (taxas de juros, valor da firma, intensidade de default).

O trabalho serve como uma ponte entre a unificação teórica estrutural/de forma reduzida e a precificação prática para estruturas de títulos realistas (cupons discretos fixos) onde soluções de forma fechada não existem.