Optimal Portfolio Choice in a Jump-Diffusion Model with Self-Exciting

Show more

1. Introduction

Empirical studies suggest that asset price encounters jumps and its volatility is stochastic. Further studies show that jumps occur in clusters, that is, a sequence of jumps occur in short time following a (big) jump which occurs after a relatively long quiet period of time. The feature of clustered jumps can be caught by a type of stochastic process known as Hawkes process. In this paper, we model occurrence of jumps by a Hawkes process hence our model is an extension of well-known jump-diffusion models, e.g. [1] . Meanwhile, we assume that such Hawkes jumps may occur in asset price itself as well as in its volatility. As a result, our model merges with the vast literature of stochastic volatility.

The contributions of the present paper are twofold. First, we solve the optimal investment problem and prove a verification result for a CRRA utility while [2] solves for a logarithm utility function. They do discuss the CRRA case but leave proofs to some references. However, our work shows that it is highly non-trivial to prove the existence and uniqueness of a classic solution to the associated Hamilton-Jacobi-Bellman (HJB) equation, which is essential for solving the optimal portfolio choice problem and implementing numerical methods. Second, our model may incorporate stochastic volatility or stochastic risk premium. Thus it shall be powerful to explain more financial phenomena, e.g. flight-to-quality, under-diversification, and possibly disclose more economic insights.

Technically, in the strand of the relevant literature, [3] introduces an Ornstein-Uhlenbeck (OU) type process of subordinator to model volatility. After that, there are several works about optimal portfolio selection in a model with OU type processes of subordinators. For example, [4] solves the portfolio choice problem in a model where volatility is a linear combination of Ornstein-Uhlenbeck type processes of subordinators. [5] solves the optimal investment and consumption problem in a similar model where a state variable (economic factor) is an Ornstein-Uhlenbeck type process of subordinator. They deal with a more general model, compared to [4] . As a result, they arrive at a nonlinear partial integro-differential equation, instead of a linear one in the latter. In both of these two papers, there are no jump components in the dynamics of asset prices. In stark contrast to the literature we incorporate jumps in the asset price in the present paper and hence confront a new/different challenge (if not more difficult) when we solve the optimal portfolio choice problem.

In our general setting, the volatility is a function of the state variable (the jump intensity) following a Hawkes process. Hence the present paper may combine two strands of research: modeling volatility by an OU type of process of subordinator and modeling jumps by Hawkes processes together. There are several papers studying the case where there are jumps in volatilities, e.g. [5] , etc., while the present paper considers jumps in both of asset price and volatility, and in jump intensity as well. This leads to a more complicated HJB equation than in [5] . As a result, our model has the features of stochastic volatility and self-exciting, while the optimal investment problem is still solvable in sense of by an iteration method, given that the iteration is proved to converge correctly to the unique classical solution of the problem.

This paper is incremental to the few aforementioned papers and dedicated to solving several technical problems related to CRRA utility functions. Nevertheless, the literature on portfolio choice is vast and is growing quickly. We would like to refer to some of them, for example, [6] , which proposes a semivariance method for diversified portfolio selection; [7] , that discusses a portfolio adjusting problem. The both assume that security returns are subject to experts’ estimations.

The organization of the paper is as follows. We formulate our model and optimal portfolio choice problem in Section 2. The Hawkes process is introduced in this section. In Section 3, we analyze the HJB equation. The existence and uniqueness of a classical solution to the equation are proved under some appropriate conditions. We also prove a verification theorem for the solution in Section 4. Two illustrative examples and one numerical example are provided in Section 5. Conclusion and further discussion are in Section 6. An extension is supplied in Appendix.

2. Problem Formulation

2.1. Hawkes Process for Self-Exciting Jumps

A Hawkes process is a counting process with self-exciting feature. Roughly speaking, it is a compounded Poisson process with stochastic intensity.

Given a complete probability space $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}{\left({\mathcal{F}}_{t}\right)}_{t\ge 0}\mathrm{,}P\right)$ , a counting process, ${\left\{{N}_{t}\right\}}_{t\ge 0}$ , satisfies

$P\left({N}_{t+\Delta t}-{N}_{t}=\mathrm{1|}{\mathcal{F}}_{t}\right)={\lambda}_{t}\Delta t+o\left(\Delta t\right)\mathrm{,}$

$P\left({N}_{t+\Delta t}-{N}_{t}>\mathrm{1|}{\mathcal{F}}_{t}\right)=o\left(\Delta t\right)\mathrm{,}$

where the intensity process ${\lambda}_{t}$ is given by the integrated form

${\lambda}_{t}={\text{e}}^{-\alpha t}{\lambda}_{0}+\left(1-{\text{e}}^{-\alpha t}\right){\lambda}_{\infty}+{\displaystyle {\int}_{0}^{t}\beta {\text{e}}^{-\alpha \left(t-s\right)}\text{d}{N}_{s}}$

or by the differentiation form

$\text{d}{\lambda}_{t}=\alpha \left({\lambda}_{\infty}-{\lambda}_{t}\right)\text{d}t+\beta \text{d}{N}_{t}\mathrm{;}$ (2.1)

here ${\lambda}_{\infty}$ is the long-run average of the jump intensity corresponding to the jump; $\alpha >0$ is the decay rate driving the jump intensity back to the long-run average before a jump occurs; $\beta \ge 0$ is a constant indicating non-negative impact of the jump occurrence on the jump intensity.

