Journal of Applied Mathematics and Physics
Vol.05 No.12(2017), Article ID:81043,9 pages
10.4236/jamp.2017.512190

Closed-Form Absolute Ruin Problems of the Risk Models with State-Dependent Switched Claims

Yuanxun Liu, Dianli Zhao*

College of Science, University of Shanghai for Science and Technology, Shanghai, China

Copyright © 2017 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: November 14, 2017; Accepted: December 11, 2017; Published: December 14, 2017

ABSTRACT

This letter mainly investigates a general risk model with the threshold dividend strategy under assumption that the claim amounts obey a state-dependent switched exponential distribution. By establishing the differential-integral equations for the Gerber-Shiu discounted penalty function, and applying the hypergeometric functions, the closed-form absolute ruin probability is derived.

Keywords:

Ruin Probability, Threshold Dividend Strategy, Switched Exponential Distribution, Hypergeometric Function

1. Introduction

Due to its practical importance, more and more people use the ruin probability, an important leading indicator, to judge whether an insurance company can survive or not. The classical risk model was introduced by Gerber [1] . The function is of the form

U ( t ) = u 0 + c t S ( t ) , S ( t ) = i = 0 N ( t ) X i , (1)

where U ( t ) is the insurance company’s surplus at time t. c represents the risk-free rate. S ( t ) is the total compensation by time t; X i is the i-th compensation; N ( t ) means the frequency of compensation. u 0 represents the initial investment for an insurance company. In order to make the classical risk model more realistic, more and more factors have been taken into account. The aggregate premium process taking a linear function of time was considered by Zou, Gao and Xie [2] , in which the explicit solutions of ruin time were obtained under assumption that the claim amount obeyed exponential distribution. Yao, Yang and Wang [3] studied a dual risk model with fixed transaction cost, and solved the optional control problem by the techniques of quasi-variational inequalities. Yuen, Zhou and Guo [4] discussed the compound Poisson risk model with debit interest and dividend payments, and the explicit expressions for the Gerber-Shiu function and the optimal barrier with exponential claim amounts were given. When the dividend strategy was considered, Yin and Yuen [5] considered a dividend policy aiming to maximize the expected discounted value of dividend until ruin. Taking the threshold dividend strategy into account, two integro-differential equations for the Gerber-Shiu discounted penalty function were obtained by Lin and Pavlova [6] . The ruin model with multiple thresholds was discussed by Lin and Sendova [7] . For more about the related risk model, we can refer to [8] [9] [10] [11] [12] and cited there in. From all the literatures mentioned above, they have a common assumption that the claim sizes obey exponential distribution, i.e. f ( x ) = η e η x . However, in fact, different dividend levels always reflect the difference of the economic environment, which results in the changes of the claim sizes. In detail, when the surplus u is more than the threshold b, the claim sizes obey an exponential distribution where the strength is η 1 ; while u < b , the claim sizes obey an exponential distribution with strength η 2 . Therefore, the exponential distribution is noted like

f ( x ) = η 1 e η 1 x I { u b } + η 2 e η 2 x I { u < b } . (2)

I { } is an indicator function, such that I { u b } = 1 and I { u < b } = 0 if u b . In order to make it more general, we set the risk model with the claim sizes obeying a switched exponential distribution as follows

f ( x ) = p η 1 e η 1 x I { u b } + q η 2 e η 2 x I { u < b } , (3)

where p and q are two suitable positive constants to be chosen later. Clearly, when we choose p = q = 1 and η 1 = η 2 , (3) is the standard exponential distribution applied in [2] .

The main aim of this letter is to compute the ruin probability under assumption that the claim sizes obey a switched exponential distribution (3).

The rest of this letter is organized as follows. In Section 2, we firstly get the Gerber-Shiu discounted penalty functions with different surplus. Then, the ruin probability is calculated by assuming the claim sizes obey a switched exponential distribution.

2. Closed Form Absolute Ruin Probability

In this letter, we consider the risk model with constant interest rate r > 0 and debit interest ρ . Besides, the studied risk model includes a general threshold dividend strategy. When the surplus is over threshold dividend b, we assume the insurance company gets insurance premium at a constant rate c 1 and earns interest at a constant rate r. If the surplus is between zero and b, it will collect insurance premium at a constant rate c 2 and earns interest at a constant rate r.

