Bond and CDS Pricing via the Stochastic Recovery Black-Cox Model

Abstract: Building on recent work incorporating recovery risk into structural models by Cohen & Costanzino (2015), we consider the Black-Cox model with an added recovery risk driver. The recovery risk driver arises naturally in the context of imperfect information implicit in the structural framework. This leads to a two-factor structural model we call the Stochastic Recovery Black-Cox model, whereby the asset risk driver At defines the default trigger and the recovery risk driver Rt defines the amount recovered in the event of default. We then price zero-coupon bonds and credit default swaps under the Stochastic Recovery Black-Cox model. Finally, we compare our results with the classic Black-Cox model, give explicit expressions for the recovery risk premium in the Stochastic Recovery Black-Cox model, and detail how the introduction of separate but correlated risk drivers leads to a decoupling of the default and recovery risk premiums in the credit spread. We conclude this work by computing the effect of adding coupons that are paid continuously until default, and price perpetual (consol bonds) in our two-factor firm value model, extending calculations in the seminal paper by Leland (1994).


Background and Motivation
Most legacy credit models assume that recovery is a constant.It is well known, however, that recovery rates are not constant and indeed are correlated to a variety of risk drivers, including default and interest rates.For example, in their study on real-world recovery rates, Altman et al. [1] show that realized recovery rates are inversely proportional to realized default rates.This phenomena has been successfully incorporated into recent economic capital models (c.f.[2][3][4][5][6][7]) but still hasn't enjoyed mainstream adoption in pricing models.One notable exception has been in CDO pricing, where stochastic recovery has been added to many legacy models.This was necessary after the most recent credit-liquidity crisis, whereby for a period of time it was not possible to calibrate standard CDO models to the complete set of CDX.IG and ITRAXX.IG tranche quotes.The inability to calibrate the standard CDO models has been attributed to the assumption that recovery at default is deterministic, and does not depend on time or state variables.This has led to adding stochastic recovery to CDO models (c.f.[8][9][10][11][12][13]), which is now standard practice in industry.
However, stochastic recovery has not yet received equally widespread interest and acceptance for other credit products such as defaultable bonds, credit default swaps and credit linked notes.In fact, there are a dearth of pricing formulas for credit products where recovery is modeled explicitly.This has become commonplace even though recovery is clearly a key component in the determination of credit spreads [14].One recent approach to stochastic recovery modeling can be found in [15].Take for example bond pricing, where the price of a zero-coupon defaultable bond is given by the risk-neutral expected present value of the payoff Π τ (where τ is the default time): The first term describes receiving the full repayment of the notional N at maturity T in the event of no default (i.e., τ > T), while the second term describes receiving a recovery R τ in the event of default before maturity (i.e., τ ≤ T).The price of a zero-coupon defaultable bond is therefore given by where ∧ is the min operator a ∧ b := min{a, b}, P t (resp.E t ) is the risk-neutral probability of default (resp.risk-neutral expectation) conditioned on information about the default and recovery process known at t, and D(t, s) is the present value of a dollar at time t received at time s.To evaluate (2) one needs a model for D, default time τ, and recovery at default R τ .As noted in [15], it is common in structural models such as the Merton and Black-Cox models to make the simplifying assumption that a single risk factor A t representing the evolution of the underlying firm assets determines both the default time τ and recovery at default R τ .Symbolically, for a fixed f ∈ [0, 1] and a filtered probability space Ω, F A t t≥0 , F , P , where the natural filtration F A t is the minimal σ−algebra containing σ(A u ) for all u ∈ [0, t].Under this assumption, the price of a defaultable zero-coupon bond simplifies to which can be computed using only the single risk factor A t .
The lack of recovery modeling is even more explicit in the case of pricing Credit Default Swaps.For instance, the value of the default protection leg of a CDS per unit notional is given by which is the expected risk-neutral discounted recovery amount given default.However, it is common to make the simplifying assumption that recovery is a constant so that R τ ≡ R in (5) and the recovery term can be taken out of the expectation to obtain the simplified expression Under this simplifying assumption, which is standard in typical CDS pricing (c.f.[16][17][18] etc.), one simply needs a model for default τ and need not concern themselves with modeling recovery in the event of default.In this paper, we explicitly model recovery and remove the constant recovery assumption in CDS pricing thereby valuing the protection leg of the CDS using (5) directly rather than (6).
As described by the above examples of bond and CDS pricing, lack of recovery modeling inherent in classical structural models such as Merton and Black-Cox stems from the fact that default and recovery are driven by the same process.This feature makes it impossible to disentangle the effects of default and recovery.For instance, in the Merton model [19], the default time is defined as This default time has the benefit that recovery at default is random in the sense that it is defined through the random variable A T .However, a deficiency in the model is that default is only possible at the maturity T of the bond (i.e., τ Merton ∈ {T, ∞}).This is in direct violation of empirical evidence showing that firms default on their bonds at times before maturity as well.Black & Cox introduced their model [20] in order to remove this deficiency by allowing default to occur at times up to and including maturity (i.e., τ BC ∈ {[0, T], ∞}).The authors are able to achieve this by introducing an additional first passage default time τ K , defined as the first time the assets A t are less than the barrier K: with the usual convention that inf{∅} = +∞.The full Black-Cox default time τ BC is then However, the very same default mechanism in the Black-Cox model that allows for default before maturity forces recovery at default to be a constant.That is, (9) implies that for default before maturity, the asset value at default is equal to the predetermined constant barrier K (i.e., A τ BC = K).Hence Black-Cox improves the default modeling of Merton at the expense of constraining recovery to a predetermined known constant K when τ BC < T. 1 The main reason for this is that the Black-Cox model (as in all one-factor structural models) intertwines default and recovery risks.In one-factor structural models, the same process that determines default, namely A t , also determines recovery.This structure results in default and recovery rates being multiplicatively linked through the credit spread, which in turn makes the separation of default risk from recovery risk impossible.To disentangle the default risk from the recovery risk in a structural model, one needs to introduce a separate recovery risk driver.This leads us to the Stochastic Recovery Black-Cox model considered herein.
This paper is organized as follows.In Section 2, we review the classic Black-Cox (BC) model as a benchmark for our later results in Section 4 where we include recovery risk.In Section 3, we motivate the need for adding a recovery risk driver into the BC model and present our correlated asset-recovery model ( 49), which we compare with Moody's PD-LGD model presented in [4].In Section 4, we define the Stochastic Recovery Black-Cox (SRBC) model and use it to price instruments based on credit quality.Specifically, in Sections 4.1 and 4.4, we price bonds and CDS with the SRBC model.The SRBC setting is essentially a two-factor Black-Cox model where the extra factor is recovery risk.We then compare the prices obtained from the two-factor SRBC model with the original Black-Cox model.In Section 5, we provide an algorithm to compute the market implied recovery rate as well as analytical formulas for the recovery risk premium.Finally, in Section 6, we investigate the addition of coupon payments to bond prices in the stochastic recovery structural model.