A jump process $\left\{{N}_{t}\right\}$ with the jump intensity described by (2.1) is called a Hawkes process [8] . It is known that the process is stationary if $\frac{\beta}{\alpha}<1$ . A Hawkes process differs from a doubly stochastic Poisson process since its increments are not independent. $\left\{{N}_{t}\right\}$ is not Markovian but $\left\{\left({N}_{t}\mathrm{,}{\lambda}_{t}\right)\right\}$ is. The compensated process ${N}_{t}-{\displaystyle {\int}_{0}^{t}{\lambda}_{s}\text{d}s}$ is a local martingale. For more information and a formal definition of Hawkes process, we refer to [8] . Figure 1 illustrates a sample path of one self-exciting (Hawkes) process.

The Hawkes process has a feature of self-exciting which is ideal to model jumps in financial markets. As one jump occurs, the jump intensity is increased by the occurred jump through the mechanics of (2.1). Hence a sequent jump happens more likely in a unit time following. In other words, as a jump happens, it impacts on the jump intensity as well as on itself. As a result, one may see a sequence of jumps in a short frame of time after one (big) jump. Thus, jump is self-triggered through the channel described by (2.1) and jumps tend to be clustered. Of course, it may not be the only channel to generate clustered jumps, but the empirical studies of [9] show an evidence that this channel is convincing.

Figure 1. A sample path of Hawkes $\left\{{N}_{t}\mathrm{,}{\lambda}_{t}\right\}$ and stock price $\left\{{S}_{t}\right\}$ .

Meanwhile, the mean-reversion property of (2.1) prevents the jump intensity from explosion given $0\le \beta <\alpha $ . Indeed, taking the expectation of (2.1) and using $\mathbb{E}\left[\text{d}{N}_{s}\right]=\mathbb{E}\left[{\lambda}_{s}\right]\text{d}s$ we obtain that