When the surplus is between c 2 ρ and zero, the insurer can borrow an amount

of money equal to the deficit at a debit force ρ . Meanwhile, the insurer will repay the debts continuously from his premium income. We denote the surplus of the insurer at time t with the credit interest r and debit interest ρ by U ( t ) which is the solution to

d U ( t ) = ( c 1 d t + r U ( t ) d t d S ( t ) ; u b , c 2 d t + r U ( t ) d t d S ( t ) ; 0 u < b , c 2 d t + ρ U ( t ) d t d S ( t ) ; c 2 ρ < u < 0 , (4)

where S ( t ) = i = 0 N ( t ) X i is the cumulative damages in time interval [ 0 , t ] , c (>0)

represents the constant rate of premium, { N ( t ) , t 0 } a Poisson process with intensity λ > 0 , which counts the claim numbers in the interval [ 0 , t ] , and { X i , i 1 } (representing the size of claims and independent of { N ( t ) , t 0 } ) is a sequence of independent and identically distributed nonnegative variables with common distribution function F ( x ) that satisfies F ( 0 ) = 0 and has a positive mean μ = 0 x d F ( x ) . Let ( Ω , F , { F t } t 0 , P ) be a filtered probability space containing all processes and random variables in this letter, satisfying the usual conditions, i.e. { F t } t 0 is right continuous and P-complete. In order to obtain the expression of the ruin probability, we should work for two things. First of all, we should get the Gerber-Shiu discounted penalty function. Then, we calculate the ruin probability. So, let’s firstly define the Gerber-Shiu discounted penalty function φ ( u , b ) .

φ ( u , b ) = E { e β T W ( U b ( T ) , | U b ( T ) | ) R ( T < ) | U b ( 0 ) = u } , (5)

where T is the time of bankrupting, R ( T < ) is an indicator function, which means R ( T < ) = 1 , if the time of bankrupting is finite value. Else if

R ( T < ) = 0 . W ( x , y ) is a non-negative measurable function satisfying

W ( x , y ) ( c 2 ρ , ) × [ c 2 ρ , ) . U b ( T ) is the surplus immediately before the

company bankrupts. Next, we give the integro-differential equations for φ i ( u , b ) ( i = 1 , 2 , 3 ) .

Theorem 2.1. Let φ i ( u , b ) ( i = 1 , 2 , 3 ) be the Gerber-Shiu discounted penalty function by the surplus u and the threshold b for model (4). Then

( r u + c 1 ) φ 1 ( u , b ) ( λ + β ) φ 1 ( u , b ) + λ 0 u b φ 1 ( u x , b ) d F ( x ) + λ u b u φ 2 ( u x , b ) d F ( x ) + λ u u + c 2 ρ φ 3 ( u x , b ) d F ( x ) + λ u + c 2 ρ W ( u , x u ) d F ( x ) = 0, ( u b ) ; (6)

( r u + c 2 ) φ 2 ( u , b ) ( λ + β ) φ 2 ( u , b ) + λ 0 u φ 2 ( u x , b ) d F ( x ) + λ u u + c 2 ρ φ 3 ( u x , b ) d F ( x ) + λ u + c 2 ρ W ( u , x u ) d F ( x ) = 0 , ( 0 u < b ) ; (7)

( ρ u + c 2 ) φ 3 ( u , b ) ( λ + β ) φ 3 ( u , b ) + λ u u + c 2 ρ φ 3 ( u x , b ) d F ( x ) + λ u + c 2 ρ W ( u , x u ) d F ( x ) = 0 , ( c 2 ρ < u < 0 ) . (8)

the boundary conditions are given like

φ 1 ( b , b ) = φ 2 ( b , b ) , φ 2 ( 0 , b ) = φ 3 ( 0 , b ) , ( r b + c 2 ) φ 2 ( b , b ) = ( r b + c 1 ) φ 1 ( b , b ) ,

φ 3 ( 0 , b ) = φ 2 ( 0 , b ) , φ 3 ( c 2 ρ , b ) = W ( c 2 ρ , c 2 ρ ) .

Proof. When u b , define the cumulative value

h = u e r t + c 1 0 t e r ( t s ) d s = u e r t + c 1 r ( e r t 1 ) . (9)

Note that the Gerber-Shiu discounted penalty functions have different expression when the surplus lies within different ranges, we discuss that with the following four cases:

1) When u b ( t ) c 2 ρ , the amount of compensation is X = h U b ( t ) . Let φ 3 ( c 2 ρ , b ) = 1 , its counterpart in the Gerber-Shiu discounted penalty function is λ t h + c 2 ρ W ( h , h x ) d F ( x ) ;