Review of Credit Risk and Pricing in the Black-Cox Model
In this section, we briefly review bond and CDS pricing in the classical Black-Cox framework.In the original Black-Cox model [20], the default boundary K was taken to be exponential in time, i.e., K(t) = Ke −β(T−t) .However to simplify the exposition and results, we consider the case of a flat-boundary corresponding to β = 0, i.e., K(t) ≡ K, although the results can be trivially extended to the exponential boundary case as well.The results of this section will mainly serve as a benchmark for which to compare the results in Sections 4.1 and 4.4, where recovery risk is incorporated into Bond and CDS pricing via the Stochastic Recovery Black-Cox model described in Section 3.
We begin by listing the assumptions that form the Black-Cox model. 1 This constraint can be improved upon slightly by introducing a curved (time-dependent) boundary τ K := arg inf{t ∈ [0, T] : A t ≤ K(t)} for a suitably chosen K(t).However, the fact remains there is still a single driver for both recovery and default.

Assumption 1. (Constant Rates)
The short rate is assumed to be a constant r so that the discount factor is given by Assumption 2. (Asset Dynamics) The P t -dynamics of the firms asset value A t at time t are assumed to follow where r, σ A are constant and W A t is a standard Brownian process.
Assumption 3. (Default Time) Default is defined as where τ BC is explicitly given in (9).
Assumption 4. (Recovery Dynamics) Recovery at default is given by the asset value at default, The collection of Assumptions 1-4 form the Black-Cox model, and can be used to price defaultable bonds as well as other products such as CDS which we do in Sections 2.1 and 2.2 below.
Assumption 4 isn't often thought of as an explicit assumption of the Black-Cox model, since, as a one-factor model, this assumption is embedded in the default time.This is another example of the intertwining of default and recovery.Furthermore, notice in Assumption 4 that if τ BC < T then by definition A τ BC ≡ K so recovery given default is constant and deterministic.The main focus of this paper is to remove Assumption 4 and allow for stochastic recovery at default by decoupling the asset risk driver from the recovery risk driver.This should not be confused with adding stochasticity to the default boundary (c.f.[21,22]) since single-factor models with stochastic boundary still cannot disentangle the credit puzzle, as again, the asset value determines both default and recovery amount at default.
The Black-Cox model can be used to price credit products and compute related credit risk measures.For instance, in the case of a zero-coupon bond the price of such a defaultable zero-coupon bond is then given by the risk-neutral pricing formula Given the price of a defaultable zero-coupon bond via (14), we define the credit spread of the bond B t,T as the spread S t,T over the risk-free rate r which reprices the bond.Thus S t,T is defined implicitly by the solution to which returns the formula The risk-neutral probability of default PD t,T at time t is defined as the probability under the risk-neutral measure P that the default event τ occurs at or before maturity T, The risk-neutral expected loss given default LGD t,T at time t is defined as The term structure of expected Recovery Given Default R t,T can easily be inferred from the term structure of Loss Given Default (18) via the relation LGD t,T = 1 − R t,T .
The other credit product we consider is a Credit Default Swap (CDS).In a CDS, one party (the buyer) pays premiums to another party (the seller) to insure against default on a bond (c.f.[16][17][18]).Pricing consists of separately modeling the present value of the fixed premiums paid by the protection buyer, and the present value of the contingent default payment leg received by the buyer.The difference between the two is then the value of the CDS.If there is no upfront fee at initiation of the contract, then the premium P is given as the value that makes the contract worthless at initiation.
To be more precise, let T be the expiry of the CDS contract and let T n := {t = t 0 , t 1 , t 2 , ...t n = T} be the premium payment dates.For i = 1...n we define ∆t i = t i − t i−1 to be the time between payments.The premium leg of the transaction is then given by the expected present value of the premium payments P t,T that the buyer pays (and seller receives), namely where A p is the accrual payment in case default occurs between two payment dates.Instead of considering premiums paid at discrete dates, we pass to the continuous limit (see [23]) and consider the continuous premium formulation where the second term is the value of the accrual, i.e., Under the same continuous premium formulation, the value of the protection (default) leg can then be written as Using usual no arbitrage principles, the CDS premium P t,T is given as the value that balances these two equations, namely To evaluate (23), we need a model for recovery R t and default τ.A standard assumption [17] in a hazard rate framework is that recovery is a constant, under which our CDS premium (23) reduces to the classical result where R = R τ /N.In Sections 2.1 and 2.2 we present Bond and CDS prices, respectively, for the Black-Cox model which does not account for recovery risk.In Sections 4.1 and 4.4, respectively, we present bond and CDS prices under the Stochastic Recovery Black-Cox model which does take into account recovery risk.