$\mathbb{E}\left[{\lambda}_{t}\right]={\lambda}_{0}{\text{e}}^{\left(\beta -\alpha \right)t}+\frac{\alpha {\lambda}_{\infty}}{\alpha -\beta}\left(1-{\text{e}}^{\left(\beta -\alpha \right)t}\right)\to \frac{{\lambda}_{\infty}}{1-\frac{\beta}{\alpha}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{as}\text{\hspace{0.17em}}\text{\hspace{0.05em}}t\to \infty .$

2.2. Asset Dynamics with Self-Exciting Jumps

We consider a market with a risk-free asset (bond) and a risky asset (stock). The price $\left\{{B}_{t}\right\}$ of the risk-free asset follows dynamics:

$\text{d}{B}_{t}=r{B}_{t}\text{d}t,$

where r is the risk-free interest rate. The price $\left\{{S}_{t}\right\}$ of the risky asset follows a stochastic differential equation

$\frac{\text{d}{S}_{t}}{{S}_{t-}}=\left[\mu \left({\lambda}_{t}\right)+r\right]\text{d}t+\sigma \left({\lambda}_{t}\right)\text{d}{W}_{t}+{Y}_{t}\text{d}{N}_{t}$ (2.2)

where $\left\{{N}_{t}\right\}$ is a Hawkes process described in the preceding section, $\left\{{W}_{t}\right\}$ is a standard Brownian motion, and ${Y}_{t}>-1$ is a random jump size independent of the jump process and the Brownian motion. To be more precise, we define

$\int}_{0}^{t}{Y}_{t}\text{d}{N}_{t}}={\displaystyle \underset{1\le i\le {N}_{t}}{\sum}{y}_{i$

where ${\left\{{y}_{i}\right\}}_{i=1}^{\infty}$ are i.i.d. random variables that are independent of either the Brownian motion or the jump process.

Note that the volatility $\sigma $ and the risk premium $\mu $ are set to be functions of the jump intensity. Given appropriate conditions about the mappings (e.g. monotonic mappings), we may write $\lambda =\lambda \left(\sigma \right)$ and $\mu =\mu \left(\sigma \right)$ , so our model is compatible with many models studied in stochastic volatility literature.

2.3. Optimal Portfolio Selection Problem

Now we turn to the Merton’s problem: An investor invests in the risky asset and the risk-free asset in a time horizon $\left[\mathrm{0,}T\right]$ . In order to maximize the expected utility of the terminal wealth, the investor needs to find an optimal investment strategy.

Let ${X}_{t}$ be the wealth of portfolio at time t and ${\pi}_{t}$ be the proportion of wealth invested into the risky asset at time t. Then

$\frac{\text{d}{X}_{t}}{{X}_{t-}}={\pi}_{t}\frac{\text{d}{S}_{t}}{{S}_{t-}}+\left(1-{\pi}_{t}\right)\frac{\text{d}{B}_{t}}{{B}_{t}}=\left({\pi}_{t}\mu +r\right)\text{d}t+{\pi}_{t}\sigma \text{d}{W}_{t}+{\pi}_{t}{Y}_{t}\text{d}{N}_{t}.$ (2.3)

An investment strategy is an adapted stochastic process
$\pi ={\left\{{\pi}_{t}\right\}}_{t\in \left[0,T\right]}$ . It is admissible if the associated wealth process is non-negative almost surely. The jump size
${Y}_{t}$ , in particular, is assumed to take a form of
${\text{e}}^{{Z}_{t}}-1$ , where
${Z}_{t}$ is a Gaussian random variable. This setting is popularly admitted in the literature of jump-diffusion model. See, for example, the seminal paper of [1] . As a result, an admissible strategy
$\pi $ shall satisfy the constraint
$0\le {\pi}_{t}\le 1$ . Thus, shorting either stock or bond is not permitted. This constraint condition may be relaxed^{1} according to distribution and support set of a specific jump size
${Y}_{t}$ . We denote all admissible strategies by
$\mathcal{A}$ .

An optimal investment strategy is a strategy that maximizes the expected utility of the terminal wealth. That is, the objective of an investor is to find V and ${\pi}^{\mathrm{*}}$ such that

$V\left(x,\lambda ,t\right)=\underset{\pi \in \mathcal{A}}{\mathrm{max}}{\mathbb{E}}_{t}^{x,\lambda}\left[U\left({X}_{T}\right)\right],\text{\hspace{1em}}{\pi}^{*}=\underset{\pi \in \mathcal{A}}{\mathrm{arg}\mathrm{max}}{\mathbb{E}}_{0}^{x,\lambda}\left[U\left({X}_{T}\right)\right]$ (2.4)

where ${\mathbb{E}}_{t}^{x\mathrm{,}\lambda}$ is the expectation conditioned on ${X}_{t}=x$ and ${\lambda}_{t}=\lambda $ . In this paper, we solve the problem for the CRRA utility $U\left(x\right)={x}^{p}/p$ . We prove the existence and uniqueness of a classical solution to the associated HJB equation when $p\in \left(\mathrm{0,1}\right)$ . By a similar approach, our framework may be extended to the case of $p<0$ regardless of an amount of efforts. The case of logarithm utility (corresponding to $p=0$ ) has been studied in [2] while assuming constant volatility and risk premium. In Section 5, we will discuss an extension case of logarithm utility with stochastic volatility and stochastic risk premium, as an application of our general results.

The Hamilton-Jacobi-Bellman (HJB) equation associated with the above stochastic optimization problem can be derived as

$\begin{array}{l}\alpha \left(\lambda -{\lambda}_{\infty}\right){V}_{\lambda}-{V}_{t}+\lambda V\\ =\underset{\pi}{\mathrm{max}}\left\{\left(r+\pi \mu \right)x{V}_{x}+\frac{1}{2}{\pi}^{2}{\sigma}^{2}{x}^{2}{V}_{xx}+\lambda \mathbb{E}\left[V\left(x\left(1+\pi Y\right)\mathrm{,}\lambda +\beta \mathrm{,}t\right)\right]\right\}\mathrm{,}\end{array}$ (2.5)

$V\left(x\mathrm{,}\lambda \mathrm{,}T\right)=U\left(x\right)\mathrm{.}$ (2.6)

^{1}If a tight bound of Y is
$-a\le Y\le b$ , where a and b are positive constants, then
$\left\{{\pi}_{t}\right\}$ is admissible if and only if
$-1/b\le {\pi}_{t}\le 1/a$ .

Here subscripts denote partial derivatives and Y is a random variable having the same distribution as ${Y}_{t}$ . In this paper, we shall make a full mathematical analysis of the HJB Equation (2.5) subject to the terminal condition (2.6), and a certain growth condition such as (2.8) below.

There are several papers studying the case where there are jumps in volatilities, e.g. [5] . Different from only jumps in volatility, our model incorporates jumps in both of asset price and volatility, and in jump intensity as well. It is worth to mention that the HJB Equation (2.5) is more complicated than that in [5] . The framework of our model is the same as that in [2] , except the setting of utility function. As well-known, CRRA utility functions generally involve much more difficult technical problems than the logarithm utility function.

Remark 2.1. The HJB Equation (2.5), together with the terminal condition (2.6), may have several solutions. From a view point of partial differential equation, it is necessary to prescribe an asymptotic behavior of the solution as $\lambda \to \infty $ . Note that, as a suboptimal strategy, investing everything into the risk-free asset gives the lower bound

$V\left(x\mathrm{,}\lambda \mathrm{,}t\right)\ge U\left(x{\text{e}}^{rt}\right)=U\left(x\right){\text{e}}^{prt}\mathrm{.}$ (2.7)

Although it is very had to estimate an upper bound, we shall prove the existence of a unique solution of (2.5)-(2.6) under the following growth condition: For some constant $C>0$ .

$\frac{1}{C}\le \frac{V\left(x,\lambda ,t\right)}{U\left(x\right)}\le C,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall \lambda \ge 0,t\in \left[0,T\right].$ (2.8)

3. Mathematical Analysis of the HJB Equation

3.1. Scaling Invariance

Note that applying the same strategy for two initial portfolios with initial condition $\left({X}_{0},{\lambda}_{0}\right)=\left(x,\lambda \right)$ and $\left({X}_{0},{\lambda}_{0}\right)=\left(1,\lambda \right)$ respectively, we find that the corresponding wealth of the two portfolios differ by a factor of x at any time $t\in \left[\mathrm{0,}T\right]$ . Hence, optimal strategies do not depend on x, and we have the scaling invariance

$V\left(x,\lambda ,t\right)=U\left(x\right)H\left(\lambda ,T-t\right),$ (3.1)

where $\tau =T-t$ is the time to expiry and $H\left(\lambda ,\tau \right):=V\left(1,x,T-\tau \right)/U\left(1\right)$ . Plugging (3.1) into (2.5) and using Remark 2.2 we obtain the equation for H:

${H}_{\tau}+\alpha \left(\lambda -{\lambda}_{\infty}\right){H}_{\lambda}=\underset{0\le \pi \le 1}{\mathrm{max}}\left\{A\left(\lambda ,\pi \right)H\left(\lambda ,\tau \right)+B\left(\lambda ,\pi \right)H\left(\lambda +\beta ,\tau \right)\right\},$ (3.2)

$H\left(\lambda ,0\right)=1,$ (3.3)

where

$A\left(\lambda ,\pi \right)=pr+p\mu \left(\lambda \right)\pi -\frac{p\left(1-p\right){\sigma}^{2}}{2}{\pi}^{2}-\lambda ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}B\left(\lambda ,\pi \right)=\lambda \mathbb{E}\left[{\left(1+\pi Y\right)}^{p}\right].$ (3.4)

^{2}We assume that the optimizer is achieved in the interior.

It is easy to check that both A and B are concave function of
$\pi $ . Hence, if
$H>0$ , then there exists a unique
${\pi}^{\mathrm{*}}$ such that^{2}

${\pi}^{*}\left(\lambda ,\tau \right)=\underset{0\le \pi \le 1}{\mathrm{arg}\mathrm{max}}\left\{A\left(\lambda ,\pi \right)H\left(\lambda ,\tau \right)+B\left(\lambda ,\pi \right)H\left(\lambda +\beta ,\tau \right)\right\}.$ (3.5)

In fact, with integrability of the jump size Y the first order condition brings us that

${\pi}^{*}=\frac{\mu \left(\lambda \right)}{\left(1-p\right)\sigma {\left(\lambda \right)}^{2}}+\frac{\lambda}{\left(1-p\right)\sigma {\left(\lambda \right)}^{2}}\mathbb{E}\left[{\left(1+{\pi}^{*}Y\right)}^{p-1}Y\right]\frac{H\left(\lambda +\beta ,\tau \right)}{H\left(\lambda ,\tau \right)}.$ (3.6)

The second term in the right hand side of the above formula is the hedging demand for the self-exciting jump risks.

In the rest of this section, we shall impose certain conditions on $\sigma (\cdot )\mathrm{,}\mu (\cdot )$ and Y and show that the HJB Equation (3.2) subject to the initial condition (3.3) and a certain growth condition admits a unique solution. In the next section we prove a verification result showing that the solution obtained solves the optimal investment problem.

For the case of CRRA utility, [2] suggest to prove the existence of a solution by verifying a contracting mapping as [5] . Our attempts show that it is not a trivial task to do that, instead, we use a different analytic method to accomplish the mission.

3.2. Basic Assumptions

First of all, we state some necessary assumptions as follow.

$p\in \left(0,1\right),\text{\hspace{0.17em}}\mu (\cdot ),{\sigma}^{2}(\cdot )\in {C}^{2}\left(\left[0,\infty \right)\right).$ (3.7)

Recall that we assume $Y~{Y}_{t}={\text{e}}^{{Z}_{t}}-1$ . Thus $Y>-1$ a.s., $\mathbb{E}\left[{Y}^{2}\right]<\infty $ .

To explain our idea in a clear manner, we assume that

$\underset{\lambda \to \infty}{\mathrm{lim}\mathrm{sup}}\left[\mu \left(\lambda \right)+\lambda \mathbb{E}\left[Y\right]\right]<0.$ (3.8)

Intuitively, the assumption (3.8) claims that the excess return of the asset will be negative if the jump frequency is high enough. In the literature of jump-diffusion models, the excess return is usually assumed with a compensation of jump risk, see, e.g. [10] . Hence the assumption (3.8) is equivalent to say that $\mu \left(\lambda \right)$ may not be enough compensated when (negative) jumps occur at a high frequency. That does not sound unreasonable.

The assumption may be relaxed to

$\underset{\lambda \to \infty}{\mathrm{lim}\mathrm{sup}}\frac{\mu \left(\lambda \right)+\lambda \mathbb{E}\left[Y\right]}{\sigma \left(\lambda \right)}<\infty .$ (3.9)

We shall discuss this later in the Appendix.

Under (3.8) we define

${\lambda}_{*}:=\mathrm{min}\left\{s\in \left[{\lambda}_{\infty},\infty \right)|\mu \left(\lambda \right)+\lambda \mathbb{E}\left[Y\right]\le 0,\forall \lambda \in \left[s,\infty \right)\right\}.$ (3.10)

3.3. An Integral Formulation

Let $\tau =T-t$ . We study the problem

$\{\begin{array}{ll}{H}_{\tau}+\alpha \left(\lambda -{\lambda}_{\infty}\right){H}_{\lambda}=\mathcal{F}\left[H\left(\cdot ,\tau \right)\right]\left(\lambda \right)\hfill & \forall \text{\hspace{0.05em}}\tau \in \left[0,T\right],\lambda \ge 0,\hfill \\ H\left(\lambda ,0\right)=1\hfill & \forall \text{\hspace{0.05em}}\lambda \ge 0,\hfill \end{array}$ (3.11)

where, with A and B as in (3.4), $\mathcal{F}$ is a non-local operator defined by

$\mathcal{F}\left[u\right]\left(\lambda \right)\mathrm{:}=\underset{0\le \pi \le 1}{\mathrm{max}}\left\{A\left(\lambda \mathrm{,}\pi \right)u\left(\lambda \right)+B\left(\lambda \mathrm{,}\pi \right)u\left(\lambda +\beta \right)\right\}\mathrm{.}$

Note that $\mathcal{F}$ satisfies, for any positive constant c and continuous functions u and v,

$\mathcal{F}\left[cu\right]=c\mathcal{F}\left[u\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathcal{F}\left[u\right]-\mathcal{F}\left[v\right]\le \mathcal{F}\left[u-v\right].$ (3.12)

We shall use characteristic curves to convert (3.11) into an equivalent integral formulation. For this we introduce

$D=\left[0,\infty \right)\times \left[0,T\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}\Lambda \left(\lambda ,\tau ,t\right):={\lambda}_{\infty}+\left(\lambda -{\lambda}_{\infty}\right){\text{e}}^{\alpha \left(t-\tau \right)}.$ (3.13)

For notational simplicity, in the sequel, we write $\Lambda \left(\lambda ,\tau ,t\right)$ simply as $\Lambda $ .

Let M be a positive constant to be determined. For fixed $\left(\lambda \mathrm{,}\tau \right)\in D$ and along the characteristic curve ${\left\{\left(\Lambda \left(\lambda \mathrm{,}\tau \mathrm{,}t\right)\mathrm{,}t\right)\right\}}_{0\le t\le \tau}$ , we obtain from (3.11) that

$\begin{array}{l}\frac{\text{d}}{\text{d}t}\left[{\text{e}}^{Mt}H\left(\Lambda \left(\lambda ,\tau ,t\right),t\right)\right]\\ ={\text{e}}^{Mt}\left[MH\left(\Lambda ,t\right)+{H}_{\tau}\left(\Lambda ,t\right)+\alpha \left(\Lambda -{\lambda}_{\infty}\right){H}_{\lambda}\left(\Lambda ,t\right)\right]\\ ={\text{e}}^{Mt}\left\{MH\left(\Lambda ,t\right)+\mathcal{F}\left[H\left(\cdot ,t\right)\right]\left(\Lambda \right)\right\},\text{\hspace{1em}}\forall t\in \left[0,\tau \right].\end{array}$

Integration over $t\in \left[\mathrm{0,}\tau \right]$ gives

$H\left(\lambda ,\tau \right)={\text{e}}^{-M\tau}+{\displaystyle {\int}_{0}^{\tau}}\left\{MH\left(\Lambda \left(\lambda ,\tau ,t\right),t\right)+\mathcal{F}\left[H\left(\cdot ,t\right)\right]\left(\Lambda \left(\lambda ,\tau ,t\right)\right)\right\}{\text{e}}^{M\left(t-\tau \right)}\text{d}t.$

Substituting the definition of $\mathcal{F}$ into the expression we obtain the integral formulation:

$H={\text{e}}^{-M\tau}+{\displaystyle {\int}_{0}^{\tau}}\underset{0\le \pi \le 1}{\mathrm{max}}\left\{\left[M+A\left(\Lambda \mathrm{,}\pi \right)\right]H\left(\Lambda \mathrm{,}t\right)+B\left(\Lambda \mathrm{,}\pi \right)H\left(\Lambda +\beta \mathrm{,}t\right)\right\}{\text{e}}^{M\left(t-\tau \right)}\text{d}t\mathrm{,}$ (3.14)

for every $\left(\lambda \mathrm{,}\tau \right)\in D$ where D and $\Lambda :=\Lambda \left(\lambda ,\tau ,t\right)$ are as in (3.13).

In the sequel we shall choose an appropriate positive constant M and solve the above integral equation by a monotonic iteration technique.

3.4. Determination of the Constant M

First we consider the function

$f\left(\lambda ,\pi \right):=A\left(\lambda ,\pi \right)+B\left(\lambda ,\pi \right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall \text{\hspace{0.05em}}\lambda \ge 0,\pi \in \left[0,1\right].$

Direct calculation yields

$f\left(\lambda ,0\right)=pr,$

${f}_{\pi}\left(\lambda ,\pi \right)=p\left[\mu \left(\lambda \right)-\left(1-p\right)\pi {\sigma}^{2}\left(\lambda \right)+\lambda \mathbb{E}\left[y{\left(1+\pi y\right)}^{p-1}\right]\right],$

${f}_{\pi \pi}\left(\lambda ,\pi \right)=-p\left(1-p\right)\left\{{\sigma}^{2}\left(\lambda \right)+\lambda \mathbb{E}\left[{y}^{2}{\left(1+\pi y\right)}^{p-2}\right]\right\}<0.$

Hence, we have the following:

Lemma 1. For each $\lambda \ge 0$ , there exists a unique ${\pi}^{\mathrm{*}}\left(\lambda \right)\in \left[\mathrm{0,1}\right]$ such that

$\zeta \left(\lambda \right):=\underset{0\le \pi \le 1}{\mathrm{max}}f\left(\lambda ,\pi \right)=f\left(\lambda ,{\pi}^{*}\left(\lambda \right)\right)\ge f\left(\lambda ,0\right)=pr.$ (3.15)

In addition, the following holds:

1) If $\lambda \in {I}_{0}:=\left\{\lambda \ge 0|\mu \left(\lambda \right)+\lambda \mathbb{E}\left[y\right]\le 0\right\}$ , then ${\pi}^{*}\left(\lambda \right)=0$ and $\zeta \left(\lambda \right)=pr$ ;

2) If $\lambda \in {I}_{1}:=\left\{\lambda \ge 0|\mu \left(\lambda \right)-\left(1-p\right){\sigma}^{2}\left(\lambda \right)+\lambda \mathbb{E}\left[y{\left(1+y\right)}^{p-1}\right]\ge 0\right\}$ , then ${\pi}^{*}\left(\lambda \right)=1$ ;

3) If $\lambda \in I:=\left[0,\infty \right)\backslash \left({I}_{1}\cup {I}_{2}\right)$ , then ${\pi}^{\mathrm{*}}\left(\lambda \right)\in \left(\mathrm{0,1}\right)$ . In addition, ${\pi}^{\mathrm{*}}(\cdot )$ is smooth on I.

Consequently, if (3.8) holds, then for ${\lambda}_{\ast}$ as in (3.10),

${\pi}^{*}\left(\lambda \right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\zeta \left(\lambda \right)=pr,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall \lambda \ge {\lambda}_{\ast}.$ (3.16)

We now define

$M=\underset{0\le \pi \le \mathrm{1,0}\le \lambda \le {\lambda}_{\infty}+\left({\lambda}_{\ast}-{\lambda}_{\infty}\right){\text{e}}^{\alpha T}}{\mathrm{max}}\left\{\left|A\left(\lambda \mathrm{,}\pi \right)\right|\right\}\mathrm{.}$

This implies that

$M+A\left(\lambda \mathrm{,}\pi \right)\ge 0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall \pi \in \left[\mathrm{0,1}\right]\mathrm{,}\lambda \in \left[\mathrm{0,}{\lambda}_{\infty}+\left({\lambda}_{\ast}-{\lambda}_{\infty}\right){\text{e}}^{\alpha T}\right]\mathrm{.}$ (3.17)

3.5. Monotonic Iteration

We now solve (3.14) by the following iteration: for each non-negative integer n and $\left(\lambda \mathrm{,}\tau \right)\in D$ , we define iteratively ${H}_{n}$ by

${H}_{0}\left(\lambda \mathrm{,}\tau \right)\mathrm{:}={\text{e}}^{pr\tau}\mathrm{,}$ (3.18)

$\begin{array}{l}{H}_{n+1}\left(\lambda \mathrm{,}\tau \right)\mathrm{:}={\text{e}}^{-M\tau}+{\displaystyle {\int}_{0}^{\tau}}\underset{0\le \pi \le 1}{\mathrm{max}}\{\left[M+A\left(\Lambda \mathrm{,}\pi \right)\right]{H}_{n}\left(\Lambda \mathrm{,}t\right)\\ \text{\hspace{1em}}\text{\hspace{1em}}\text{\hspace{1em}}\text{\hspace{1em}}\text{\hspace{1em}}\text{\hspace{0.17em}}+B\left(\Lambda \mathrm{,}\pi \right){H}_{n}\left(\Lambda +\beta \mathrm{,}t\right)\}{\text{e}}^{M\left(t-\tau \right)}\text{d}t\mathrm{.}\end{array}$ (3.19)

Since for every and, we see that is a well-defined family of continuous functions on D. We introduce

Lemma 2. For each integer and every,.

Proof. We use a mathematical induction. Assume that for all. Then when and, we have

Thus, and for all. Consequently,. Also, by (3.16),. It then follows from (3.19) that

Hence, by mathematical induction, for every and every non-negative integer n. This completes the proof.

In the sequel, we focus on the case. Note that

Hence,

(3.20)

Lemma 3. For each integer, on D.

Proof. We need only consider the case. When, we obtain from (3.18), (3.19), and the definition of in (3.15) that

Thus,. Next using (3.20) we can show by a mathematical induction that on. This completes the proof.

Now we establish an upper bound. We define

(3.21)

Lemma 4. For each integer and,.

Proof. We use an induction argument. Assume that the assertion holds for some. Then for, by (3.19)-(3.20) we obtain

Thus, the assertion of the Lemma holds for. This completes the proof.

3.6. Solution of the Integral Equation

The family is a bounded monotonic family, so exists. We wish to prove a uniform convergence. For this, we introduce a norm by

Note that when,

Hence, using we obtain, when and,

By a mathematical induction, one can derive that

Hence, we have the following:

Lemma 5. There exists a function such that

In addition, H is a solution of (3.14) and has the following properties:

3.7. Lipschitz Continuity

For each, we define

For, we consider the function

Lemma 6. Let and. Then for each there exist a unique such that

In addition, is Lipschitz continuous on, , and

(3.22)

Proof. When, we have so

It remains to consider the case. Note that for,

Thus, there exists a unique such that on attains its maximum at. Note that

Here we define if; in this case we have for any. Hence, as, we see that is continuous and that I is an open set. When, by the implicitly function theorem,

Since for almost every, by the continuity of, we hence know that and. This completes the proof.

We now calculate the norm of. Notice that

When, we have so. When and,

where

Hence, we have

Lemma 7. Under the condition of Lemma 6, we have

Now, applying this estimate for (3.19) and using and we obtain

(3.23)

(3.24)

Thus,

Lemma 8. For every non-negative integer n, and for each,

(3.25)

Proof. Clearly, the assertion (3.25) holds when. Assume that (3.25) holds for some. Then we have

Thus, by mathematical induction, the assertion of the Lemma holds.

Similarly, using (3.24), we can obtain an estimate for. We omit the details.

3.8. Estimate

Now assume that

Differentiating with respect to we obtain

, a.e.

Hence, using estimate (3.22) for, we find that

We can calculate, when and,

Hence,

where

Next, we calculate for and,

Note that, for,

Hence,

It then follows that

Thus, we obtain

Hence, we have the following lemma.

Lemma 9. There exists a constant such that for each non-negative integer n,

(3.26)

Proof. The estimate for follows from a mathematical induction. The estimate for and follows by differentiating (3.24).

Now sending, obtain the following:

Theorem 1. Assume (3.7) and (3.8). Then the problem (3.11) admits a classical solution H that has the following properties:

1) For each and,;

2) For each,

