Next Article in Journal
Value-at-Risk Effectiveness: A High-Frequency Data Approach with Semi-Heavy Tails
Next Article in Special Issue
Intangible Assets and Analysts’ Overreaction and Underreaction to Earnings Information: Empirical Evidence from Saudi Arabia
Previous Article in Journal
The Regime-Switching Structural Default Risk Model
Previous Article in Special Issue
Market Equilibrium and the Cost of Capital with Heterogeneous Investment Horizons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Longevity-Indexed Bond in Risk Management of Aggregated Defined Benefit Pension Scheme

1
School of Economics and Management, Hebei University of Technology, Tianjin 300401, China
2
School of Finance, Capital University of Economics and Business, Beijing 100070, China
3
School of Mathematical Sciences, Nankai University, Tianjin 300071, China
*
Author to whom correspondence should be addressed.
Risks 2024, 12(3), 49; https://doi.org/10.3390/risks12030049
Submission received: 27 January 2024 / Revised: 22 February 2024 / Accepted: 4 March 2024 / Published: 6 March 2024
(This article belongs to the Special Issue Optimal Investment and Risk Management)

Abstract

:
Defined benefit (DB) pension plans are a primary type of pension schemes with the sponsor assuming most of the risks. Longevity-indexed bonds have been used to hedge or transfer risks in pension plans. Our objective is to study an aggregated DB pension plan’s optimal risk management problem focusing on minimizing the solvency risk over a finite time horizon and to investigate the investment strategies in a market, comprising a longevity-indexed bond and a risk-free asset, under stochastic nominal interest rates. Using the dynamic programming technique in the stochastic control problem, we obtain the closed-form optimal investment strategy by solving the corresponding Hamilton–Jacobi–Bellman (HJB) equation. In addition, a comparative analysis implicates that longevity-indexed bonds significantly reduce solvency risk compared to zero-coupon bonds, offering a strategic advantage in pension fund management. Besides the closed-form solution and the comparative study, another novelty of this study is the extension of actuarial liability (AL) and normal cost (NC) definitions, and we introduce the risk neutral valuation of liabilities in DB pension scheme with the consideration of mortality rate.

1. Introduction and Motivation

There are two primary types of pension schemes: defined contribution (DC) plans and DB plans (Recent works on DC and DB pension plans can be found in Ng and Chong (2023); Guan et al. (2022), respectively). In a DC pension scheme, contributions are deterministic or a percentage of the member’s salary (which may be stochastic), with benefits depending on the accumulation of contributions and investment returns. In this scenario, the individual bears the majority of the inherent risk. Conversely, in a DB pension scheme, the benefits received upon retirement are clearly defined by the plan’s rules (typically related to salary and years of service), with the sponsor assuming most of the risks, including investment, inflation, longevity and solvency risks.
Longevity risk, the concern of policyholders or affiliates outliving insurers’ or pension sponsors’ expectations, has garnered significant academic and practical interest recently. The work of Cox et al. (2013) proposes hedging strategies for managing this risk, yet natural hedging methods cannot mitigate all aspects of longevity risk within nonparametric mortality models; see Zhu and Bauer (2014). Further studies on longevity risk are detailed in subsequent research: Blake and Cairns (2021); Broeders et al. (2021); Rong et al. (2023), etc. To transfer or hedge against longevity risk, the literature extensively discusses products related to longevity. The work of Chen et al. (2023) investigates optimal longevity risk transfer under asymmetric information using indemnity longevity swaps. Moreover, the market offers longevity caps, longevity floors and longevity-indexed bonds; see Bravo and Nunes (2021).
Longevity-indexed bonds are financial instruments linked to mortality rates, bridging financial and actuarial risks. Blake et al. (2006) thoroughly analyzed such assets, with Barbarin (2008) applying the Heath–Jarrow–Morton (HJM) methodology for a more realistic market model of longevity bond prices. Wills and Sherris (2010) focused on the structuring and pricing of longevity securities using mortgage and credit risk techniques. Other works on longevity bond can be found in Bauer et al. (2010); Leung et al. (2018).
Longevity-indexed bonds and other longevity-linked securities are widely studied in stochastic optimal consumption/portfolio decision problems. In order to maximize the terminal utility function for an economic agency, Menoncin (2008) investigated the role of a longevity-indexed bond in optimal portfolios and found that the value function in the market with a longevity bond is equal or greater than that without a longevity bond. By trading a longevity-linked security with continuous rate payments, Biagini et al. (2013) examined the mean-variance hedging for insurance products through longevity-linked securities. Further research on optimal strategies for longevity bonds is documented in various studies; see Kort and Vellekoop (2017); Menoncin and Regis (2017).
However, for pension plans, longevity-indexed bonds and other longevity-linked securities are not comprehensively studied in stochastic optimization problems. Usually, the pension sponsor allocates the fund in financial market consisting of several kinds of risky and risk-less assets; thus, it is important to consider the kind of assets to hold and the number of shares to buy. There are a few works that study the role of longevity bonds and longevity insurances in pension schemes. In DC pension schemes, Zhang (2018) studied optimal investment decisions in a longevity bond, implying the optimal expected utility in the market with a longevity bond is higher than that without it. Berstein and Morales (2021) examined the role of longevity insurances and found that longevity insurances could finance higher pension efficiently. In hybrid pension plans, Liu et al. (2023) examined a longevity bond and revealed that the optimal benefit payment with a longevity-indexed bond is higher than that without it. In DB pension schemes, Cox et al. (2013) investigated optimal investment decisions numerically, but not analytically, and found that longevity bonds can be used to hedge longevity risks, within a discrete time framework.
Our motivation is to find the closed-form investment decisions on longevity bonds within a continuous time framework and examine whether longevity bonds play a positive role or not in a DB pension plan. Theoretically, with the longevity risk increasing, the solvency risk also increases in most cases. As longevity bonds play a positive role and can hedge longevity risks in Zhang (2018); Liu et al. (2023); Cox et al. (2013), we think it is reasonable to suppose that investment in longevity bond may reduce solvency risk in DB pension plans.
This paper considers an optimization problem in an aggregated DB pension plan which pays close attention to solvency risk. It involves continuously adjusting investment weights within a market comprising a longevity bond and a bank account, aiming to minimize terminal solvency risk within a finite time frame.
Our contribution is three-fold. Theoretically, our work offers closed-form investment decisions on a longevity bond. Moreover, this research advances previous work (see Josa-Fombellida and Rincón-Zapatero 2010) by providing an explicit form of optimal solvency risk. Empirically, this work investigates the role of longevity bond in an aggregated DB pension plan. Practically, we extend the definitions of AL and NC in Josa-Fombellida and Rincón-Zapatero (2004, 2010) with the consideration of mortality rate, thereby deriving a new actualization rate.
The remainder of this study is organized as follows: Section 2 introduces the pension model and three assets within a stochastic interest rate framework. Section 3 identifies the optimal investment strategy and the minimum solvency risk with investment in a complete market consisting of a longevity-indexed bond and a risk-free asset. We derive the HJB equation using the dynamic programming approach. Section 4 examines the case replacing the indexed bond with a zero-coupon bond. Section 5 conducts a sensitivity analysis and compares the solvency risks in Section 3 and Section 4. Finally, Section 6 concludes the paper and suggests directions for further research.