Bond Pricing with the Black-Cox Model
In this section we price a zero coupon bond with default risk using the Black-Cox model.The structure of the price depends on whether the boundary K is larger or smaller than the notional N.

Definition 1. (Weak and Strong
Covenants in Black-Cox Model) Let K be the default barrier and N be the notional of the zero-coupon bond in the Black-Cox framework.We say the bond has a weak (resp.strong) covenant if K ≤ N (resp.K ≥ N).Definition 2. (Distances to Default) Define the following distance to default, which we will use in bond and CDS pricing under both weak and strong covenants. • Notice that the weak covenant distances to default d w 0 , d w 1 are the usual Merton distances to default (c.f.[15]).Lemma 1. (Probability of Hitting the Barrier) Let τ K be the first passage time (8).Then for any constant default barrier K > 0 and A where Φ is the cumulative normal distribution function and d s 0 , x s 0 are given by (26).
Proof.We begin this proof by noting that the value of a bond in the Black-Cox model requires a barrier at covenant value K.One can therefore calculate directly the probability of no-default using barrier option theory.To carry out this computation, we define the function W corresponding to a digital option: t From this definition (28) of W, if we add a lower barrier at K (See for example Ch.10 in [24]) then the value W of a digital option with a lower barrier at K is Therefore, by setting N = K in ( 25) and ( 29), we obtain the result and ( 27) above follows.

