It is well-known that decision-making problems from stochastic control can be formulated by means of a forward–backward stochastic differential equation (FBSDE). Recently, the authors of Ji et al. (2022) proposed an efficient deep learning algorithm based on the stochastic maximum principle (SMP). In this paper, we provide a convergence result for this deep SMP-BSDE algorithm and compare its performance with other existing methods. In particular, by adopting a strategy as in Han and Long (2020), we derive a-posteriori estimate, and show that the total approximation error can be bounded by the value of the loss functional and the discretization error. We present numerical examples for high-dimensional stochastic control problems, both in the cases of drift- and diffusion control, which showcase superior performance compared to existing algorithms.

The pressure-free two-fluid model (PFTFM) is a recent reformulation of the one-dimensional two-fluid model (TFM) for stratified incompressible flow in ducts (including pipes and channels), in which the pressure is eliminated through intricate use of the volume constraint. The disadvantage of the PFTFM was that the volumetric flow rate had to be specified as an input, even though it is an unknown quantity in case of periodic boundary conditions. In this work, we derive an expression for the volumetric flow rate that is based on the demand for energy (and momentum) conservation. This leads to PFTFM solutions that match those of the TFM, justifying the validity and necessity of the derived choice of volumetric flow rate. Furthermore, we extend an energy-conserving spatial discretization of the TFM, in the form of a finite volume scheme, to the PFTFM. We propose a discretization of the volumetric flow rate that yields discrete momentum and energy conservation. The discretization is extended with an energy-conserving discretization of the source terms related to gravity acting in the streamwise direction. Our numerical experiments confirm that the discrete energy is conserved for different problem settings, including sloshing in an inclined closed tank, and a traveling wave in a periodic domain. The PFTFM solutions and the volumetric flow rates match the TFM solutions, with reduced computation time, and with exact momentum and energy conservation.

Generative adversarial networks (GANs) have shown promising results when applied on partial differential equations and financial time series generation. We investigate if GANs can also be used to approximate one-dimensional Ito ^ stochastic differential equations (SDEs). We propose a scheme that approximates the path-wise conditional distribution of SDEs for large time steps. Standard GANs are only able to approximate processes in distribution, yielding a weak approximation to the SDE. A conditional GAN architecture is proposed that enables strong approximation. We inform the discriminator of this GAN with the map between the prior input to the generator and the corresponding output samples, i.e. we introduce a ‘supervised GAN’. We compare the input-output map obtained with the standard GAN and supervised GAN and show experimentally that the standard GAN may fail to provide a path-wise approximation. The GAN is trained on a dataset obtained with exact simulation. The architecture was tested on geometric Brownian motion (GBM) and the Cox–Ingersoll–Ross (CIR) process. The supervised GAN outperformed the Euler and Milstein schemes in strong error on a discretisation with large time steps. It also outperformed the standard conditional GAN when approximating the conditional distribution. We also demonstrate how standard GANs may give rise to non-parsimonious input-output maps that are sensitive to perturbations, which motivates the need for constraints and regularisation on GAN generators.

In this work, we consider rule-based investment strategies for managing a defined contribution pension savings scheme, under the Dutch pension fund testing model. We find that dynamic, rule-based investment strategies can outperform traditional static strategies, by which we mean that the investor may achieve the target retirement income with a higher probability or limit the shortfall when the target is not met. In comparison with dynamic programming-based strategies, the rule-based strategies have more stable asset allocations throughout time and avoid excessive transactions that may be hard to explain to an investor. We also study a combined strategy of a rule-based target with dynamic programming. A key feature of our setting is that there is no risk-free asset, instead, a matching portfolio is introduced for the investor to avoid unnecessary risk.

In this paper, we propose third-order semi-discretized schemes in space based on the tempered weighted and shifted Grunwald difference (tempered-WSGD) operators for the tempered fractional diffusion equation. We also show stability and convergence analysis for the fully discrete scheme based a Crank–Nicolson scheme in time. A third-order scheme for the tempered Black–Scholes equation is also proposed and tested numerically. Some numerical experiments are carried out to confirm accuracy and effectiveness of these proposed methods.

Since the 2007/2008 financial crisis, the total value adjustment (XVA) should be included when pricing financial derivatives. In the present paper, the derivative values of European and American options have been priced where we take into account counterparty risk. Whereas European and American options considering counterparty risk have already been priced under Black-Scholes dynamics in [1], here the novel contribution is the introduction of stochastic volatility resulting in a Heston stochastic volatility type partial differential equation to be solved. We derive the partial differential equation modeling the XVA when stochastic volatility is assumed. For both European and American options, a linear and a nonlinear problem have been deduced. In order to obtain a numerical solution, suitable and appropriate boundary conditions have been considered. In addition, a method of characteristics for the time discretization combined with a finite element method in the spatial discretization has been implemented. The expected exposure and potential future exposure are also computed to compare the current model with the associated Black–Scholes model.

