21st Winter school on Mathematical Finance
Abstracts



Minicourses

Christa Cuchiero: Signature methods in finance

Signature methods represent a non-parametric way for extracting characteristic features from time series data which is essential in machine learning tasks. This explains why these techniques become more and more popular in Econometrics and Mathematical Finance. Indeed, signature based approaches allow for data-driven and thus more robust model selection mechanisms, while first principles like no arbitrage can still be easily guaranteed.
In this course we shall focus on the use of signature as universal linear regression basis of continuous functionals of paths for financial applications. We first give a brief introduction to continuous rough paths and show how to embed continuous semimartingales into the rough path setting. Indeed our main focus lies on signature of semimartingales, one of the main modeling tools in finance. By relying on the Stone-Weierstrass theorem we show how to prove the universal approximation property of linear functions of the signature in appropriate topologies on path space.
In the financial applications that we have in mind one key quantity that one needs to compute is the expected signature of some underlying process. Surprisingly this can be achieved for generic classes of diffusions, called signature-SDEs (with possibly path dependent characteristics), via techniques from affine and polynomial processes. More precisely, we show how the signature process of these diffusions can be embedded in the framework of affine and polynomial processes. These classes of processes have been -- due to their tractability -- the dominating process class prior to the new era of highly over-parametrized dynamic models. Following this line we obtain that the infinite dimensional Feynman Kac PDE of the signature process can generically be reduced to an infinite dimensional ODE either of Riccati or linear type.
In terms of financial applications, we shall treat two main topics: stochastic portfolio theory and signature based asset price models.
In the context of stochastic portfolio theory we introduce a novel class of portfolios which we call linear path-functional portfolios. These are portfolios which are determined by certain transformations of linear functions of a collections of feature maps that are non-anticipative path functionals of an underlying semimartingale. As main example for such feature maps we consider signature of the (ranked) market weights. Relying on the universal approximation theorem we show that every continuous (possibly path-dependent) portfolio function of the market weights can be uniformly approximated by signature portfolios. Besides these universality features, the main numerical advantage lies in the fact that several optimization tasks like maximizing expected logarithmic utility or mean-variance optimization within the class of linear path-functional portfolios reduces to a convex quadratic optimization problem, thus making it computationally highly tractable. We apply our method to real market data, indicating out-performance on the considered out-of-sample data even under transaction costs.
In view of asset price models we consider a stochastic volatility model where the dynamics of the volatility are described by linear functions of the (time extended) signature of a primary underlying process, which is supposed to be some multidimensional continuous semimartingale. Under the additional assumption that this primary process is of polynomial type, we obtain closed form expressions for the VIX squared, exploiting the fact that the truncated signature of a polynomial process is again a polynomial process. Adding to such a primary process the Brownian motion driving the stock price, allows then to express both the log-price and the VIX squared as linear functions of the signature of the corresponding augmented process. This feature can then be efficiently used for pricing and calibration purposes. Indeed, as the signature samples can be easily precomputed, the calibration task can be split into an offline sampling and a standard optimization. For both the SPX and VIX options we obtain highly accurate calibration results, showing that this model class allows to solve the joint calibration problem without adding jumps or rough volatility. Slides: Slides 1, Slides 2, Slides 3, Slides 4
References
[1] P. Friz, Peter and N. Victoir. Multidimensional stochastic processes as rough paths: theory and applications. Vol. 120. Cambridge University Press, 2010, in particular chapters 7-9.
[2] I. Chevyrev and A. Kormilitzin. A primer on the signature method in machine learning. arXiv:1603.03788 (2016).
[3] C. Cuchiero, S. Svaluto-Ferro, and J.Teichmann. Signature SDEs from an affine and polynomial perspective. arXiv:2302.01362 (2023).
[4] C. Cuchiero, F. Primavera, and S. Svaluto-Ferro. Universal approximation theorems for continuous functions of càdlàg paths and Lévy-type signature models. arXiv:2208.02293 (2022).
[5] C. Cuchiero and J. Möller. Signature Methods in Stochastic Portfolio Theory. arXiv:2310.02322 (2023).
[6] C. Cuchiero, G. Gazzani, and S. Svaluto-Ferro. Signature-based models: Theory and calibration. SIAM journal on financial mathematics 14.3 (2023): 910-957.
[7] C. Cuchiero, G. Gazzani, J. Möller and S. Svaluto-Ferro. Joint calibration to SPX and VIX options with signature-based models. arXiv:2301.13235 (2023).

Carole Bernard: Robust risk management

Making sound decisions under uncertainty generally requires quantitative analysis and the use of models. However, a “perfect” model does not exist since some divergence between the model and the reality it attempts to describe cannot be avoided. In a broad sense, model risk is about the extent to which the quality of model-based decisions is sensitive to underlying model deviations and data issues. An example concerns the establishment of the capital buffers banks need to put aside to absorb unforeseen losses for a portfolio of risky loans. Doing so requires accurate estimates of the likelihood that various obligors default together, which is very difficult due to a scarcity of data. We will discuss various approaches on how to plan for worst-case scenarios in a financial and insurance context. More specifically, we will discuss the problem of quantifying model uncertainty, such as departures from assumed independence, incomplete dependence information, factor models that are only partially specified, or portfolio information that is only available on an aggregate level (e.g., mean and variance of the portfolio loss), desired properties (e.g., unimodality, symmetry, non-negativity). We will develop necessary tools to quantify this model uncertainty and in particular to determine the best upper and lower risk bounds for various risk aggregation functionals of interest including Value-at-Risk, Range Value-at-Risk and distorted risk measures.
Slides: Slides 1, Slides 2, Slides 3
References
[1] L. Rüschendorf , S. Vanduffel, C. Bernard. Model Risk Management: Risk Bounds under Uncertainty. Cambridge University Press; 2024. Link.
[2] C. Bernard, S. Pesenti , S. Vanduffel. Robust Distortions Measures. 2023, Mathematical Finance, forthcoming,
available at SSRN, published in an open access.
[3] C. Bernard, A. Müller (University of Siegen) and M. Oesting. \(L_p\)-norm spherical copulas, Journal of Multivariate Analysis, 2023, forthcoming,
arXiv:2206.10180.
[4] C. Bernard, C. De Vecchi, S. Vanduffel. Impact of Correlation on the (Range) Value-at-Risk. 2023, Scandinavian Actuarial Journal , (6), 531–564, published version.
[5] C. Bernard, R. Kazzi , S. Vanduffel. Range Value-at-Risk Bounds for Unimodal Distributions under Partial Information. 2020, Insurance: Mathematics and Economics, 94, 9-24,
published version, available at SSRN.
[6] C. Bernard, O. Bondarenko (UIC) and S. Vanduffel. Rearrangement Algorithm and Maximum Entropy. 2018, Annals of Operational Research, 261(1-2), 107-134,
PDF in open access, available at SSRN.
[7] C. Bernard, L. Rüschendorf , S. Vanduffel, R. Wang. Risk bounds for factor models. 2017, Finance and Stochastics , 21(3), 631-659,
published version, available at SSRN.
[8] C. Bernard, L. Rüschendorf , S. Vanduffel. VaR Bounds with Variance Constraint. 2017, Journal of Risk and Insurance , 84(3), 923-959,
published version, available at SSRN.
[9] C. Bernard, X. Jiang and R. Wang. Risk Aggregation with Dependence Uncertainty. 2014, Insurance: Mathematics and Economics , 54, 93-108,
published version.

Special invited lectures

Griselda Deelstra: Some topics related to stochastic mortality and/or interest rates in the valuation of life insurance products

In the main part of this talk, we focus on death-linked contingent claims (GMDBs) paying a random financial return at a random time of death in the general case where financial returns follow a regime switching model with two-sided phase-type jumps. We approximate the distribution of the remaining lifetime by either a series of Erlang distributions or a Laguerre series expansion. More precisely, we develop a Laurent series expansion of the discounted Laplace transform of the subordinated process at an Erlang distributed time, which leads to explicit formulae for European-type GMDB as well as related risk measures such as the Value-at-Risk (VaR) and the Conditional-Tail-Expectation (CTE). We further concentrate upon path-dependent GMDBs with lookback features like dynamic fund protection or dynamic withdrawal benefits, by relying on a Sylvester equation approach. The advantage of our approaches is that our results are of semi-closed form, avoiding numerical Fourier inversion or Monte-Carlo simulation, leading to fast evaluation. These results have implications beyond life-insurance and GMDBs, namely in all situations where randomization or Erlangization replaces known quantities. In Finance, it is for example well-known that a random maturity time leads to much more convenient valuation formulas which well approximate its non-random counterpart. If time allows, we focus on the relaxation of the traditional assumption of independence between mortality risk and interest rate risk. We investigate the impact of the inclusion of correlation on the best estimate of usual life insurance contracts and study in which cases the non-taking into account of correlation can lead to underestimations of the best estimate, and where one could have hedging possibilities in the presence of correlation. Therefore, we first focus upon a framework of two correlated Hull and White processes for modelling the stochastic mortality and interest rates. (Slides)
This talk is based on joint work with Peter Hieber (University of Lausanne); and with Pierre Devolder (UCLouvain) & Benjamin Roelants du Vivier (ULB).
References
[1] Deelstra, G., Hieber, P., (2023), Randomization and the valuation of guaranteed minimum death benefits, European journal of operational research, Vol. 309, Issue 3, 1218-1236, open access link.
[2] Deelstra, G., Devolder, P., Roelants du Vivier, B. (2024), Impact of correlation between interest rates and mortality rates on the valuation of various life insurance products, Working paper.

Claudio Fontana: Term structure modelling beyond stochastic continuity

Overnight rates, such as the Secured Overnight Financing Rate (SOFR), are central to the current reform of interest rate benchmarks. A striking feature of overnight rates is the presence of jumps and spikes occurring at predetermined dates due to monetary policy interventions and liquidity constraints. This corresponds to stochastic discontinuities (i.e., discontinuities occurring at ex ante known points in time) in their dynamics. In this work, we propose a generalized Heath-Jarrow-Morton (HJM) setup allowing for stochastic discontinuities and characterize absence of arbitrage. We extend the classical short-rate approach to accommodate stochastic discontinuities, developing a tractable setup driven by affine semimartingales. In a Gaussian setting, we provide explicit valuation formulas for bonds and caplets. Furthermore, we investigate hedging in the sense of local risk-minimization when the underlying term structures feature stochastic discontinuities. Time permitting, we shall explore applications beyond the interest rate setting.
Based on joint work with Z. Grbac and T. Schmidt. (Slides)
References
[1] C. Fontana, Z. Grbac, S. Gümbel, T. Schmidt (2020), Term structure modeling for multiple curves with stochastic discontinuities. Finance Stoch. 24, 465 - 511.
[2] C. Fontana, Z. Grbac, T. Schmidt (2023), Term structure modeling with overnight rates beyond stochastic continuity Mathematical Finance, forthcoming.

Antonis Papapantoleon: Model-free and data driven methods in mathematical finance

Academics, practitioners and regulators have understood that the classical paradigm in mathematical finance, where all computations are based on a single "correct" model, is flawed. Model-free methods, were computations are based on a variety of models, offer an alternative. More recently, these methods are driven by information available in financial markets. In this talk, we will discuss model-free and data driven methods and bounds and present how ideas from probability, statistics, optimal transport and optimization can be applied in this field. (Slides)
References
[1] A. Neufeld, A. Papapantoleon, Q. Xiang: Model-free bounds for multi-asset options using option-implied information and their exact computation. Management Science 69, 2051–2068, 2023. arXiv 2006.14288.
[2] A. Papapantoleon, P. Yanez Sarmiento: Detection of arbitrage opportunities in multi-asset derivatives markets. Dependence Modeling 9, 439–459, 2021. arXiv 2002.06227.
[3] D. Bartl, M. Kupper, T. Lux, A. Papapantoleon, S. Eckstein (appendix): Marginal and dependence uncertainty: bounds, optimal transport, and sharpness. SIAM Journal on Control and Optimization 60, 410–434, 2022. arXiv 1709.00641.
[4] T. Lux, A. Papapantoleon: Model-free bounds on Value-at-Risk using extreme value information and statistical distances. Insurance: Mathematics and Economics 86, 73–83, 2019. arXiv 1610.09734.
[5] T. Lux, A. Papapantoleon: Improved Fréchet–Hoeffding bounds for d-copulas and applications in model-free finance. Annals of Applied Probability 27, 3633–3671, 2017. arXiv:1602.08894.

Short lectures

Kristoffer Andersson: Convergence of a robust deep FBSDE method for stochastic control"

In this talk, a deep learning based numerical scheme for strongly coupled FBSDEs is presented. More specifically, problems stemming from stochastic control problems are considered and the method is a modification of the deep BSDE method in which the initial value of the backward equation is not a free parameter, and with a new loss function being the weighted sum of the cost of the control problem, and a variance term which coincides with the mean squared error in the terminal condition. It is demonstrated by a numerical example that a direct extension of the classical deep BSDE method to FBSDEs, fails for a simple linear-quadratic control problem, and motivate why the new method works. Under regularity and boundedness assumptions on the exact controls of time continuous and time discrete control problems, an error analysis for our method is provided. Finally, empirical experiments demonstrate that the method converges for three different problems, one being the one that failed for a direct extension of the deep BSDE method.

