1. Background and Motivation
In his seminal paper,
Wang (
2002) proposed a transform method to price both liabilities and contingent claims, whether traded or not. One application of the Wang Transform is in option pricing, which for example can be used in the structural modeling of defaultable bonds. In the practice established by actuaries, the time value of money is another important aspect of bond pricing, and in the same issue that Wang’s paper appeared,
Bühlmann (
2002) proposed a paradigm shift in thinking about this technology. Indeed, in his editorial, Bühlmann advocated that actuaries take a
numeraire based approach to financial risk.
This approach to pricing lines up with that of quantitative finance, where
Merton (
1974) has laid the groundwork for structural models of corporate bond pricing that provide internally consistent methods to price credit-linked instruments. One major benefit of Merton’s model is in its elegant simplicity, where he assumes that the default event, as well as post-default recoverable values, are driven by a common stochastic factor: the firm’s asset value. Merton defines equity as a call option on the firm value, and so a change of numeraire also appears in his solution to the bond price via the Miller–Modigliani framework of capital structure. However, a single factor for pricing such bonds assumes that bond and equity returns are perfectly correlated, which need not be true. This also leads to perfect correlation between probability of default and recovery-given-default, which is also not empirically observed. Adding a coupled stochastic recovery driver addresses these issues.
1.1. Recent Extensions
New methods, developed in
Cohen and Costanzino (
2015,
2017), address the entanglement of default risk (i.e., the risk borne by the holder of the financial instrument that the counterparty will not fulfill its obligations to pay the contract holder) with recovery risk (i.e., the risk of not fully recovering the principal in the event of default) in the pricing of defaultable bonds. The authors show that this mixture can lead to a mis-estimation of the total risk that investors in these bonds have undertaken. We also note here that recovery risk has not received as much attention as default risk, and empirical studies have shown that investors are taking on a significant amount of recovery risk for which they are not being properly compensated (see, for example,
Schläfer and Uhrig-Homburg (
2014) and references within.)
As we will see, the effect of this extra uncertainty about post-default recovery can be observed through the distortion of the original distance to default in the Merton model, and a distortion of the probability of default in the general structural model of default, when calculating the recoverable value of a bond post-default. In our short note, we develop a transform approach for a multifactor bond-pricing model that represents the partial information available to firm managers about recoverable value. In the spirit of
Bühlmann (
2002), this approach is dependent on the choice of numeraire as well as the conditions for default. We calculate this probability distortion for the two-factor stochastic recovery Merton model to be a version of the Wang Transform, and present multiple examples to illuminate the extended transform method, including models that incorporate jumps in the asset-recovery values and enforced bound on recovered value. This uncertainty also factors into the calculation of a term structure of expected recovery, which we develop below.
Finally, we establish here that our approach is firmly in the structural model framework pioneered by
Merton (
1974) , and extended by others such as
Black and Cox (
1976) and
Leland (
1994) to name just a few. The reader is directed to
Cohen and Costanzino (
2015,
2017) for a more robust reference list on the structural approach, as well as the reduced form approach, to credit risk modeling and their differences in application to data.
1.2. Paper Outline
The rest of the paper is developed as follows: in
Section 2, we set upon the modeling framework. In
Section 3, we revisit the bond pricing techniques presented in
Cohen and Costanzino (
2015,
2017) and begin their extension via the definition of a Partial Information Transform of the risk-neutral, dollar denominated probability of default. With this definition, we can recast the results in
Cohen and Costanzino (
2015,
2017) via this technique.
Section 4 sees the extension of the models presented in
Cohen and Costanzino (
2015,
2017) to ones that also incorporate bounded recovery in the the pricing of the bond. In
Section 5, we allow for discontinuous information about the firm and so non-predictable defaults are incorporated into our approach to stochastic recovery. Results on short-term credit spreads are also developed, with alternate derivations of these results presented in the
Appendix A.
2. Modeling Framework
In order to incorporate recovery risk into classical structural models and understand its effect on the credit spreads, a coupled stochastic recovery risk driver
R was added in
Cohen and Costanzino (
2015,
2017) to the one-dimensional structural models of Merton and Black–Cox that model firm value (assets)
A. To enable this extension, the asset process
A in these classical models is extended to the joint process of
. This joint process lives in a filtered probability space
satisfying the usual conditions, where
is a risk-neutral measure. In this setting, the joint dynamics of
under
are given by
where
is the firm’s asset value at time t.
is the firm’s recoverable value at time t, if default were to happen at t.
is the rate of return of firm value at time t.
is the rate of return of recoverable value at time t.
is the filtration generated by A, and
is the filtration generated by .
Under the risk-neutral measure, we observe that
is proportional to the non-diversifiable portion
of recovery risk defined via
Remark 1. The constraint of continuous evolution (1) is relaxed in Section 5, where a jump-to-0 asset-recovery model is utilized. For a general reference on such processes in credit risk, please consult Zhou (1997). In classical structural models of credit risk, recovery is indeed stochastic by the modeling choice
. However, this assumes that there is only one driver
A for credit risk and so (under the Miller–Modigliani framework assumed in
Merton (
1974) and subsequent work), equity and bond prices are constrained to be perfectly correlated. We make two further assumptions in the modeling of recovery and firm assets:
Assumption 1. We assume a constant risk-free rate r, although our results can be generalized to time-varying rate .
Assumption 2. Under , we assume that for all , almost surely for risk-free rate r. This allows for us to develop discounted assets and discounted recovery as local martingales.
Remark 2. This second assumption is further explored in Cohen and Costanzino (2015), where equity is modeled as a call option on recoverable assets R, rather than pre-default firm value A. The reader is also directed to the excellent resource Korn and Kraft (2004) and references within for further treatment of this assumption of proxies on assets. 2.1. Classical Bond Pricing
Consider now the case of a zero-coupon bond issued by a firm with notional
N and redemption time
T. In the classical one-dimensional setting of firm value
A and information
at time
T, default is modeled as a stopping time
. We can also study default and its role in pricing via the set
. In this setting, the price
of a bond is the exchange option on underlying asset
A for notional
N:
where stochastic recovery is simply the firm’s assets upon default (NR).
Remark 3. Although this set D is most commonly generated in the literature as the set of a default time τ, defining default via a set allows for a slightly more general approach to default modeling.
3. Bond Pricing in Stochastic Recovery Models
Under the stochastic recovery framework (SR), specifically for Merton and Black–Cox models but extended via the optional sampling theorem, the price of a defaultable bond is shown in
Cohen and Costanzino (
2015,
2017) to be similar in nature to the traditional structural models without stochastic recovery, where
A remains the sole trigger for default:
:
We establish here that Equation (
4) is true in a general sense when discounted recovery is a local martingale, where the optional sampling theorem allows us to calculate for
and so
Depending on the default set
D, computation of the bond price in Equation (
4) requires the solution of a two-factor PDE with initial and perhaps boundary conditions. See, for example,
Cohen and Costanzino (
2015,
2017) for futher details on this approach. However, with the strong influence of probabilistic methods in solving such PDEs, it is reasonable to consider parallel solution methods. One such approach is seen in the classical one factor models presented in
Ishizaka and Takaoka (
2003). By appealing to barrier option theory, the authors in
Ishizaka and Takaoka (
2003) utilize the measure
with the Radon–Nikodym density process defined by the firm value
A
with the notational conventions of
that the bond price without stochastic recovery can be written as
This line of reasoning, that of changing numeraire, can be carried over into the case where
recovery is used as the numeraire, to build a new measure
to estimate probability of default, with Radon–Nikodym density process and conditional probability defined by
Armed with this measure, we use recovery
R to estimate the probability of default triggered by asset value
A. The authors in
Cohen and Costanzino (
2015,
2017) present closed form solutions in the Merton and Black–Cox cases, which can now be shown to be equal to
Using the idea of changing the base-numeraire to the two-factor model recovery
R, from the classical one-factor recovery model numeraire
A, we can look to solve the two-factor bond price (
11) by solving for the quantity
:
Definition 1. We define the Partial Information Transform (PIT) of for a model with assets A and default set , to be the transform which maps .
For the two-factor model, with Asset–Recovery as a pair of correlated Geometric Brownian Motions, we can compute the
PIT via the following Lemma:
Lemma 1. For classical structural models that are driven by our Asset-Recovery pair (1), with constant volatilities , the PIT has the form Proof. Under Assumptions 1 and 2, Ito’s formula applied to (
1) results in
Substitution into our pricing formula (
11) leads to
☐
Before using the above for classical structural models, a few points to consider:
Remark 4. Market Data and Classical Limits:
Remark 5. Partial Information and the Term Structure of Recovery Rates:
We point out here that the distortion also finds its way into the decoupling of recovery at default and probability of default. The result is that the recovery rate is the product of recovery modeled at t (per discounted notional) and the ratio of distorted to risk-neutral probability of default. Indeed, we have This decomposition also highlights the effect of distortion of the numeraire-based probability of default on the term structure of recovery rates, and allows for greater flexibility in calibrating to empirical data. This form is a generalization of the one-factor approach to implied recoveries undertaken in the important work of Das and Hanouna (2009). The addition of a term structure of recovery rates addresses the standard assumption of industrial models that pin expected recovery at a constant , only to find that subsequent calculations of quantities such as spreads are far from observed values. For more on the term structure of expected recovery rates, please refer to Doshi et al. (2015) and references within. Two examples of the Partial Information Transform are presented in this section, where volatilites
in (
1) are constant. In the first case, that of the Merton Model, we recover a version of the Wang Transform as a special case of the Partial Information Transform.
3.1. Merton Model with Stochastic Recovery (SRM)
In the classical Merton case, the probability of default is known to be
where
and distances-to-default
are defined via
This leads to the Partial Information Transform for the Merton model:
Proof. Using the result
, shown explicitly in
Cohen and Costanzino (
2015), we can compute the asset and recovery numeraire defined probability of default as
. Furthermore,
We remark that the distance-to-default under has shifted to a new value under . As mentioned above, this new metric requires estimation of only one new parameter, that of the in our recovery R compared to asset A. From this interpretation, it follows that there are multiple definitions of distance-to-default, depending on the numeraire used. ☐
3.2. Black–Cox Model with Stochastic Recovery (SRBC)
A similar effect is noted when a level covenant
on firm assets is introduced, such as in the Black–Cox case, and the probability of default in the Stochastic Recovery Black–Cox case via both asset as well as recovery numeraire are
Cohen and Costanzino (
2017):
where it is assumed
and we define
In adding a level barrier K, we are removed from the class of Wang Transforms on standard normal probabilities, but still observe a deformation of the default probability as we pass from firm assets A to recoverable value R as our choice for numeraire.
4. Bounded Stochastic Recovery
In the previous section, we revisited the models developed in
Cohen and Costanzino (
2015,
2017), where recovery may indeed surpass notional. As expressed in those papers, this approach trades bounded recovery for more tractable calculations, with estimates of the probability of the event that recovery exceeds notional. In this section, we present two models that extend the previous approach, while enforcing bounded recovery and, in the second case, an inclusion of stochastic volatility. We should point out here that bounded recovery is part of the original Merton default model, and, in
Cohen and Costanzino (
2015,
2017), the recovery is potentially unbounded, but can be rectified by calibrating the model to observed (positive) spreads.
4.1. Merton Model with Stochastic Recovery and Enforced Positive Credit Spread (SRM-CR)
One of the features of the partial information model proposed in
Cohen and Costanzino (
2015,
2017) is the non-zero probability attached to the event
In this setting, with
, there is correspondingly a chance that, in default, the recovery is greater than notional and spread can be negative. The size of the negative spread depends on model parameters, and the calibrated model will only produce a negative credit spread under very unique circumstances. While negative credit spreads have been observed in the market (c.f.
Bhanot and Guo 2011), one can eliminate this possibility with a slight modification of the model. Consider instead the modification of Equation (
4) via
It follows that
and so the corresponding credit spread is non-negative. With a little algebra, we can show the relationship
In the Merton default setting, with
, we have
Hence, the partial information approach in Equation (
4) that uses
R for recovery, under Merton default timing, also allows for robust modeling of credit spreads. In this section, we cap stochastic recovery in the Merton model by notional to restrict to non-negative credit spreads. We present below a closed-formula for both bond and (risk-neutral) credit-spread:
Theorem 1. The two-factor structural model (1) results in bond price and credit spread as follows:where andis a bivariate normal probability, and we define the distances-to-default for asset and recovery asand Proof. The structural approach to modeling leads to the initial formulation
for bond price. By definition (
30), and via the fact that
, direct computation (also carried out in full detail in
Cohen and Costanzino (
2015)) leads to the formula for bond price
Correspondingly, we compute the credit spread using its definition
☐
Remark 6. This extended Stochastic Recovery Merton model that enforces positive credit spreads by strictly capping the payout to bondholders to notional N is not an example of the Partial Information Transform approach, as the triggering event now depends on the joint value of , and not just . It is tempting, of course, to expand the definition of the default event D to a trigger event for the joint process asand so In this set-up, the bondholder is paid N in full if either firm assets at T are valued to be at least N, or if firm assets are viewed to be less than N, but recovery is valued to be at least N. It follows that 4.2. A Bounded Discounted Recovery Model with Stochastic Volatility
Consider a recovery R that lives on our probability space, with the same recovery-driving Brownian motion (with same correlation to ) as in the previous section, but now the recovery process has a volatility that is driven by R as well: This recovery process (
36) can be rewritten in terms of another process
:
where
and so
is a local martingale with values residing in the interval
almost surely, and so is in fact a martingale. We point out that this process
π has a deep history in Mathematical Statistics (
Shiryaev 2007) and Actuarial Science (
Gerber 1997), to name but two important historical resources, and more recent applications of bounded stochastic processes have found their way into election modeling (
Taleb 2017) and financial instrument modeling (
Carr 2017).
In this setting, the bond price reduces to
However, to find a closed form for this transformation , say under a Black–Cox default time with barrier K, would require the solution of a PDE with barrier and initial conditions. A middle path would be to restrict the original Geometric Brownian Motion approach to Recovery to be capped at notional N when used in pricing, which we presented in the previous subsection.
4.3. Further Extensions
It should be noted that the examples presented above can be generalized to add extra conditions on the default. Computation of the PIT could involve the joint distribution of
with quantities such as the running maximum, minimum, and average of
A over the interval
, among other quantities, depending on the new default set
D. For example, one could consider the default time corresponding to the cumulative Parisian option, as detailed in
Ishizaka and Takaoka (
2003) and references within.
5. Stochastic Recovery Merton Model with Jumps
In this last example, we extend the two-factor model in Equation (
1) to allow for
abnormal price variations. These variations are different than those that
arise from new information that causes marginal changes in the stock’s value (
Merton 1976). Such a model was proposed in
Samuelson (
1973) and used to solve for call option prices in
Merton (
1976).
In the
Samuelson (
1973) model, the underlying is allowed to suddenly jump to zero value, due to information that relates to the firm specifically or perhaps the larger industry to which it belongs. The jump-to-zero model is reasonable for equity prices as typically equity-holders get very little in the event of default. However, a jump-to-zero model for credit is less intuitive. Indeed, in a one-factor model where the assets
drive both default and recovery, a jump-to-zero model is extremely restrictive since it implies that we have zero recovery at default. In a two-factor model, such as Equation (
1) whereby there is a separate risk driver for assets and for recovery, a jump-to-zero model for assets is reasonable since the assets define the default risk and not recovery risk. Therefore, even though
may jump-to-zero indicating default, the amount recovered at default is driven by
, which is not necessarily zero. This is an added benefit of having an additional recovery risk driver.
As we shall see below, the Partial Information Transform in this model is computed directly via integration, and is a modified version of Equation (
12) that incorporates jumps as well as diffusion in modeling
A. 5.1. The Jump-to-Zero Diffusion Model
The jump-to-zero diffusion model we define allows for a jump to zero in the asset value , and therefore unpredictable default, but not necessarily a jump-to-zero in the recovery value . Unlike in the one-factor model, the decoupling of the default risk driver from the recovery risk driver allows us to treat jump-to-zero of the asset value as a trigger for default, which does not affect the recovery at default.
To define the model, suppose that the asset and its recovery follow the jump-diffusion, modified from Equation (
1):
where
and
are
independent driving Poisson processes with constant rates
and
, respectively. The process
models an unpredictable default event with rate
whereby the recovery value upon default is non-zero, whereas the process
models a (catastrophic) unpredictable default event with rate
whereby both the asset and recovery value upon default are exactly zero.
Consider a filtration
generated by the process
, and the larger filtration
generated by the coupled jump-diffusion in Equation (
40). This evolution is similar to the one proposed in
Martin (
2007). We also note here that firm accounted asset’s jump to 0 may or not signal that recoverable assets jump to 0.
In this setting, the solution to Equation (
40), conditioned on information at time 0, is
Here,
is the solution pair of Equation (
1), as before.
5.2. Effect of Jumps on Zero-Recovery Bonds
If we return to the simplest case where
, then we can see for the zero-recovery, zero-coupon bond
which has the solution (
Merton 1976)
This extra jump due to unexpected, abnormal information, returns a significant modeling advantage, as it adds to the no-jump asset model credit spread:
Thus, in the short-term as
, formula (
44) returns a short-term credit spread of
.
5.3. Effect of Jumps on Bond Prices with Stochastic Recovery
If we allow for stochastic recovery in this setting, then
Notice that default is still only driven by
, which now allows for the surprise jump-to-0 value in asset. The correlated jump-diffusion in Equation (
40) is still a martingale, in that
and so, using Equation (
40), we can rewrite Equation (
45) as
Closed form solutions for the components in the last line of Equation (
47) are retained from Lemma 2, but adapted with the parameter change of
and corresponding distances-to-default
and
under
denoted by
and
, respectively. We use these to write the bond price with stochastic recovery and default-to-0 as
Here,
is the zero-coupon bond price under stochastic recovery, no jumps in asset
or recovery value, and the new risk-free rate
, which is also used to calculate the new distance to default
under
instead of
r.
5.4. Pricing via Change of Numeraire
From the calculations above in
Section 5.3, we have the following result:
Proof. Once again, consider that bond pricing will occur via change of measure. Using the form (
40) as well as the results in Equations (
43) and (
49) above, we compute that, for the Stochastic Recovery Merton Model with short-term credit spread built in, the pair of numeraire measures are
where an adjusted risk-free rate
is used to calculate
☐
Remark 7. In analyzing these results, we consider the following points:
The total default probability transformation upon change of numeraire now results in a weighted average of
- 1.
a Wang Transform of the classical Merton default probability , and
- 2.
a probability of default based on a jump-to-0, all under a new distortion of the risk-free rate from r to .
This is due to the nature of default now being extended to allow for external triggers to force assets to 0.
Note that the recovery value when , and so recovery, under default triggered by , does not factor into our bond price. If we distort , then does indeed factor into our calculations. Please see the subsections below for more details, and the effect of jumps in this model on short-term credit spreads.
To conclude this section, we present an analysis of the short-term credit risk in this Merton Model that we have extended to allow for the effect of jumps on the credit spread. We begin by looking at zero-coupon, zero-recovery bonds.
5.5. Widening the Credit Spread
The effect of jumps in the asset value directly carries over into the short-term credit spread:
Theorem 2. The effect upon short-term credit spreads of adding jumps into both assets and recovery of a firm, in a classical asset model used for structural modeling such as Merton, is reflected in the relationship Proof. Using the result in Equation (
49), we repeat the calculation in Equation (
44) to compute our credit spread as
If we let
, holding
fixed and assuming
is also fixed and known, we denote by
the distance to default with only
varying. By L’Hopital’s rule,
To complete this limit calculation, we compute that
and, since
we have the limit
A similar argument is used to show that
and thus our main result:
☐
Remark 8. The structure of short-term credit spreads, obtained in this result, suggests another reason for incorporating recovery into structural modeling: to allow for short-term credit spreads that are stochastic in nature. A related viewpoint is that one could start from observed stochastic credit spreads to model recovery in bond pricing. Please see (A8) for the comparison with classical Merton (with jump-to-0) credit spreads. Finally, we offer that the effect of short term credit spreads can also be seen if jump-to-0 is combined with bounded recovery in the pricing of the bond. These calculations follow the same format as those carried out in the sections above, and the interested reader is encouraged to compute the relevant quantities for their own needs. 6. Conclusions
By incorporating recovery risk into pricing of defaultable bonds, we have shown that the effect mimics default risk by the distortion of default probability estimation. This distortion reflects the extra risk investors take on when bonds recover a non-constant amount, or non-constant fraction of assets, upon default. This technology can be further incorporated into credit default swaps
Cohen and Costanzino (
2015,
2017) and other financial instruments, as well as used to perhaps better calibrate risk measures arising from observed data. From an industry standpoint, the extra flexibility of allowing for stochastic recovery is also reflected in the ability to model a term structure for recovery rates by understanding the effects the model choice for recovery
R has on the distortion probability
that factors directly into these calculations. Furthermore, by allowing jumps to 0 in the asset as a trigger for default, we are able to provide a hybrid model that incorporates stochastic recovery as well as non-zero short term credit spread that could possibly be implied from observed market values.