Free boundary problem pricing defaultable corporate bonds with multiple credit rating migration risk and stochastic interest rate

: In this paper, valuation for a defaultable corporate bond subject to multiple credit rating migration risk and stochastic volatility of interest rate is addressed in the structure framework through a free boundary problem, which is derived by a series of transformations. The existence, uniqueness and regularity of solution to the free boundary problem are obtained to verify the rationality of the bond pricing model. Furthermore, we show that the solution of the free boundary problem is convergent to a close form steady status, which may provide some information on the developing characteristics of the bond price. As the coexistence of stochastic interest rate and defaulting boundary, this convergence is achieved through an auxiliary free boundary problem and a Lyapunov argument. Interestingly, the converged steady status can be explicitly solved, which is not the case in the existing literatures on multiple credit rating migration. Finally, we present an explicit formula for valuating this defaultable bond with multiple credit rating migration risk and stochastic interest rate.


Introduction
The globalization of financial markets has been developing rapidly, which requires more on credit risk management. Credit risk, which is attracting more and more attention from people in academics and practices, refers to not only default risk but also credit rating migration risk. The credit rating migration risk is playing a more and more significant role in financial markets and risk management, especially after the outbreak and spread of 2008 financial crisis. In particular, the credit rating migration risk makes difference to the corporate bond pricing.
There are two traditional models for default risk, involving the structural models and reduced form models. In the reduced form models, the default event is depicted and captured by introducing an exogenous variant. The default time is modeled by a stochastic default intensity in this approach, see [7,14,17] and so forth. The assumption on the structural models is that if the below bound of some insolvency threshold is met by corporate value, default occurs. In the model proposed by Merton [28], a default event may only occur at the maturity. Subsequently, Black and Cox [2] extended Merton's model to a first-passage-time model, where default may occur at any time before the debt maturity, see also [1,3,18,25] and so forth.
With regard to literatures working on credit rating migration, a commonly adopted approach is the Markov chain, which is captured by transition intensity matrix coming from general statistic data, see [5,8,9] and so forth. The framework of reduced form then can be directly developed for dynamic credit rating migration process, see [6,15,30] and so forth. However, the Markov chain ignores the role played by the corporate value when modelling credit rating migration. In fact, the corporate value is an important factor in the credit rating migration and should be taken into consideration. Accordingly, from the corporate perspective, Liang et al. [20] started to model and analyze credit rating migration risk by structural model based on Merton's model. They set a predetermined migration threshold to divide the corporate value into high and low rating regions, where the corporate value follows different stochastic processes. However in practice, the threshold dividing credit ratings is usually not predetermined. To solve this problem, Hu et al. [11] improved the model proposed by Liang et al. [20]. They determined the migration boundary by the dynamic proportion between corporate debt and corporate value, which results in a free boundary problem. Subsequently, Liang et al. [21] incorporated a risk discount factor, which measures the sensibility of credit rating migration to the proportion, into the model and showed that an asymptotic traveling wave solution exists in the free boundary problem. Problem on credit rating migration in switching macro regions can also be referred to Wu and Liang [36], while credit contingent interest rate swap with credit rating migration can be seen in Liang and Zou [23].
In particular, the aforementioned works [11,20,21,23,36] only take two credit ratings into consideration in credit rating migration problem. The credit region is divided into the high rating region and the low rating region, which results in only one free boundary in the corresponding free boundary problem. However, in practice, we should notice that there are usually more than two credit ratings used when accessing the corporate credit levels. The Standard & Poor's, an international rating agency, downgraded the long-term sovereign credit rating of Greece from A-to BBB+ on the evening of December 16, 2009. This verifies the fact and inspires us to consider multiple credit ratings in migration problems. Wu and Liang [34] provided some numeric results for multiple credit rating migration problem. Wang et al. [32] presented some theoretical results by showing that the asymptotic traveling wave solution obtained in Liang et al. [21] persists in the free boundary problem with multiple free boundaries. By considering stochastic interest rate in reality [19,24,29], Yin et al. [37] improved the model by replacing the constant interest rate with a stochastic version. This improved model covers the previous works where only two credit ratings are involved with an interest-dependent volatility [22] or multiple credit ratings migration with constant interest rate [32,34]. Then Huang et al. [12] continued to study the bond pricing model with multiple credit rating migration and stochastic interest rate. They contributed to establishing the asymptotic traveling wave solution in the time-heterogeneous free boundary problem with multiple free boundaries.
The aforementioned models for credit rating migration in the structure framework are based on the Merton's model, i.e., default may only occur at the maturity. However, in practice, default may occur at any time up to maturity [4,26,31]. Wu et al. [35] relaxed the default restriction in the credit rating migration model by setting a predetermined threshold capturing the first-passage time when the default occurs. Once the corporate value falls below the threshold at any time, the default occurs. This results in a free boundary problem subject to a new boundary condition. Again we have to notice that the credit region in their credit rating migration model pricing a defaultable corporate bond is still divided into the high and low rating regions and meanwhile, a constant interest rate is also considered. Hence, motivated by these existing works on the effect of credit rating migration risk when valuating a corporate bond, in this paper, we devote to studying a pricing model for a defaultable corporate bond with both multiple credit rating migration risk and stochastic interest rate. Our work extends the existing works [12,37], where the multiplicity of credit rating and stochasticity of interest rate are involved in their models, by inserting the default risk. Meanwhile, we improve the results of Wu et al. [35], who considered default risk in pricing the corporate bond with only one credit rating migration boundary and constant interest rate, to fit the effects of multiple credit ratings and stochastic interest rate.
The difficulties in analysis are generated from the joint effects and mutual restrictions among the multiplicity of credit ratings, stochastic volatility of interest rate and presence of default boundary. Besides that the free boundary problem turns into a initial-boundary problem, in particular, it is perplexed by not only a time-dependent and discontinuous coefficient but also a time-dependent process arisen in the problem. In addition to an irreducible barrier boundary, these indeed cause some troubles in deriving necessary estimates and then proving the existence and uniqueness of solution to the free boundary problem. Another contribution is the asymptotic behavior of solution to the free boundary problem. We prove that the solution converges to some steady status, which is the spatially homogeneous solution of an auxiliary free boundary problem, whose coefficients are the long time limits of the time-dependent coefficients in the original free boundary problem. This convergence, which shows us the developing tendency of solution, is established by two steps. In the first step, it is shown that the solution of the original free boundary problem converges to the solution of the auxiliary free boundary problem with time tending to infinity, while in the second step, it is shown that the solution of the auxiliary free boundary problem converges to the steady status. Moreover, the steady status can be solved explicitly. Thus, we present an explicit formula to valuate the defaultable corporate bond with multiple credit rating migration risk and stochastic volatility of interest rate.
The paper is organized as follows. In Section 2, the pricing model is constructed. In Section 3, the model is reduced into a free boundary problem with initial condition and boundary conditions involving default boundary and migrating boundaries. In Section 4, an approximated free boundary problem is analyzed and some preliminary lemmas for uniform estimates are collected. In Section 5, through the approximated free boundary problem, the existence and uniqueness of solution to the free boundary problem are obtained. In Section 6, it is proved that the solution is convergent to a steady status by a Lyapunov argument. Then we conclude the paper by presenting an explicit pricing formula for the defaultable corporate bond in Section 7.