Artificial neural networks (ANNs) have recently also been applied to solve partial differential equations (PDEs). The classical problem of pricing European and American financial options, based on the corresponding PDE formulations, is studied here. Instead of using numerical techniques based on finite element or difference methods, we address the problem using ANNs in the context of unsupervised learning. As a result, the ANN learns the option values for all possible underlying stock values at future time points, based on the minimization of a suitable loss function. For the European option, we solve the linear Black–Scholes equation, whereas for the American option we solve the linear complementarity problem formulation. Two-asset exotic option values are also computed, since ANNs enable the accurate valuation of high-dimensional options. The resulting errors of the ANN approach are assessed by comparing to the analytic option values or to numerical reference solutions (for American options, computed by finite elements). In the short note, previously published, a brief introduction to this work was given, where some ideas to price vanilla options by ANNs were presented, and only European options were addressed. In the current work, the methodology is introduced in much more detail.

This study contributes to understanding Valuation Adjustments (xVA) by focussing on the dynamic hedging of Credit Valuation Adjustment (CVA), corresponding Profit & Loss (P&L) and the P&L explain. This is done in a Monte Carlo simulation setting, based on a theoretical hedging framework discussed in existing literature. We look at hedging CVA market risk for a portfolio with European options on a stock, first in a Black-Scholes setting, then in a Merton jump-diffusion setting. Furthermore, we analyze the trading business at a bank after including xVAs in pricing. We provide insights into the hedging of derivatives and their xVAs by analyzing and visualizing the cash-flows of a portfolio from a desk structure perspective. The case study shows that not charging CVA at trade inception results in an expected loss. Furthermore, hedging CVA market risk is crucial to end up with a stable trading strategy. In the Black-Scholes setting this can be done using the underlying stock, whereas in the Merton jump-diffusion setting we need to add extra options to the hedge portfolio to properly hedge the jump risk. In addition to the simulation, we derive analytical results that explain our observations from the numerical experiments. Understanding the hedging of CVA helps to deal with xVAs in a practical setting.

Since the 2007/2008 financial crisis, the total value adjustment (XVA) should be included when pricing financial derivatives. In the present paper, the derivative values of European and American options have been priced where we take into account counterparty risk. Whereas European and American options considering counterparty risk have already been priced under Black-Scholes dynamics in [2], here the novel contribution is the introduction of stochastic volatility resulting in a Heston stochastic volatility type partial differential equation to be solved. We derive the partial differential equation modeling the XVA when stochastic volatility is assumed. For both European and American options, a linear and a nonlinear problem have been deduced. In order to obtain a numerical solution, suitable and appropriate boundary conditions have been considered. In addition, a method of characteristics for the time discretization combined with a finite element method in the spatial discretization has been implemented. The expected exposure and potential future exposure are also computed to compare the current model with the associated Black–Scholes model.

In this paper, we propose a neural network-based method for CVA computations of a portfolio of derivatives. In particular, we focus on portfolios consisting of a combination of derivatives, with and without true optionality, e.g., a portfolio of a mix of European- and Bermudan-type derivatives. CVA is computed, with and without netting, for different levels of WWR and for different levels of credit quality of the counterparty. We show that the CVA is overestimated with up to 25% by using the standard procedure of not adjusting the exercise strategy for the default-risk of the counterparty. For the Expected Shortfall of the CVA dynamics, the overestimation was found to be more than 100% in some non-extreme cases.

In this work, we developed a Python demonstrator for pricing total valuation adjustment (XVA) based on the stochastic grid bundling method (SGBM). XVA is an advanced risk management concept which became relevant after the recent financial crisis. This work is a follow-up work on Chau and Oosterlee in (Int J Comput Math 96(11):2272–2301, 2019), in which we extended SGBM to numerically solving backward stochastic differential equations (BSDEs). The motivation for this work is basically two-fold. On the application side, by focusing on a particular financial application of BSDEs, we can show the potential of using SGBM on a real-world risk management problem. On the implementation side, we explore the potential of developing a simple yet highly efficient code with SGBM by incorporating CUDA Python into our program.

