Viscosity Solutions

Goal-based portfolio selection with mental accounting ArXiv ID: 2506.06654 “View on arXiv” Authors: Erhan Bayraktar, Bingyan Han Abstract We present a continuous-time portfolio selection framework that reflects goal-based investment principles and mental accounting behavior. In this framework, an investor with multiple investment goals constructs separate portfolios, each corresponding to a specific goal, with penalties imposed on fund transfers between these goals, referred to as mental costs. By applying the stochastic Perron’s method, we demonstrate that the value function is the unique constrained viscosity solution of a Hamilton-Jacobi-Bellman equation system. Numerical analysis reveals several key features: the free boundaries exhibit complex shapes with bulges and notches; the optimal strategy for one portfolio depends on the wealth level of another; investors must diversify both among stocks and across portfolios; and they may postpone reallocating surplus from an important goal to a less important one until the former’s deadline approaches. ...

Approximation and regularity results for the Heston model and related processes ArXiv ID: 2504.21658 “View on arXiv” Authors: Edoardo Lombardo Abstract This Ph.D. thesis explores approximations and regularity for the Heston stochastic volatility model through three interconnected works. The first work focuses on developing high-order weak approximations for the Cox-Ingersoll-Ross (CIR) process, essential for financial modelling but challenging due to the square root diffusion term preventing standard methods. By employing the random grid technique (Alfonsi & Bally, 2021) built upon Alfonsi’s (2010) second-order scheme, the work proves that weak approximations of any order can be achieved for smooth test functions. This holds under a condition that is less restrictive than the famous Feller’s one. Numerical results confirm convergence for both CIR and Heston models and show significant computational time improvements. The second work extends the random grid technique to the log-Heston process. Two second-order schemes are introduced (one using exact volatility simulation, another using Ninomiya-Victoir splitting under a the same restriction used above). Convergence to any desired order is rigorously proven. Numerical experiments validate the schemes’ effectiveness for pricing European and Asian options and suggest potential applicability to multifactor/rough Heston models. The third work investigates the partial differential equation (PDE) associated with the log-Heston model. It extends classical solution results and establishes the existence and uniqueness of viscosity solutions without relying on the Feller condition. Uniqueness is proven even for certain discontinuous initial data, relevant for pricing instruments like digital options. Furthermore, the convergence of a hybrid numerical scheme to the viscosity solution is shown under relaxed regularity (continuity) for the initial data. An appendix includes supplementary results for the CIR process. ...

Numerical analysis of American option pricing in a two-asset jump-diffusion model ArXiv ID: 2410.04745 “View on arXiv” Authors: Unknown Abstract This paper addresses an important gap in rigorous numerical treatments for pricing American options under correlated two-asset jump-diffusion models using the viscosity solution framework, with a particular focus on the Merton model. The pricing of these options is governed by complex two-dimensional (2-D) variational inequalities that incorporate cross-derivative terms and nonlocal integro-differential terms due to the presence of jumps. Existing numerical methods, primarily based on finite differences, often struggle with preserving monotonicity in the approximation of cross-derivatives, a key requirement for ensuring convergence to the viscosity solution. In addition, these methods face challenges in accurately discretizing 2-D jump integrals. We introduce a novel approach to effectively tackle the aforementioned variational inequalities while seamlessly handling cross-derivative terms and nonlocal integro-differential terms through an efficient and straightforward-to-implement monotone integration scheme. Within each timestep, our approach explicitly enforces the inequality constraint, resulting in a 2-D Partial Integro-Differential Equation (PIDE) to solve. Its solution is expressed as a 2-D convolution integral involving the Green’s function of the PIDE. We derive an infinite series representation of this Green’s function, where each term is non-negative and computable. This facilitates the numerical approximation of the PIDE solution through a monotone integration method. To enhance efficiency, we develop an implementation of this monotone scheme via FFTs, exploiting the Toeplitz matrix structure. The proposed method is proved to be both $\ell_{"\infty"} $-stable and consistent in the viscosity sense, ensuring its convergence to the viscosity solution of the variational inequality. Extensive numerical results validate the effectiveness and robustness of our approach. ...