The baseline pricing model
In this section, we set up the baseline pricing model for a defaultable corporate bond subject to multiple credit rating migration risk and stochastic interest rate. Some necessary assumptions are put forward as follows.

Assumptions
Let (Ω, F , P) be a complete probability space. Suppose that the corporation issues a defaultable bond, which is a contingent claim of its value on the space (Ω, F , P).
Let S denote the corporate value in the risk neutral world. It satisfies the Black-Scholes model where r is the time-varying interest rate and σ is the heterogeneous volatility with respect to credit ratings. W t is the Brownian motion generating the filtration {F t }. The credit is divided into n ratings. In different credit ratings, the corporation shows different volatilities of its value. We denote the volatility in the i'th rating by σ(t) = σ i , i = 1, 2, · · · , n, and in addition, they satisfy 0 < σ 1 < σ 2 < · · · < σ n−1 < σ n < ∞, which means that in the highest credit rating, the corporation shows the smallest volatility σ 1 and in the lowest credit rating, it shows the largest volatility σ n . The stochastic interest rate r is supposed to satisfy the Vasicek model [27,33] dr = a(t)(θ(t) − r)dt + σ r (t)dW r t , which is widely popular in financial application, where the parameters a, θ, σ r are supposed to be positive constants in this paper. σ r is the volatility of the interest rate. θ is considered as the central location or the long-term value. a determines the speed of adjustment.
We suppose that the corporation issues only one defaultable bond with face value F. The effect of corporate value on the bond value is focused on and the discount value of bond is considered. Denote by φ t the discount value of bond at time t. The corporation exhibits two risks, the default risk and credit rating migration risk. The corporation can default before maturity time T . The default time τ d is the first moment when the corporate value falls below the threshold K, namely that where K < F · D(t, T ), where 0 < D(t, T ) < 1 is the discount function. Once the corporation defaults, the investors will get what is left. Hence, φ t (K) = K and at the maturity time T , the investors can get φ T = min{S T , F}. On the other hand, the credit regions are determined by the leverage γ(t) = φ t /S t . Denote the thresholds of leverage γ(t) by γ i , i = 1, 2, · · · , n − 1, and they satisfy 0 < γ 1 < γ 2 < · · · < γ n−2 < γ n−1 < 1.
The credit rating migration times are the first moments when the corporate credit rating is upgraded or downgraded. They are defined as follows: τ 1 is the first moment that the corporation degrades from the highest credit rating. τ n is the first moment that the corporation upgrades from the lowest credit rating. τ i,i+1 is the first moment that the corporation jumps up into the i + 1'th credit rating from the i'th credit rating, while τ i,i−1 is the first moment that the corporation jumps down into the i − 1'th credit rating from the i'th credit rating.