2. Model Assumptions and Notations

We let ( Ω , F , F , P ) be a filtered probability space. Filtration F = F t t [ 0 , T 0 ] satisfies the usual conditions, i.e., F t t [ 0 , T 0 ] is right continuous and F 0 contains all the P -negligible events in F . Here, T 0 is a constant terminal time (which is shown later in Section 2.2). We let E be the expectation with respect to probability measure P , and E ( · | F t ) denotes conditional expectation with respect to probability measure P given by information available up to time t.

2.1. The Pension Model

In an aggregated DB pension plan, with the purpose of delivery of benefits to the affiliates after retirement time, the plan sponsor withdraws time-varying funds continuously. We use the following notations of elements for the pension plan:
  • A L ( t ) : Actuarial liability at time t, that is, the level of reserve to be detained by the pension fund and fixed by the chosen funding method, which is the liability part of the pension plan;
  • F ( t ) : Fund at time t, which is the market value of assets covering the actuarial liability (asset part);
  • U A L ( t ) : Unfunded actuarial liability at time t, which is the difference between the actuarial liability (reserve that we should have) and the fund (reserve that we really have). This amount can be positive (deficit or underfunding) or negative (surplus or overfunding);
  • P ( t ) : Benefits promised to affiliates at time t;
  • C ( t ) : Contribution rate at time t;
  • N C ( t ) : Normal cost at time t;
  • S C ( t ) : Supplementary contribution rate, which satisfies C ( t ) = N C ( t ) + S C ( t ) ;
  • M ( x ) : Percentage of promised benefits accumulated up to age x [ a , d ] , a , d R + ; hence, we have M ( a ) = 0 and M ( d ) = 1 . We assume that M is differentiable;
  • m ( x ) : The probability density of M;
  • δ ( t ) : Technical rate of actualization;
  • λ ( t ) : Force of mortality, which is defined by λ ( t ) = lim Δ t 0 P ( τ x < t + Δ t | τ x > t ) Δ t , where τ x is the future residual lifetime of a life who survives to age x.
Remark 1. 
Force of mortality λ is deterministic in the context. There are lots of different assumptions on the form of λ. For simplicity, λ has a form of ω ω t in Guan and Liang (2014), where ω is the largest survival age.
For analytical tractability, we have two assumptions following in the work of Josa-Fombellida and Rincón-Zapatero (2010) and the work of Liang et al. (2014):
Assumption 1. 
Benefit P follows a drifted Brownian Motion:
d P ( t ) = μ P P ( t ) d t + σ P d W P ( t ) ,
with initial liability P ( 0 ) = P 0 . μ P , σ P R + and W P is a standard Brownian Motion under probability measure P .
Assumption 2. 
Supplementary contribution is proportional to the unfunded actuarial liability, so the total contribution satisfies
C ( t ) = N C ( t ) + k ( A L ( t ) F ( t ) ) ,
where k [ 0 , 1 ] is a constant.
We extend the definitions of the actuarial functions A L and N C in Josa-Fombellida and Rincón-Zapatero (2004) (with deterministic interest rate) and Josa-Fombellida and Rincón-Zapatero (2010) (with stochastic interest rate), with consideration of mortality rate in Definition 1.
Definition 1. 
For every t 0 , A L ( t ) and N C ( t ) are defined as follows:
A L ( t ) = E a d e t t + d x δ ( s ) d s M ( x ) P ( t + d x ) d x p t d x | F t = E a d e t t + d x [ δ ( s ) + λ ( s ) ] d s M ( x ) P ( t + d x ) d x | F t ;
N C ( t ) = E a d e t t + d x δ ( s ) d s m ( x ) P ( t + d x ) d x p t d x | F t = E a d e t t + d x [ δ ( s ) + λ ( s ) ] d s m ( x ) P ( t + d x ) d x | F t ,
where   t p x is the probability that a life aged x will survive at least t years. A L ( t ) is the total expected benefits accumulated by percentage M ( x ) with consideration of mortality rate λ ( t ) , discounted at technical rate δ ( t ) . Analogous explanation can be given to N C ( t ) according to probability density function m ( x ) .
Remark 2. 
If we set λ = 0 in A L and N C in Definition 1, it degenerates into the case given by Josa-Fombellida and Rincón-Zapatero (2010) without consideration of mortality rate.
Using basic properties of conditional expectation, A L and N C in Definition 1 can be rewritten as
A L ( t ) = a d E e t t + d x [ δ ( s ) + λ ( s ) ] d s P ( t + d x ) | F t M ( x ) d x ;
N C ( t ) = a d E e t t + d x [ δ ( s ) + λ ( s ) ] d s P ( t + d x ) | F t m ( x ) d x .
The identities of A L and N C and the connection between them are given in Proposition 1.
Proposition 1. 
Actuarial functions A L and N C satisfy A L = ψ A L P and N C = ψ N C P , and the connection of A L and N C is given by
μ P δ ( t ) λ ( t ) + ξ A L ( t ) ψ A L ( t ) A L ( t ) N C ( t ) + P ( t ) = 0 ,
where ψ A L , ψ N C and ξ A L are the following stochastic processes:
ψ A L ( t ) = a d e t t + d x [ μ P δ ( u ) λ ( u ) ] d u M ( x ) d x ,
ψ N C ( t ) = a d e t t + d x [ μ P δ ( u ) λ ( u ) ] d u m ( x ) d x ,
ξ A L ( t ) = a d e t t + d x μ P δ ( u ) λ ( u ) d u μ P δ ( t + d x ) λ ( t + d x ) M ( x ) d x μ P δ ( t ) λ ( t ) ψ A L ( t ) .
Moreover, the actuarial liability satisfies
d A L ( t ) = μ P + ξ A L ( t ) ψ A L ( t ) A L ( t ) d t + σ P A L ( t ) d W P ( t ) ,
with A L ( 0 ) = ψ A L ( 0 ) P 0 A L 0 .
Proof of Proposition 1. 
The proof is similar to that in the work of Josa-Fombellida and Rincón-Zapatero (2010), so it is omitted here. □