3.9. Uniqueness

Theorem 2. There exists a unique solution of (3.11) in the following class

(3.27)

Proof. Let H and be two solutions of in (3.11) in. Suppose is not true. Then there exists such that. Set

Then. Now for each we define

Then. Since, we see that. Denote by the characteristic curve passing through. Then

Hence, denoting we obtain

This implies that

Now define. Then we derive that

This implies that, for,

which contradicts the boundedness of. Thus we must have. Similarly, we can show that, so. This completes the proof.

4. Verification Theorem and Optimality of the Solution

That the optimal investment problem can be solved through the preceding results is based on a verification result which guarantees that the solution to the HJB equation is the value function corresponding to the optimal investment problem.

Theorem 3. Suppose. Assume that solves (2)--(3). Set and set by (5). Then, the value function defined in (2.4); also the optimal investment strategy is given by.

4.1. Preliminary Results

Lemma 10. Let be positive constants and be stochastic process satisfying

(4.1)

(4.2)

Then: 1) With, for each constant and,

where is the solutions of the o.d.e.

(4.3)

2) Assume that. Denote by the unique positive root of. Then where a) if, then for every, b) if, then at time

Proof. The assertion of the Lemma follows from a general result for the joint characteristic function of and (see e.g. [9] ):

(4.4)

