1. Introduction
A pension fund is an important financial instrument for individuals to reallocate incomes and sustain consumption after retirement. Generally, according to the determination of benefits, there are two typical types of pension plans: defined benefit (DB) and defined contribution (DC) pension plans. In DB plans, benefits are fixed in advance, while, in the DC case, contributions are fixed by the trustee. There are two phases in a pension scheme: the accumulation phase, which is the period from the entry time to the retirement time, and the decumulation phase, which is the period from the retirement time to death.
Next, we review the literature relevant to our paper. In the accumulation phase of a DC pension scheme, the contributor contributes part of his/her salary to the fund. Since the salary is related to the profitability of the company, the works of Bodie et al. [
1] and Dybvig and Liu [
2] assume that the salary process is spanned by the stock price. In addition, Guan and Liang [
3] and Li and Wang [
4] describe the salary process using a Heston stochastic volatility model, i.e., the salary is correlated with the volatility of the stock. Furthermore, Zeng et al. [
5] assume that the salary process is related to stochastic volatility. Based on [
3,
4,
5], we add an independent random process to the stochastic salary process to be closer to the reality.
It is appropriate to insert a jump process in the stochastic salary due to a promotion and job-hopping. Moreover, it is realistic to introduce jumps in volatility, which represent some unexpected events, such as an economic crisis or policy adjustments by the government. In our model, the contribution rate of the pension scheme is proportional to the salary of the individual, the dynamics of which follows a Heston stochastic volatility model with jumps. In addition, salary and variance are correlated by means of a common shock. In reality, a common component may depict an event that has an impact on both salary and volatility. Common shock models are widely used in the area of actuarial science. For instance, in Liang et al. [
6], the insurance risk model is modulated by a compound Poisson process, and the two-jump-number processes are correlated through a common shock. Liang et al. [
7] assume that jumps in both the risky asset and insurance risk process are correlated through common shock dependence.
Since the period of a pension scheme is usually long, the inflation risk should be considered during the optimization phase. There are various studies focused on the stochastic optimization problem for DC pension plans under the inflation risk. For example, ref. [
8] explores the optimal asset allocation problem with downside protection and stochastic inflation risk. Yao et al. [
9] solve an optimal portfolio decision problem under the inflation risk and mean-variance criterion. Other relevant works on optimal control under inflation risk can be found in [
10,
11]. The inflation risk is also involved in our model. Instead of a pure diffusion process, we introduce a jump diffusion process to model the nominal price level of a representative bundle of commodity goods in the market. In other words, the dynamics of the price index given by Zhang et al. [
12] and Zhang and Ewald [
13] are extended in our model, and a Poisson jump is included in the evolution of the index price.
Stochastic processes are commonly used to model the uncertainty in the financial market. For further study of stochastic processes, we refer interested readers to [
14,
15,
16,
17,
18,
19,
20]. In our paper, we assume that the pension manager is allowed to invest in two types of assets: the bond and the stock. The dynamics of the bond price follows a geometric Brownian process. The stock price is driven by a drifted Brownian motion and a Poisson jump. The similar asset model is widely used in various asset allocation problems. For example, Merton [
21] considers Poisson jumps in an optimal dynamic portfolio decision problem. In a DC pension funding framework, Sun et al. [
22] deal with the pre-commitment and equilibrium investment strategies by incorporating jumps into the risky asset process. More relevant works on jump diffusion asset allocation problems in pension management can be found in [
23,
24,
25,
26,
27,
28].
The aim of pension management is to find the optimal investment and minimize the expected distance between the terminal wealth and two given targets. To find the optimal policy, the dynamic programming principle is used to derive the Hamilton–Jacobi–Bellman (HJB) equation. From the classical optimal control theory, once a continuously differentiable solution of the HJB equation is explicitly solved, the optimal value function and the optimal policy can be derived. In our paper, by solving the HJB equation, we show the explicit form of the optimal investment policy and the optimal value function. The dynamic programming principle and HJB equation are applicable to various optimization problems. However, the drawback is that if there is no explicit solution for the HJB equation, then the dynamic programming principle fails to solve the problem. In our paper, there is a continuously differentiable solution for the HJB equation, since there is only one specific boundary condition in the HJB equation. If there are more than two boundary conditions for the optimization problem, it will be more difficult to find the explicit expression for the value function as well as the optimal policy. When there is no explicit solution for the HJB equation, it is more appropriate to adopt the maximum principle, Martingale approach or viscosity solution.
The rest of the paper is structured as follows.
Section 2 describes the financial market with the jump diffusion price index, as well as two tradable assets that are of interest for pension management. This section also gives the pension model.
Section 3 deals with a stochastic optimal control problem in order to minimize the fluctuation in the final real wealth over a finite horizon. The closed form of the investment strategy is given by solving the HJB equation. Finally,
Section 4 gives the sensitivity analysis and
Section 5 establishes the conclusions.
2. Model Assumptions and Notations
Consider a probability space , with as the real-world probability measure on and . The filtration is generated by a five-dimensional standard Brownian motion , i.e., , which represents the risk sources of the interest rate, inflation, stock price, salary and its volatility, respectively. and are correlated, which is captured by the coefficient . and are also correlated, which is captured by the coefficient . Let the filtration be generated by a five-dimensional Poisson process with intensity , where , i.e., , which represent the jumps in inflation, stock price, salary, volatility and common shocks between two jumps, respectively. Suppose that Poisson processes are mutually independent. Moreover, Brownian motions are independent of Poisson processes on .
2.1. The Financial Market
Following the work of Eisenberg [
29], we assume that the discount factor is a geometric Brownian motion:
where
, and
.
The financial market consists of two underlying instruments that are traded continuously over time and perfectly divisible. Suppose that there are no transaction costs or taxes in the given context. The bond
evolves according to the the following dynamics:
with initial price
.
Besides the cash account, the trustee also has the opportunity to invest the fund into a stock with the dynamics
where
is the appreciation rate for the stock.
is the volatility associated with the diffusion component of the stock price.
denotes the magnitude of a jump. We state that
to prevent the process from jumping to a value below zero.
describes the fluctuation, and
describes the jump in the stock price. For simplicity, it is assumed that
and
are independent stochastic processes.
2.2. The Pension Model
This paper considers the accumulation phase of a DC-type pension plan. Assume that the entry time of a pensioner is the initial time 0, and his/her retirement time is the terminal time in our model. We denote the pensioner’s death time as
, which is a positive random variable defined on the probability space
. The mortality rate
is defined as
In a general pension plan, the pensioner pays contributions before the retirement time
T, where
. The level of the contribution rate is usually defined as a proportion
of the pensioner’s salary. In previous works, such as refs. [
3,
4,
5], it is assumed that the stochastic salary is driven by the Heston stochastic volatility model, i.e., the salary process has stochastic volatility and the salary return variance is governed by a mean-reverting process. To be more realistic, we add a Brownian motion
in the salary process to describe the fluctuation in the salary itself. In addition, we also assume that there are possible common Poisson jumps between the salary process and the stochastic variance. The dynamics of the investor’s salary
follow
while the stochastic volatility
is governed by
Regarding the parameters in (
3),
denotes the instantaneous expected rate of the salary, and
is the instantaneous volatility scale factor measuring how the risk source of the stock price affects the salary.
and
denote the magnitude of the jumps associated with Poisson processes
and
, respectively.
describes the jump in the salary itself, and
describes the possible common jumps between the salary (given by Equation (
3)) and the stochastic volatility
(given by Equation (
4)). As introduced earlier, we state that the Brownian motion
in Equation (
3) and the Brownian motion
in Equation (
4) are correlated with the coefficient
.
Regarding the parameters in (
4),
denotes the mean-reversion rate,
denotes the long-run mean, and
is the volatility coefficient. The assumption
is proposed to guarantee the volatility process
. The Brownian motion
describes the fluctuation in the volatility, and the Poisson process
describes the jump in the volatility. We state that
and
to prevent the process
from jumping to a value below zero. It should be noted that all Poisson processes are mutually independent. Moreover, Brownian motions are independent of Poisson processes.
3. The Optimal Portfolio
The aim of the stochastic control problem is to find the optimal investment decision. The pension trustee continuously decides on the weights invested into the cash account and the stock. We denote the nominal wealth at time
t as
. Under the investment policy chosen, it is easy to obtain the following stochastic differential equation, which describes the evolution of the wealth:
with
.
denotes the weight invested into the stock at time
t. The remainder,
, is the proportion invested into the cash account. Borrowing and short-selling are permitted in the given context. A negative value of
means that the pension trustee takes a short position in the stock, while a negative value of
reflects that the trustee borrows money from the bank to purchase the risky asset.
By substituting Equations (
1) and (
2) into Equation (
5), we obtain that
As mentioned in
Section 1, the time horizon for the accumulation phase of a pension fund (in our model, from time 0 to
T) is usually long; hence, the influence of inflation is considered in the given context.
The price index at time
t is denoted by
, which refers to the purchase power per unit of money. The dynamics are driven by a jump diffusion process of the following type:
with initial value
.
is the instantaneous expected inflation rate.
is the instantaneous volatility associated with the diffusion component, and
denotes the magnitude of a jump with the condition
to ensure that the price index remains strictly positive. In Equation (
6),
and
are independent stochastic processes.
Next, we define the corresponding real salary process as the following:
Definition 1. The real salary process is defined by Applying the quotient rule of Itô’s formula,
is given by
with initial value .
Then, the real wealth process with the consideration of inflation follows
with initial condition .
Next, we restrict the strategies in order to fulfil some technical conditions. We call a strategy
an admissible control process if it is
-measurable, Markovian and stationary and satisfies the condition
Denote
the set of all admissible controls, i.e., it is the set of all measurable processes
, which satisfies Equation (
9). Next, we try to find the optimal investment strategy for the DC pension plan manager under
.
Assume that the pension trustee has a preference to minimize the expected value of the fluctuations in the terminal wealth until time
, where
T is the terminal time of the control problem. The objective is to minimize
with
as the conditional expectation given the filtration
.
and
are two positive constants representing the target funds of the plan at time
T and
, respectively. The deviation between the actual fund and the target fund is called the discontinuity risk; see Wang et al. [
30].
In quadratic loss functions (
10), any deviations between
and
(or
) are penalized. To be more specific, we assumed that
and
in Equation (
10) to characterize that under-funding is more penalized than over-funding. A similar setting can be seen in Devolder, Janssen and Manca [
31] and Zhang and Guo [
32].
For a better understanding, we use a flowchart to describe the whole research process in
Figure 1.
According to
and
the objective function with an uncertain lifetime can be converted into the following deterministic horizontal function:
The dynamic programming approach is used to solve the stochastic optimization problem. Define the value function as
In stochastic optimal control theory, the HJB equation accomplishes the connection between the value function and the optimal control; see, for instance, the books [
33,
34,
35,
36,
37] and the papers [
38,
39,
40]. The HJB equation is
where
with terminal condition . , , , , , , , and denote the first- and second-order partial derivatives of the value function with respect to t, , and V, respectively.
If there exists a twice continuously differentiable solution of Equation (
12), strictly convex, then the minimizer of the investment strategy is obtained by the optimal functional
, which satisfies the following necessary conditions:
We shall frequently use the following notations. Define
By using the first-order condition and solving the HJB equation, the explicit form of the optimal investment decision is given by the following theorem.
Theorem 1. (Main result) The optimal investment strategy on the stock is given by where , , , , and are given by Equation (A6), Equation (A7), Equation (A28), Equation (A27), Equation (A8) and Equation (A9), respectively. 4. Sensitivity Analysis
In order to investigate the influence of the parameters on the optimal investment decision, we provide a sensitivity analysis in this section. Unless otherwise stated, the employed parameters of the model are based on the following annualized benchmark values presented in
Table 1. In what follows, we mainly explore the impacts of the volatility
, the jump magnitude
and the jump intensity
,
on optimal dividend policy
.
Figure 2 gives a possible path simulation of the optimal investment decision
. Denoting the optimal wealth process by
, a possible path simulation is given in
Figure 3. As we can see, the path simulation of the wealth of the pension fund is increasing as time passes, which highly coincides with reality.
Figure 4 analyzes the relationship between the optimal investment decision
and the stock volatility
. First, we can see that the optimal investment decision
is always negative under the given parameters in
Table 1, which means that the pension manager prefers to be a short-seller and invests more money in less risky bonds. From
Figure 4, we can also see that when
increases, the optimal investment policy
increases. This can be explained by the fact that as
increases, the risk and return of the pension fund simultaneously increase. To achieve the given target value
and
, the manager prefers to take more risks as well as achieve more profits.
In
Figure 5, as
increases, the optimal investment policy
also increases. This phenomenon can be explained as follows. The parameter
measures the jump magnitude of the salary. When
is positive and increasing, which means that there is a higher jump in the pensioner’s salary surplus, or, in other words, the pensioner will input more money into the pension fund and eventually the manager of the pension fund will increase the investment amount for the risky stocks to achieve higher profits.
Figure 6 and
Figure 7 show the impact of the Poisson intensity
,
on optimal policy
, respectively. As we can see, as
,
increase, the optimal investment amount
increases. This can be explained as follows.
,
represent the intensity of a positive jump in salary surplus. If
(or
) increases, then the pensioner will be more positive about the future and will make more contributions to the pension fund. Eventually, the pension manager is able to invest more money in the risky asset.
Figure 6 depicts the impact of
on the optimal policy under different parameters
.
Figure 7 depicts the impact of
on the optimal policy under different parameters
. Both figures show the same conclusion, i.e., that the higher the salary jump magnitude, the more money should be invested in the stock market. The impact of
on
is similar to that of
since
and
both measure the jump magnitude of the salary.
Figure 8 depicts the impact of the jump magnitude of the salary on the optimal policy
with different
. As we can see, as
increases, the investment policy
also increases. This can be explained as follows: if the average wage of society increases, then the investment enthusiasm will increase.
Figure 9 depicts the impact of
with different
. The increasing
leads to an increase in investment policy
. This shows that if the volatility of the salary is high, then the optimal choice is to increase the investment in the risky asset to reach the desired target as soon as possible, to avoid possible losses. From the above analysis, from the perspective of the government, an increase in salary will increase the investment in risky assets. On the other hand, to encourage the pension manager to invest more money in a risk-less bond, the government should be reduce the intensity of wage growth.