2) When c 2 ρ < u b ( t ) < 0 , λ t h h + c 2 ρ φ 3 ( h x , b ) d F ( x ) ;

3) When 0 u b ( t ) < b , λ t h b h φ 2 ( h x , b ) d F ( x ) ;

4) When u b ( t ) b , λ t 0 h b φ 1 ( h x , b ) d F ( x ) .

Therefore, the Gerber-Shiu discounted penalty function is given like,

φ 1 ( u , b ) = λ t h + c 2 ρ W ( h , h x ) d F ( x ) + λ t h h + c 2 ρ φ 3 ( h x , b ) d F ( x ) + λ t h b h φ 2 ( h x , b ) d F ( x ) + λ t 0 h b φ 1 ( h x , b ) d F ( x ) + ( 1 λ t ) φ 1 ( h , b ) + o ( t ) . (10)

For φ 1 ( h , b ) , by using Taylor expansion, we get

φ 1 ( h , b ) = φ 1 ( u , b ) + ( 2 φ 1 ( h , b ) u 2 2 u t 2 ) | ( u , b ) t 2 + ( φ 1 ( h , b ) u u t ) | ( u , b ) t + o ( t ) = φ 1 ( u , b ) + ( r u + c 1 ) t φ 1 ( u , b ) u + o ( t ) . (11)

Insert (11) into (10) and divide it by t. Let t 0 , we get the desired result. Since (7) and (8) can be proved by using the same method, the detail is omitted. □

Next, we prove the main result on the absolute ruin probability. For convenience, we denote

M ˜ = [ 0 0 0 0 L 5 ( c 2 ρ ) L 6 ( c 2 ρ ) L 1 ( b ) L 2 ( b ) L 3 ( b ) L 4 ( b ) 0 0 0 0 L 3 ( 0 ) L 4 ( 0 ) L 5 ( 0 ) L 6 ( 0 ) q 1 ( b ) q 2 ( b ) p 3 ( b ) p 4 ( b ) 0 0 0 0 L 3 ( 0 ) L 4 ( 0 ) L 5 ( 0 ) L 6 ( 0 ) 0 0 0 0 L 5 ( c 2 ρ ) L 6 ( c 2 ρ ) ] ,

p i ( u ) = ( r u + c 2 ) L i ( u ) , ( i = 3 , 4 ) , q i ( u ) = ( r u + c 1 ) L i ( u ) , ( i = 1 , 2 ) , C = ( C 3 , C 4 , C 5 , C 6 , C 7 , C 8 ) T , B = ( 1 , 0 , 0 , 0 , 0 , 0 ) T . The hyper-geometric function used later is introduced as follows

M ( a , b , x ) = Γ ( b ) Γ ( b a ) Γ ( a ) 0 1 e x t t a 1 ( 1 t ) b a 1 d t , b > a > 0 , (12)

where the Gamma function is Γ ( x ) = 0 e t t ( x 1 ) d t .

In view of (6), let β = 0 and W ( x , y ) = 1 in the Gerber-Shiu discounted penalty function φ i ( u , b ) ( i = 1 , 2 , 3 ) , we can get the following theorem on the ruin probability.

Theorem 2.2. In model (4), the closed-form ruin probability φ i ( u ) ( i = 1 , 2 , 3 ) are given like

φ 1 ( u ) = C 3 e η 1 u M ( r λ p r , 1 λ r , ( r u + c 1 ) η 1 r ) + C 4 ( ( r u + c 1 ) η 1 r ) λ r e η 1 u M ( r + ( 1 p ) λ r , 1 + λ r , ( r u + c 1 ) η 1 r ) ; (13)

φ 2 ( u ) = C 5 e η 2 u M ( r λ q r , 1 λ r , ( r u + c 2 ) η 2 r ) + C 6 ( ( r u + c 2 ) η 2 r ) λ r e η 2 u M ( r + ( 1 q ) λ r , 1 + λ r , ( r u + c 2 ) η 2 r ) ; (14)