Cash flow
Once the credit rating migrates before the maturity T , a virtual substitute termination happens, namely that the bond is virtually terminated and substituted by a new one with a new credit rating [35]. Thus, there is a virtual cash flow of the bond. Denote the bond values in different credit ratings by φ i (t, S ), i = 1, 2, · · · , n. Then they are the conditional expectations as follows and for i = 2, 3, · · · , n − 1, where χ is the indicative function, satisfying χ = 1 if the event happens and otherwise, χ = 0.

PDE problem
Suppose that the correlation between the interest rate and the corporate value is given by dW r t ·dW t = ρt, −1 ≤ ρ ≤ 1. By the Feynman-Kac formula, we can derive that φ i , i = 1, 2, · · · , n, are functions of time t, interest rate r and value S . They satisfy the following PDE in their regions and for i = 2, 3, · · · , n − 1, and boundary condition φ n (t, r, K) = K.
The bond value is continuous when it passes a rating threshold, i.e., φ i = φ i+1 on the rating migration boundaries, where i = 1, 2, · · · , n − 1. Also, if we construct a risk free portfolio π by longing a bond and shorting ∆ amount asset value S , i.e., π = φ − ∆S and such that dπ = rπ, this portfolio is also continuous when it passes the rating migration boundaries, namely that π i = π i+1 or ∆ i = ∆ i+1 on the rating migration boundaries, where i = 1, 2, · · · , n − 1. By Black-Scholes theory [16], it is equivalent to Denote by P(t, r) the value of a guaranteed zero-coupon bond with face value 1 at the maturity t = T , where the interest rate follows the Vasicek model. By the Feynman-Kac formula, P(t, r) satisfies the following PDE ∂P ∂t with terminal condition P(T, r) = 1, whose explicit solution is solved as P(t, r) = e A(T −t) [13], where and Then ψ i , i = 1, 2, · · · , n satisfy and for i = 2, 3, · · · , n − 1, and ∂ψ n ∂t +σ 2 n y 2 2 ∂ 2 ψ n ∂y 2 = 0, ψ n > γ n−1 y, (2.9) The terminal conditions are given as and the boundary conditions are and ∂y , on the rating migration boundary (2.12) for i = 1, 2, · · · , n − 1.
The domain will be divided into n rating regions Q i , i = 1, 2, · · · , n. We will prove that the domain can be separated by n − 1 free boundaries x = λ i (t), i = 1, 2, · · · , n − 1. These boundaries are a prior unknown since they should be solved by equations where u is also a priori unknown. Since we have assumed that (3.1) is valid cross the free boundaries, (2.12) implies that for i = 1, 2, · · · , n − 1, In the work [35] and [12], where the former model is subject to constant interest rate and the later one is subject to stochastic interest rate but without default boundary, the processṡ(t) in (3.6) is replaced by a constant. The presence ofṡ(t) indeed leads to some technical differences in deriving the estimates in the following argument. We can rewrite the formula ofṡ(t) aṡ It is not difficult to analyze and derive that one of the following conditions holds, then there holdṡ 4. Approximated problem with some uniform estimates
We approximate H(ξ) by a C ∞ function H (ξ) satisfying Consider the approximated free boundary problem and Problem (4.1)-(4.3) admits a unique classical solution u . Now we proceed to derive some estimates for u .