2.2. The Financial Market

We assume that the nominal instantaneous interest rate follows the Vasicek model, which satisfies the stochastic differential equation (henceforth SDE):
d r ( t ) = ( a b r ( t ) ) d t + σ r d W r ( t ) ,
with initial value r ( 0 ) = r 0 , where W r is a standard Brownian Motion under P , and a, b and σ r are all positive constants. In the context of this work, we assume that W r is correlated with W P and the correlation coefficient is ρ , since inflation has effects on both interest rate and salary, and P ( t ) is usually related to the salary at the moment of retirement.
Now, we define W 0 ( t ) such that
d W P ( t ) = ρ ( t ) d W r ( t ) + 1 ρ 2 ( t ) d W 0 ( t ) ,
with W 0 ( 0 ) = 0 . Thus, W r and W 0 are independent Brownian motions transforming Equations (1) and (3) into the following forms:
d P ( t ) = μ P P ( t ) d t + σ P ρ ( t ) d W r ( t ) + 1 ρ 2 ( t ) d W 0 ( t ) ,
d A L ( t ) = μ P + ξ A L ( t ) ψ A L ( t ) A L ( t ) d t + σ P A L ( t ) ρ ( t ) d W r ( t ) + 1 ρ 2 ( t ) d W 0 ( t ) .
We assume that there are three continuously tradable and perfectly divisible underlying instruments in the financial market. Moreover, we suppose that borrowing and short-selling are permitted.
a. A risk-less asset S 0 ( t ) is modeled by
d S 0 ( t ) S 0 ( t ) = r ( t ) d t ,
with S 0 ( 0 ) = 1 .
b. A zero-coupon bond (or a T-bond) pays one dollar at its expiration time T 0 which can be thought of as a derivative of interest rate. Its value B ( t , T 0 ) at time t can be written as the expected present value of its future cash flow under the equivalent martingale measure Q , under which the present value process, e 0 t r ( u ) d u B ( t , T 0 ) , is a local martingale:
B ( t , T 0 ) = E Q e t T 0 r ( u ) d u | F t ,
where E Q denotes the expectation with respect to measure Q .
By the martingale property of the discounted process, e 0 t r ( u ) d u B ( t , T 0 ) , and following the work of Menoncin (2008), we obtain
d B ( t , T 0 ) = B ( t , T 0 ) r ( t ) d t + B ( t , T 0 ) r ( t ) σ r d W r Q ( t ) ,
with terminal condition B ( T 0 , T 0 ) = 1 , where W r Q ( t ) is a standard Brownian motion under measure Q . By the Girsanov theorem, we have d W r Q ( t ) = ζ r d t + d W r ( t ) , where ζ r is the market price for interest rate risk. The corresponding Radon–Nikodym derivative is Λ ( T 0 ) = d Q / d P = exp [ ζ r W r ( T 0 ) 1 2 ζ r 2 T 0 ] . Consequently, we have the SDE of B ( t , T 0 ) under physical measure P :
d B ( t , T 0 ) B ( t , T 0 ) = ( r ( t ) + r B σ r ζ r ) d t + r B σ r d W r ( t ) ,
where r B = B ( t , T 0 ) r ( t ) 1 B ( t , T 0 ) is the semi-elasticity of bond price B ( t , T 0 ) with respect to interest rate r. r B is negative because bond price usually negatively reacts to fluctuates in interest rate. It has an explicit form of e b ( T t ) 1 b under the Vasicek model given by Equation (4). Thus, ζ r is negative, because the bond has positive premium compared with the bank account.
Remark 3. 
The explicit formula of the zero-coupon bond process ( B ( t , T 0 ) ) t [ 0 , T 0 ] with Vasicek’s interest rates is B ( t , T 0 ) = A 1 ( t , T 0 ) e A 2 ( t , T 0 ) r ( t ) , with
A 1 ( t , T 0 ) = exp a b σ r 2 2 b 2 A 2 ( T 0 t ) ( T 0 t ) σ r 2 4 b A 2 ( t , T 0 ) 2 ,
and
A 2 ( t , T 0 ) = 1 b 1 e b ( T 0 t ) .
c. A longevity-indexed bond is defined with the same terminal time T 0 , from which the investor can obtain l T 0 / l t at the maturity time, where l t is the amount of people of a given population (whose ages are all x at initial time 0) who have survived from time 0 to time t. By the fundamental theorem of asset pricing, we obtain the price of the longevity bond at time t:
L ( t , T 0 ) = E Q l T 0 l t e t T 0 r ( u ) d u | F t = E Q T 0 t p x + t e t T 0 r ( u ) d u | F t = E Q e t T 0 λ ( u ) d u e t T 0 r ( u ) d u | F t = e t T 0 λ ( u ) d u B ( t , T 0 ) .
Following Zhang (2018), we have
d L ( t , T 0 ) L ( t , T 0 ) = λ ( t ) d t + d B ( t , T 0 ) B ( t , T 0 ) = ( λ ( t ) + r ( t ) + r B σ r ζ r ) d t + r B σ r d W r ( t )
with terminal condition L ( T 0 , T 0 ) = 1 . Thus, we obtain a correlation between the prices of a longevity bond and a zero-coupon bond.
Remark 4. 
If we set λ = 0 , the SDE of the longevity bond degenerates to the SDE of the zero-coupon bond; thus, the latter can serve as a special case of the former.
Remark 5. 
The market may also consist of stocks or other kind of assets, but the purpose of this work is to investigate the behavior of the longevity bond in solvency risk management, so it makes little sense for our problem.

3. DB Pension Fund Management with Investment of the Longevity Bond

In this section, we attempt to find the optimal proportion of fund assets put into a risk-less asset and a longevity-indexed bond in an aggregated DB pension scheme in order to minimize the solvency risk at terminal time T ( T T 0 ) .

3.1. Risk-Neutral Valuation of Liabilities

As mentioned in Josa-Fombellida and Rincón-Zapatero (2010), it is significant to derive a fair valuation of the liability by a technical rate of actualization, which can be served as a modification of the short interest rate. It is given in Proposition 2.
Proposition 2. 
The instantaneously technical rate of actualization δ L satisfies:
δ L ( t ) = r ( t ) + μ P ν L ( t ) λ ( t ) ,
where ν L ( t ) = μ P ζ r σ P ρ ( t ) λ ( t ) σ P ρ ( t ) r B σ r .
Proof of Proposition 2. 
The proof is similar to that of Josa-Fombellida and Rincón-Zapatero (2010), so it is omitted. □

