Sovereign Default Model with Long-Term Debt

Computed with taste shocks (discrete choice methods) for $B'$ and $d$. Fortran parallelized with OpenMP.

Value and default option:

V\left(y, B\right) = \rho_D \log \left\\{ \exp\left[ \dfrac{V^d(y)}{\rho_D} \right] + \exp \left[ \dfrac{V^r(y, B)}{\rho_D} \right] \right\\}

$\Pr\left(d = 1 \middle| y, B\right) = \dfrac{ \exp\left[ \dfrac{V^d(y)}{\rho_D} \right] }{\exp\left[ \dfrac{V^d(y)}{\rho_D} \right] + \exp \left[ \dfrac{V^r(y, B)}{\rho_D} \right]} = \dfrac{1}{1 + \exp \left[ \dfrac{V^r(y, B) - V^d(y)}{\rho_D} \right]}$

Default value:

V^d\left(y\right) = u\left[h\left(y\right)\right] + \beta \mathbb{E}_{y'|y} \left\\{ \gamma V\left( y', 0 \right) + (1-\gamma) V^d\left(y'\right) \right\\}

Repayment values:

$W\left(y, B, B'\right) = u\left[ y - \kappa B + q\left(y, B'\right) \left( B' - (1-\delta) B \right) \right] + \beta \mathbb{E}_{y'|y} V\left(y', B'\right)$

$V^r\left(y, B\right) = \rho_B \log \sum_{B'} \exp \left[ \dfrac{W\left(y, B, B'\right)}{\rho_B} \right]$

Choice probabilities for $B'$:

$\Pr\left(B' = x \middle| y, B \right) = \dfrac{ \exp \left[ \dfrac{W\left(y, B, x\right)}{\rho_B} \right] }{\sum_{i} \exp \left[ \dfrac{W\left(y, B, i\right)}{\rho_B} \right]}$

Bond price schedule:

$q\left(y, B'\right) = \dfrac{1}{1+r} \mathbb{E}_{y'|y} \Pr\left(d=0 \middle| y', B' \right) \left[ \kappa + (1-\delta) \mathcal{Q}(y', B') \right]$

with

$\mathcal{Q}(y', B') = \sum_{B''} \Pr\left( B'' \middle| y', B' \right) q\left(y', B''\right)$

Functional forms and shocks:

$u\left(c\right) = \dfrac{c^{1-\sigma} - 1}{1-\sigma}$

$h\left(y\right) = y - \max\\{ 0, \lambda_0 y + \lambda_1 y^2 \\}$

$\log y' = - (1-\rho) \dfrac{\sigma_y^2}{2 \left( 1 - \rho^2\right)} + \rho \log y + \sigma_y \varepsilon, \quad \varepsilon \sim \mathcal{N}(0, 1)$