where satisfy the Riccati equation:

(4.5)

(4.6)

4.2. Proof of the Verification Theorem

Let be an admissible strategy and be the resulting wealth of the portfolio. Let be the solution of HJB equation obtained in the previous section. Subject to, by Itô Lemma, we have

(4.7)

Under the condition, since and H is bounded, we can verify that

Therefore, taking expectation of (4.7) we obtain

Thus,

After taking the supreme over all admissible strategies, we then obtain

Finally, one can check that if we take, the above gives us an identity. Thus,

This completes the verification theorem.

The above theorem provides a verification result for a general framework of stochastic volatility and double jump models.

5. Two Illustrative Cases and a Numerical Example

In the following, we present two illustrative cases. The first case of logarithm utility function is studied in [2] , where and are constants. In the present paper, we extend the case to allow both of them stochastic. In the second case of, the Equation (3.2) has a solution in form of. This type of model is used to study effects of rare events in [10] . They do not show the verification of this solution though.

5.1. Logarithm Utility

When the utility function is logarithm, i.e., the solution is more explicit. Applying the same strategy for initial condition and one finds that the optimal strategy does not depend on x. Hence, setting we have. Then we obtain a simpler first order condition for:

(5.1)

Note that, so there exists a unique solution such that. The PDE for becomes