3.2. The Wealth Process and Optimization Problem

The pension plan manager invests the fund wealth in a portfolio consisting of risk-less asset S 0 ( t ) , which is given in Equation (6), and longevity-indexed bond L ( t , T 0 ) , which is given in Equation (7), in which case the market is complete. The dynamics of returns on the fund process is
d F ( t ) = ω L ( t ) L ( t , T 0 ) d L ( t , T 0 ) L ( t , T 0 ) + F ( t ) ω L ( t ) L ( t , T 0 ) d S 0 ( t ) S 0 ( t ) + C ( t ) P ( t ) d t ,
with initial value F ( 0 ) = F 0 > 0 , where ω L ( t ) is the number of the longevity bond held in the portfolio at time t.
Substituting Equations (2), (6) and (7) into Equation (8), we have
d F ( t ) = [ F ( t ) r ( t ) + ω L ( t ) L ( t , T 0 ) λ ( t ) + r B σ r ζ r + ξ A L ( t ) ψ A L ( t ) + μ P δ L ( t ) λ ( t ) A L ( t )   + k A L ( t ) F ( t ) ] d t + ω L ( t ) L ( t , T 0 ) r B σ r d W r ( t ) ,
with consideration of Proposition 2.
Next, we restrict investment strategy ω L to fulfill some technical conditions:
Definition 2. 
Strategy ω L is called admissible if the following three conditions are satisfied:
a. ω L ( t ) is progressively measurable with respect to F t t [ 0 , T 0 ] ;
b. E 0 T ω L ( t ) L ( t , T 0 ) r B σ r 2 d t < ;
c. Equation (9) has a unique strong solution for the initial value 0 , F 0 , A L 0 , r 0 .
The set of all admissible controls is denoted by Π.
By setting X = F A L and inserting Equation (5) into d X , we have
d X ( t ) = ( r k ) X ( t ) + r δ L ( t ) λ ( t ) A L ( t ) + ω L ( t ) L ( t , T 0 ) λ ( t ) + r B σ r ζ r d t   + ω L ( t ) L ( t , T 0 ) r B σ r σ P ρ ( t ) A L ( t ) d W r ( t ) σ P 1 ρ 2 ( t ) A L ( t ) d W 0 ( t ) ,
with initial condition X ( 0 ) = F 0 A L 0 X 0 .
The preference of the pension fund sponsor is to minimize the terminal solvency risk. The objective function is given by
J ( ( X 0 , A L 0 , r 0 ) ; ω L ) = E X 0 , A L 0 , r 0 X 2 ( T ) .
The dynamic programming approach is used to solve the problem. We denote
V ( t , X , A L , r ) = min ω L Π J ( ( t , X , A L , r ) ; ω L ) :   subject   to ( 10 ) , ( 5 ) , ( 4 )
as the value function of this optimization problem, which is non-negative and strictly convex (which is shown in Appendix A). In optimal control problem, value function and optimal investment strategy can be obtained by solving the corresponding HJB equation; see Fleming and Soner (1993) or Yong and Zhou (1999):
min ω L Π Ψ ( ω L L ) = 0 ,
where
Ψ ( ω L L ) = V t + ( r k ) X + ( r δ L λ ) A L + ω L L ( λ + r B σ r ζ r ) V X + ξ A L ψ A L + μ P A L V A L + ( a b r ) V r + 1 2 ( ω L L r B σ r σ P ρ A L ) 2 + σ P 2 ( 1 ρ 2 ) A L 2 V X X + 1 2 σ P 2 A L 2 V A L , A L + 1 2 σ r 2 V r r + ω L L r B σ r σ P ρ A L σ P 2 A L 2 V X , A L + σ r ( ω L L r B σ r σ P ρ A L ) V X , r + σ r σ P ρ A L V A L , r = 0 ,
with
V ( T , X , A L , r ) = X 2 ,
where V t , V X , V A L , V r , V X X , V A L , A L , V r r , V X , A L , V X , r and V A L , r are the first- and second-order partial derivatives of V with respect to t, X, A L and r, respectively.
ω L L is the amount invested in the longevity-indexed bond. Then, min ω L Π Ψ ( ω L L ) can be obtained by utilizing the following two necessary conditions satisfied by the optimal amount ω L * L :
Ψ ( ω L * L ) = 0 , d Ψ d ( ω L L ) ( ω L * L ) = 0 .
The optimal investment amount ω L * L in terms of derivatives of the value function is given by
ω L * L = λ r B 2 σ r 2 + ζ r r B σ r V X V X X σ P ρ A L r B σ r V X , A L V X X 1 r B V X , r V X X + σ P ρ A L r B σ r .
Substituting Equation (14) into Equations (11) and (12), we finally obtain the explicit form of the optimal investment strategy and the minimized solvency risk. The results are shown in the following theorem.
Theorem 1. 
We denote the value function by V ( t , X , A L , r ) V L ( t , X , A L , r ) and let τ T . The minimum terminal solvency risk is given by
V L ( t , X , A L , r ) = g L ( t ) e γ L ( t ) r · X 2 + t T ε L ( t ; τ ) e η L ( t ; τ ) r d τ A L 2 .
The optimal investment strategy is
ω L * L ( t ) = λ r B 2 σ r 2 + ζ r r B σ r 2 b 1 e b ( T t ) 1 r B X + σ P ρ r B σ r A L ,
where
γ L ( t ) = 2 b ( 1 e b ( T t ) ) ; g L ( t ) = e t T h L ( s ) d s ; h L ( t ) = 2 k ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 + a 2 σ r λ r B σ r + ζ r γ L ( t ) 1 2 σ r 2 γ L 2 ( t ) ; η L ( t ; τ ) = γ L ( τ ) e b ( τ t ) ; ε L ( t ; τ ) = σ P 2 ( 1 ρ 2 ) g L ( τ ) e t τ y L ( s ; τ ) d s ; y L ( t ; τ ) = 2 ξ A L ψ A L + μ P + σ P 2 + a + 2 σ r σ P ρ η L ( t ; τ ) + 1 2 σ r 2 η L 2 ( t ; τ ) .
Proof. 
See Appendix A. □
Remark 6. 
Note that short-selling is permitted in the context. From Equation (16), ω L * L ( t ) might be negative. This may occur when X ( t ) is positive, i.e., the pension is overfunded. On one hand, longevity bonds can be used to hedge longevity risks according to Cox et al. (2013), but perhaps it is unnecessary to hedge this kind of risk when it is overfunded. On the other hand, a longevity bond is a kind of risky asset itself. When funds are much larger than liabilities, the pension sponsor may attempt to hold more bank account, say, more that 100% since there is no constraint in the market, thus the amount held in the longevity bond might be negative.
Substituting the optimal investment strategy ω L * L ( t ) into X ( t ) , we obtain the corresponding X * ( t ) . Here, we offer a possible path of ω L * L ( t ) and X * ( t ) in Figure 1 and Figure 2, respectively.