Jori Hoencamp: A static replication approach for callable interest rate derivatives: Efficient estimation of forward prices and SIMM-MVA

The computation of credit risk measures such as exposure and Credit Valuation Adjustments (CVA) requires the simulation of future portfolio prices. Recent metrics, such as dynamic Initial Margin (IM) and Margin Valuation Adjustments (MVA) additionally require the simulation of future conditional sensitivities. For portfolios with non-linear instruments that do not admit closed-form valuation formulas, this poses a significant computational challenge. This problem is addressed by proposing a static replication algorithm for interest rate options with early-exercise features under an affine term-structure model. Under the appropriate conditions we can find an equivalent portfolio of vanilla options that replicate these products. Specifically we decompose the product into a portfolio of European swaptions. The weights and strikes of the portfolio are obtained by regressing the target option value with interpretable, feed-forward neural networks. Once an equivalent portfolio of European swaptions is determined, we can leverage on closed-form expressions to obtain the conditional prices and sensitivities, which serve as an input to exposure and SIMM- driven MVA quantification. For a consistent forward sensitivity estimation, this in- volves the differentiation of the portfolio-weights. The accuracy and convergence of the method is demonstrated through several representative numerical examples, benchmarked against the established least-square Monte Carlo method.

Thomas van der Zwaard: Valuation Adjustments with an Affine-Diffusion-based Interest Rate Smile

Affine Diffusion (AD) dynamics are frequently used for Valuation Adjustments (xVA) calculations due to their analytic tractability. However, these models cannot capture the market-implied skew and smile, which are relevant when computing the xVA metrics. Hence, additional degrees of freedom are required to capture these market features. In this paper, we address this through an SDE with state-dependent coefficients. The SDE is consistent with the convex combination of a finite number of different AD dynamics. In this paper, we combine Hull-White one-factor (HW) models where one model parameter is varied. We use the Randomized AD (RAnD) technique to parameterize the combination of HW models. We refer to our SDE with state-dependent coefficients and the RAnD parametrization of the HW models as the rHW model. The rHW model allows for fast semi-analytic calibration to European swaptions through the analytic tractability of the HW dynamics. A regression-based Monte Carlo simulation is used to calculate exposures. In this setting, we assess the effect of skew and smile on exposures of interest rate derivatives. (slides)

Evgenii Vladimirov: Estimating Option Pricing Models Using a Characteristic Function-Based Linear State Space Representation.

We develop a novel filtering and estimation procedure for parametric option pricing models driven by general affine jump-diffusions. Our procedure is based on the comparison between an option-implied, model-free representation of the conditional log-characteristic function and the model-implied conditional log-characteristic function, which is functionally affine in the model’s state vector. We formally derive an associated linear state space representation and establish the asymptotic properties of the corresponding measurement errors. The state space representation allows us to use a suitably modified Kalman filtering technique to learn about the latent state vector and a quasi-maximum likelihood estimator of the model parameters, which brings important computational advantages. We analyze the finite-sample behavior of our procedure in Monte Carlo simulations. The applicability of our procedure is illustrated in two case studies that analyze S&P 500 option prices and the impact of exogenous state variables capturing Covid-19 reproduction and economic policy uncertainty (joint work with Peter Boswijk and Roger Laeven)
Reference
H. P. Boswijk, R. J. A. Laeven, E. Vladimirov: Estimating Option Pricing Models Using a Characteristic Function-Based Linear State Space Representation. 2022. arXiv:2210.06217.

Poster presentations

  • Purba Banerjee (Indian Institute of Science): Multiperiod static hedging of European options
  • Younhee Lee (Department of Mathematics, Chungnam National University): Real option pricing under the regime-switching model with jumps on a finite-time horizon
  • Ivo Richert (Kiel University): Quasi-maximum likelihood estimation of partially observed affine and polynomial processes
  • Benjamin Robinson (Universität Wien): Bicausal optimal transport for SDEs with irregular coefficients
  • Henrik Valett (Kiel University): Parameter estimation for polynomial processes
  • Mark van den Bosch (Leiden University): Multidimensional stability of planar traveling waves for stochastically forced reaction-diffusion systems

To the homepage of the Winter School on Mathematical Finance