Estimates for approximated problem
Some uniform estimates are presented in this section, which are sufficient to obtain the existence and uniqueness of solution to problem (3.6)-(3.10).
Proof. It is easy to verify that 0 is the lower solution of u and meanwhile, e −(x+log K) and 1 are upper solutions. The result is a direct application of comparison principle.
Proof. It is easy to see thatσ 2 can be written aŝ Differentiating (4.1) with respect to x gives It is known that ∂u then letting x → 0, it holds that ∂u ∂x (t, 0) ≤ 0. Thus it follows by maximum principle that there holds ∂u ∂x ≤ 0.
On the other hand, sinceṡ(t) is uniformly bounded, take an appropriate value C > 0, such that It follows by the comparison principle that there holds ∂u ∂x ≥ −C. Lemma 4.3. Let u be the solution of problem (4.1)-(4.3). Suppose that one of the conditions (3.11)-(3.13) holds. Then there exist constants C 1 , C 2 , C 3 and C 4 , independent of , such that Proof. Differentiating (4.1) with respect to t gives Since and
we can write According to the formulas ofσ 2 On the other hand, there exists λ i (t) such that H (u (t, λ i (t)) − γ i ) attains its maximum and Denote by y(t) the solution of the following ODE where the constant C is given as the one in Lemma 4.2. At x = − log K, ∂ 2 u ∂x 2 (0, x) produces a Dirac measure of density −1. Thus ∂u ∂t (0, x) ≤ 0 in the distribution sense. In addition, since the second order compatibility condition is satisfied at (0, 0), we have ∂u ∂t is continuous at (0, 0). Meanwhile, ∂u ∂t (0, t) = 0. By further approximating the initial data with smooth function if necessary, there holds by the comparison principle that if we set y 0 = 0. The ODE (4.5) can be solved formally as Since ∂σ 2 n ∂t → 0 as t → ∞ and H converges to the Dirac measure as → 0, this implies that exp( t 0 p(s)ds) is uniformly bounded with respect to t. With regard to q(t), in addition tos(t) → 0 as t → ∞, it is known that q(t) is also uniformly bounded. Hence, we conclude that y(t) is uniformly bounded with respect to t.
for |x + log K| < ρ 2 , 0 < t ≤ ρ 2 4 . Note that (4.4) can be rewritten as Asṡ(t) is uniformly bounded, we take a sufficiently large constant C 3 such that C 3 ≥ sup t≥0 |ṡ(t)|C, where constant C is the one given in Lemma 4.2. Then if ∂u ∂t < −C 3 , there holds Denote by z(t) the solution of the following ODE which can be solved as As ∂σ 2 n ∂t ands(t) tends to 0 as t → ∞, z(t) is uniformly bounded. Moreover, z(t) is also decreasing if z 0 ≤ 0. Take the initial data |z 0 | and constant C 3 sufficiently large such that We claim that the region is an empty set. If not, on the parabolic boundary of this region, we clearly have ∂u ∂t ≥ −C 3 , which implies by the comparison principle that in this region. This is a contradiction.
Remark 4.4. In the work [35], it was proved that ∂u ∂t ≤ 0, which is different from the result shown in Lemma 4.3. However, although we get a similar result to Lemma 5.4 in [12], the proof is very different and more technical. This is due to the joint effect of stochastic interest rate and default boundary.
Lemma 4.6. Let λ i , i = 1, 2, · · · , n − 1, be the approximated free boundaries defined in (4.6). Suppose that one of the conditions (3.11)-(3.13) holds. Then there exist constants C 1 , C 2 , independent of , such that This means that region {x > − log γ 1 K} is in the highest rating region and hence Denote by m = sup t≥0ṡ (t) and v x .