4. Special Case: DB Pension Fund Management with Investment of Zero-Coupon Bond

This section considers a similar optimal problem with an absent mortality rate, i.e., λ 0 , which makes the dynamic of the longevity-indexed bond degenerates to the dynamic of a normal zero-coupon bond. Compared with the work of Josa-Fombellida and Rincón-Zapatero (2010), we find the explicit form of the value function. Here, we abuse the notation F ( t ) to represent the pension wealth at time t. The market consisting of a zero-coupon bond and a bank account is again complete. As the counterpart of Proposition 2, we state the technical rate in this special case in Proposition 3 to define a fair liability.
Proposition 3. 
The technical rate of actualization is
δ B ( t ) = r ( t ) + μ P ν B ( t ) ,
where ν B ( t ) = μ P ζ r σ P ρ ( t ) .
Similarly, the dynamic of the wealth process is
d F ( t ) = ω B ( t ) B ( t , T 0 ) d B ( t , T 0 ) B ( t , T 0 ) + F ( t ) ω B ( t ) B ( t , T 0 ) d S 0 ( t ) S 0 ( t ) + C ( t ) P ( t ) d t = r ( t ) k F ( t ) A L ( t ) + r ( t ) + ξ A L B ( t ) ψ A L B ( t ) + μ P δ ( t ) A L ( t ) + ω B ( t ) B ( t , T 0 ) r B σ r ζ r d t   + ω B ( t ) B ( t , T 0 ) r B σ r d W r ( t ) ,
with initial value F ( 0 ) = F 0 > 0 , where ω B is the number of the zero-coupon bond held in the portfolio. Here,
ψ A L B ( t ) = a d e t t + d x [ μ P δ ( u ) ] d u M ( x ) d x , ξ A L B ( t ) = a d e t t + d x μ P δ ( u ) d u μ P δ ( t + d x ) M ( x ) d x μ P δ ( t ) ψ A L B ( t ) .
A strategy is called admissible if it satisfies similar conditions in Definition 2, and the set of all admissible controls is again denoted by Π .
Setting X = F A L , we have
d X ( t ) = ( r k ) X ( t ) + r δ ( t ) A L ( t ) + ω B ( t ) B ( t , T 0 ) r B σ r ζ r d t   + ω B ( t ) B ( t , T 0 ) r B σ r σ P ρ ( t ) A L ( t ) d W r ( t ) σ P 1 ρ 2 ( t ) A L ( t ) d W 0 ( t ) ,
with X ( 0 ) = X 0 .
Here, the objective function and the optimal problem are the same as those in Section 3, except that the longevity-indexed bond is changed into a normal zero-coupon bond. Analogously, denoting V B ( t , X , A L , r ) as the value function, we offer the optimal results in Theorem 2 without any proof.
Theorem 2. 
We set τ T and denote the value function in Section 4 by V B ( t , X , A L , r ) . The minimum terminal solvency risk is given by
V B ( t , X , A L , r ) = g B ( t ) e γ B ( t ) r · X 2 + t T ε B ( t ; τ ) e η B ( t ; τ ) r d τ A L 2 .
The optimal investment strategy is
ω B * B ( t ) = ζ r r B σ r 2 b 1 e b ( T t ) 1 r B X + σ P ρ r B σ r A L ,
where
γ B ( t ) = 2 b ( 1 e b ( T t ) ) = γ L ( t ) γ ( t ) ; g B ( t ) = e t T h B ( s ) d s ; h B ( t ) = 2 k ζ r 2 + ( a 2 σ r ζ r ) γ ( t ) 1 2 σ r 2 γ 2 ( t ) ; η B ( t ; τ ) = γ B ( τ ) e b ( τ t ) = η L ( t ; τ ) η ( t ; τ ) ; ε B ( t ; τ ) = σ P 2 ( 1 ρ 2 ) g B ( τ ) e t τ y B ( s ; τ ) d s ; y B ( t ; τ ) = 2 ξ A L B ψ A L B + μ P + σ P 2 + ( a + 2 σ r σ P ρ ) η ( t ; τ ) + 1 2 σ r 2 η 2 ( t ; τ ) .

5. Sensitivity Analysis