A novel generating mechanism for non-strict bivariate Archimedean copulas via the Lorenz curve of a non-negative random variable is proposed. Lorenz curves have been extensively studied in economics and statistics to characterize wealth inequality and tail risk. In this paper, these curves are seen as integral transforms generating increasing convex functions in the unit square. Many of the properties of these "Lorenz copulas", from tail dependence and stochastic ordering, to their Kendall distribution function and the size of the singular part, depend on simple features of the random variable associated to the generating Lorenz curve. For instance, by selecting random variables with a lower bound at zero it is possible to create copulas with asymptotic upper tail dependence.An"alchemy" of Lorenz curves that can be used as general framework to build multiparametric families of copulas is also discussed.

We discuss a competitive alternative to stochastic local volatility models, namely the Collocating Volatility (CV) framework, introduced in [L. A. Grzelak (2019) The CLV framework-A fresh look at efficient pricing with smile, International Journal of Computer Mathematics 96 (11), 2209-2228]. The CV framework consists of two elements, a "kernel process"that can be efficiently evaluated and a local volatility function. The latter, based on stochastic collocation-e.g. [I. Babuška, F. Nobile & R. Tempone (2007) A stochastic collocation method for elliptic partial differential equations with random input data, SIAM Journal on Numerical Analysis 45 (3), 1005-1034; B. Ganapathysubramanian & N. Zabaras (2007) Sparse grid collocation schemes for stochastic natural convection problems, Journal of Computational Physics 225 (1), 652-685; J. A. S. Witteveen & G. Iaccarino (2012) Simplex stochastic collocation with random sampling and extrapolation for nonhypercube probability spaces, SIAM Journal on Scientific Computing 34 (2), A814-A838; D. Xiu & J. S. Hesthaven (2005) High-order collocation methods for differential equations with random inputs, SIAM Journal on Scientific Computing 27 (3), 1118-1139]-connects the kernel process to the market and allows the CV framework to be perfectly calibrated to European-type options. In this paper, we consider three different kernel process choices: The Ornstein-Uhlenbeck (OU) and Cox-Ingersoll-Ross (CIR) processes and the Heston model. The kernel process controls the forward smile and allows for an accurate and efficient calibration to exotic options, while the perfect calibration to liquid market quotes is preserved. We confirm this by numerical experiments, in which we calibrate the OU-CV, CIR-CV and Heston-CV frameworks to FX barrier options.

In the collocating volatility (CLV) model, the stochastic collocation technique is used as a convenient representation of the terminal distribution of the market option prices. A specific dynamic is added in the form of a stochastic driver process, which allows more control over the prices of forward starting options. This is reminiscent of the Markov functional models. (Grzelak uses a single-factor Ornstein–Uhlenbeck process as the driver for the CLV model, and Fries uses a single-factor Wiener process with time-dependent volatility in his equity Markov functional model. Van der Stoep et al consider a Heston stochastic volatility driver process and show that it offers more flexibility to capture the forward smile in the context of foreign exchange options.) In this paper, we discuss all aspects of derivative pricing under the Heston– CLV model: calibration with an efficient Fourier method; a Monte Carlo simulation with second-order convergence; and accurate partial differential equation pricing through implicit and explicit finite-difference methods.

In this work, we apply the Stochastic Grid Bundling Method (SGBM) to numerically solve backward stochastic differential equations (BSDEs). The SGBM algorithm is based on conditional expectations approximation by means of bundling of Monte Carlo sample paths and a local regress-later regression within each bundle. The basic algorithm for solving the backward stochastic differential equations will be introduced and an upper error bound is established for the local regression. A full error analysis is also conducted for the explicit version of our algorithm and numerical experiments are performed to demonstrate various properties of our algorithm.

We present a method for conditional time series forecasting based on an adaptation of the recent deep convolutional WaveNet architecture. The proposed network contains stacks of dilated convolutions that allow it to access a broad range of historical data when forecasting. It also uses a rectified linear unit (ReLU) activation function, and conditioning is performed by applying multiple convolutional filters in parallel to separate time series, which allows for the fast processing of data and the exploitation of the correlation structure between the multivariate time series. We test and analyze the performance of the convolutional network both unconditionally and conditionally for financial time series forecasting using the Standard & Poor’s 500 index, the volatility index, the Chicago Board Options Exchange interest rate and several exchange rates, and we extensively compare its performance with those of the well-known autoregressive model and a long short-term memory network. We show that a convolutional network is well suited to regression-type problems and is able to effectively learn dependencies in and between the series without the need for long historical time series, that it is a time-efficient and easy-to-implement alternative to recurrent-type networks, and that it tends to outperform linear and recurrent models.

We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that is modeled using the Richards equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared with the original MC estimator.