BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christian Bayer

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

ResearcherUnclaimed source record
math.NAmath.PRq-fin.CPNumerical Analysismath.OCmath-phmath.MPq-fin.MFComputational Engineering, Finance, and ScienceMachine Learningmath.STq-fin.PMq-fin.PRq-fin.TRStatistics Theory

Catalog footprint

What is connected

21works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

21 published item(s)

preprint2026arXiv

Dimension reduction for path signatures

This paper focuses on the mathematical framework for reducing the complexity of models using path signatures. The structure of these signatures, which can be interpreted as collections of iterated integrals along paths, is discussed and their applications in areas such as stochastic differential equations (SDEs) and financial modeling are pointed out. In particular, exploiting the rough paths view, solutions of SDEs continuously depend on the lift of the driver. Such continuous mappings can be approximated using (truncated) signatures, which are solutions of high-dimensional linear systems. In order to lower the complexity of these models, this paper presents methods for reducing the order of high-dimensional truncated signature models while retaining essential characteristics. The derivation of reduced models and the universal approximation property of (truncated) signatures are treated in detail. Numerical examples, including applications to the (rough) Bergomi model in financial markets, illustrate the proposed reduction techniques and highlight their effectiveness.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

A Pontryagin Maximum Principle on the Belief Space for Continuous-Time Optimal Control with Discrete Observations

We study a continuous time stochastic optimal control problem under partial observations that are available only at discrete time instants. This hybrid setting, with continuous dynamics and intermittent noisy measurements, arises in applications ranging from robotic exploration and target tracking to epidemic control. We formulate the problem on the space of beliefs (information states), treating the controller's posterior distribution of the state as the state variable for decision making. On this belief space we derive a Pontryagin maximum principle that provides necessary conditions for optimality. The analysis carefully tracks both the continuous evolution of the state between observation times and the Bayesian jump updates of the belief at observation instants. A key insight is a relationship between the adjoint process in our maximum principle and the gradient of the value functional on the belief space, which links the optimality conditions to the dynamic programming approach on the space of probability measures. The resulting optimality system has a prediction and update structure that is closely related to the unnormalised Zakai equation and the normalised Kushner-Stratonovich equation in nonlinear filtering. Building on this analysis, we design a particle based numerical scheme to approximate the coupled forward (filter) and backward (adjoint) system. The scheme uses particle filtering to represent the evolving belief and regression techniques to approximate the adjoint, which yields a practical algorithm for computing near optimal controls under partial information. The effectiveness of the approach is illustrated on both linear and nonlinear examples and highlights in particular the benefits of actively controlling the observation process.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Markovian approximations of stochastic Volterra equations with the fractional kernel

We consider rough stochastic volatility models where the variance process satisfies a stochastic Volterra equation with the fractional kernel, as in the rough Bergomi and the rough Heston model. In particular, the variance process is therefore not a Markov process or semimartingale, and has quite low Hölder-regularity. In practice, simulating such rough processes thus often results in high computational cost. To remedy this, we study approximations of stochastic Volterra equations using an $N$-dimensional diffusion process defined as solution to a system of ordinary stochastic differential equation. If the coefficients of the stochastic Volterra equation are Lipschitz continuous, we show that these approximations converge strongly with superpolynomial rate in $N$. Finally, we apply this approximation to compute the implied volatility smile of a European call option under the rough Bergomi and the rough Heston model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reinforced optimal control

Least squares Monte Carlo methods are a popular numerical approximation method for solving stochastic control problems. Based on dynamic programming, their key feature is the approximation of the conditional expectation of future rewards by linear least squares regression. Hence, the choice of basis functions is crucial for the accuracy of the method. Earlier work by some of us [Belomestny, Schoenmakers, Spokoiny, Zharkynbay. Commun.~Math.~Sci., 18(1):109-121, 2020](arXiv:1808.02341) proposes to reinforce the basis functions in the case of optimal stopping problems by already computed value functions for later times, thereby considerably improving the accuracy with limited additional computational cost. We extend the reinforced regression method to a general class of stochastic control problems, while considerably improving the method's efficiency, as demonstrated by substantial numerical examples as well as theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pricing high-dimensional Bermudan options with hierarchical tensor formats

An efficient compression technique based on hierarchical tensors for popular option pricing methods is presented. It is shown that the "curse of dimensionality" can be alleviated for the computation of Bermudan option prices with the Monte Carlo least-squares approach as well as the dual martingale method, both using high-dimensional tensorized polynomial expansions. This discretization allows for a simple and computationally cheap evaluation of conditional expectations. Complexity estimates are provided as well as a description of the optimization procedures in the tensor train format. Numerical experiments illustrate the favourable accuracy of the proposed methods. The dynamical programming method yields results comparable to recent Neural Network based methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hierarchical adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model

The rough Bergomi (rBergomi) model, introduced recently in [5], is a promising rough volatility model in quantitative finance. It is a parsimonious model depending on only three parameters, and yet remarkably fits with empirical implied volatility surfaces. In the absence of analytical European option pricing methods for the model, and due to the non-Markovian nature of the fractional driver, the prevalent option is to use the Monte Carlo (MC) simulation for pricing. Despite recent advances in the MC method in this context, pricing under the rBergomi model is still a time-consuming task. To overcome this issue, we have designed a novel, hierarchical approach, based on i) adaptive sparse grids quadrature (ASGQ), and ii) quasi-Monte Carlo (QMC). Both techniques are coupled with a Brownian bridge construction and a Richardson extrapolation on the weak error. By uncovering the available regularity, our hierarchical methods demonstrate substantial computational gains with respect to the standard MC method, when reaching a sufficiently small relative error tolerance in the price estimates across different parameter constellations, even for very small values of the Hurst parameter. Our work opens a new research direction in this field, i.e., to investigate the performance of methods other than Monte Carlo for pricing and calibrating under the rBergomi model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pricing Options Under Rough Volatility with Backward SPDEs

In this paper, we study the option pricing problems for rough volatility models. As the framework is non-Markovian, the value function for a European option is not deterministic; rather, it is random and satisfies a backward stochastic partial differential equation (BSPDE). The existence and uniqueness of weak solution is proved for general nonlinear BSPDEs with unbounded random leading coefficients whose connections with certain forward-backward stochastic differential equations are derived as well. These BSPDEs are then used to approximate American option prices. A deep leaning-based method is also investigated for the numerical approximations to such BSPDEs and associated non-Markovian pricing problems. Finally, the examples of rough Bergomi type are numerically computed for both European and American options.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Functional Limit Theorem for Limit Order Books with State Dependent Price Dynamics

We consider a stochastic model for the dynamics of the two-sided limit order book (LOB). Our model is flexible enough to allow for a dependence of the price dynamics on volumes. For the joint dynamics of best bid and ask prices and the standing buy and sell volume densities, we derive a functional limit theorem, which states that our LOB model converges in distribution to a fully coupled SDE-SPDE system when the order arrival rates tend to infinity and the impact of an individual order arrival on the book as well as the tick size tends to zero. The SDE describes the bid/ask price dynamics while the SPDE describes the volume dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From rough path estimates to multilevel Monte Carlo

New classes of stochastic differential equations can now be studied using rough path theory (e.g. Lyons et al. [LCL07] or Friz--Hairer [FH14]). In this paper we investigate, from a numerical analysis point of view, stochastic differential equations driven by Gaussian noise in the aforementioned sense. Our focus lies on numerical implementations, and more specifically on the saving possible via multilevel methods. Our analysis relies on a subtle combination of pathwise estimates, Gaussian concentration, and multilevel ideas. Numerical examples are given which both illustrate and confirm our findings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the probability density function of baskets

The state price density of a basket, even under uncorrelated Black-Scholes dynamics, does not allow for a closed from density. (This may be rephrased as statement on the sum of lognormals and is especially annoying for such are used most frequently in Financial and Actuarial Mathematics.) In this note we discuss short time and small volatility expansions, respectively. The method works for general multi-factor models with correlations and leads to the analysis of a system of ordinary (Hamiltonian) differential equations. Surprisingly perhaps, even in two asset Black-Scholes situation (with its flat geometry), the expansion can degenerate at a critical (basket) strike level; a phenomena which seems to have gone unnoticed in the literature to date. Explicit computations relate this to a phase transition from a unique to more than one "most-likely" paths (along which the diffusion, if suitably conditioned, concentrates in the afore-mentioned regimes). This also provides a (quantifiable) understanding of how precisely a presently out-of-money basket option may still end up in-the-money.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Efficient Forward-Reverse Expectation-Maximization Algorithm for Statistical Inference in Stochastic Reaction Networks

In this work, we present an extension to the context of Stochastic Reaction Networks (SRNs) of the forward-reverse representation introduced in "Simulation of forward-reverse stochastic representations for conditional diffusions", a 2014 paper by Bayer and Schoenmakers. We apply this stochastic representation in the computation of efficient approximations of expected values of functionals of SNR bridges, i.e., SRNs conditioned to its values in the extremes of given time-intervals. We then employ this SNR bridge-generation technique to the statistical inference problem of approximating the reaction propensities based on discretely observed data. To this end, we introduce a two-phase iterative inference method in which, during phase I, we solve a set of deterministic optimization problems where the SRNs are replaced by their reaction-rate Ordinary Differential Equations (ODEs) approximation; then, during phase II, we apply the Monte Carlo version of the Expectation-Maximization (EM) algorithm starting from the phase I output. By selecting a set of over dispersed seeds as initial points for phase I, the output of parallel runs from our two-phase method is a cluster of approximate maximum likelihood estimates. Our results are illustrated by numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational error estimates for Born-Oppenheimer molecular dynamics with nearly crossing potential surfaces

The difference of the values of observables for the time-independent Schroedinger equation, with matrix valued potentials, and the values of observables for ab initio Born-Oppenheimer molecular dynamics, of the ground state, depends on the probability to be in excited states and the electron/nuclei mass ratio. The paper first proves an error estimate (depending on the electron/nuclei mass ratio and the probability to be in excited states) for this difference of microcanonical observables, assuming that molecular dynamics space-time averages converge, with a rate related to the maximal Lyapunov exponent. The error estimate is uniform in the number of particles and the analysis does not assume a uniform lower bound on the spectral gap of the electron operator and consequently the probability to be in excited states can be large. A numerical method to determine the probability to be in excited states is then presented, based on Ehrenfest molecular dynamics and stability analysis of a perturbed eigenvalue problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forward-reverse EM algorithm for Markov chains: convergence and numerical analysis

We develop a forward-reverse EM (FREM) algorithm for estimating parameters that determine the dynamics of a discrete time Markov chain evolving through a certain measurable state space. As a key tool for the construction of the FREM method we develop forward-reverse representations for Markov chains conditioned on a certain terminal state. These representations may be considered as an extension of the earlier work Bayer and Schoenmakers [2013] on conditional diffusions. We proof almost sure convergence of our algorithm for a Markov chain model with curved exponential family structure. On the numerical side we give a complexity analysis of the forward-reverse algorithm by deriving its expected cost. Two application examples are discuss to demonstrate the scope of possible applications ranging from models based on continuous time processes to discrete time Markov chain models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Option pricing in affine generalized Merton models

In this article we consider affine generalizations of the Merton jump diffusion model [Merton, J. Fin. Econ., 1976] and the respective pricing of European options. On the one hand, the Brownian motion part in the Merton model may be generalized to a log-Heston model, and on the other hand, the jump part may be generalized to an affine process with possibly state dependent jumps. While the characteristic function of the log-Heston component is known in closed form, the characteristic function of the second component may be unknown explicitly. For the latter component we propose an approximation procedure based on the method introduced in [Belomestny et al., J. Func. Anal., 2009]. We conclude with some numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulation of forward-reverse stochastic representations for conditional diffusions

In this paper we derive stochastic representations for the finite dimensional distributions of a multidimensional diffusion on a fixed time interval, conditioned on the terminal state. The conditioning can be with respect to a fixed point or more generally with respect to some subset. The representations rely on a reverse process connected with the given (forward) diffusion as introduced in Milstein, Schoenmakers and Spokoiny [Bernoulli 10 (2004) 281-312] in the context of a forward-reverse transition density estimator. The corresponding Monte Carlo estimators have essentially root-$N$ accuracy, and hence they do not suffer from the curse of dimensionality. We provide a detailed convergence analysis and give a numerical example involving the realized variance in a stochastic volatility asset model conditioned on a fixed terminal value of the asset.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cubature on Wiener space: pathwise convergence

Cubature on Wiener space [Lyons, T.; Victoir, N.; Proc. R. Soc. Lond. A 8 January 2004 vol. 460 no. 2041 169-198] provides a powerful alternative to Monte Carlo simulation for the integration of certain functionals on Wiener space. More specifically, and in the language of mathematical finance, cubature allows for fast computation of European option prices in generic diffusion models. We give a random walk interpretation of cubature and similar (e.g. the Ninomiya--Victoir) weak approximation schemes. By using rough path analysis, we are able to establish weak convergence for general path-dependent option prices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Utility Maximization in a Binomial Model with transaction costs: a Duality Approach Based on the Shadow Price Process

We consider the problem of optimizing the expected logarithmic utility of the value of a portfolio in a binomial model with proportional transaction costs with a long time horizon. By duality methods, we can find expressions for the boundaries of the no-trade-region and the asymptotic optimal growth rate, which can be made explicit for small transaction costs. Here we find that, contrary to the classical results in continuous time, the size of the no-trade-region as well as the asymptotic growth rate depend analytically on the level of transaction costs, implying a linear first order effect of perturbations of (small) transaction costs. We obtain the asymptotic expansion by an almost explicit construction of the shadow price process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How accurate is molecular dynamics?

Born-Oppenheimer dynamics is shown to provide an accurate approximation of time-independent Schrödinger observables for a molecular system with an electron spectral gap, in the limit of large ratio of nuclei and electron masses, without assuming that the nuclei are localized to vanishing domains. The derivation, based on a Hamiltonian system interpretation of the Schrödinger equation and stability of the corresponding Hamilton-Jacobi equation, bypasses the usual separation of nuclei and electron wave functions, includes caustic states and gives a different perspective on the Born-Oppenheimer approximation, Schrödinger Hamiltonian systems and numerical simulation in molecular dynamics modeling at constant energy microcanonical ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Semi-Closed Form Cubature and Applications to Financial Diffusion Models

Cubature methods, a powerful alternative to Monte Carlo due to Kusuoka~[Adv.~Math.~Econ.~6, 69--83, 2004] and Lyons--Victoir~[Proc.~R.~Soc.\\Lond.~Ser.~A 460, 169--198, 2004], involve the solution to numerous auxiliary ordinary differential equations. With focus on the Ninomiya-Victoir algorithm~[Appl.~Math.~Fin.~15, 107--121, 2008], which corresponds to a concrete level $5$ cubature method, we study some parametric diffusion models motivated from financial applications, and exhibit structural conditions under which all involved ODEs can be solved explicitly and efficiently. We then enlarge the class of models for which this technique applies, by introducing a (model-dependent) variation of the Ninomiya-Victoir method. Our method remains easy to implement; numerical examples illustrate the savings in computation time.

0 citations0 reviewsNo code linkedOpen access