This section investigates the influence of the employed parameters of the model on optimal control and value function. We mainly explore the impacts of mortality rate λ and terminal time T on optimal investment decision ω L * L and value function V L .
In the pension model, k [ 0 , 1 ] represents the surplus (or deficit) amortized rate. Here, we choose 0.5 . Parameter σ P ( > 0 ) is the volatility of benefit, which is 8 % in Josa-Fombellida and Rincón-Zapatero (2010). Here, σ P = 30 % , since we try to investigate the influence of the employed parameters under relative high risks. Due to the matching of risks and returns, i.e., appreciation rate growing with volatility increasing, we set μ P = 0.1 , which is also a bit higher than 0.04 in Josa-Fombellida and Rincón-Zapatero (2010). Parameter ρ is the correlation coefficient between W P and W r . In the market, if inflation rate grows, both the nominal interest rate and salary grow, and the promised benefits also grow, since benefits are usually positively correlated with the salary at the moment of retirement. Thus, it is rational to assume that ρ [ 0 , 1 ] and here we choose ρ = 0.5 .
In the interest rate model, σ r > 0 is the volatility of interest rate; here, we choose 30 % since we attempt to conduct the sensitivity analysis under high risks. Parameter b represents the mean reversion coefficient, while a / b is the long-run mean of interest rate. In the work of Han and Hung (2017), the speed of mean-reversion is 0.0395 and the long-run mean is 0.0369. Here, we choose almost the same speed of mean-reversion, which is 0.03, but a bit higher long-run mean, which is 0.05, since our volatility is also higher than that in Han and Hung (2017). Parameter ζ r is the market price of interest rate risk, i.e., the bond’s excess return divided by volatility, which should be negative. Theoretically, bonds’ return rate is usually higher than risk-less return rate, for example, 5 % , and bonds’ volatility should not excess stocks’ volatility. So it is reasonable to assume that, for instance, ζ r [ 0.01 , 2 ] . The value of ζ r is 0.46 in Menoncin (2008). Here, we choose −0.5 without loss of generality. The parameters are given in Table 1.
When the mortality rate is high, the longevity risk is usually low. Likewise, ω L * L might be low because the sponsor may not need a lot of longevity bonds to hedge longevity risk. From another aspect, when the mortality rate is high, the price of longevity bond is low, thus ω L * L might be low, too. To test this point of view, Figure 3 demonstrates the effect of mortality rate λ on optimal investment strategy ω L * L . When the mortality rate becomes larger, the investment amount becomes smaller. Figure 3 also shows that with the same λ , ω L * L becomes larger with a longer T. A rational explanation is that the DB pension scheme faces higher risks in the long run, so investment amount in the risky bond (in a short position) becomes smaller in order to reduce risk, which is more clear in Figure 4.
If longevity risk is hedged by the longevity-indexed bond, we may have a lower solvency risk with a lower mortality rate. Figure 5 demonstrates the effect of mortality rate λ on value function V L , and it shows that the pension scheme has a lower value function (which means lower solvency risk) with a lower mortality rate (which may cause higher longevity risk). It seems that investments in the longevity bond reduce the solvency risk in a DB pension plan. In addition, Figure 5 also shows that the solvency risk is lower with a smaller λ in a longer period. It shows that the hedging strategy is more efficient in the long run.
In addition, if longevity risk is reduced by the longevity bond, we may have a lower solvency risk in the market with the longevity bond than that with a normal zero-coupon bond. V L and V B denote the value function in Section 3 (with investment in the longevity bond) and Section 4 (with investment in the zero-coupon bond), respectively. To test this, Figure 6 demonstrates the performance of V L / V B with different λ values. The ratio V L / V B is lower than 1 with λ from 0.1 to 0.8 , implying the solvency risk in Section 3 is lower than that in Section 4.

6. Conclusions

This paper investigates, by means of dynamic programming technique, an optimal control problem for an aggregated DB pension plan in a stochastic interest rate framework. The pension sponsor’s objective is to minimize the terminal solvency risk by managing the fund wealth in finite time horizon through a portfolio comprising a longevity-indexed bond and a risk-less asset. The explicit investment policy is derived by solving the corresponding HJB equation. By replacing the longevity-indexed bond with a normal zero-coupon bond, we find that the value function in the market with a longevity bond is lower than that with a zero-coupon bond, i.e., a longevity bond plays a positive role and has a strategic advantage in hedging solvency risk in such DB pension scheme.
It may be significant to consider the benefit process with the influence of mortality rate in the future, which may generate more realistic sample paths of the benefit. Mortality rate may be stochastic. It is also interesting to consider constrained control variable ω 0 , i.e., short-selling is prohibited, which is consistent with some regulations and risk management policies. This makes the value function unsmooth; thus, the technique of viscosity solutions in the dynamic programming approach might be applied. We leave these two points for future study.

Author Contributions

Conceptualization, J.G.; methodology, J.G. and X.Z.; software, Y.L.; writing—original draft preparation, X.Z.; writing—review and editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 12201174, No. 11931018, No. 12271274 and No. 12371468) and the Science and Technology Foundation of Universities in Hebei Province (No. QN2021215).

Data Availability Statement

No new data were created or analyzed in this study.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DCDefined contribution
DBDefined benefit
HJBHamilton–Jacobi–Bellman
SDEStochastic differential equation
PDEPartial differential equation

Appendix A