φ 3 ( u ) = C 7 e η 2 u M ( ρ λ q ρ , 1 λ ρ , ( ρ u + c 2 ) η 2 ρ ) + C 8 ( ( ρ u + c 2 ) η 2 ρ ) λ ρ e η 2 u M ( ρ + ( 1 q ) λ ρ , 1 + λ ρ , ( ρ u + c 2 ) η 2 ρ ) ; (15)

where C = ( C 3 , C 4 , , C 8 ) T is solution of M ˜ C = B .

Proof. By Theorem 2.1, we get

( r u + c 1 ) φ 1 ( u , b ) + ( r λ + ( r u + c 1 ) η 1 ) φ 1 ( u , b ) ( 1 p ) λ η 1 φ 1 ( u , b ) = 0 , u b ; (16)

( r u + c 2 ) φ 2 ( u , b ) + ( r λ + ( r u + c 2 ) η 2 ) φ 2 ( u , b ) ( 1 q ) λ η 2 φ 2 ( u , b ) = 0 , 0 u < b ; (17)

( ρ u + c 2 ) φ 3 ( u , b ) + ( ρ λ + ( ρ u + c 2 ) η 2 ) φ 3 ( u , b ) ( 1 q ) λ η 2 φ 3 ( u , b ) = 0 , c 2 ρ < u < 0. (18)

The boundary conditions are

φ 3 ( c 2 ρ , b ) = 1 , lim u φ 1 ( u , b ) = 0 , φ 1 ( b , b ) = φ 2 ( b , b ) , φ 2 ( 0 , b ) = φ 3 ( 0 , b ) , ( r b + c 2 ) φ 2 ( b , b ) = ( r b + c 1 ) φ 1 ( b , b ) , φ 3 ( 0 , b ) = φ 2 ( 0 , b ) , φ 3 ( c 2 ρ , b ) = 0.

In order to obtain the expression of ruin probabilities, we employ the confluent hyper-geometric function. In (16), make auxiliary function φ 1 ( u ) = X 1 ( u ) e a u where a is a constant to be determined later and then

( r u + c 1 ) X 1 ( u ) + [ ( 2 a + η 1 ) ( r u + c 1 ) + r λ ] X 1 ( u ) + [ ( 1 p ) λ η 1 + a 2 ( r u + c 1 ) + a ( r λ + ( r u + c 1 ) η 1 ) ] X 1 ( u ) = 0 , (19)

Define r u + c 1 = z , we have

z X 1 ( z c 1 r ) + ( 1 λ r + 2 a + η 1 r z ) X 1 ( z c 1 r ) + 1 r 2 [ ( 1 p ) λ η 1 + a ( r λ ) + ( a 2 + a η 1 ) z ] X 1 ( z c 1 r ) = 0. (20)

By choosing a = η 1 ,

z X 1 ( z c 1 r ) + ( 1 λ r η 1 r z ) X 1 ( z c 1 r ) + ( λ p r r η 1 r ) X 1 ( z c 1 r ) = 0 , (21)

Define η 1 r z = v , we have by (21)

v X 1 ( r v η 1 c 1 r η 1 ) + ( 1 λ r v ) X 1 ( r v η 1 c 1 r η 1 ) + ( λ p r r ) X 1 ( r v η 1 c 1 r η 1 ) = 0 , (22)

Let X 1 ( r v η 1 c 1 r η 1 ) = Y 1 ( v ) , (22) can be rewritten as

v Y 1 ( v ) + ( 1 λ r v ) Y 1 ( v ) ( r λ p r ) Y 1 ( v ) = 0, (23)

By using the first kind of confluent hypergeometric function, we obtain that

Y 1 ( v ) = C 3 M ( r λ p r , 1 λ r ; v ) + C 4 ( v ) λ r M ( r + ( 1 p ) λ r , 1 + λ r ; v ) , (24)

where C 3 and C 4 are constants. Without extra claims, C i are constants for i 3 in the following.

Next, the same method can be used to calculate (17) and (18). In (17), set

φ 2 ( u ) = X 2 ( u ) e a u , r u + c 2 = z , a = η 2 , η 2 r z = v and X 2 ( r v η 2 c 2 r η 2 ) = Y 2 ( v ) , then we have