where is the differential-difference operator associated with the self-exciting Hawkes process and F is a function associated with the asset dynamics defined by

By Feynman-Kac formula, can be written as an expectation of a function of the process:

(5.2)

When and are constants, [2] study the optimal investment and consumption problem, and prove a verification result for the logarithm utility case. In the present, we may allow both and be stochastic.

For this logarithm case, it is straightforward to show that if all jumps are negative, i.e., then, hence we find the fact consistent with the phenomena known as flight-to-quality: as a market crash happens, all positions in risky assets are reduced.

5.2. The Case of Stochastic Volatility

We let and. Substituting the conditions into the model and setting, we arrive at

This model is close to the one used in [10] to study effects of rare events, except that there is no diffusion term in the variance dynamics here. The dynamics of the variance becomes a type of OU process of subordinations if we extensively replace by a pure jump Lévy process with no drift and positive increments (subordination). Such a process is used to model volatility in [3] , [4] , or [5] , where no simultaneous jump is assumed in their asset price.

Given the objective problem and the utility function, we let and derive the HJB equation:

The solution can be expressed as

(5.3)

and solve

(5.4)

with. The optimal strategy is given by

(5.5)

When the assumption (3.8) or holds, we know a classical solution to the HJB equation exists by the results in the preceding sections. In such a case, classical solutions to the ordinary differential equations (5.4) exist. We show below that could be explosive in a finite time horizon if the assumption is not satisfied.