Consider the region
R } are completely and uniformly within the highest or lowest rating region (independent of ). Thus by compactness and the strong maximum principle, on these four boundaries, it holds that − ∂u ∂x ≥C > 0 for someC independent of . It follows that which completes the proof of the lemma.

By the classical parabolic theory, it is also clear that the solution is in
where and for i = 2, 3, · · · , n − 1 and Now we prove the uniqueness of solution to the problem (3.6)-(3.10).
Besides, at t = 0, λ i (0) =λ i (0) = − log K − log γ i , i = 1, 2, · · · , n − 1. As shown in Lemma 4.7, there exists a constant C * > 0, such that u x < −C * andũ x < −C * on the region Ω * defined by (4.7). Then by the implicit function theorem, there exists ρ > 0, such that when 0 < t < ρ, where C is a positive constant, whose value may change line on line but makes no difference. Let w = u −ũ and denote byσ andσ the corresponding coefficients, then w satisfies u,ũ and their derivatives decay exponentially fast to 0 as x → ∞. Multiplying (5.2) by w on both sides and integrating x from 0 to ∞ gives and ∂ũ ∂t and ∂ũ ∂x are uniformly bounded outside the region Q t 0 , we conclude that they are bounded for . Since w decays exponentially to 0 as x → ∞, for any t > 0, there exists x 0 < ∞ such that max 0<x<∞ w 2 (t, x) = w 2 (t, x 0 ).
for some x * ∈ (x 0 , x 0 + ). Then there holds We now proceed to estimate the left side of (5.3). First, we have asṡ(t) is uniformly bounded. Combining the above inequalities (5.5)-(5.8), taking into account (5.3), we drive
It is easy to see that there exists a constant C 0 > 0 such that Then for sufficiently small , by applying the Gronwall's inequality, we conclude that w ≡ 0. This proves the uniqueness for 0 ≤ t ≤ ρ. A close examination of the proof indicates that the uniqueness result can be extended to any time interval, where ∂u ∂x is strictly negative, which is already verified in (4.8).