v Y 2 ( v ) + ( 1 λ r v ) Y 2 ( v ) ( r λ q r ) Y 2 ( v ) = 0 , (25)

The solution of (25) has the form

Y 2 ( v ) = C 5 M ( r λ q r , 1 λ r ; v ) + C 6 ( v ) λ r M ( r + ( 1 q ) λ r , 1 + λ r ; v ) . (26)

Similarly, the solution of (18) is

Y 3 ( v ) = C 7 M ( ρ λ q ρ , 1 λ ρ ; v ) + C 8 ( v ) λ ρ M ( ρ + ( 1 q ) λ ρ , 1 + λ ρ ; v ) , (27)

In the sequel, we turn to compute the expression of ruin probability. Inserting

v = η 1 r z and X 1 ( r v η 1 c 1 r η 1 ) = Y 1 ( v ) into (24) yields that

X 1 ( z c 1 r ) = C 3 M ( r λ p r , 1 λ r ; η 1 z r ) + C 4 ( η 1 z r ) λ r M ( r + ( 1 p ) λ r , 1 + λ r ; η 1 z r ) . (28)

Denote r u + c 1 = z , then (28) yields that

X 1 ( u ) = C 3 M ( r λ p r ,1 λ r ; ( r u + c 1 ) η 1 r ) + C 4 ( ( r u + c 1 ) η 1 r ) λ r M ( r + ( 1 p ) λ r ,1 + λ r ; ( r u + c 1 ) η 1 r ) . (29)

Define φ 1 ( u ) = X 1 ( u ) e η 1 u , we have by (29)

φ 1 ( u ) = C 3 e η 1 u M ( r λ p r , 1 λ r ; ( r u + c 1 ) η 1 r ) + C 4 ( ( r u + c 1 ) η 1 r ) λ r e η 1 u M ( r + ( 1 p ) λ r , 1 + λ r ; ( r u + c 1 ) η 1 r ) . (30)

Then we can use the same method to solve (26) and (27) like that

φ 2 ( u ) = C 5 e η 2 u M ( r λ q r , 1 λ r , ( r u + c 2 ) η 2 r ) + C 6 ( ( r u + c 2 ) η 2 r ) λ r e η 2 u M ( r + ( 1 q ) λ r , 1 + λ r , ( r u + c 2 ) η 2 r ) ; (31)

φ 3 ( u ) = C 7 e η 2 u M ( ρ λ q ρ , 1 λ ρ , ( ρ u + c 2 ) η 2 ρ ) + C 8 ( ( ρ u + c 2 ) η 2 ρ ) λ ρ e η 2 u M ( ρ + ( 1 q ) λ ρ , 1 + λ ρ , ( ρ u + c 2 ) η 2 ρ ) . (32)

To calculate the coefficients of the ruin probability C = ( C 3 , C 4 , , C 8 ) T , we firstly rewrite (30)-(32) to be

φ 1 ( u , b ) = C 3 L 1 ( u ) + C 4 L 2 ( u ) ; φ 2 ( u , b ) = C 5 L 3 ( u ) + C 6 L 4 ( u ) and φ 3 ( u , b ) = C 7 L 5 ( u ) + C 8 L 6 ( u ) . (33)

By inserting boundary conditions into (30), (31) and (32), we get