We consider the special case:. Assume which conflicts with the assumption (3.8). In this case, (5.4) becomes

where. Note that given, we have

We suppose and choose such that

then we have

Hence blows up during where is given by

As a conclusion, the optimal portfolio choice problem in the stochastic volatility model may meet the issue of exploding in a finite time horizon. This paper provides a sufficient condition under which a classical solution exists and exploding does not happen. Our results shall be useful in studies of a stochastic volatility model or alike.

5.3. A Numerical Example

In our proof, the monotonic iteration (3.18) and (3.19) actually suggest an iteration algorithm to find the solution numerically. Figure 2 gives a numerical example. The upper edge of curves stands for the limit of iteration corresponding to and, respectively. The figure indicates that the iteration converges monotonically. The following parameters are assumed for this example.

where in makes the total excess return goes to negative as jump occurs at a high frequency. A constant, e.g. satisfying the condition (3.8) as well results in a figure with the same pattern. The time horizon T is 0.2, and the iteration procedure is stopped as soon as

where ‘s are points of a partition of.

Figure 2. This figure shows the numerical iteration algorithm generates sequences converging monotonically to the limits that correspond to and, respectively.

6. Conclusions and Discussions

We study the optimal portfolio choice problem in a jump-diffusion model where the jump likelihood is increased by jump itself. We establish the existence and uniqueness of a classical solution to the corresponding HJB equation. A verification theorem which guarantees the optimality of the solution is proved. Our approach relies on a monotonic iteration procedure which naturally hints a numerical algorithm. We consider CRRA utilities and a stochastic investment opportunity set which may nest interesting models in the literature.

