Valuation of GLWB-LTC Annuities with L\'evy Equity Dynamics, Stochastic Interest Rates and Health-State Transitions
Pith reviewed 2026-06-28 23:51 UTC · model grok-4.3
The pith
A hybrid tree-IMEX scheme prices GLWB-LTC annuities under Lévy equity jumps and Hull-White rates, showing these risks materially alter fair fees and surrender incentives.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The hybrid tree-IMEX method delivers stable long-maturity prices consistent with simulation benchmarks. Lévy equity dynamics and stochastic interest rates have a material impact on fair fees and surrender incentives, and affect the decomposition of contract value. The findings highlight the importance of modelling financial tail risk and interest-rate risk jointly when pricing long-term insurance guarantees with LTC-contingent benefits.
What carries the argument
The hybrid tree-IMEX scheme: a recombining Hull-White trinomial tree for stochastic rates coupled to an implicit-explicit finite-difference discretization of the Lévy-driven fund value on a seven-state health model.
If this is right
- Fair fees computed under Lévy dynamics differ materially from those under geometric Brownian motion.
- Adding stochastic interest rates changes the optimal surrender strategy and the value split between guarantee and optionality.
- The seven-state health model can be embedded directly into the same tree-FD grid without losing recombining structure.
- Long-maturity contract values remain computable in seconds rather than hours once the hybrid scheme is calibrated.
- Joint modeling of tail risk and rate risk is required to avoid systematic mispricing of LTC-contingent withdrawal guarantees.
Where Pith is reading between the lines
- The same hybrid scheme could be reused for other health-contingent guarantees such as critical-illness or disability income riders.
- Calibration of Lévy parameters to equity option surfaces would be needed before the method could be used for live pricing.
- If the health-transition matrix is estimated from limited claims data, small changes in intensities could shift the fair fee by several basis points.
- The method's stability at 30-plus years suggests it may also handle variable-annuity guarantees with similar maturity and optionality.
Load-bearing premise
The seven-state health model and its transition intensities accurately represent real disability dynamics, and the numerical scheme remains stable and accurate for long contract maturities without material discretization bias.
What would settle it
Monte Carlo prices for a 30-year GLWB-LTC contract computed under the same Lévy and Hull-White parameters deviate by more than a few percent from the hybrid-tree prices at the reported fair-fee level.
Figures
read the original abstract
This paper develops a valuation framework for guaranteed lifetime withdrawal benefit (GLWB) contracts with long-term care (LTC) features when the reference fund follows exponential Levy dynamics and the short rate follows the Hull-White model. The contract combines financial guarantees, longevity protection, health-contingent LTC payments, and surrender optionality, requiring the joint treatment of jump risk, stochastic discounting, and disability risk. The numerical method couples a recombining Hull-White trinomial tree with an implicit-explicit (IMEX) finite difference scheme. The framework incorporates a seven-state health model, annual fees, LTC payments, guaranteed withdrawals, and bang-bang policyholder actions, and is benchmarked against Monte Carlo simulation. Numerical results show that the hybrid tree-IMEX method delivers stable long-maturity prices consistent with simulation benchmarks. They also show that Levy equity dynamics and stochastic interest rates have a material impact on fair fees and surrender incentives, and affect the decomposition of contract value. The findings highlight the importance of modelling financial tail risk and interest-rate risk jointly when pricing long-term insurance guarantees with LTC-contingent benefits.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a hybrid numerical valuation framework for GLWB-LTC annuities under exponential Lévy equity dynamics and Hull-White stochastic interest rates, coupled to a seven-state Markov health model with LTC payments, guaranteed withdrawals, annual fees, and bang-bang surrender options. The method combines a recombining Hull-White trinomial tree for the short rate with an IMEX finite-difference scheme for the remaining state variables and is benchmarked against Monte Carlo simulation. Numerical experiments are used to demonstrate stability over long maturities and to quantify the material effects of Lévy jumps and stochastic rates on fair fees, surrender incentives, and value decomposition.
Significance. If the reported numerical stability and Monte Carlo agreement hold under scrutiny, the work provides a practical route to pricing high-dimensional insurance guarantees that jointly incorporate jump risk, interest-rate risk, and health-state transitions. The explicit benchmarking against simulation is a positive feature that strengthens the central claim of reliable long-maturity valuation.
major comments (2)
- [Numerical implementation and results sections] The central claim of stable long-maturity prices consistent with Monte Carlo benchmarks rests on the hybrid tree-IMEX scheme, yet the manuscript provides no quantitative convergence study (grid sizes, time-step sizes, or error tables) for the non-local Lévy operator over 30-year horizons; without such evidence it is impossible to rule out material discretization bias in the reported fair fees and surrender values.
- [Model formulation and numerical scheme] The coupling of the seven-state health transitions to the jump-diffusion and short-rate processes inside the IMEX scheme is described at a high level; the paper does not specify how the non-local integral term arising from the Lévy measure is discretized or whether any truncation or approximation is introduced, which directly affects the accuracy of the LTC-contingent cash-flow valuation.
minor comments (2)
- Notation for the health-state intensities and the decomposition of contract value into guarantee, LTC, and surrender components could be made more uniform across tables and figures.
- A brief statement on the parameter calibration procedure for the Lévy process and Hull-White model would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on the numerical validation and scheme details. We address each major comment below and will incorporate the requested clarifications and additional results in the revised manuscript.
read point-by-point responses
-
Referee: [Numerical implementation and results sections] The central claim of stable long-maturity prices consistent with Monte Carlo benchmarks rests on the hybrid tree-IMEX scheme, yet the manuscript provides no quantitative convergence study (grid sizes, time-step sizes, or error tables) for the non-local Lévy operator over 30-year horizons; without such evidence it is impossible to rule out material discretization bias in the reported fair fees and surrender values.
Authors: The referee correctly notes that the manuscript does not contain a dedicated quantitative convergence study with explicit grid sizes, time-step sizes, and error tables for the non-local Lévy operator over 30-year horizons. While Monte Carlo benchmarking is used to demonstrate consistency of the reported prices, this does not substitute for a systematic discretization-error analysis. In the revision we will add a new subsection (or appendix) presenting convergence tables for spatial and temporal refinements of the IMEX scheme applied to the Lévy integral term, covering representative 30-year contracts and reporting both absolute and relative errors against a fine-grid reference. revision: yes
-
Referee: [Model formulation and numerical scheme] The coupling of the seven-state health transitions to the jump-diffusion and short-rate processes inside the IMEX scheme is described at a high level; the paper does not specify how the non-local integral term arising from the Lévy measure is discretized or whether any truncation or approximation is introduced, which directly affects the accuracy of the LTC-contingent cash-flow valuation.
Authors: We agree that the current description of the IMEX discretization of the non-local Lévy integral and its coupling to the seven-state health Markov chain is at a high level. The revised manuscript will include an expanded subsection that (i) states the quadrature rule or truncation used for the Lévy integral (including the cutoff radius and any singularity handling), (ii) details how the health-state transition intensities enter the implicit/explicit splitting, and (iii) specifies the boundary conditions and interpolation between the Hull-White tree nodes and the finite-difference grid for the remaining variables. revision: yes
Circularity Check
Numerical valuation framework benchmarked against independent Monte Carlo simulation exhibits no circularity
full rationale
The paper constructs a hybrid numerical scheme (recombining Hull-White trinomial tree coupled to IMEX finite differences) to value GLWB-LTC contracts under Lévy equity dynamics and stochastic rates, then directly compares the resulting prices to separate Monte Carlo benchmarks. All load-bearing claims (stability for long maturities, material impact of jumps and rates on fees/surrender) are presented as outputs verified against this external simulation rather than being defined by or fitted to the same quantities inside the method. No self-definitional steps, fitted-input predictions, or load-bearing self-citations appear in the derivation chain; the seven-state health model and policyholder actions are exogenous inputs whose effects are measured, not presupposed.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Apicella, G., Gaudenzi, M., and Molent, A. (2025). The life care annuity: Enhancing product features and refining pricing methods.Decisions in Economics and Finance, 48(2):873–911
2025
-
[2]
R., Maggistro, R., and Zoccolan, I
Bacinello, A. R., Maggistro, R., and Zoccolan, I. (2024). Risk-neutral valuation of GLWB riders in variable annuities.Insurance: Mathematics and Economics, 114:1–14
2024
-
[3]
R., Millossovich, P., and Montealegre, A
Bacinello, A. R., Millossovich, P., and Montealegre, A. (2016). The valuation of GMWB variable annuit- ies under alternative fund distributions and policyholder behaviours.Scandinavian Actuarial Journal, 2016(5):446–465
2016
-
[4]
R., Millossovich, P., Olivieri, A., and Pitacco, E
Bacinello, A. R., Millossovich, P., Olivieri, A., and Pitacco, E. (2011). Variable annuities: A unifying valuation approach.Insurance: Mathematics and Economics, 49(3):285–297
2011
-
[5]
Briani, M., Caramellino, L., Terenzi, G., and Zanette, A. (2019). Numerical stability of a hybrid method for pricing options.International Journal of Theoretical and Applied Finance, 22(7):1–46
2019
-
[6]
and Mercurio, F
Brigo, D. and Mercurio, F. (2006).Interest Rate Models: Theory and Practice. Springer, 2 edition
2006
-
[7]
and Warshawsky, M
Brown, J. and Warshawsky, M. (2013). The life care annuity: A new empirical examination of an insurance innovation that addresses problems in the markets for life annuities and long-term care insurance.Journal of Risk and Insurance, 80(3):677–703. 31
2013
-
[8]
Chen, S., Cui, Z., Zhang, Z., and Zhong, W. (2026). Valuation of GLWB annuities with optional conversion to combo products providing LTC benefits.Scandinavian Actuarial Journal, 2026(3):263–294
2026
-
[9]
and Tankov, P
Cont, R. and Tankov, P. (2004).Financial Modelling with Jump Processes. Chapman & Hall/CRC, Boca Raton
2004
-
[10]
and Voltchkova, E
Cont, R. and Voltchkova, E. (2005). A finite difference scheme for option pricing in jump diffusion and exponential Lévy models.SIAM Journal on Numerical Analysis, 43(4):1596–1626
2005
-
[11]
Forsyth, P. A. and Vetzal, K. R. (2014). An optimal stochastic control framework for determining the cost of hedging of variable annuities.Journal of Economic Dynamics and Control, 44:29–53. Goudenège, L., Molent, A., and Zanette, A. (2021). Gaussian process regression for pricing variable annuities with stochastic volatility and interest rate.Decisions i...
2014
-
[12]
L., Chiu, Y.-F., and Chen, Y.-C
Hsieh, M.-h., Wang, J. L., Chiu, Y.-F., and Chen, Y.-C. (2018). Valuation of variable long-term care annuities with guaranteed lifetime withdrawal benefits: A variance reduction approach.Insurance: Mathematics and Economics, 78:246–254
2018
-
[13]
Molent, A. (2020). Taxation of a GMWB variable annuity in a stochastic interest rate model.ASTIN Bulletin: The Journal of the IAA, 50(3):1001–1035
2020
-
[14]
M., Spillman, B
Murtaugh, C. M., Spillman, B. C., and Warshawsky, M. J. (2001). In sickness and in health: An annuity approach to financing long-term care and retirement income.Journal of Risk and Insurance, 68(2):225–254
2001
-
[15]
Pritchard, D. J. (2006). Modeling disability in long-term care insurance.North American Actuarial Journal, 10(4):48–75. A Log-return moment calculation with Hull–White interest rates This appendix summarises the moment calculation used in Table 4. Throughout the appendix,T >0denotes a generic time horizon over which the log-return moments are computed. Le...
2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.