Convergence
The coefficients in the model of [35] are time-homogeneous and moreover, the solution is decreasing in time. However, both of these are not the cases in our model. We cannot take advantage of the decreasing property of solution to obtain the convergence. In this paper, we are motivated by the idea of [12] and obtain the convergence by two steps. The first step is to show that the solution of problem (3.6)-(3.10) converges to the solution of some auxiliary problem defined below, while in the second step, we show that the solution of auxiliary problem converges to the solution of static problem (6.1)-(6.3). Now define an auxiliary free boundary problem as follows with initial conditionū 0 and boundary conditionū(t, 0) = 1, All the results derived above involving the existence, uniqueness and some properties of solution presented in Lemmas 4.1-4.6 hold for solution of problem (6.18)- (6.19). We have to notice that as the presence of default boundary, although we follow the idea of [12], the technical proofs are different, especially in the step, i.e., the convergence from the original solution to the auxiliary solution.

Convergence from u toū
Sinceσ i (t) →σ i , i = 1, 2, · · · , n, andṡ(t) → θ − σ 2 r 2a 2 as t → ∞, then for any > 0, there exists a T > 0 such that for t ≥ T , Letū 0 (x) = u(T, x) and denote by u T (t, x) = u(t + T, x) for t ≥ 0. We have u T (t, λ i (t)) =ū(t,λ i (t)) = γ i , i = 1, 2, · · · , n − 1, and Similarly to the proof of Theorem 5.2, by the implicit function theorem, there exists a ρ > 0, such that when 0 < t < ρ, where C is a positive constant, whose value may change line on line but makes no difference. Let w = u T −ū. Then w satisfies whereṡ T (t) =ṡ(t + T ), and As u,ū and their derivatives decay exponentially fast to 0 as x → ∞, multiplying (6.22) by w on both sides and integrating x from 0 to ∞ gives By Lemmas 4.2 and 4.3, we know that ∂ū ∂t and ∂ū ∂x are uniformly bounded as the initial data is set to bē u 0 (x) = u(T, x) for a sufficiently large T . It follows that

AIMS Mathematics
Volume 5, Issue 6, 7746-7775. where hold by (6.20) and the exponential decay of w. Using (5.4), there holds On the other hand, there holds We now proceed to estimate the left side of (6.23) by The remaining terms in (6.23) can be estimated by Combining the inequalities (6.24)-(6.28), we have
It follows that Then taking a sufficiently small , which means a sufficiently large T , by applying the Gronwall's inequality, we conclude that w = 0, namely that u T (t, x) =ū(t, x) for 0 ≤ t ≤ ρ. The result can be extended to any interval of t.

Convergence fromū to ψ
The convergence from the solution of auxiliary problem (6.18)- (6.19) to the solution of problem (6.1)-(6.3) is proved by a Lyapunov argument, which is similar to the procedure in [12,21,32], but with some necessary modifications to fit the model in this paper. For instance, we have extend the solution with domain [0, ∞) to the whole real line. The first step is to present the formal construction of a Lyapunov function, ignoring the integrability of any arisen integral. The second step is to verify the integrability of those integrals arisen in the formal construction. The third step is to complete the proof of convergence.
Denote byŪ the extension ofū, who is the solution of auxiliary problem (6.18)- (6.19), from x ∈ [0, ∞) to the real line R, namely that Following [12], let V(x, u, q) be a undetermined function and set Formally, assuming the integrability, we also have Denote by ρ = V qq . Suppose that V(x, u, 0) = V q (x, u, 0) = 0 as in [12]. Then we have It follows that and Then (6.29) can be written as To ensure (6.29) holds, it should be Formally, denote by v the solution of the following equation Then by (6.30) and (6.31), there holds which can be solved as ρ(x 0 , u 0 , q 0 ) = exp Replacing x 0 by x, u 0 by u and q 0 by q, we have dz . (6.32) Integrating the Lyapunov function and assuming E(t) ≥ 0, we have Following the formal construction, we proceed to the second step, namely that verify the integrability of those integrals arisen in the formal construction of the Lyapunov function. As indicated by Liang et al. [21], there are two problems, where the first one is that ρ grows exponentially as x → ±∞, while the second one is that the coefficient in (6.31) is discontinuous and the theory of ODE cannot be applied directly. These two problems will be overcame through the approximated solution of problem (4.1)-(4.3) with all the uniform estimates. We begin this step by defining for R > 0, whereŪ is the approximation ofŪ with the coefficient of the approximated problem given asσ and ρ is defined by (6.32) withσ 2 replaced byσ 2 . Thus, (6.31) can be solved on the real line x ∈ R and V is well defined. Meanwhile, it also satisfies Then we have Lemma 6.1. Letū be the approximated solution of problem (6.18)-(6.19). Then for any K 1 > 0, there exist constants C 0 , K 2 > 0 such that for x > C 2 , there holds where C 2 is the constant given in Lemma 4.5.
Proof. We know that for x > C 2 where the operator By for x > C 2 . Then for any K 1 > 0, there exist constants C 0 , K 2 > 0 such that From the formula of ρ ρ (x, u, q) = exp Then there holds where the constant C is independent of , which implies that according to (6.28) and the setting ofŪ.

