Citation: M Toufiq Reza, Jessica Nover, Benjamin Wirth, Charles J Coronella. Hydrothermal carbonization of glucose in saline solution: sequestration of nutrients on carbonaceous materials[J]. AIMS Energy, 2016, 4(1): 173-189. doi: 10.3934/energy.2016.1.173
[1] | Longfei Wei, Lu Liu, Jialong Hou . Pricing hybrid-triggered catastrophe bonds based on copula-EVT model. Quantitative Finance and Economics, 2022, 6(2): 223-243. doi: 10.3934/QFE.2022010 |
[2] | Lianzhang Bao, Guangliang Zhao, Zhuo Jin . A new equilibrium trading model with asymmetric information. Quantitative Finance and Economics, 2018, 2(1): 217-229. doi: 10.3934/QFE.2018.1.217 |
[3] | Manuela Larguinho, José Carlos Dias, Carlos A. Braumann . Pricing and hedging bond options and sinking-fund bonds under the CIR model. Quantitative Finance and Economics, 2022, 6(1): 1-34. doi: 10.3934/QFE.2022001 |
[4] | Jie Xiong, Geng Deng, Xindong Wang . Extension of SABR Libor Market Model to handle negative interest rates. Quantitative Finance and Economics, 2020, 4(1): 148-171. doi: 10.3934/QFE.2020007 |
[5] | Lichao Tao, Yuefu Lai, Yanting Ji, Xiangxing Tao . Asian option pricing under sub-fractional vasicek model. Quantitative Finance and Economics, 2023, 7(3): 403-419. doi: 10.3934/QFE.2023020 |
[6] | Richard J. Cebula, Don Capener . The impact of federal income tax rate cuts on the municipal bond market in the U.S.: A brief exploratory empirical note. Quantitative Finance and Economics, 2018, 2(2): 407-412. doi: 10.3934/QFE.2018.2.407 |
[7] | Serçin ŞAHİN, Serkan ÇİÇEK . Interest rate pass-through in Turkey during the period of unconventional interest rate corridor. Quantitative Finance and Economics, 2018, 2(4): 837-859. doi: 10.3934/QFE.2018.4.837 |
[8] | Guillermo Peña . Interest rates affect public expenditure growth. Quantitative Finance and Economics, 2023, 7(4): 622-645. doi: 10.3934/QFE.2023030 |
[9] | Bonang N. Seoela . Efficacy of monetary policy in a currency union? Evidence from Southern Africa's Common Monetary Area. Quantitative Finance and Economics, 2022, 6(1): 35-53. doi: 10.3934/QFE.2022002 |
[10] | Mustafa Tevfik Kartal, Özer Depren, Serpil Kılıç Depren . The determinants of main stock exchange index changes in emerging countries: evidence from Turkey in COVID-19 pandemic age. Quantitative Finance and Economics, 2020, 4(4): 526-541. doi: 10.3934/QFE.2020025 |
In global financial markets, the pricing of a bond becomes more complicated due to different credit rating and different interest rate. A change of credit rating or a change of interest rate may cause drastic swing of bond prices. To manage the complication, one needs to understand how the credit risk impacts on the bond markets and how it can be measured. There are two types of credit risks: defaults and credit rating migrations. In the previous researches, the defaults are paid more attentions. However, in the financial crisis in 2008 and European debt crisis in 2010, the credit rating migrations played a much more important role in the finance and economics than what people expected. Therefore, the study of credit rating migrations becomes significant. Another observed factor in a financial market is that bond prices strongly depend on the interest rate. This leads to a crucial question of how to price a corporate bond with credit rating migration risks and interest rate change.
There are many existing research papers on bond pricing models, see examples, Briys and de Varenne (Briys and Varenne, 1997), Tsiveriotis and Fernandes (Tsiveriotis and Fernandes, 1998).To measure default risks, there are two basic frameworks: a structure model, and a reduced form one. It is well-known that the first structure model for credit default proposed by Merton (Merton, 1974) is under the assumption that the firm's value follows a geometric Brownian motion process and the corporate bonds are considered as contingent claims. The default event may only occur at the maturity of the claim. In 1976, Black and Cox extended Merton's model to a first-passage-time model by introducing safety covenants. This extends the original model which gives bondholders the right to recognize a firm's value if its asset falls below some exogenously given threshold. In the modified model(Black and Cox, 1976), the default is allowed to occur any time up to the debt's maturity. By the reduced form framework, the credit events are considered to be governed by exogenous reason stochastically, see (Duffie and Singleton, 1999; Jarrow and Turnbull, 1995; Lando, 1998; Leland and Toft, 1996).
In study of credit rating migration, a major tool is the Markov chain. Most researchers mainly adopted a transfer intensity matrix. By this approach, the reduced form framework is naturally developed for the dynamic process of credit rating migrations, see Jarrow et al. (1997), Das and Tufano (1996), Hurd and Kuznetsov (2007), Lando (1998), Thomas et al. (2002) and so forth. The transfer intensity matrix usually comes from general statistical data, which do not include any particulary firm's information. However, a firm's own value plays a key role of its credit rating migration. In this situation, the Markov chain alone cannot fully capture the credit-rating mitigation for a firm. From the existing literatures, few research took this into account. Liang et al. (2015) in 2015 used a structure model for pricing bonds with credit rating migrations at first time. They gave a predetermined migration threshold where a firm's value is divided into high and low rating regions, under the assumption that the value of the firm follows a stochastic process. By using Feynman-Kac theorem, these models can be reduced to certain boundary value problems of partial differential equations. This PDE method is very different from the traditional approaches in the research field of corporative bond modelling. However, the rating migration boundary is usually not predetermined in the real financial world. It depends on the proportion of the debt and the value of the firm. From this point of view, the migration boundary should be a free boundary. Recently, Hu et al. (2015) and liang et al. (2016) proposed a new model to reflect this fact under the fixed-interest rate. They obtained the existence and uniqueness for the free boundary problem. Also, they proved that there is a travelling wave in the problem.
For pricing a bond, interest rate is a sensitive factor. Many stochastic models are used for the model, see Longstarff and Schwartz (1995) as an example. However, when a stochastic interest rate is considered, the bond model becomes two dimensions, which is much harder to analyze.
So far, all valuations for credit rating migration via pricing a corporate bond by structure model are assumed a constant interest rate. However, the stochastic interest rate for a bond is essential. So the extension for the model to stochastic interest rate is necessary. However, the extension meets a great difficult that the original free boundary model from one dimension to a higher one. That is the main contribution and solved problem in this paper.
In the present paper, we follow the ideas from (Hu et al., 2015) to derive a new model for a bond price in which credit-rating may change and the interest rate follows a stochastic process. To be more precise, we assume that the interest rate follows a Vasicek's stochastic process and the bond price follows a geometric Brownian motion with a different volatility in different credit-rating region. A coupled system of PDEs is derived. By employing a dimension-reduction technique we are able to simplify the coupled system to a single equation with a nonlinear jump coefficient in the leading term of the equation. Under certain conditions, the existence and uniqueness are established for the new model. We would like to indicate that our approach is different from that in (Hu et al., 2015; Liang et al., 2016). Moreover, the uniqueness proof in Section 4 is new to our knowledge. The method is also suitable for n-dimensional problems. Furthermore, the regularity of the free boundary is also considered. By transforming the problem to a Stefan-like problem, we prove that the free boundary is smooth. Some numerical results are also illustrated in the present paper. Those numerical results are well matching the theoretical analysis.
The paper is organized as follows: In Section 2, a new model with credit-rating migration and stochastic interest rate is rigorously derived under usual assumptions. In Section 3, we use the dimension reduction technique to simplify the model to a free boundary problem. In Section 4, we first use various PDE techniques to derive uniform estimates for the approximate problem. The existence and uniqueness are then established by using a compactness argument. Uniqueness is proved via an energy method. In Section 5, further properties of the solution are shown. In Section 6, some numerical results are presented. Some concluding remarks are given in section 7.
Let $(\Omega, \cal{F}, P)$ be a complete probability space. We assume that the firm issues a corporate bond, which is a contingent claim of its value on the space $\cal{F}$.
Assumption 2.1. Interest rate $r_t$ follows Vasicek process:
$ dr_t=a(t)(\theta(t)-r_t)dt+\sigma_r(t)dW_t^r, $ |
where $a(t)$, $\theta(t)$ and $\sigma_r(t)$ are given positive functions, which are supposed to be positive constants $a, \theta, \sigma _r$ for the simplification in this paper. $W_t^r$ is the Brownian motion which generates the filtration $\{{\cal F}_t\}$.
We would like to point out that our model is a big improvement from the assumption of a constant interest rate to the stochastic one. There are many stochastic models for the interest rate, Vesicek model is popular one. The advantage point of this model is simply and easily calculable, while the weak point is the possibility of negative interest rate, which is proved to be very low. We assume that the interest rate follows Vesicek model also because that in our dimension reducing technique, only Vesicek model works.
Assumption 2.2 (the firm asset with credit rating migration). Let $S_t$ denote the firm's value in the risk neutral world. It satisfies
$
dSt={rtStdt+σH(t)StdWt,inhighratingregion,rtStdt+σL(t)StdWt,inlowratingregion,
$
|
where $r_t$ is the risk free interest rate, and
$ \sigma_{H}(t)<\sigma_{L}(t) $ | (1) |
represent volatilities of the firm under the high and low credit grades, respectively. They are assumed to be deterministic and differentiable on $[0, T]$ and have a positive lower bound:
$ 0<\tilde\sigma_0\leq \sigma_{H}(t)<\sigma_{L}(t)\leq \tilde\sigma_1<\infty, 0\leq t\leq T. $ | (2) |
$W_t$ is the Brownian motion which generates the filtration $\{{\cal F}_t\}$, and
$Cov(dW_t, dW_t^r)=\rho dt, $ |
where $-1 < \rho < 1$ is a constant.
It is reasonable to assume (1), namely, that the volatility in high rating region is lower than the one in the low rating region. This implies that there is a jump for the volatility when the rating changes.
Assumption 2.3 (the corporate bond). The firm issues only one corporate zero-coupon bond with face value $F$. Denote $\Phi_t$ the discount value of the bond at time $t$, which is also the function of $r_t$. Therefore, on the maturity time $T$, an investor can get $\Phi_T=\min\{S_T, F\}$.
Assumption 2.4 (the credit rating migration time). High and low rating regions are determined by the proportion of the debt and value. The credit rating migration time $\tau_1 $ and $\tau_2 $ are the first moments when the firm's grade is downgraded and upgraded, respectively, as follows:
$ \tau_1=\inf\{t>0|\Phi_0/S_0<\gamma, \Phi_t/S_t\geq \gamma\}, \tau_2=\inf\{t>0|\Phi_0/S_0>\gamma, \Phi_t/S_t\leq \gamma\}, $ |
where $\Phi_t=\Phi_t(S_t, r, t)$ is a contingent claim with respect to $S_t$ and
$ 0<\gamma<1 $ | (3) |
is a positive constant representing the threshold proportion of the debt and value of the firm's rating.
Once the credit rating migrates before the maturity $T$, a virtual substitute termination happens, i.e., the bond is virtually terminated and substituted by a new one with a new credit rating. There is a virtual cash flow of the bond. Denoted by $\Phi_H(y, r, t)$ and $\Phi_L(y, r, t)$ the values of the bond in high and low grades respectively. Then, they are the conditional expectations as follows:
$
ΦH(y,r,t)=Ey,t[e−∫Ttrsdsmin(ST,F)⋅1τ1≥T+e−∫τ1trsdsΦL(Sτ1,rt,τ1)⋅1t<τ1<T|St=y>1γΦH(y,t),rt=r],
$
|
(4) |
$
ΦL(y,r,t)=Ey,t[e−∫Ttrsdsmin(ST,F)⋅1τ2≥T+e−∫τ2trsdsΦH(Sτ2,rt,τ2)⋅1t<τ2<T|St=y<1γΦL(y,t),rt=r],
$
|
(5) |
where ${\bf 1}_{event}=\left\{1, if "event" happens ,0, otherwise.
By Feynman-Kac formula (Dixit and Pindyck, 1994), it is not difficult to drive that $\Phi_H$ and $\Phi_L$ are the function of the firm value $S$, interest rate $r$ and time $t$. They satisfy the following partial differential equations in their regions:
$
∂ΦH∂t+12σ2HS2∂2ΦH∂S2+σrσHρS∂2ΦH∂S∂r+12σ2r∂2ΦH∂r2+rS∂ΦH∂S+a(θ−r)∂ΦH∂r−rΦH=0,S>1γΦH,0<t<T,
$
|
(6) |
$
∂ΦL∂t+12σ2LS2∂2ΦL∂S2+σrσLρS∂2ΦL∂S∂r+12σ2r∂2ΦL∂r2+rS∂ΦL∂S+a(θ−r)∂ΦL∂r−rΦL=0,S<1γΦL,0<t<T,
$
|
(7) |
with the terminal condition:
$ \Phi_H(S, T)=\Phi_L(S, T)=\min\{S, F\}. $ | (8) |
(4) and (5) imply that the value of the bond is continuous when it passes the rating threshold, i.e.,
$ \Phi_H=\Phi_L =\gamma S\quad \mbox{on the rating migration boundary}. $ | (9) |
Also, as in (Hu et al., 2015; Jin et al., 2016), we have
$ \frac {\partial \Phi_H}{\partial S} =\frac{\partial \Phi_L}{\partial S} \quad \mbox{on the rating migration boundary}. $ | (10) |
This is a two-dimensional free boundary problem of an terminate value parabolic problem with discontinuous leading coefficient.
The equations (6) and (7) are two-dimensional. For the purpose of simplifying the problem, we could use a technique to reduce the dimension.
Denote $P(r_t, t:T)$ to be the value of a zero-coupon bond faced 1 at time $T$. So that, it satisfies
$ \frac{\partial P}{\partial t} +\frac 12 \sigma _r^2 \frac{\partial ^2 P}{\partial r^2} +a(\theta-r)\frac{\partial P}{\partial r}-rP = 0, $ | (11) |
which admits a unique solution (Hall, 1989)
$ P(r, t:T)= A(t)e^{-rB(t)}, $ | (12) |
where
$ A(t)=e^{\{\frac 1{a^2}[B^2(t)-(T-t)](a^2\theta-\frac {\sigma _r^2}2)-\frac{\sigma ^2}{4a}B^2(t)\}}, \quad B(t)=\frac 1a(1-e^{-a(T-t)}).$ |
Now make a transform by
$y=\frac S{P(r, t;T)}, \quad V_H(y, t)=\frac{\Phi_H(S, r, t)}{P(r, t;T)}, \quad V_L(y, t)=\frac{\Phi_L(S, r, t)}{P(r, t;T)}.$ |
Then $V(y, t)$ satisfies
$ \frac{\partial V_H}{\partial t} +\frac 12\hat{\sigma }_H^2y^2 \frac{\partial ^2 V_H}{\partial y^2} = 0, \quad y>\frac 1\gamma V_H, \ 0<t<T, $ | (13) |
$ \frac{\partial V_L}{\partial t} +\frac 12\hat{\sigma }_L^2y^2 \frac{\partial ^2 V_L}{\partial y^2} = 0, \quad y<\frac 1\gamma V_L, \ 0<t<T, $ | (14) |
$ V_H(y, T)=V_L(y, T)=\min\{y, F\}, $ | (15) |
where
$ \hat{\sigma }_i=\sqrt{\sigma _i^2+2\rho\sigma _i\sigma _rB(t)+\sigma _r^2B^2(t)}, \quad i=H, L. $ | (16) |
And on credit rating migration boundary, there holds
$ V_H(y, t)=V_L(y, t)=\gamma y, \qquad \frac{\partial V_H}{\partial y}= \frac{\partial V_L}{\partial y}. $ | (17) |
Using the standard change of variables $x=\log y$ and renaming $T-t$ as $t$, and defining
$
\phi(x, t)=\left\{VH(ex,T−t) in high rating region,VL(ex,T−t) in low rating region, \right.
$
|
using also (6) and (7), we then derive the following equation from (6), (7):
$ \frac{\partial \phi}{\partial t} -\frac12 \sigma ^2\frac{\partial ^2\phi}{\partial x^2} +\frac12\sigma ^2\frac{\partial \phi}{\partial x} = 0, \quad (x, t)\in Q_T\backslash \Gamma_T $ | (18) |
where
$ Q_T=(-\infty, +\infty)\times(0, T], \quad \Gamma_T=\{(x, t)|\phi(x, t)=\gamma e^{x}\} $ | (19) |
and $\sigma $ is a function of $\phi$ and $(x, t)$, i.e.,
$
\sigma = \sigma (\phi, x, t) = \left\{ˆσHifϕ<γex,ˆσLifϕ≥γex. \right.
$
|
(20) |
The constant $\gamma $ is defined in (3), and $\hat\sigma _H, \hat\sigma _L$ are defined in (16). There exist constants $\sigma _0, \sigma _1$, such that
$ 0<\sigma_0\le\hat{\sigma }_H, \hat{\sigma }_L\le\sigma _1<\infty. $ | (21) |
Without loss of generality, we assume $F=1$. Equation (18) is supplemented with the initial condition (derived from (8))
$ \phi(x, 0) = \phi_0(x)=\min\{e^x, 1\}, \quad -\infty<x<\infty. $ | (22) |
The domain will be divided into the high rating region $Q_{T}^H$ where $\phi < \gamma e^{x}$ and a low rating region $Q_{T}^L$ where $\phi > \gamma e^{x}$. We shall prove that these two domains will be separated by a free boundary $x=s(t)$, and
$ Q_T^H = \{x>s(t), 0<t<T\}, \qquad Q_T^L = \{x<s(t), 0<t<T\}. $ | (23) |
In another word, $s(t)$ is apriorily unknown since it should be the solved by the equation
$ \phi(s(t), t)=\gamma e^{s(t)}, $ |
where the solution $\phi$ is apriorily unknown.
Let $H(\xi)$ be the Heaviside function, i.e., $H(\xi)=0$ for $\xi < 0$ and $H(\xi)=1$ for $\xi>0$. Then we can rewrite (20) as
$ \sigma (t, \phi) = \hat\sigma _H + (\hat\sigma _L-\hat\sigma _H) H(\phi - \gamma e^{x}). $ |
We use a standard approximation for $H(\xi)$ and $\phi_0(x)$ such that
$ H_\varepsilon (\xi)\in C^{\infty}(\mathbb{R}^1)\quad \mbox{and}\quad H_\varepsilon (\xi)=H(\xi) \;\;\text{for}\;\; |\xi|\ge\varepsilon , \\ \phi_{0\varepsilon }(x)\in C^3(\mathbb{R}^1)\cap W^{1, \infty}(\mathbb{R}^1), \\ \phi'_{0\varepsilon }(x)-\gamma e^x<0, \quad \;\;\text{for}\;\;0<x<\ln\frac1\gamma . $ |
Now we consider the following approximated problem
$ \frac{\partial \phi_\varepsilon }{\partial t} -\frac12 \sigma ^2_\varepsilon \Big(\frac{\partial ^2\phi_\varepsilon }{\partial x^2} -\frac{\partial \phi_\varepsilon }{\partial x}\Big) = 0, \quad (x, t)\in Q_T, $ | (24) |
$ \phi_\varepsilon (x, 0) = \phi_{0\varepsilon }(x), \quad -\infty<x<\infty, $ | (25) |
where
$ \sigma _\varepsilon = \hat\sigma _H + (\hat\sigma _L-\hat\sigma _H) H_\varepsilon (\phi_\varepsilon - \gamma e^{x}). $ | (26) |
Lemma 4.1. The problem (24)-(25) admits a unique solution
$ \phi_\varepsilon (x, t)\in C^{2, 1}(Q_T)\cap C^{\alpha, \frac\alpha2}(\overline{Q}_T), \quad \;\;\text{for}\;\; \;\;\text{for}\;\;all\alpha\in(0, 1) $ |
Moreover, the following estimates hold uniformly with respect to $\varepsilon $:
a) $0\le\phi_\varepsilon (x, t)\le 1$, $0\le\phi_{\varepsilon x}(x, t)\le C_1$
b) $\sup\limits_{0\le t\le T}\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2dx+\iint\limits_{Q_T}(\phi_{\varepsilon xx}^2+\phi_{\varepsilon t}^2) dxdt\le C_2$
where $C_1$ and $C_2$ depend only on known data.
Proof. By the theory of parabolic equation, we know that the problem (24)-(25) admits a unique classic solution $\phi_\varepsilon \in C^{2, 1}(Q_T)\cap C^{\alpha, \frac\alpha2}(\overline{Q}_T)$.
Now we derive the uniform estimates. First of all, by the maximum principle for Cauchy problem, we see
$ 0\le\phi_\varepsilon (x, t)\le 1, \quad (x, t)\in\overline{Q}_T. $ |
Let $\psi_\varepsilon (x, t)=\phi_{\varepsilon x}(x, t)$, then
$ \psi_{\varepsilon t}-\frac{\partial }{\partial x}\big[\frac12\sigma _\varepsilon ^2\big(\psi_{\varepsilon x}-\psi_\varepsilon \big)\big], \quad (x, t)\in\overline{Q}_T, \\ \psi_\varepsilon (x, 0)=\phi'_{0\varepsilon }(x)\ge 0, \quad -\infty<x<+\infty $ |
Again, since $\phi'_{0\varepsilon }(x)\ge 0$, by maximum principle, we obtain
$ \psi_\varepsilon (x, t)\ge 0. $ |
On the other hand,
$ \psi_{\varepsilon t}-\frac{\partial }{\partial x}\big(\frac12\sigma _\varepsilon ^2 \psi_{\varepsilon x}\big)=-\frac{\partial }{\partial x}\big(\frac12\sigma _\varepsilon ^2 \psi_{\varepsilon }\big) $ | (27) |
We multiply equation (27) $\psi_\varepsilon ^p$ for any $p\ge 1$, to obtain
$ \frac{d}{dt}\frac{1}{p+1}\int_{-\infty}^{+\infty}\psi_\varepsilon ^{p+1}dx+p\int_{-\infty}^{+\infty}\frac12\sigma _\varepsilon ^2 \psi^2_{\varepsilon x}\psi_\varepsilon ^{p-1}dx\\ =p\int_{-\infty}^{+\infty}\frac12\sigma _\varepsilon ^2 \psi_\varepsilon ^p\psi_{\varepsilon x}dx\\ \le\delta \int_{-\infty}^{+\infty}p\cdot\frac12\sigma _\varepsilon ^2\psi_\varepsilon ^{p-1}\psi_{\varepsilon x}^2dx+C(\delta )p\int_{-\infty}^{+\infty}\frac12\sigma _\varepsilon ^2\psi_\varepsilon ^{p+1}dx. $ |
If we choose $\delta $ sufficiently small, we have
$ \frac{d}{dt}\frac{1}{p+1}\int_{-\infty}^{+\infty}\psi_\varepsilon ^{p+1}dx+\frac12 p\int_{-\infty}^{+\infty}\frac12\sigma _\varepsilon ^2\psi^2_{\varepsilon x}\psi_\varepsilon ^{p-1}dx \le C\int_{-\infty}^{+\infty}\frac12\sigma _\varepsilon ^2 \psi_\varepsilon ^{p+1}dx $ |
Since $0 < \sigma _0\le\sigma _\varepsilon \le\sigma _1 < +\infty$, we have
$ \int_{-\infty}^{+\infty}\psi_\varepsilon ^{p+1}dx+\int_{-\infty}^{+\infty}\int_{0}^{t}\psi^2_{\varepsilon x}\psi_\varepsilon ^{p-1}dxdt\\ \le\int_{-\infty}^{+\infty}\phi'_{0\varepsilon }(x)^2dx+C\int_{0}^{t}\int_{-\infty}^{+\infty}\psi_\varepsilon ^{p+1}dxdt $ |
Gronwall's inequality yields,
$ \int_{-\infty}^{+\infty}\psi_\varepsilon ^{p+1}dx+\int_{-\infty}^{+\infty}\int_{0}^{T}\psi^2_{\varepsilon x}\psi_\varepsilon ^{p-1}dxdt\le C, $ |
where $C$ depends only on known data and $p$, but not on $\varepsilon $.
Now from equation (27),
$ \psi_{\varepsilon t}-\frac{\partial }{\partial x}\big(\frac12\sigma _\varepsilon ^2 \psi_{\varepsilon x}\big)=-\frac{\partial }{\partial x}\big(\frac12\sigma _\varepsilon ^2 \psi_{\varepsilon }\big), \quad (x, t)\in Q_T\\ \psi_\varepsilon (x, 0)=\phi'_{0\varepsilon }(x), \quad \quad -\infty<x<+\infty $ |
where $\sigma _\varepsilon ^2 \psi_{\varepsilon }\in L^p(\mathbb{R}^1)$ for any $p>1$. Then $|| \psi_\varepsilon || _{L^{\infty}(Q_T)}\le C$. To derive the $W^{2, 1}_2(Q_T)-$estimate, we multiply equation (24) by $\phi_{\varepsilon xx}$ and integrate over $Q_t$ to obtain
$ \sigma _0\iint\limits_{Q_t}\phi_{\varepsilon xx}^2dxdt+\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2dx-\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2(x, 0)dx\\ =\iint\limits_{Q_t}\frac12\sigma _\varepsilon ^2\phi_{\varepsilon x}\phi_{\varepsilon xx}dxdt\\ \le \delta \iint\limits_{Q_t}\phi_{\varepsilon xx}^2dxdt+C(\delta )\iint\limits_{Q_t}\phi_{\varepsilon x}^2dxdt $ |
By choosing $\delta $ sufficiently small, we have
$ \sigma _0\iint\limits_{Q_t}\phi_{\varepsilon xx}^2dxdt+\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2(x, t)dx \le\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2(x, 0)dx+C\iint\limits_{Q_t}\phi_{\varepsilon x}^2dxdt $ |
Gronwall's inequality yields
$ \sup\limits_{0\le t\le T}\int_{-\infty}^{+\infty}\phi_{\varepsilon x}^2dx+\iint\limits_{Q_t}\phi_{\varepsilon xx}^2dxdt\le C. $ |
From the equation (24), we see
$ \iint\limits_{Q_t}\phi_{\varepsilon t}^2dxdt\le C $ |
Lemma 4.2. The solution $\phi_\varepsilon (x, t)\in C^{1+\alpha, \frac{1+\alpha}{2}}({Q}_{T})$ uniformly with respect to $\varepsilon $, for some $\alpha\in (0, 1)$. Moreover, for any $\rho>0$ there exists a constant $C(\rho)$ such that
$ ||\psi_{\varepsilon }||_{C^{\alpha, \frac{\alpha}{2}}(\bar{Q}_{T}\backslash B_{\rho}(0, 0))}\leq C(\rho), $ |
where $B_{\rho}(0, 0)=\{(x, t): -\rho < x < \rho, 0 < t < \rho\}$ and $C(\rho)$ depends only on known data and $\rho$.
Proof. Indeed, we see $\psi_\varepsilon (x, t)$ is a solution of the following problem:
$ \psi_{\varepsilon t}-\frac{\partial }{\partial x}\big[\frac12\sigma ^2\big(\psi_{\varepsilon x}-\psi_{\varepsilon }\big)\big]=0, \quad (x, t)\in {Q}_T, \\ \psi_{\varepsilon }(x, 0)=\phi'_{\varepsilon 0}(x)\ge 0, \quad -\infty<x<+\infty. $ |
By the standard theory for parabolic equations, we see that $\psi_{\varepsilon }(x, t)$ is Hölder continuous in $Q_{T}$. Moreover, the Hölder continuity of $\psi_\varepsilon $ holds up to the initial time $t=0$ except a neighborhood of the singular point $(0, 0)$. Hence, for any $\rho>0$, we have
$ ||\psi_{\varepsilon }||_{C^{\alpha, \frac{\alpha}{2}}(\bar{Q}_{T}\backslash B_{\rho}(0, 0))}\leq C(\rho), $ |
where $B_{\rho}(0, 0)$ defined above and $C(\rho)$ depends only on known data and $\rho$, but not on $\varepsilon $.
Theorem 4.1. The problem (17), (21) admits a solution $\phi(x, t)\in W^{1, 0}_\infty(Q_T)\cap W^{2, 1}_2(Q_T)\cap C^{\alpha, \frac{\alpha}{2}}(\overline{Q}_T)$, for any $\alpha\in(0, 1)$. Furthermore, the solution satisfies
a) $0\le\phi(x, t)\le 1$, $0\le\phi_{x}(x, t)\le C_1$,
b) $\int_{-\infty}^{+\infty}\phi_{x}^2dx+\iint\limits_{Q_T}(\phi_{xx}^2+\phi_{t}^2) dxdt\le C_2$,
where $C_1$ and $C_2$ are constants which depend only on known data.
Proof. By the standard compactness argument, we see that there exists a function $\phi(x, t)\in C^{\alpha, \frac\alpha2}(\overline{Q}_T)\cap W^{2, 1}_2(Q_T)$ such that
$ \phi_{\varepsilon }(x, t)\to \phi(x, t) \quad \mbox{uniformly in } C^{\alpha, \frac\alpha2}(\overline{Q}_T), \\ \phi_{\varepsilon }(x, t)\to \phi(x, t) \quad \mbox{weakly in } W^{2, 1}_2(Q_T). $ |
Moreover, $0\le\phi(x, t)\le 1$, $0\le\phi_{x}(x, t)\le C_1$. For the limit function $\phi$, and $(x_0, t_0)$ with $(x_0, t_0)\notin\Gamma_T$ and $t_0>0$, where $\Gamma_T$ is defined in (19), there exists a small neighborhood $P_\delta (x_0, t_0)=\{(x, t)|(x-x_0)^2+(t-t_0)^2 < \delta \}$ such that $P_\delta (x_0, t_0)\cap\Gamma_T=\emptyset$. In $P_\delta (x_0, t_0)$, $\sigma _\varepsilon (x, t, \phi)=\sigma (x, t, \phi).$ By regularity theory of parabolic equation, $\phi(x, t)\in C^{2, 1}(P_\delta (x_0, t_0))$, and
$ \phi_t-\frac12\sigma ^2(\phi_{xx}-\phi_{x})=0. $ |
It follows that for any $(x, t)\notin\Gamma_T$
$ \phi_t-\frac12\sigma ^2(\phi_{xx}-\phi_{x})=0, \quad \mbox{in the classical sense in $Q_T\backslash\Gamma_T$}, \\ \phi(x, 0)=\phi_0(x), \quad -\infty<x<+\infty. $ |
Notice also by Lemma 4.2, $\phi, \phi_x$ are continuous across $\Gamma_T.$ That is, $\phi(x, t)$ is a solution in $Q_T$.
We can obtain a better regularity for the strong solution obtained in Theorem 4.1.
Theorem 4.2. The solution $\phi(x, t)\in C^{1+\alpha, \frac{1+\alpha}{2}}({Q}_{T})$ for some $\alpha\in (0, 1)$. Moreover, $\phi_x(x, t)$ is Hölder continuous up to the initial time $t=0$ except a neighborhood of $(0, 0)$.
Proof. Indeed, we see $\psi(x, t)=\phi_x(x, t)$ is a weak solution of the following problem:
$ \psi_{ t}=\frac{\partial }{\partial x}\big[\frac12\sigma ^2\big(\psi_{ x}-\psi\big)\big], \quad (x, t)\in {Q}_T, \\ \psi(x, 0)=\phi'_{0}(x)\ge 0, \quad -\infty<x<+\infty. $ |
By the standard regularity theory of parabolic equations, we see that $\psi(x, t)$ is Hölder continuous in $Q_{T}$. Moreover, the Hölder continuity of $\psi$ holds up to the initial time $t=0$ except a neighborhood of the singular point $(0, 0)$.
In this section we use a novel approach to prove the uniqueness. This method is also suitable for a general n-dimensional problem.
Theorem 4.3. The solution obtained in Theorem 4.1 is unique.
Proof. We rewrite Eq.(18) as follows:
$ k(\phi, x, t)\phi_t=\phi_{xx}-\phi_x, $ | (28) |
where
$k(\phi, x, t)=\Big\{2/ˆσ2L(t), if ϕ>γex,2/ˆσ2H(t), if ϕ<γex. $
|
Define
$
e(ϕ,x,t)=(1/ˆσ2L(t)+1/ˆσ2H(t))(ϕ−γex)+(1/ˆσ2L(t)−1/ˆσ2H(t))|ϕ−γex|={2(ϕ−γex)/ˆσ2L(t), if ϕ>γex,2(ϕ−γex)/ˆσ2H(t), if ϕ<γex.
$
|
(29) |
which satisfies that
$
k(\phi, x, t)= \frac{\partial e}{\partial
\phi}=\Big\{2/ˆσ2L(t), if ϕ>γex,2/ˆσ2H(t), if ϕ<γex.
$
|
(30) |
Therefore,
$ \frac{\partial }{\partial t}[e(\phi, x, t)]-[\phi_{xx}-\phi_x]=\frac{\partial e}{\partial t}\equiv f(\phi, x, t), $ | (31) |
where
$
f(ϕ,x,t):=(1/ˆσ2L(t)+1/ˆσ2H(t))′(ϕ−γex)+(1/ˆσ2L(t)−1/ˆσ2H(t))′|ϕ−γex|.
$
|
We also have
$ e(\phi, x, t)|_{t=0}=e_0(x), $ | (32) |
where
$
e0(x)=(1/ˆσ2L(0)+1/ˆσ2H(0))(ϕ0(x)−γex)+(1/ˆσ2L(0)−1/ˆσ2H(0))|ϕ0(x)−γex|,
$
|
To prove the uniqueness, we assume that there are two solutions $\phi_i(x, t), \ i=1, 2$ of the problem (17), (21), which satisfy the estimates in Theorem 4.1. Set
$ \phi(x, t)=\phi_1(x, t)-\phi_2(x, t), \hspace{1cm}(x, t)\in Q_T. $ |
Now we multiply $\int_t^\tau\phi(x, s) ds$ in the both sides of (31), where $\tau\in [0, T]$, and then integrate it in $\mathbb{R}$, we have
$
∂∂t∫R[(e(ϕ1,x,t)−e(ϕ2,x,t))∫τtϕds]dx+∫R[(e(ϕ1,x,t)−e(ϕ2,x,t))ϕ]dx=∫R[(ϕxx−ϕx)∫τtϕds]dx+∫R[(f(ϕ1,x,t)−f(ϕ2,x,t))∫τtϕds]dx:=A1(x,t,τ)+A2(x,t,τ),
$
|
where
$
A1(t,τ)=∫R[(ϕxx−ϕx)∫τtϕds]dx,A2(t,τ)=∫R[(f(ϕ1,x,t)−f(ϕ2,x,t))∫τtϕds]dx.
$
|
Now,
$
A1(t,τ)=−∫R[(ϕx−ϕ)∫τtϕxds]dx=12∂∂t∫R[∫τtϕxds]2dx+∫R[ϕ∫τtϕxds]dx.
$
|
$ A_2(t, \tau) = \int_{\mathbb{R}}\Big[\Big(f_\phi(\xi, t)\phi\Big)\int_t^\tau\phi ds\Big]dx, $ |
where $\xi$ is the mean-value between $\phi_1$ and $\phi_2$.
We take integration to the above equation with respect to $t$ from $0$ to $\tau$ to obtain:
$
∫τ0∫R[e(ϕ1,x,t)−e(ϕ2,x,t)]ϕdxdt+12∫R[∫τ0ϕxds]2dx=∫τ0∫R+[ϕ∫τtϕxds]dxdt+∫τ0∫R(fϕ(ξ,x,t)ϕ)∫τtϕdsdxdt≤2ε∫τ0∫Rϕ2dxdt+C(ε)∫τ0∫R([∫τ0ϕds]2+[∫τ0ϕxds]2)dxdt.
$
|
Here at the final step we have used the uniform boundedness of $f_{u}(u, t)$ and Cauchy-Schwarz's inequality with small parameter $\varepsilon $.
Notice (29), we have
$[e(\phi_1, x, t)-e(\phi_2, x, t)]\ge c_0(\phi_1-\phi_2)=(2/\sigma _1^2)\phi, $ |
where $\sigma _1$ defined in (21) has a positive lower bound. We choose $\varepsilon$ sufficiently small to obtain
$ \int_0^\tau \int_{\mathbb{R}}\phi^2 dxdt+\int_{\mathbb{R}}\Big[\int_{0}^{\tau}\phi _xds\Big]^2 dx\\ \leq C\Big[\int_0^\tau [\int_0^t \int_{\mathbb{R}}\phi^2 dxds]dt+\int_{0}^{\tau}\int_{\mathbb{R}} \Big[\int_0^\tau\phi_x ds\Big]^2 dxdt\Big]. $ |
Consequently, Gronwall's inequality yields $\phi(x, t)\equiv 0$ on $Q_{\tau}$. As $\tau$ is arbitrary over $[0, T]$, we have $\phi(x, t)\equiv 0$ on $Q_{T}$.
Thin concludes the proof of uniqueness.
In this section, we show that $\Gamma_T$ is the graph of a smooth free boundary $x=s(t)$.
Lemma 5.1. There exists a number $T_0>0$, such that $\Gamma_{T_0}$ defined in (19) is the graph of a function $x=s(t), 0\le t\le T_0$ with $s(0)=s_0=\ln \frac{1}{\gamma}$.
Proof. Note that $\phi_0(x)=\min\{e^x, 1\}$, then $s_0:=s(0)=\ln 1/\gamma >0$, it follows that $\phi'_0(x)=0$, for $x>0$; $\phi'_0(x)=e^x$, for $x < 0$. It follows that if $x>0$,
$ \phi_{0}^{'}(x)-\gamma e^x=-\gamma e^x<0. $ |
Since $\phi_x(x, t)$ is Hölder continuous except $(0, 0)$, there exists small numbers $t_0>0$ and $\delta_0>0$ such that
$ \phi_x(x, t)-\gamma e^x\le m_0<0, \quad s_0-\delta_0<x<s_0+\delta_0, \quad 0\le t\le t_0, $ |
where $m_0$ depends on $t_0$ and $\delta_0$.
By the implicit function theorem, there exists a function, denoted by $x=s(t)$ solves $\phi(s(t), t)=\gamma e^{s(t)}$ such that
$ \Gamma_{t_0}=\{(x, t)|x=s(t), 0\le t\le t_0\}. $ |
It will be seen in the next theorem that we can extend $t_0$ as long as $\phi_x^-(s(t), t)-\gamma e^{s(t)} < 0$.
Theorem 5.1. There exists a function $x=s(t) \in C^{\beta}[0, T]$ with $\beta\in (\frac{1}{2}, 1)$ such that
$ \Gamma_T = \{(s(t), t): 0\leq t\leq T\}. $ |
Moreover, there exists a constant $\varepsilon _0>0$ such that
$ \phi_{x}(s(t), t)-\gamma e^{s(t)}\leq -\varepsilon _0<0, 0\leq t\leq T. $ |
Proof. From Lemma 5.1, we know that there exists at least a small internal $[0, T_{0}]$ such that
$ \Gamma_{T_{0}}=\{(s(t), t): 0\leq t\leq T_0.\} $ |
Moreover,
$ \phi_x(s(t), t)-\gamma e^{s(t)}<0, t\in [0, T_{0}]. $ |
Since $F(x, t):=\phi(x, t)-\gamma e^x$ is differentiable in $Q_{T}$ with respect to $x$, by the implicit function theorem, the set
$ \Gamma_{t_0}=\{(x, t); F(x, t)=0, 0\leq t\leq t_0\} $ |
can be extended as the graph of a curve $x=s(t)$ as long as
$ F_x(x, t)<0, \mbox{for all $(x, t)\in \Gamma_{T}$}. $ |
Moreover, $x=s(t)$ is decreasing on $[0, T]$.
Define
$ T^*=\sup\{ T_{0}: \phi_x(s(t), t)-\gamma e^{s(t)}<0, 0\leq t\leq T_{0}\}. $ |
We claim $T^*=T$. Assume the contrary with $T^* < T$. It is clear by the definition of $T^*$ and the continuity of $\phi_x(x, t)$ that
$ \lim\limits_{t\to T^{*}}(\phi_x^-(s(t), t)-\gamma e^{s(t)})=0 . $ |
We define $W(x, t)=\phi(x, t)-\gamma e^{x}$, then
$ W_t=\phi_t, \quad W_x=\phi_x-\gamma e^{x}, \quad W_{xx}=\phi_{xx}-\gamma e^{x} $ |
Therefore, $W(x, t)$ satisfies
$ W_t-\frac12\sigma _{H}^2(W_{xx}-W_x)= 0, \, -\infty<x<s(t), 0<t<T^*, \\ W(s(t), t)=0, \hspace{1cm} 0<t<T^*, \\ W(x, 0)=W_0(x), \hspace{1cm}-\infty<x<s_0, $ |
where $W_0(x):=\phi_0(x)-\gamma e^x\geq 0$ on $(-\infty, s_0]$.
Since $s(t)$ is decreasing on $[0, T^*)$, Hopf's lemma holds at $(s(T^*), T^*)$. Hence, by Hopf's lemma, there exists a number $\varepsilon _0>0$ such that $W_x(s(T^*), T^*)\leq -\varepsilon _0 < 0$, which contradicts with the definition of $T^*$. It follows that $T=T^*$. Consequently, there exists a curve, denoted by $x=s(t)$, such that
$ \Gamma_{T}=\{(s(t), t): 0\leq t\leq T.\}. $ |
Moreover, there exists a constant $\varepsilon _0>0$ such that
$ \phi_{x}(s(t), t)-\gamma e^{s(t)}\leq -\varepsilon _0<0, 0\leq t\leq T. $ | (33) |
To prove the Hölder continuity of $s(t)$ on $[0, T]$, we take any $t_1 < t_2$ in $[0, T]$ with $|t_1-t_2|$ sufficiently small,
$ \phi(s(t_2), t_2)-\phi(s(t_1), t_1)=\gamma (e^{s(t_{2})}-e^{s(t_{1})}). $ |
Then,
$ \phi(s(t_2), t_1)-\phi(s(t_1), t_1)-\gamma (e^{s(t_{2})}-e^{s(t_{1})})\\ = \phi(s(t_2), t_1)-\phi(s(t_2), t_2). $ |
We use the mean-value theorem to obtain
$ |s(t_2)-s(t_1)|\leq \frac{|\phi(s(t_{2}), t_{2})-\phi(s(t_{2}), t_{1})|}{|F_{x}(\theta, t_{1})|}, $ |
where $\theta$ is the mean-value between $s(t_{1})$ and $s(t_{2})$.
Since $|F_x(x, t)|$ is bounded and $|F_x(x, t)|$ has a positive lower bound in a neighborhood of $x=s(t), 0\leq t\leq T$, we find that there exists a constant $C$ such that
$ |s(t_2)-s(t_1)|\leq C|t_1-t_2|^{\beta}, $ |
where $\beta=\frac{1+\alpha}{2}$.
It follows that $s(t)$ is Hölder continuous on $[0, T]$ with a Hölder exponent $\beta\in (\frac{1}{2}, 1)$.
Theorem 5.2. The free boundary $s(t)$ satisfies
$ s'(t)=-\frac{ \sigma _{H}(t)(\phi_{xx}(s(t)-, t)-\phi_x(s(t)-, t))}{\phi_{x}(s(t), t)-\gamma e^{s(t)}}, 0\leq t\leq T. $ |
Moreover, $s(t)\in C^{\infty}[0, T]$, provided that $\sigma_{L}(t)$ and $\sigma_{H}(t)$ are smooth on $[0, T]$.
Proof. By Theorem 5.1, we can define $Q_{T}^{L}$ and $Q_{T}^{H}$ as in Section 3. From the theory of parabolic equations (Ladyženskaja et al., 1968), we know that $ \phi(x, t)$ is smooth in $Q_{T}^L$ and $Q_{T}^H$, respectively. Moreover, $\phi_x(x, t)$ is continuous up the boundary $x=s(t)$ from the left-hand side and the right-hand side on $[0, T]$.
To prove the further regularity of $s(t)$, we note that $\psi(x, t)=\phi_x(x, t)$ is a classical solution of the following free boundary problem
$ \psi_t=\frac{\partial }{\partial x}[\frac{1}{2}\sigma ^2 (\psi_{xx}-\psi)], \hspace{1cm} (x, t)\in Q_{T}\backslash \Gamma_{T}, \\ \psi(s(t)-, t)=\psi(s(t)+, t), \hspace{1cm} 0<t<T, \\ \sigma _{H}(\psi_x(s(t)-, t)-\psi(s(t), t))=\sigma _{L}(\psi_x(s(t)+, t)-\psi(s(t)+, t)), 0\leq t\leq T, \\ s'(t)=-\frac{\sigma _{H}(t) (\psi_{x}(s(t)-, t)-\psi(s(t), t))}{\psi(s(t), t)-\gamma e^{s(t)}}, \hspace{1cm} 0\leq t\leq T, \\ \psi(x, 0)=\phi_{0}^{'}(x), \hspace{1cm} -\infty<x<\infty. $ |
We can apply the classical technique for the Musket problem Evans (1978) to obtain $s(t)\in C^{\infty}[0, \infty]$.
Theorem 5.3. $s(t)$ is non-increasing on $[0.T]$.
Proof. We go back the approximate problem (24)-(25). Note that $ \phi_{0}^{'}(x)-\phi_{0}(x)=0 $ if $x < 0$ and $\phi_{0}^{'}(x)-\phi_{0}(x)=-1$ if $x>0$. We can construct a smooth approximation for $\phi_0(x)$ such that
$ \phi_{\varepsilon 0}^{''}-\phi_{0}^{'}(x)<0, -\infty<x<\infty. $ |
Let $w(x, t)=\phi_{\varepsilon t}(x, t)$. Then $w(x, t)$ satisfies
$ w_t=\frac{1}{2}\sigma ^2[w_{xx}-w_{x}]+a(x, t)w, (x, t)\in Q_T. $ |
where
$ a(x, t):=\frac{ (\sigma ^{2}(t, \phi))_{t} }{\sigma ^{2}(t, \phi)}. $ |
Moreover,
$ w(x, 0)=\frac{1}{2}\sigma _{\varepsilon }(0, \phi_0)^2[\phi_{0\varepsilon }^{''}(x)-\phi_{0\varepsilon }^{'}(x)]<0, -\infty<x<\infty. $ |
By the maximum principle, we see $ w(x, t)\leq 0$ on $Q_T$. After passing the limit, we see
$ \phi_{xx}(x, t)-\phi_x(x, t)\leq 0, (x, t)\in Q_{T}, $ |
in the $L^2$-sense. On the other hand,
$ \phi_x(s(t), t)-\gamma e^{s(t)}<0, 0\leq t\leq T, $ |
we see from the representation of $s(t)$ that
$ s'(t)\leq 0, \hspace{1cm} t\in [0, T]. $ |
The numerical solution of the Eq.(17), (21) can be calculated by using a standard PDE method. However, we are more interested in the. For convenience, the free boundary is denoted by $s(r, t)$, which depends on the interest rate $r$.
In order to draw the graph $s(r, t)$, we use the explicit difference scheme (Jiang, 2005), which is equivalent to the well-known the Binomial-Tree method. The steps of the computation are as follows:
1. Discrete $t\in(0, T)$ by $t_i=T -i\Delta t$, where $T=t_0, ..., t_I=0$, $I=10, $ $\Delta t =\frac{t_I-t_0}{I}; $ $r\in(0.01, 0.04)$ by $r_k=0.01 + k\Delta r$, where $0.01=r_0, ..., r_K=0.04$, $\Delta r =\frac{r_K-r_0}{K}; $ in space using a mesh with difference $\Delta x$: $0=x_0, ..., x_j=2$. And at every point $(t_i, x_j, r_k)$, denote the function $\phi(t_i, x_j, r_k)$ by $\phi_{i, j, k}$.
2. For $t_0=T$, $\phi_{0, j, k}=\min(x_j, F)$, then solve $\phi_{i, j, k}$ by finite difference scheme:
$
ϕi,j,k=ϕi−1,j,k+σ2Δt2Δx2x2j[ϕi−1,j−1,k+ϕi−1,j+1,k−2ϕi−1,j,k]+σ2rΔt2Δr2r2k[ϕi−1,j,k−1+ϕi−1,j,k+1−2ϕi−1,j,k]+rkxjΔt2Δx[ϕi−1,j+1,k−ϕi−1,j−1,k]+a(θ−rk)Δt2Δr[ϕi−1,j,k+1−ϕi−1,j,k−1]−rkϕi−1,j,kΔt+σσrΔtρxj4ΔrΔx[ϕi−1,j+1,k+11−ϕi−1,j+1,k−1−ϕi−1,j−1,k+1+ϕi−1,j−1,k−1],
$
|
where
$
\sigma = \left\{σL, if ϕi−1,j,k>γxk,σH, otherwise, \right.
$
|
Denote the joint points to be $s_{j, k}$, which represents the approximated credit rating migration joints.
3. Repeat the above process.
4. Draw the graph based on $\{s_{j, k}\}$ in space $r.S, t.$
Now we have the 2-d free boundary interface.
The numerical solution of the free boundary is shown in Figure 1, where the parameters and the ranges of the variables are chosen as follows
$ a=1, \, \theta = 0.03, \, F=1, \, \gamma=0.8, \, \\ \sigma_L=0.4, \, \sigma_H=0.2, \, \sigma_r=0.3, \, \rho=0.5, T=5, \\ r\in(0.01, 0.04), \quad \ t\in(0, 5). $ |
From the graphs, we see that a free boundary $y=s(r, t)$ divides the region into two parts: high and low rating regions. Moreover, the free boundary is decreasing with respect to $r$ and increasing with respect to $t$ as expected (see Theorem 4.1 and Theorem 5.2), where $t$ is the original life time of the bond. From the graph Fig. 1, we also see that the change of free boundary is more sensitive with respect to the interest rate than that of time.
In this paper, we derived a new model for pricing corporate bonds with credit-rating migration and a stochastic interest rate. The new model consists of a coupled system of partial differential equations. By employing a dimension-deduction technique, we obtain a new free boundary problem. By using various analysis techniques for partial differential equations, we established the existence, uniqueness and regularities of the solution under certain assumptions. Moreover, some properties of the solution are derived. A $C^\infty$-regularity of the free boundary is proved. Their financial implications are explained. The model is new to our knowledge and is a big improvement for the current model. The analysis is delicate. Finally, numerical simulations are also presented.
This work is supported by National Natural Science Foundation of China (No.11671301).The authors would like to thank Professor Bei Hu for some discussions during the preparation of the paper. Many thanks also go to anonymous referees for their comments and suggestions. Chen and Yin would like to thank the hospitalities of Tongji University during their visiting where the part of this work was completed.
All authors declare no conflicts of interest in this paper.
[1] |
Reza MT, Lynam JG, Vasquez VR, et al. (2012) Pelletization of Biochar from Hydrothermally Carbonized Wood. Environ Prog Sustain 31: 225-234. doi: 10.1002/ep.11615
![]() |
[2] |
Reza MT, Lynam JG, Uddin MH, et al. (2013) Hydrothermal carbonization: Fate of inorganics. Biomass Bioenerg 49: 86-94. doi: 10.1016/j.biombioe.2012.12.004
![]() |
[3] | Reza MT, Borrego AG, Wirth B (2014) Optical texture of hydrochar from maize silage and maize silage digestate. Int J of Coal Geology 134–135: 74-79. |
[4] |
Reza MT, Becker W, Sachsenheimer K, et al. (2014) Hydrothermal carbonization (HTC): Near infrared spectroscopy and partial least-squares regression for determination of selective components in HTC solid and liquid products derived from maize silage. Bioresource Technol 161: 91-101. doi: 10.1016/j.biortech.2014.03.008
![]() |
[5] | Lynam J, Reza MT, Yan W, et al. (2014) Hydrothermal carbonization of various lignocellulosic biomass. Biomass Conv Bioref : 1-9. |
[6] | Coronella C, Lynam J, Reza MT, et al. (2014) Hydrothermal Carbonization of Lignocellulosic Biomass. In: Jin F, editor. Application of Hydrothermal Reactions to Biomass Conversion: Springer Berlin Heidelberg: 275-311. |
[7] |
Liu ZG, Quek A, Parshetti G, et al. (2013) A study of nitrogen conversion and polycyclic aromatic hydrocarbon (PAH) emissions during hydrochar-lignite co-pyrolysis. Appl Energ 108: 74-81. doi: 10.1016/j.apenergy.2013.03.012
![]() |
[8] | Reza MT, Andert J, Wirth B, et al. (2014) Hydrothermal Carbonization of Biomass for Energy and Crop Production. Applied Bioenergy 1: 11-29. |
[9] |
Demir-Cakan R, Makowski P, Antonietti M, et al. (2010) Hydrothermal synthesis of imidazole functionalized carbon spheres and their application in catalysis. Catal Today 150: 115-118. doi: 10.1016/j.cattod.2009.05.003
![]() |
[10] |
Bandura AV, Lvov SN (2006) The ionization constant of water over wide ranges of temperature and density. J Phys Chem Ref Data 35: 15-30. doi: 10.1063/1.1928231
![]() |
[11] |
Reza MT, Yan W, Uddin MH, et al. (2013) Reaction kinetics of hydrothermal carbonization of loblolly pine. Bioresource Technol 139: 161-169. doi: 10.1016/j.biortech.2013.04.028
![]() |
[12] |
Reza MT, Uddin MH, Lynam J, et al. (2014) Hydrothermal carbonization of loblolly pine: reaction chemistry and water balance. Biomass Conv Bioref 4: 311-321. doi: 10.1007/s13399-014-0115-9
![]() |
[13] |
Wiedner K, Naisse C, Rumpel C, et al. (2013) Chemical modification of biomass residues during hydrothermal carbonization - What makes the difference, temperature or feedstock? Org Geochem 54: 91-100. doi: 10.1016/j.orggeochem.2012.10.006
![]() |
[14] |
Funke A (2015) Fate of Plant Available Nutrients during Hydrothermal Carbonization of Digestate. Chemie Ingenieur Technik 87: 1713-1719. doi: 10.1002/cite.201400182
![]() |
[15] |
Heilmann SM, Molde JS, Timler JG, et al. (2014) Phosphorus Reclamation through Hydrothermal Carbonization of Animal Manures. Environ Sci Technol 48: 10323-10329. doi: 10.1021/es501872k
![]() |
[16] | Funke A, Mumme J, Koon M, et al. (2013) Cascaded production of biogas and hydrochar from wheat straw: Energetic potential and recovery of carbon and plant nutrients. Biomass Bioenerg 58: 229-237. |
[17] |
Heilmann SM, Davis HT, Jader LR, et al. (2010) Hydrothermal carbonization of microalgae. Biomass Bioenerg 34: 875-882. doi: 10.1016/j.biombioe.2010.01.032
![]() |
[18] |
Sevilla M, Fuertes AB (2009) The production of carbon materials by hydrothermal carbonization of cellulose. Carbon 47: 2281-2289. doi: 10.1016/j.carbon.2009.04.026
![]() |
[19] |
Funke A, Ziegler F (2010) Hydrothermal carbonization of biomass: A summary and discussion of chemical mechanisms for process engineering. Biofuel Bioprod Bior 4: 160-1677. doi: 10.1002/bbb.198
![]() |
[20] |
Kruse A, Funke A, Titirici MM (2013) Hydrothermal conversion of biomass to fuels and energetic materials. Curr Opin Chem Biol 17: 515-521. doi: 10.1016/j.cbpa.2013.05.004
![]() |
[21] | Zhu XD, Liu YC, Luo G, et al. (2014) Facile Fabrication of Magnetic Carbon Composites from Hydrochar via Simultaneous Activation and Magnetization for Triclosan Adsorption. Environ Sci Technol 48: 5840-5848. |
[22] |
Kumar S, Loganathan VA, Gupta RB, et al. (2011) An Assessment of U(VI) removal from groundwater using biochar produced from hydrothermal carbonization. J Environ Manage 92: 2504-2512. doi: 10.1016/j.jenvman.2011.05.013
![]() |
[23] |
Regmi P, Moscoso JLG, Kumar S, et al. (2012) Removal of copper and cadmium from aqueous solution using switchgrass biochar produced via hydrothermal carbonization process. J Environ Manage 109: 61-69. doi: 10.1016/j.jenvman.2012.04.047
![]() |
[24] |
Hu B, Wang K, Wu LH, et al. (2010) Engineering Carbon Materials from the Hydrothermal Carbonization Process of Biomass. Adv Mater 22: 813-828. doi: 10.1002/adma.200902812
![]() |
[25] |
Demir-Cakan R, Baccile N, Antonietti M, et al. (2009) Carboxylate-Rich Carbonaceous Materials via One-Step Hydrothermal Carbonization of Glucose in the Presence of Acrylic Acid. Chem Mater 21: 484-490. doi: 10.1021/cm802141h
![]() |
[26] | Baccile N, Laurent G, Babonneau F, et al. (2009) Structural Characterization of Hydrothermal Carbon Spheres by Advanced Solid-State MAS C-13 NMR Investigations. J Phys Chem C 113: 9644-9654. |
[27] |
Baccile N, Antonietti M, Titirici MM (2010) One-Step Hydrothermal Synthesis of Nitrogen-Doped Nanocarbons: Albumine Directing the Carbonization of Glucose. Chemsuschem 3: 246-253. doi: 10.1002/cssc.200900124
![]() |
[28] |
Zhao L, Bacsik Z, Hedin N, et al. (2010) Carbon Dioxide Capture on Amine-Rich Carbonaceous Materials Derived from Glucose. Chemsuschem 3: 840-845. doi: 10.1002/cssc.201000044
![]() |
[29] |
Fechler N, Wohlgemuth SA, Jaker P, et al. (2013) Salt and sugar: direct synthesis of high surface area carbon materials at low temperatures via hydrothermal carbonization of glucose under hypersaline conditions. J Mater Chem A 1: 9418-9421. doi: 10.1039/c3ta10674h
![]() |
[30] |
Reza MT, Rottler E, Tolle R, et al. (2015) Production, characterization, and biogas application of magnetic hydrochar from cellulose. Bioresource Technol 186: 34-43. doi: 10.1016/j.biortech.2015.03.044
![]() |
[31] |
Sevilla M, Fuertes AB (2009) Chemical and Structural Properties of Carbonaceous Products Obtained by Hydrothermal Carbonization of Saccharides. Chem-Eur J 15: 4195-4203. doi: 10.1002/chem.200802097
![]() |
[32] |
Titirici MM, Antonietti M (2010) Chemistry and materials options of sustainable carbon materials made by hydrothermal carbonization. Chem Soc Rev 39: 103-116. doi: 10.1039/B819318P
![]() |
[33] |
Sevilla M, Fuertes AB (2009) Chemical and structural properties of carbonaceous products obtained by hydrothermal carbonization of saccharides. Chemistry 15: 4195-4203. doi: 10.1002/chem.200802097
![]() |
[34] |
Cao XD, Ma LN, Gao B, et al. (2009) Dairy-Manure Derived Biochar Effectively Sorbs Lead and Atrazine. Environ Sci Technol 43: 3285-3291. doi: 10.1021/es803092k
![]() |
[35] |
Diakite M, Paul A, Jager C, et al. (2013) Chemical and morphological changes in hydrochars derived from microcrystalline cellulose and investigated by chromatographic, spectroscopic and adsorption techniques. Bioresource Technol 150: 98-105. doi: 10.1016/j.biortech.2013.09.129
![]() |
[36] | Si Y, Ren T, Li Y, et al. (2012) Fabrication of magnetic polybenzoxazine-based carbon nanofibers with Fe3O4 inclusions with a hierarchical porous structure for water treatment. Carbon 50: 5176-5185. |
[37] |
Lynam JG, Reza MT, Vasquez VR, et al. (2012) Effect of salt addition on hydrothermal carbonization of lignocellulosic biomass. Fuel 99: 271-273. doi: 10.1016/j.fuel.2012.04.035
![]() |
[38] |
Ryu J, Suh YW, Suh DJ, et al. (2010) Hydrothermal preparation of carbon microspheres from mono-saccharides and phenolic compounds. Carbon 48: 1990-1998. doi: 10.1016/j.carbon.2010.02.006
![]() |
[39] | Peterson AA, Vogel F, Lachance RP, et al. (2008) Thermochemical biofuel production in hydrothermal media: A review of sub- and supercritical water technologies. Energy Env Sci 1: 32-65. |
[40] | Danso-Boateng E, Shama G, Wheatley AD, et al. (2008) Hydrothermal carbonisation of sewage sludge: Effect of process conditions on product characteristics and methane production. Bioresource technology 177: 318-327. |
1. | Zhehao Huang, Tianpei Jiang, Zhenzhen Wang, On a multiple credit rating migration model with stochastic interest rate, 2020, 43, 0170-4214, 7106, 10.1002/mma.6435 | |
2. | Hong-Ming Yin, Jin Liang, Yuan Wu, On a New Corporate Bond Pricing Model with Potential Credit Rating Change and Stochastic Interest Rate, 2018, 11, 1911-8074, 87, 10.3390/jrfm11040087 | |
3. | Shilong Li, Xia Zhao, Chuancun Yin, Zhiyue Huang, Stochastic interest model driven by compound Poisson process and Brownian motion with applications in life contingencies, 2018, 2, 2573-0134, 246, 10.3934/QFE.2018.1.246 | |
4. | Juri Kandilarov, Lubin Vulkov, 2019, Chapter 47, 978-3-030-10691-1, 416, 10.1007/978-3-030-10692-8_47 | |
5. | Shilong Li, Xia Zhao, Chuancun Yin, Zhiyue Huang, Stochastic interest model driven by compound Poisson process and Brownian motion with applications in life contingencies, 2018, 2, 2573-0134, 731, 10.3934/QFE.2018.1.731 | |
6. | Zhehao Huang, Zhenghui Li, Zhenzhen Wang, Utility Indifference Valuation for Defaultable Corporate Bond with Credit Rating Migration, 2020, 8, 2227-7390, 2033, 10.3390/math8112033 | |
7. | Stephanie Heck, Corporate bond yields and returns: a survey, 2022, 36, 1934-4554, 179, 10.1007/s11408-021-00394-4 | |
8. | Jin Liang, Yang Lin, A traveling wave with a buffer zone for asymptotic behavior of an asymmetric fixed credit migration model, 2022, 0170-4214, 10.1002/mma.8974 | |
9. | Zhongdi Cen, Jian Huang, Aimin Xu, Anbo Le, An iterative method for solving a PDE with free boundary arising from pricing corporate bond with credit rating migration, 2023, 8, 2473-6988, 3286, 10.3934/math.2023169 | |
10. | Xinfu Chen, Jin Liang, A Free Boundary Problem for Corporate Bond Pricing and Credit Rating Under Different Upgrade and Downgrade Thresholds, 2021, 12, 1945-497X, 941, 10.1137/20M1343592 | |
11. | Chen Xinfu, Hu Bei, Liang Jin, Yin Hongming, The free boundary problem for measuring credit rating migration risks, 2023, 1674-7216, 10.1360/SSM-2022-0151 | |
12. | Jin Liang, Yang Lin, A Bond Pricing Model with Credit Migration Risk: Different Upgrade and Downgrade Thresholds, 2023, 39, 0168-9673, 765, 10.1007/s10255-023-1082-3 | |
13. | Yuchao Dong, Jin Liang, Claude-Michel Brauner, Double free boundary problem for defaultable corporate bond with credit rating migration risks and their asymptotic behaviors, 2023, 372, 00220396, 505, 10.1016/j.jde.2023.06.050 | |
14. | Xinfu Chen, Jin Liang, Variational inequalities arising from credit rating migration with buffer zone, 2023, 0956-7925, 1, 10.1017/S095679252300030X | |
15. | Yang Lin, Jin Liang, A Multi-Credit-Rating Migration Model with Asymmetric Migration Boundaries, 2024, 29, 2297-8747, 7, 10.3390/mca29010007 | |
16. | Jin Liang, Bei Hu, 2024, Chapter 6, 978-981-97-2178-8, 109, 10.1007/978-981-97-2179-5_6 | |
17. | Jin Liang, Bei Hu, 2024, Chapter 8, 978-981-97-2178-8, 167, 10.1007/978-981-97-2179-5_8 | |
18. | Alicia Vidler, Toby Walsh, 2025, Decoding OTC Government Bond Market Liquidity: An ABM Model for Market Dynamics, 979-8-3315-0831-9, 1, 10.1109/CiFer64978.2025.10975734 |