( C 7 L 5 ( c 2 ρ ) + C 8 L 6 ( c 2 ρ ) = 1 ; C 3 L 1 ( b ) + C 4 L 2 ( b ) C 5 L 3 ( b ) C 6 L 4 ( b ) = 0 ; C 5 L 3 ( 0 ) + C 6 L 4 ( 0 ) C 7 L 5 ( 0 ) C 8 L 6 ( 0 ) = 0 ; ( r b + c 2 ) C 5 L 3 ( b ) + ( r b + c 2 ) C 6 L 4 ( b ) ( r b + c 1 ) C 3 L 1 ( b ) ( r b + c 1 ) C 4 L 2 ( b ) = 0 ; C 7 L 5 ( 0 ) + C 8 L 6 ( 0 ) C 5 L 3 ( 0 ) C 6 L 4 ( 0 ) = 0 ; C 7 L 5 ( c 2 ρ ) + C 8 L 6 ( c 2 ρ ) = 0 ,

(34)

which can be simplified as M ˜ C = B . The proof is complete. □

Remark 1. Let p = q = 1 and η 1 = η 2 , Theorem 2.1 reduces to the related results given in [2] . In this sense, the previously-known results are generalized. Further, Theorem 2.2 shows us the closed-form solution of the ruin probability which is valuable in theoretical applications.

3. Conclusion

In this letter, we discuss the ruin probability of the risk model under assumption that the claim amounts obey the switched exponential distribution. The closed-form expressions of the ruin probability are obtained by using the first kind of hypergeometric function. The obtained results may give us guidance in facing up to the economic crisis.

Cite this paper

Liu, Y.X. and Zhao, D.L. (2017) Closed-Form Absolute Ruin Problems of the Risk Models with State-Dependent Switched Claims. Journal of Applied Mathematics and Physics, 5, 2326-2334. https://doi.org/10.4236/jamp.2017.512190

References

  1. 1. Gerber, H.U. (1981) On the Probability of Ruin in the Presence of a Linear Dividend Barrier. Scandinavian Actuarial Journal, 1981, 105-115. https://doi.org/10.1080/03461238.1981.10413735

  2. 2. Zou, W., Gao, J. and Xie, J. (2014) On the Expected Discounted Penalty Function and Optimal Dividend Strategy for a Risk Model with Random Incomes and Interclaim-Dependent Claim Sizes. Journal of Computational and Applied Mathematics, 255, 270-281. https://doi.org/10.1016/j.cam.2013.05.004

  3. 3. Yao, D., Yang, H. and Wang, R. (2011) Optimal Dividend and Capital Injection Problem in the Dual Model with Proportional and Fixed Transaction Costs. European Journal of Operational Research, 211, 568-576. https://doi.org/10.1016/j.ejor.2011.01.015

  4. 4. Yuen, K., Zhou, M. and Guo, J. (2008) On a Risk Model with Debit Interest and Dividend Payments. Statistics and Probability Letters, 78, 2426-2432. https://doi.org/10.1016/j.spl.2008.02.021

  5. 5. Yin, C. and Yuen, K. (2011) Optimality of the Threshold Dividend Strategy for the Compound Poisson Model. Statistics and Probability Letters, 81, 1841-1846. https://doi.org/10.1016/j.spl.2011.07.022

  6. 6. Lin, X. and Pavlova, K.P. (2006) The Compound Poisson Risk Model with a Threshold Dividend Strategy. Insurance: Mathematics and Economics, 38, 57-80. https://doi.org/10.1016/j.insmatheco.2005.08.001

  7. 7. Lin, X. and Sendova, K.P. (2008) The Compound Poisson Risk Model with Multiple Thresholds. Insurance: Mathematics and Economics, 42, 617-627. https://doi.org/10.1016/j.insmatheco.2007.06.008

  8. 8. Fu, K.-A. (2016) On Joint Ruin Probability for a Bidimensional Levy-Driven Risk Model with Stochastic Returns and Heavy-Tailed Claims. Journal of Mathematical Analysis and Applications, 442, 17-30. https://doi.org/10.1016/j.jmaa.2016.04.042

  9. 9. Zhao, Y. and Yao, J. (2005) Actuarial Module of Enterprise-Complemented Pension Scheme While Taking Vasicek as Stochastic Interest Rate Module. Journal of University of Shanghai for Science and Technology, 27, 268-270. (In Chinese)

  10. 10. Wang, Q., Xu, Z. and Gao, Y. (2009) Research on the Pricing Problem of Whole-Life-Equity-Linked Policies. Journal of University of Shanghai for Science and Technology, 31, 80-84. (In Chinese)

  11. 11. Chadjiconstantinidis, S. and Papaioannou, A.D. (2013) On a Perturbed by Diffusion Compound Poisson Risk Model with Delayed Claims and Multi-Layer Dividend Strategy. Journal of Computational and Applied Mathematics, 253, 26-50. https://doi.org/10.1016/j.cam.2013.02.014

  12. 12. Lu, D. and Zhang, B. (2016) Some Asymptotic Results of the Ruin Probabilities in a Two-Dimensional Renewal Risk Model with Some Strongly Subexponential Claims. Statistics and Probability Letters, 114, 20-29. https://doi.org/10.1016/j.spl.2016.03.005