The proof of Theorem 1. 
Substituting Equation (14) into Equations (11) and (12), we determine that the value function V satisfies
V t + ( r k ) X + ( r δ L λ ) A L + σ P ρ r B σ r ( λ + r B σ r ζ r ) A L V X + ξ A L ψ A L + μ P A L V A L + ( a b r ) V r + 1 2 σ P 2 ( 1 ρ 2 ) A L 2 V X X + 1 2 σ P 2 A L 2 V A L , A L + 1 2 σ r 2 V r r σ P 2 ( 1 ρ 2 ) A L 2 V X , A L + σ r σ P ρ A L V A L , r 1 2 ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 V X 2 V X X σ P ρ λ r B σ r + ζ r A L V X V X , A L V X X 1 2 σ P 2 ρ 2 A L 2 V X , A L 2 V X X σ r λ r B σ r + ζ r V X V X r V X X 1 2 σ r 2 V X r 2 V X X σ r σ P ρ A L V X , A L V X r V X X = 0 ,
with the final value given in Equation (13). The structure of the HJB equation suggests a quadratic form as follows:
V ( t , X , A L , r ) = f X X ( t , r ) X 2 + f A L , A L ( t , r ) A L 2 + f X , A L ( t , r ) X A L ,
where f X X ( · , · ) , f A L , A L ( · , · ) and f X , A L ( · , · ) are undetermined functions with terminal values f X X ( T , r ) = 1 , f A L , A L ( T , r ) = 0 and f X , A L ( T , r ) = 0 .
Differentiating Equation (A2), we obtain
V t = ( f X X ) t X 2 + ( f A L , A L ) t A L 2 + ( f X , A L ) t X A L , V X = 2 f X X X + f X , A L A L , V A L = 2 f A L , A L A L + f X , A L X , V r = ( f X X ) r X 2 + ( f A L , A L ) r A L 2 + ( f X , A L ) r X A L , V X X = 2 f X X , V A L , A L = 2 f A L , A L , V X , A L = f X , A L , V X , r = 2 ( f X X ) r X + ( f X , A L ) r A L , V A L , r = 2 ( f A L , A L ) r A L + ( f X , A L ) r X , V r r = ( f X X ) r r X 2 + ( f A L , A L ) r r A L 2 + ( f X , A L ) r r X A L ,
where ( f X X ) t , ( f X X ) r , ( f X X ) r r are the first- and second-order partial derivatives of f X X with respect to t and r, respectively. Other partial derivatives in Equation (A3) are defined similarly.
Substituting Equation (A3) into Equation (A1) and rearranging it by the order of X 2 , A L 2 and X A L , three PDEs are obtained together with boundary conditions:
{ ( f X X ) t + 2 ( r k ) ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 f X X + ( a b r ) 2 σ r λ r B σ r + ζ r ( f X X ) r + 1 2 σ r 2 ( f X X ) r r σ r 2 ( f X X ) r 2 f X X = 0 , f X X ( T , r ) = 1 ,
{ ( f A L , A L ) t + 2 ξ A L ψ A L + μ P + σ P 2 f A L , A L + ( a b r ) + 2 σ r σ P ρ ( f A L , A L ) r + 1 2 σ r 2 ( f A L , A L ) r r + σ P 2 ( 1 ρ 2 ) f X X = 0 , f A L , A L ( T , r ) = 0 ,
{ ( f X , A L ) t + ( r k ) ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 + ξ A L ψ A L + μ P σ P ρ r B σ r ( λ + r B σ r ζ r ) f X , A L + [ ( a b r ) σ r λ r B σ r + ζ r σ r σ P ρ ] ( f X , A L ) r + 1 2 σ r 2 ( f X , A L ) r r σ r λ r B σ r + ζ r σ r σ P ρ f X , A L ( f X X ) r f X X σ r 2 ( f X , A L ) r ( f X X ) r f X X + 2 ( r δ L λ ) + σ P ρ r B σ r ( λ + r B σ r ζ r ) f X X = 0 , f X , A L ( T , r ) = 0 ,
since Equation (A1) satisfies for every X 2 , A L 2 and X A L .
To solve Equation (A4), we conjecture that f X X ( t , r ) = g ( t ) e γ ( t ) r , with terminal conditions g ( T ) = 1 and γ ( T ) = 0 . So we have ( f X X ) t = [ g ( t ) / g ( t ) + r γ ( t ) ] f X X , ( f X X ) r = γ ( t ) f X X and ( f X X ) r r = γ 2 ( t ) f X X . Substituting ( f X X ) t , ( f X X ) r and ( f X X ) r r into Equation (A4), we obtain the following equation:
r γ ( t ) b γ ( t ) + 2 + g ( t ) g ( t ) 2 k ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 + a 2 σ r λ r B σ r + ζ r γ ( t ) 1 2 σ r 2 γ 2 ( t ) = 0 .
From Equation (A7), we have
γ ( t ) b γ ( t ) + 2 = 0 ,
and
g ( t ) g ( t ) 2 k ( λ + r B σ r ζ r ) 2 r B 2 σ r 2 + a 2 σ r λ r B σ r + ζ r γ ( t ) 1 2 σ r 2 γ 2 ( t ) = 0 .
For convenience, we denote g = g L and γ = γ L in Theorem 1 in order to draw comparison with the results given in Section 4. g L ( t ) and γ L ( t ) are given in Equation (17).
The method solving Equation (A4) cannot be used to find an analytical solution of Equation (A5) due to the nonhomogeneous term σ P 2 ( 1 ρ 2 ) f X X . By Proposition 2 in the work of Yao et al. (2013), we can derive the solution of Equation (A5) by solving the associated homogeneous PDE in the following:
{ υ t ( t , r ; τ ) + 2 ξ A L ψ A L + μ P + σ P 2 υ ( t , r ; τ ) + ( ( a b r ) + 2 σ r σ P ρ ) υ r ( t , r ; τ ) + 1 2 σ r 2 υ r r ( t , r ; τ ) = 0 , υ ( τ , r ; τ ) = σ P 2 ( 1 ρ 2 ) g ( τ ) e γ ( τ ) r ,
with τ T ; then, the analytical solution of Equation (A5) is
f A L , A L ( t , r ) = t T υ ( t , r ; τ ) d τ .
The analytical solution of Equation (A8) can be derived easily by using the same method as in solving f X X ( t , r ) . We try solution υ ( t , r ; τ ) = ε ( t ; τ ) e η ( t ; τ ) r , with terminal values ε ( τ ; τ ) = σ P 2 ( 1 ρ 2 ) g ( τ ) and η ( τ ; τ ) = γ ( τ ) . ε and η are obtained after substituting the first- and second-order derivatives of υ with respect to t and r into Equation (A8). We denote ε = ε L and η = η L , in order to make comparison with results given in Section 4. ε L ( t ) and η L ( t ) are given in Equation (17).
Using Proposition 2, we find that the first equation in Equation (A6) is homogeneous. So it is easy to prove that f X , A L = 0 , since its terminal condition is 0.
It is obvious that V X X = 2 f X X ( t , r ) = g L ( t ) e γ L ( t ) r > 0 . The optimal investment strategy ω L * L given in Equation (16) is obtained by substituting the derivatives of the value function into Equation (14). □

References

  1. Barbarin, Jérôme. 2008. Heath-Jarrow-Morton modelling of longevity bonds and the risk minimization of life insurance portfolios. Insurance: Mathematics and Economics 43: 41–55. [Google Scholar] [CrossRef]
  2. Bauer, Daniel, Matthias Börger, and Jochen Ruß. 2010. On the pricing of longevity-linked securities. Insurance: Mathematics and Economics 46: 139–49. [Google Scholar] [CrossRef]
  3. Berstein, Solange, and Marco Morales. 2021. The role of a longevity insurance for defined contribution pension systems. Insurance: Mathematics and Economics 76: 233–40. [Google Scholar] [CrossRef]
  4. Biagini, Francesca, Thorsten Rheinländer, and Jan Widenmann. 2013. Hedging mortality claims with longevity bonds. ASTIN Bulletin 43: 123–57. [Google Scholar] [CrossRef]
  5. Blake, David, and Andrew J. G. Cairns. 2021. Longevity risk and capital markets: The 2019–2020 update. Insurance: Mathematics and Economics 99: 395–439. [Google Scholar]
  6. Blake, David, Andrew J. G. Cairns, and Kevin Dowd. 2006. Living with mortality: Longevity bonds and other mortality-linked securities. British Actuarial Journal 12: 153–97. [Google Scholar] [CrossRef]
  7. Bravo, Jorge M., and João Pedro Vidal Nunes. 2021. Pricing longevity derivatives via Fourier transforms. Insurance: Mathematics and Economics 96: 81–97. [Google Scholar] [CrossRef]
  8. Broeders, Dirk, Roel Mehlkopf, and Annick van Ool. 2021. The economics of sharing macro-longevity risk. Insurance: Mathematics and Economics 99: 440–58. [Google Scholar] [CrossRef]
  9. Chen, An, Hong Li, and Mark B. Schultze. 2023. Optimal longevity risk transfer under asymmetric information. Economic Modelling 120: 106179. [Google Scholar] [CrossRef]
  10. Cox, Samuel H., Yijia Lin, Ruilin Tian, and Jifeng Yu. 2013. Managing capital market and longevity risks in a defined benefit pension plan. The Journal of Risk and Insurance 80: 585–619. [Google Scholar] [CrossRef]
  11. Fleming, Wendell H., and Halil Mete Soner. 1993. Controlled Markov Processed and Viscosity Solutions. New York: Springer. [Google Scholar]
  12. Guan, Guohui, and Zongxia Liang. 2014. Optimal management of DC pension plan in a stochastic interest rate and stochastic volatility framework. Insurance: Mathematics and Economics 57: 58–66. [Google Scholar] [CrossRef]
  13. Guan, Guohui, Jiaqi Hu, and Zongxia Liang. 2022. Robust equilibrium strategies in a defined benefit pension plan game. Insurance: Mathematics and Economics 106: 193–217. [Google Scholar] [CrossRef]
  14. Han, Nan-Wei, and Mao-Wei Hung. 2017. Optimal consumption, portfolio, and life insurance policies under interest rate and inflation risks. Insurance: Mathematics and Economics 73: 54–67. [Google Scholar] [CrossRef]
  15. Josa-Fombellida, Ricardo, and Juan Pablo Rincón-Zapatero. 2004. Optimal risk management in defined benefit stochastic pension funds. Insurance: Mathematics and Economics 34: 489–503. [Google Scholar] [CrossRef]
  16. Josa-Fombellida, Ricardo, and Juan Pablo Rincón-Zapatero. 2010. Optimal asset allocation for aggregated defined benefit pension funds with stochastic interest rates. European Journal of Operational Research 201: 211–21. [Google Scholar] [CrossRef]
  17. Kort, Jan, and Michel H. Vellekoop. 2017. Existence of optimal consumption strategies in markets with longevity risk. Insurance: Mathematics and Economics 72: 107–21. [Google Scholar]
  18. Leung, Melvern, Man Chung Fung, and Colin O’Hare. 2018. A comparative study of pricing approaches for longevity instruments. Insurance: Mathematics and Economics 82: 95–116. [Google Scholar] [CrossRef]
  19. Liang, Xiaoqing, Lihua Bai, and Junyi Guo. 2014. Optimal time-consistent portfolio and contribution selection for defined benefit pension schemes under mean-variance criterion. The ANZIAM Journal 56: 66–90. [Google Scholar] [CrossRef]
  20. Liu, Zilan, Huanying Zhang, and Lei He. 2023. Optimal assets allocation and benefit adjustment strategy with longevity risk for target benefit pension plans. Journal of Industrial and Management Optimization 19: 3931–51. [Google Scholar] [CrossRef]
  21. Menoncin, Francesco. 2008. The role of longevity bonds in optimal portfolios. Insurance: Mathematics and Economics 42: 343–58. [Google Scholar] [CrossRef]
  22. Menoncin, Francesco, and Luca Regis. 2017. Longevity-linked assets and pre-retirement consumption/portfolio decisions. Insurance: Mathematics and Economics 76: 75–86. [Google Scholar] [CrossRef]
  23. Ng, Kenneth Tsz Hin, and Wing Fung Chong. 2023. Optimal investment in defined contribution pension schemes with forward utility preferences. Insurance: Mathematics and Economics 114: 192–211. [Google Scholar] [CrossRef]
  24. Rong, Ximin, Cheng Tao, and Hui Zhao. 2023. Target benefit pension plan with longevity risk and intergenerational equity. ASTIN Bulletin 53: 84–103. [Google Scholar] [CrossRef]
  25. Wills, Samuel, and Michael Sherris. 2010. Securitization, structuring and pricing of longevity risk. Insurance: Mathematics and Economics 46: 173–85. [Google Scholar] [CrossRef]
  26. Yao, Haixiang, Zhou Yang, and Ping Chen. 2013. Markowitz’s mean-variance defined contribution pension fund management under inflation: A continuous-time model. Insurance: Mathematics and Economics 53: 851–63. [Google Scholar] [CrossRef]
  27. Yong, Jiongmin, and Xun Yu Zhou. 1999. Stochastic Controls: Hamiltonian Systems and HJB Equations. New York: Springer. [Google Scholar]
  28. Zhang, Xiaoyi. 2018. The role of longevity bond in DC pension plan during both accumulation and decumulation phases. Chinese Journal of Engineering Mathematics 37: 347–69. [Google Scholar]
  29. Zhu, Nan, and Daniel Bauer. 2014. A cautionary note on natural hedging of longevity risk. North American Actuarial Journal 18: 104–15. [Google Scholar] [CrossRef]
Figure 1. A possible path of the optimal investment strategy ω L * L ( t ) .
Figure 1. A possible path of the optimal investment strategy ω L * L ( t ) .
Risks 12 00049 g001
Figure 2. A possible path of X * ( t ) .
Figure 2. A possible path of X * ( t ) .
Risks 12 00049 g002
Figure 3. This figure plots the relationship between mortality rate λ and optimal investment strategy ω L * L with different T.
Figure 3. This figure plots the relationship between mortality rate λ and optimal investment strategy ω L * L with different T.
Risks 12 00049 g003
Figure 4. This figure plots the relationship between terminal time T and optimal investment strategy ω L * L with λ = 0.1 .
Figure 4. This figure plots the relationship between terminal time T and optimal investment strategy ω L * L with λ = 0.1 .
Risks 12 00049 g004
Figure 5. This figure plots the relationship between mortality rate λ and value function V L with different terminal time T.
Figure 5. This figure plots the relationship between mortality rate λ and value function V L with different terminal time T.
Risks 12 00049 g005
Figure 6. This figure plots the relationship between mortality rate λ and ratio V L / V B with different terminal time T.
Figure 6. This figure plots the relationship between mortality rate λ and ratio V L / V B with different terminal time T.
Risks 12 00049 g006
Table 1. Parameters of the model.
Table 1. Parameters of the model.
X 0 AL 0 κ μ P σ P ρ ζ r σ r ab r 0
110.50.10.30.5−0.50.30.00150.030.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, X.; Li, Y.; Guo, J. The Role of Longevity-Indexed Bond in Risk Management of Aggregated Defined Benefit Pension Scheme. Risks 2024, 12, 49. https://doi.org/10.3390/risks12030049

AMA Style

Zhang X, Li Y, Guo J. The Role of Longevity-Indexed Bond in Risk Management of Aggregated Defined Benefit Pension Scheme. Risks. 2024; 12(3):49. https://doi.org/10.3390/risks12030049

Chicago/Turabian Style

Zhang, Xiaoyi, Yanan Li, and Junyi Guo. 2024. "The Role of Longevity-Indexed Bond in Risk Management of Aggregated Defined Benefit Pension Scheme" Risks 12, no. 3: 49. https://doi.org/10.3390/risks12030049

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop