Source author record

Kathrin Glau

Kathrin Glau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

Stability and convergence of Galerkin schemes for parabolic equations with application to Kolmogorov pricing equations in time-inhomogeneous Lévy models

Two essential quantities for the analysis of approximation schemes of evolution equations are stability and convergence. We derive stability and convergence of fully discrete approximation schemes of solutions to linear parabolic evolution equations governed by time dependent coercive operators. We consider abstract Galerkin approximations in space combined with theta-schemes in time. The level of generality of our analysis comprises both a large class of time-dependent operators and a large choice of approximating Galerkin spaces. In particular the results apply to partial integro differential equations for option pricing in time-inhomogeneous Lévy models and allows for a large variety of option types and models. The derivation builds on the strong foundation laid out by von Petersdorff and Schwab (2003) who provide the respective results for the time-homogeneous case. We discuss the assumptions in the context of option pricing.

preprint2016arXiv

A Flexible Galerkin Scheme for Option Pricing in Lévy Models

One popular approach to option pricing in Lévy models is through solving the related partial integro differential equation (PIDE). For the numerical solution of such equations powerful Galerkin methods have been put forward e.g. by Hilber et al. (2013). As in practice large classes of models are maintained simultaneously, flexibility in the driving Lévy model is crucial for the implementation of these powerful tools. In this article we provide such a flexible finite element Galerkin method. To this end we exploit the Fourier representation of the infinitesimal generator, i.e. the related symbol, which is explicitly available for the most relevant Lévy models. Empirical studies for the Merton, NIG and CGMY model confirm the numerical feasibility of the method.

preprint2016arXiv

A unified view of LIBOR models

We provide a unified framework for modeling LIBOR rates using general semimartingales as driving processes and generic functional forms to describe the evolution of the dynamics. We derive sufficient conditions for the model to be arbitrage-free which are easily verifiable, and for the LIBOR rates to be true martingales under the respective forward measures. We discuss when the conditions are also necessary and comment on further desirable properties such as those leading to analytical tractability and positivity of rates. This framework allows to consider several popular models in the literature, such as LIBOR market models driven by Brownian motion or jump processes, the Lévy forward price model as well as the affine LIBOR model, under one umbrella. Moreover, we derive structural results about LIBOR models and show, in particular, that only models where the forward price is an exponentially affine function of the driving process preserve their structure under different forward measures.

preprint2016arXiv

Calibration to American Options: Numerical Investigation of the de-Americanization

American options are the reference instruments for the model calibration of a large and important class of single stocks. For this task, a fast and accurate pricing algorithm is indispensable. The literature mainly discusses pricing methods for American options that are based on Monte Carlo, tree and partial differential equation methods. We present an alternative approach that has become popular under the name de-Americanization in the financial industry. The method is easy to implement and enjoys fast run-times. Since it is based on ad hoc simplifications, however, theoretical results guaranteeing reliability are not available. To quantify the resulting methodological risk, we empirically test the performance of the de-Americanization method for calibration. We classify the scenarios in which de-Americanization performs very well. However, we also identify the cases where de-Americanization oversimplifies and can result in large errors.

preprint2016arXiv

Chebyshev Interpolation for Parametric Option Pricing

Recurrent tasks such as pricing, calibration and risk assessment need to be executed accurately and in real-time. Simultaneously we observe an increase in model sophistication on the one hand and growing demands on the quality of risk management on the other. To address the resulting computational challenges, it is natural to exploit the recurrent nature of these tasks. We concentrate on Parametric Option Pricing (POP) and show that polynomial interpolation in the parameter space promises to reduce run-times while maintaining accuracy. The attractive properties of Chebyshev interpolation and its tensorized extension enable us to identify criteria for (sub)exponential convergence and explicit error bounds. We show that these results apply to a variety of European (basket) options and affine asset models. Numerical experiments confirm our findings. Exploring the potential of the method further, we empirically investigate the efficiency of the Chebyshev method for multivariate and path-dependent options.

preprint2016arXiv

Improved error bound for multivariate Chebyshev polynomial interpolation

Chebyshev interpolation is a highly effective, intensively studied method and enjoys excellent numerical properties. The interpolation nodes are known beforehand, implementation is straightforward and the method is numerically stable. For efficiency, a sharp error bound is essential, in particular for high-dimensional applications. For tensorized Chebyshev interpolation, we present an error bound that improves existing results significantly.

preprint2016arXiv

Magic points in finance: Empirical integration for parametric option pricing

We propose an offline-online procedure for Fourier transform based option pricing. The method supports the acceleration of such essential tasks of mathematical finance as model calibration, real-time pricing, and, more generally, risk assessment and parameter risk estimation. We adapt the empirical magic point interpolation method of Barrault, Nguyen, Maday and Patera (2004) to parametric Fourier pricing. In the offline phase, a quadrature rule is tailored to the family of integrands of the parametric pricing problem. In the online phase, the quadrature rule then yields fast and accurate approximations of the option prices. Under analyticity assumptions the pricing error decays exponentially. Numerical experiments in one dimension confirm our theoretical findings and show a significant gain in efficiency, even for examples beyond the scope of the theoretical results.

preprint2016arXiv

Martingale Property in Terms of Semimartingale Problems

Starting from the seventies mathematicians face the question whether a non-negative local martingale is a true or a strict local martingale. In this article we answer this question from a semimartingale perspective. We connect the martingale property to existence, uniqueness and topological properties of semimartingale problems. This not only leads to valuable characterizations of the martingale property, but also reveals new existence and uniqueness results for semimartingale problems. As a case study we derive explicit conditions for the martingale property of stochastic exponentials driven by infinite-dimensional Brownian motion.

preprint2016arXiv

Martingale property of exponential semimartingales: a note on explicit conditions and applications to financial models

We give a collection of explicit sufficient conditions for the true martingale property of a wide class of exponentials of semimartingales. We express the conditions in terms of semimartingale characteristics. This turns out to be very convenient in financial modeling in general. Especially it allows us to carefully discuss the question of well-definedness of semimartingale Libor models, whose construction crucially relies on a sequence of measure changes.

preprint2016arXiv

Model reduction for calibration of American options

American put options are among the most frequently traded single stock options, and their calibration is computationally challenging since no closed-form expression is available. Due to the higher flexibility in comparison to European options, the mathematical model involves additional constraints, and a variational inequality is obtained. We use the Heston stochastic volatility model to describe the price of a single stock option. In order to speed up the calibration process, we apply two model reduction strategies. Firstly, a reduced basis method (RBM) is used to define a suitable low-dimensional basis for the numerical approximation of the parameter-dependent partial differential equation ($μ$PDE) model. By doing so the computational complexity for solving the $μ$PDE is drastically reduced, and applications of standard minimization algorithms for the calibration are significantly faster than working with a high-dimensional finite element basis. Secondly, so-called de-Americanization strategies are applied. Here, the main idea is to reformulate the calibration problem for American options as a problem for European options and to exploit closed-form solutions. Both reduction techniques are systematically compared and tested for both synthetic and market data sets.

preprint2015arXiv

Feynman-Kac formula for Lévy processes with discontinuous killing rate

The challenge to fruitfully merge state-of-the-art techniques from mathematical finance and numerical analysis has inspired researchers to develop fast deterministic option pricing methods. As a result, highly efficient algorithms to compute option prices in Lévy models by solving partial integro differential equations have been developed. In order to provide a solid mathematical foundation for these methods, we derive a Feynman-Kac representation of variational solutions to partial integro differential equations that characterize conditional expectations of functionals of killed time-inhomogeneous Lévy processes. We allow for a wide range of underlying stochastic processes, comprising processes with Brownian part, and a broad class of pure jump processes such as generalized hyperbolic, multivariate normal inverse Gaussian, tempered stable, and $α$-semi stable Lévy processes. By virtue of our mild regularity assumptions as to the killing rate and the initial condition of the partial differential equation, our results provide a rigorous basis for numerous applications, not only in financial mathematics but also in probability theory and relativistic quantum mechanics.

preprint2015arXiv

Parametric Integration by Magic Point Empirical Interpolation

We derive analyticity criteria for explicit error bounds and an exponential rate of convergence of the magic point empirical interpolation method introduced by Barrault et al. (2004). Furthermore, we investigate its application to parametric integration. We find that the method is well-suited to Fourier transforms and has a wide range of applications in such diverse fields as probability and statistics, signal and image processing, physics, chemistry and mathematical finance. To illustrate the method, we apply it to the evaluation of probability densities by parametric Fourier inversion. Our numerical experiments display convergence of exponential order, even in cases where the theoretical results do not apply.

preprint2012arXiv

Sobolev index: A classification of Lévy processes

We classify Lévy processes according to the solution spaces of the associated parabolic PIDEs. This classification reveals structural characteristics of the processes and is relevant for applications such as for solving PIDEs numerically for pricing options in Lévy models. The classification is done via the Fourier transform i.e. via the symbol of the process. We define the Sobolev index of a Lévy process by a certain growth condition on the symbol. It follows that for Lévy processes with Sobolev index $α$ the corresponding evolution problem has a unique weak solution in the Sobolev-Slobodeckii space $H^{α/2}$. We show that this classification applies to a wide range of processes. Examples are the Brownian motion with or without drift, generalised hyperbolic (GH), CGMY and (semi) stable Lévy processes. A comparison of the Sobolev index with the Blumenthal-Getoor index sheds light on the structural implication of the classification. More precisely, we discuss the Sobolev index as an indicator of the smoothness of the distribution and of the variation of the paths of the process. This highlights the relation between the $p$-variation of the paths and the degree of smoothing effect that stems from the distribution.

preprint2010arXiv

Analyticity of the Wiener-Hopf factors and valuation of exotic options in Lévy models

This paper considers the valuation of exotic path-dependent options in Lévy models, in particular options on the supremum and the infimum of the asset price process. Using the Wiener--Hopf factorization, we derive expressions for the analytically extended characteristic function of the supremum and the infimum of a Lévy process. Combined with general results on Fourier methods for option pricing, we provide formulas for the valuation of one-touch options, lookback options and equity default swaps in Lévy models.

preprint2009arXiv

Analysis of Fourier transform valuation formulas and applications

The aim of this article is to provide a systematic analysis of the conditions such that Fourier transform valuation formulas are valid in a general framework; i.e. when the option has an arbitrary payoff function and depends on the path of the asset price process. An interplay between the conditions on the payoff function and the process arises naturally. We also extend these results to the multi-dimensional case, and discuss the calculation of Greeks by Fourier transform methods. As an application, we price options on the minimum of two assets in Lévy and stochastic volatility models.