Note:
The result ( 27) can also be found in, for instance, Appendix B of [25].We provide another direct proof within the constructive proof of Theorem 3 below, by computing in the strong covenant case.
Theorem 1. (Bond Pricing under the Black-Cox Model) Suppose that Assumptions 1-4 are satisfied.Then the price of a defaultable zero-coupon bond ( 14) is given by: and the risk-neutral PD and LGD in the case of a weak covenant are and II.Strong Covenant Case.If K ≥ N the price of a zero-coupon bond is given by and the risk-neutral PD and LGD in the case of a strong covenant are and Proof.The proof uses the same integral techniques employed in Theorem 3 below.Using the barrier option characterization of the addition of a covenant, the price of a bond is expressed as This computation requires the joint density Pt [A T ∈ da, τ K > T], which is found in [26].We compute the associated integrals explicitly in our constructive proof of Theorem 3 for the case of stochastic recovery in the Black-Cox framework.
To carry out these integrals for the standard Black-Cox model, we set γ = ρ A,R = A 0 R 0 = 1 in the integrals constructed for the proof of Theorem 3. The difference between the weak and strong covenant cases (with or without stochastic recovery) is that in the strong covenant case, This fact is reflected in the appropriate substitution in the lower limits of the integrals against the joint density Pt [A T ∈ da, τ K > T].
Remark 1. (Special Case where K = N) Note that the special case where K = N is included in both the strong and weak covenant formulas in that

CDS Pricing with the Black-Cox Model
We now price CDS premiums using the Black-Cox model.In keeping with the continuous CDS pricing model, we extend the premium rate defined in [23] to include recovery at τ BC while protection is paid as long as τ K has not occurred.
In the strong covenant case, where K ≥ N, this reduces to the closed formula In the weak covenant case, where K ≥ N, this reduces to the closed formula

Note
• The numerator in (40) follows from direct substitution of the bond price calculated in Theorem 1 into the numerator in the general CDS formula (23).This direct substitution of the risky bond price also reflects the flexibility in assigning a weak or strong covenant, and will in fact be the only change observed when stochastic recovery is included in Section 5. • To complete the proof of Theorem 2, we will need to compute the denominator in (23).In the BC, and SRBC model forthcoming, this reduces to computing We are able to compute (43) by reducing to something more manageable via the fact that the discounted asset price is a local martingale.The reduced form follows from the optional sampling theorem: Proof.From Assumption 4, the general CDS premium (23) reduces to The main quantity to solve for now is E t e −r(τ BC −t) 1 {τ BC ≤T} , and upon it's computation and insertion into (45), we have the formula for the premium.By the structural definition of τ BC , it follows that and so In the case K ≥ N, the event {τ BC > T} reduces to the event {τ K > T}.Consequently, we utilize (44), (45), and the integral in (81) below to return the value (41).In the case that K ≤ N, similar calculations lead to (42).