AIMS Mathematics
Volume 5, Issue 6, 7746-7775. Now we arrive at the third step, completing the proof of convergence fromū to ψ. Denote bȳ u m (t, x) =ū(t + m, x) and considerū m as a sequence of functions on [0, 1] × (0, ∞). Sinceū m is a bounded sequence in W 1,2 ∞ ([0, 1] × (0, ∞)), we derive by the embedding theorem that there exists a subsequence m j of m and a functionψ such that as m j → ∞, there holds u m j →ψ in C Clearly, the second term on the right side converges to the corresponding integral ofψ as m = m j → ∞.
With regard to the first term on the right side of (6.36), there holds By the weak-star convergence, the first term on the right side of (6.32) converges to the corresponding integral ofψ. The second term on the right side of (6.37) is bounded by C ∞ 0 |σ 2 (ū m ) −σ 2 (ψ)| f dx, which converges to 0 by the dominated convergence theorem. By the convergence (6.35), we have i (t + m j ), i = 1, 2, · · · , n − 1. (6.38) We choose t min i, j , t max i, j ∈ [0, 1] such that inf 0≤t≤1λ i (t + m j ) =λ i (t min i, j + m j ), sup 0≤t≤1λ i (t + m j ) =λ i (t max i, j + m j ).
However, by the uniqueness of solution to static problem (6.1)-(6.3), there should hold λ min i =λ max i = λ * i , i = 1, 2, · · · , n − 1, and ψ ≡ψ. In addition, the uniqueness implies that all subsequences limit should be uniform and thus the full sequence must converge as m → ∞.

Conclusion
In this paper, we study a free boundary problem for pricing a defaultable corporate bond with multiple credit rating migration risk and stochastic interest rate. By using PDE techniques, the existence, uniqueness some regularities of solution are obtained to support the rationality of the model to pricing a defaultable corporate bond. In [35], it is shown that the solution and rating boundaries of the free boundary problem pricing a defaultable corporate bond with constant interest rate are all decreasing with respect to time, which do not hold any more in our model with stochastic interest rate, as the coefficients of model are all time heterogeneous. Furthermore, we present the asymptotic behavior of solution to this pricing model. Asymptotic behaviors of solution to free boundary problems pricing corporate bonds with credit rating migration risk have been analyzed in [12,21,32]. The asymptotic solution of the model with only one migration boundary can be solved explicitly [21], while in the works [12,32], where the models are subject to multiple migration boundaries, it is not the case. However, interestingly, in this paper, although our model is also subject to multiple migration boundaries, it is proved that the asymptotic solution can be solved explicitly. We conclude that if the maturity T is sufficiently large, we can valuate the defaultable corporate bond with multiple credit rating migration risk and stochastic interest rate by an explicit pricing formula as follows: where S (t) is the corporate value, K is default threshold and ψ is the steady status given in Section 6, whose explicit form can be obtained by solving the linear algebraic equation set.