By similar steps, our work may be extended to the case of the risk aversion () greater than 1. It is also possible to extend our framework to a multi-dimension case like [2] . With multi-assets, jumps may be not only self-excited but also mutually excited. The latter feature may be suitable to study financial contagions (see, e.g. [9] ). At last, as mentioned earlier, the assumption (3.8) can be relaxed to be

(6.1)

We provide a proof of existence of a continuous solution under this relaxed condition in the appendix. However, we need certain structural conditions on derivatives of and, in order to follow a similar line of proofs presented in Section 3 and show the regularity and uniqueness of the solution. We will investigate these conditions in the future.

Acknowledgements

The authors acknowledge support from NSFC (No. 11371280 and No. 71771142) and NSF (DMS-1516344).

Appendix. Existence of a Continuous Solution under a Relaxed Condition

In this appendix, we establish the existence of a solution H of (3.11) under the condition (3.9) instead of the condition (3.8). For each, we define

We write the partial differential equation in (3.11) as

where

We define a family iteratively as follows:. If is defined, we define as the unique solution of the linear problem

In term of characteristic curves, the linear problem for is well-posed.

1) Note that

It then follows by a comparison principle for the equation for that.

2) Next, assume that is an integer and. Then

Here we use the fact that and. Thus, by comparison principle, we have. Consequently, by mathematical induction, is an increasing family.

3) Note that, since is a concave function,

Thus,

Note that if we have. Thus, under the assumption

(A.1)

we have

Now assume that. Then we have

It then follows by comparison principle that.

Thus we have, for each integer,

Then is a bounded monotonic sequence of continuous functions. It follows from the Dini theorem that the sequence converges uniformly on compact subsets of D; so is continuous. Under certain structural conditions on derivatives of and we may follow a similar line of proofs presented in Section 3 to show the regularity and uniqueness of the solution.

References

[1] Merton, R. (1976) Option Pricing When Underlying Stock Returns Are Discontinous. Journal of Financial Economics, 3, 125-144.

https://doi.org/10.1016/0304-405X(76)90022-2

[2] Aït-Sahalia, Y. and Hurd, T.R. (2015) Portfolio Choice in Markets with Contagion. Journal of Financial Econometrics, 14, 1-28.

https://doi.org/10.1093/jjfinec/nbv024

[3] Barndorff-Nielsen, O.E. and Shephard, N. (2001) Non-Gaussian Ornstein-Uhlenbeck-Based Models and Some of Their Uses in Financial Mathematics. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 63, 167-241.

https://doi.org/10.1111/1467-9868.00282

[4] Benth, F.E., Karlsen, K.H. and Reikvam, K. (2003) Merton’s Portfolio Optimization Problem in Black-Scholes Market with Non-Gaussian Stochastic Volatility of Ornstein-Uhlenbeck Type. Mathematical Finance, 13, 215-244.

https://doi.org/10.1111/1467-9965.00015

[5] Delong, L. and Klüppelberg, C. (2008) Optimal Investment and Consumption in a Black-Scholes Market with Lévy-Driven Stochastic Coefficients. Annals of Applied Probability, 18, 879-908.

https://doi.org/10.1214/07-AAP475

[6] Chen, L., Peng, J., Zhang, B. and Rosyida, I. (2017) Diversified Models for Portfolio Selection Based on Uncertain Semivariance. International Journal of Systems Science, 48, 637-648.

https://doi.org/10.1080/00207721.2016.1206985

[7] Qin, Z., Kar, S. and Zheng, H. (2016) Uncertain Portfolio Adjusting Model Using Semiabsolute Deviation. Soft Computing, 20, 717-725.

https://doi.org/10.1007/s00500-014-1535-y

[8] Hawkes, A.G. (1971) Spectra of Some Self-Exciting and Mutually Exciting Point Processes. Biometrika, 58, 83-90.

https://doi.org/10.1093/biomet/58.1.83

[9] Aït-Sahalia, Y., Cacho-Diaz, J. and Laeven, R.J.A. (2015) Modeling Financial Contagion Using Mutually Exciting Jump Processes. Journal of Financial Economics, 117, 585-606.

https://doi.org/10.1016/j.jfineco.2015.03.002

[10] Liu, J., Longstaff, F. and Pan, J. (2003) Dynamic Asset Allocation with Event Risk. Journal of Finance, 58, 231-259.

https://doi.org/10.1111/1540-6261.00523