Modeling Recovery Risk within a Structural Framework
The original Merton and Black-Cox structural models have been extended in several directions by many different authors, including stochastic interest rates, bankruptcy costs, taxes, debt subordination, strategic default, time dependent and stochastic default barrier, jumps in the asset value process, etc.While these extensions are by no means exhaustive (c.f.[16] for a more thorough discussion), they relax several of the main assumptions in the original models.However, none consider recovery as a risk factor, and assume the only risk driver is the asset value itself, or perhaps interest rates if stochastic rates are modeled.Empirical research, however, has shown that recoveries need not be constant in time and that typically the time-series of default probabilities and recoveries are inversely correlated [1,27].To make matters even more complicated, Hillebrand [28] shows that this correlation does not necessarily have to vary co-monotonously over the whole economic cycle.The correlation between PD and LGD has been investigated by several researchers in the context of credit capital.For instance, Giese [3] incorporates PG-LGD correlations into a single-factor Vasicek framework and finds that capital increases by up to 35% at the 99.9% confidence interval for high-yield credit portfolios.While investigating stressed LGDs, Miu and Ozdemir [6] find that in order to compensate for neglecting the PD-LGD correlation in credit capital modeling, the mean LGD must be increased by about 37% from its unbiased estimate in order to compensate for the lack of correlations.
We incorporate recovery risk into the Black-Cox model for three main reasons.First, as discussed in Section 2, the Black-Cox model has a constant recovery K in the event of default before redemption time.However, there is a large body of empirical evidence showing that recovery rates are inversely correlated with probabilities of default [1,14,27,[29][30][31][32][33][34] and so recovery risk is a real financial phenomena that should be modeled.Second, this correlation can have a very large effect on credit capital [5,6].Ignoring this effect in a pricing model could potentially lead to large mispricings or significant misestimation of risk.Finally, we show that including recovery risk in a structural framework is a mechanism which allows for larger spreads.This is important because empirical literature suggests that structural models tend to underestimate observed credit spreads by 10-15% on average (c.f.[35][36][37][38]).In [39], Gemmill argues that the extra observed spread is explained by other factors such as liquidity risk.However, just as in the Merton case [15], we show in Lemma 9 that adding recovery risk can lead to an increase in the credit spread over constant recovery models, and suggest that recovery risk is another possible mechanism to explain the additional observed spread.

The Correlated Asset-Recovery Model
Let A t denote the asset price at time t > 0 and let R t denote the recovery amount at time t > 0. The unobservable process R t is interpreted as the amount that would be recovered if default were to occur at t.The asset and recovery processes are modeled as two correlated geometric Brownian motions on (Ω, F , P) given by where (W A t , W R t ) are correlated, standard Brownian motions on our probability space.

Some Preliminary Results
Lemma 2. (Existence of Risk Neutral Measure) Let (A t , R t ) be the coupled measurable stochastic processes on (Ω, F , P) given by (48).Then there exists a risk-neutral measure P t such that (X t , F X t ) is martingale under P t , where X is the coupled process X t := e −rt (A t , R t ) .Furthermore, (A t , R t ) satisfy dA t = rA t dt + σ A A t dW A t dR t = rR t dt + σ R R t dW R t . (49) Proof.The proof follows from the results in [40] as the pair (A t , R t ) is a two-dimensional diffusion.
Lemma 3. (Solution to the PD-LGD Equations) Let (A t , R t ) be given by (48).Then, under the physical measure, (A t , R t ) is given by and under the risk-neutral measure, As an application of Lemma 3, we give a simple proof of Theorems 1-3 in [4].This result is the first comparison of recovery to asset values when the underlying asset is in state K. Then where Proof.We consider first the behavior of X under the physical measure, and prove (56).By Lemmas 2 and 3, and by (48), we can consider now the solution form under two standard Brownian motions W A and W on the probability space (Ω, F , P) where Using (58) the recovery process R t under the physical measure can be explicitly written as where γ and δ P are defined by ( 52) and (53) respectively.Then, using (59) we have To prove (57) we simply set µ A = µ R = r in (59) to recognize R under the risk-neutral measure as a shift from δ P to δ P, and (55) follows directly by dividing the closed form solution (57) by (56).
Under this framework, the price of a zero-coupon bond with face N is where default τ depends on A t and recovery upon default R τ depends on R t through (49), thus correlating default and recovery.
In Sections 4.1 and 4.4 we compute the bond and CDS prices (and related credit measures such as probability of default, loss-given-default and credit spreads) in the Stochastic Recovery Black-Cox default model where default time is defined via ( 9), but recovery is stochastic and correlated to asset value via (48).

Comparison of Recovery to Asset Upon Default
One of the features of the stochastic recovery model (48) we work with is that it is possible for the recovery value modeled, R t , to surpass the asset value A t at some time t ≤ T. This is due to the fact that A is the manager's estimated asset value instead of the actual asset value upon recovery.Upon default, it may come to pass that the actual value is higher than the estimated asset value A. Hence, it is possible that when default occurs, the debt holders are paid in full, and there is still some remaining capital.So, R t > N > A t is possible, for example.In fact, this is entirely possible when default occurs due to liquidity issues.In this context, A t is the default driving process, and R t the actual firm value.Hence, this model we work with proposes that there is less information about the asset than in the classic single-name models of Merton and Black and Cox.We also note here that in the end, it is the firm manager that decides firm default.Whether due to capital structure reasons or more pressingly due to day-to-day operational costs, it is the manager who has final say in this matter.This reinforces the fact that A is used to determine default, but R determines post-default recovery.We note that this allows for the possibility that credit spreads can be negative, at times.Given the correlated asset-recovery process (48), it is natural to estimate the probability that assets priced in recovery exceed the barrier K or notional N .The following Lemma provides such an estimate.Lemma 5.At the first passage time τ K , we have the estimates (64) Proof.By the Markov Inequality and Optional Sampling Theorem, it follows that after once again setting (wlog) t = 0, The numerator in (65) is computed using the same integral techniques employed in the proof of Theorem 3, except with a lower limit of K instead of N (i.e., strong covenant.)The second inequality, for the probability that recovered value is above notional, is proved in the same fashion: (66)

Connection between Recovery Risk and PD-LGD Correlation
In related work, Moody's KMV has recently proposed a PD-LGD Correlation model in the context of credit capital [4].While the structural model they propose is the same two-factor structural model in (49), the motivation for the model is different.The motivation for Moody's model is the empirically observed PD-LGD correlation presented in Altman et.al.[1], which is justified economically in [27,34,41] among other works.However, the study in [1] was conducted in the physical measure (realized-post-ante default rates) rather than in the risk neutral measure used in pricing, and the economic considerations are different.Nevertheless, in [4], the authors attempt to price a bond with this model by integration.Indeed, we are able price bonds and CDS's in this model via a martingale analysis, completing the analysis initially suggested in [4] by integrating against a joint density for process and first passage time.In particular, in Sections 4.1 and 4.4 below, we explicitly compute the bond price using stochastic calculus and the Optional Sampling Theorem, returning closed form solutions.

The Black-Cox Model with Recovery Risk
In this section we introduce the assumptions that define the Stochastic Recovery Black-Cox model (SRBC) and use it to price bonds and CDS in Sections 4.1 and 4.4 respectively.The SRBC model essentially relaxes the recovery assumption, Assumption 4, by replacing it with a weaker assumption on the dynamics of the recovery value, Assumption 6, allowing for randomness in recovery.In particular, the SRBC Model assumes the following: Assumption 5. (Constant Rates) The short rate is assumed to be a constant r so that the discount factor is given by D(t, s) = e −r(s−t) .
(67) Assumption 6. (Correlated Asset-Recovery Dynamics) The P t -dynamics of the firms asset value A t and recovery value R t at time t are assumed to follow where r, σ A , σ R , ρ A,R are constants and W A t , W R t are standard Brownian processes.

Assumption 7. (Default Time)
Default is given by the standard Black-Cox (BC) default time: where τ BC is also explicitly given in (9).
The collection of Assumptions 5-7 form the Stochastic Recovery Black-Cox model, and can be used to price defaultable bonds as well as other products such as CDS.which we do in Sections 4.1 and 4.4 below.Note that, just as in the original Black-Cox model, the constant interest-rate assumption can be relaxed to include time varying deterministic rates r = r(t) with little effort.

Bond Pricing with Recovery Risk
If default can occur before the maturity T, say if the bond issuer is forced into default if assets A ever fall below a default point K, then the bond price must reflect this extra possibility.Our main result in this section is the computation of a closed formula for such a price, where a recovery that is correlated to the asset is substituted at default.This incorporates the model first presented by the authors in [4].Theorem 3. (Bond Price under Stochastic Recovery Black Cox Model) Suppose Assumptions 5,6,7 hold.Then the general price of a zero-coupon bond and related risk metrics are: For sake of consistency, we point out that the SRBC bond prices reduce to the Stochastic Recovery Merton (SRM) bond prices computed in [15] as K → 0. To enable comparison with the BC model, we once again present the result for both weak and strong covenants: The risk-neutral PD and LGD in the case of a weak covenant are II.Strong Covenant Case.If K ≥ N the price of a zero-coupon bond is given by The risk-neutral PD and LGD in the case of a strong covenant are (74) Remark 2. We remark that the risk-adjusted SRBC distances-to-default d γ , x γ in (75) reduce to the standard distances-to-default (d 0 , x 0 ) and (d 1 , x 1 ) of the BC model if γ = 0 or γ = 1, respectively.This adjustment for gamma reflects the uncertainty of the firm manager in the partial information setting of what the recoverable value of the firm's assets truly are, and affects only the recovery term.It should be noted that the probability of default is the same as in the case of no stochastic recovery.It is only the Loss-Given-Default that is affected by adding R as a recovery driver.
The price for the zero-coupon bond under the SRBC setting is computed by employing the optional sampling theorem to the local martingale e −rt R t .We use τ := min {τ K , T}, a bounded stopping time adapted to the filtration generated by the joint process (A, R).The same technique is also used in the proof of the strong covenant case.
Proof.Without loss of generality, set t = 0.Under the risk-neutral measure, the corresponding bond price at issue is We are able to calculate these two integrals using a result found in [26]: Begin with a standard Brownian motion W on a probability space and define (77) Proof.See the proof of Theorem 2.1(i) in [26].
We employ the notation of Theorem 2.1(i) in [26] and focus on standard Brownian motions by setting From these definitions, it follows that Using this notation, we now compute the remaining integrals in (76): It follows from completing the square in the exponent of the normal density that Substitution of the results of (82) into (81) and leads to Similiarily, substituting γ = 0 into (81) results in Assembling the computations in (83) and (84), along with substitution of w, y in terms of K, N and A 0 , and using the form (76) leads to the closed-form solution for the weak covenant case.Note that we have w = ln N A 0 as our lower limit of integration in the weak covenant case.To switch to the strong covenant case, we substitute y = ln K A 0 for our lower limit and the result for bond price and associated spreads and loss given default follows.

Consistency and Reduction to Black-Cox Model
As a result of the closed form above for B SRBC t,T , the consistency of the model as T → t follows quickly for both weak and strong covenants: Lemma 6.If A t > N in the weak covenant case, or if A t > K in the strong covenant case, then Proof.Consider that for both weak and strong covenants, the explicit formulae for distances-to-default lead to lim Hence, lim T→t + B SRBC t,T (•) = N follows from the closed-formula given in Theorem 3, as does the reduction to the Black-Cox price when γ → 1 and R t → A t .

Greeks and Comparison with Standard Black-Cox model
Using our closed formula for bond price, we also compute the Greeks assuming stochastic recovery:

CDS Pricing with Recovery Risk
In this section, we price CDS using the SRBC model.We again consider both the weak and strong-covenant cases, and the results below depend on the bond prices in both covenant settings.Because of the general definition of a CDS premium under a structural model in Section 2, the premium here is the same as that computed for the BC model, with the exception that recovery is now defined by the correlated process R: Theorem 5. (CDS Premium under Stochastic Recovery Black-Cox Model) Suppose Assumptions 5-7 hold.Then the CDS premium (23) is given by: In the strong covenant case, where K ≥ N, this reduces to the closed formula In the weak covenant case, where K ≥ N, this reduces to the closed formula Proof.Changing to the SRBC model requires only a change to the corresponding bond price in (23), as the default trigger (time) remains τ BC .
Because of this closed form, we have consistency in our model as our parameters return to the standard Black-Cox model without stochastic recovery: Lemma 8. (Premiums Are Model-Consistent) The SRBC CDS premiums (88) are consistent with the BC CDS premiums (40), in that lim Proof.This result follows directly from the closed form for CDS premiums in Theorems 2 and 5, and the result lim in Lemma 6.

Implied Recovery Rates From Observed CDS Premia
We introduce notation in this section only for the function of A and t that defines the probability of default.Specifically, define for any α ∈ R : Here, κ α is the quantity defined in (87).Furthermore, define Using the results of Theorem

The Effect of Coupons
If coupons are paid at rate C per unit time, in addition to notional N, then (conditioned on the information (A t , R t ) = (A, R)) the total value of a coupon bond in our model is now derived from the expected present value of coupons paid until τ K ∧ T: By using moment-generating functions similar to those employed in calculating the premium-leg of our CDS, we compute the expected present value of coupon income to be An interesting limit in these calculation is the passage to an infinite horizon for redemption time.In the next section, we calculate the value of perpetual bonds and show consistency with the value (102) computed for coupon bonds under finite horizon.

Perpetual Bonds
For bondholders interested in long-term bonds, say of the 100-year variety [42,43], an approximation of this long-term bond as a perpetual bond with stochastic recovery is a useful measure of the risk involved with such an instrument.The addition of stochastic recovery, even when the firm-value level that triggers default is known, reflects the uncertainty of what the true value of a firm's assets may be worth upon default when bankruptcy occurs far off into the future.
We begin with some modeling assumptions: • Managers set a bankruptcy level K and accordingly we recall the definition (8) of τ K as the first passage time of assets A to level K. • Coupons are paid continuously at rate C, there is a risk-free interest rate r, and recovery at bankruptcy is R τ K .
Based on the above assumptions, and conditioned on the information (A 0 , R 0 ) = (A, R), the value of debt at time 0 (wlog) is (103) Lemma 10.For market parameters (γ, r, σ A ) such that the value of a perpetual bond under stochastic recovery is Remark 4. In the Leland model [44] of perpetual bond issuance under a structural model for firm value, the recovery upon default is modeled as for fractional bankruptcy cost ratio α.If we substitute R = (1 − α)A and γ = 1 into (105), then we retain the classical Leland [44] formula Proof.We begin by rewriting (103) as and orthogonal decomposition Ito-Calculus leads to, ∀t ≤ s, Substituting t = 0 and s = τ K returns If we define f K (t) as the density of our first passage time to level K, then we see that Assuming that at default, the bondholder receives instead the stochastically varying amount R τ K , it follows that, conditioned on (A 0 , R 0 ) = (A, R), (114) When rewritten, with r * := γr − 1 2 γ(1 − γ)σ 2 A , the bond value is Using well known results for moment-generating functions of first-passage times of Brownian motions with drift, assuming r * > 0, we obtain the closed form which results in the closed-form expression for a perpetual bond with stochastic recovery By utilizing (102), (116), and Theorem 3, we see that if in addition to r * > 0, we also have r > 1 2 − γ σ 2 A , then it follows that ∀t ≥ 0 lim and so we achieve consistency in the passage to infinite horizon.

Conclusions
In this work we introduced the Stochastic Recovery Black-Cox (SRBC) model which is essentially a Black-Cox model with an extra recovery risk driver.We then consider pricing with recovery risk and in particular explicitly compute closed form prices for both bonds and CDS under this framework, as well as associated risk metrics.This framework allows us to compute the Recovery implies by bond and CDS prices as well as compute the recovery risk premium.

Theorem 2 .
(CDS Premium under Black-Cox Model) Suppose that Assumptions 1-4 are satisfied.Then the CDS premium (23) is given by

Lemma 4 .
(Expected Recovery at Default) Consider the Asset-Recovery model (48) and define γ

Theorem 4 .
(Joint Density for Location and Maximum)

Lemma 7 .
(Greeks for SRBC).Let B SRBC t,T be the zero-coupon bond price from the Black-Cox model with stochastic recovery (70) under a weak covenant.Then the Recovery Greeks are given by + E 0 [PV 0,T [Cpns]] .(106)