Optimal Transport: A Semi-Discrete Approach

Assume that an agent models a financial asset through a measure Q with the goal to price / hedge some derivative or optimize some expected utility. Even if the model Q is chosen in the most skilful and sophisticated way, she is left with the possibility that Q does not provide an "exact" description of reality. This leads us to the following question: will the hedge still be somewhat meaningful for models in the proximity of Q? 
If we measure proximity with the usual Wasserstein distance (say), the answer is NO. Models which are similar w.r.t. Wasserstein distance may provide dramatically different information on which to base a hedging strategy. 
Remarkably, this can be overcome by considering a suitable "adapted" version of the Wasserstein distance which takes the temporal structure of pricing models into account. This adapted Wasserstein distance is most closely related to the nested distance as pioneered by Pflug and Pichler \cite{Pf09,PfPi12,PfPi14}. It allows us to establish Lipschitz properties of hedging strategies for semimartingale models in discrete and continuous time. Notably, these abstract results are sharp already for Brownian motion and European call options.

Adapted Wasserstein Distances and Stability in Mathematical Finance

Assume that an agent models a financial asset through a measure ℚ with the goal to price/hedge some derivative or optimise some expected utility. Even if the model ℚ is chosen in the most skilful and sophisticated way, the agent is left with the possibility that ℚ does not provide an exact description of reality. This leads us to the following question: will the hedge still be somewhat meaningful for models in the proximity of ℚ? If we measure proximity with the usual Wasserstein distance (say), the answer is No. Models which are similar with respect to the Wasserstein distance may provide dramatically different information on which to base a hedging strategy. Remarkably, this can be overcome by considering a suitable adapted version of the Wasserstein distance which takes the temporal structure of pricing models into account. This adapted Wasserstein distance is most closely related to the nested distance as pioneered by Pflug and Pichler (SIAM J. Optim. 20:1406–1420, 2009, SIAM J. Optim. 22:1–23, 2012, Multistage Stochastic Optimization, 2014). It allows us to establish Lipschitz properties of hedging strategies for semimartingale models in discrete and continuous time. Notably, these abstract results are sharp already for Brownian motion and European call options.

/pdf/adapted-wasserstein-distances-and-stability-in-mathematical-1c82e5l40n.pdf

Adapted Wasserstein distances and stability in mathematical finance

The optimal weak transport problem has recently been introduced by Gozlan et al. (J Funct Anal 273(11):3327–3405, 2017). We provide general existence and duality results for these problems on arbitrary Polish spaces, as well as a necessary and sufficient optimality criterion in the spirit of cyclical monotonicity. As an application we extend the Brenier–Strassen Theorem of Gozlan and Juillet (On a mixture of brenier and strassen theorems. arXiv:1808.02681, 2018) to general probability measures on $$\mathbb {R}^d$$ under minimal assumptions. A driving idea behind our proofs is to consider the set of transport plans with a new (‘adapted’) topology which seems better suited for the weak transport problem and allows to carry out arguments which are close to the proofs in the classical setup.

/pdf/existence-duality-and-cyclical-monotonicity-for-weak-4ow9axrsug.pdf

Existence, duality, and cyclical monotonicity for weak transport costs

Motivated by the approximation of Martingale Optimal Transport problems, we
study sampling methods preserving the convex order for two probability measures
$\mu$ and $\nu$ on $\mathbb{R}^d$, with $\nu$ dominating $\mu$. When
$(X_i)_{1\le i\le I}$ (resp. $(Y_j)_{1\le j\le J}$) are i.i.d. according $\mu$
(resp. $\nu$), the empirical measures $\mu_I$ and $\nu_J$ are not in the convex
order. We investigate modifications of $\mu_I$ (resp. $\nu_J$) smaller than
$\nu_J$ (resp. greater than $\mu_I$) in the convex order and weakly converging
to $\mu$ (resp. $\nu$) as $I,J\to\infty$. In dimension 1, according to Kertz
and R\"osler (1992), the set of probability measures with a finite first order
moment is a lattice for the increasing and the decreasing convex orders. From
this result, we can define $\mu\vee\nu$ (resp. $\mu\wedge\nu$) that is greater
than $\mu$ (resp. smaller than $\nu$) in the convex order. We give efficient
algorithms permitting to compute $\mu\vee\nu$ and $\mu\wedge\nu$ when $\mu$ and
$\nu$ are convex combinations of Dirac masses. In general dimension, when $\mu$
and $\nu$ have finite moments of order $\rho\ge 1$, we define the projection
$\mu\curlywedge_\rho \nu$ (resp. $\mu\curlyvee_\rho\nu$) of $\mu$ (resp. $\nu$)
on the set of probability measures dominated by $\nu$ (resp. larger than $\mu$)
in the convex order for the Wasserstein distance with index $\rho$. When
$\rho=2$, $\mu_I\curlywedge_2 \nu_J$ can be computed efficiently by solving a
quadratic optimization problem with linear constraints. It turns out that, in
dimension 1, the projections do not depend on $\rho$ and their quantile
functions are explicit, which leads to efficient algorithms for convex
combinations of Dirac masses. Last, we illustrate by numerical experiments the
resulting sampling methods that preserve the convex order and their application
to approximate Martingale Optimal Transport problems.

Sampling of probability measures in the convex order by Wasserstein projection

Weak optimal transport has been recently introduced by Gozlan et al. The original motivation stems from the theory of geometric inequalities; further applications concern numerics of martingale optimal transport and stability in mathematical finance. 
In this note we provide a complete geometric characterization of the 'weak' version of the classical monotone rearrangement between measures on the real line, complementing earlier results of Alfonsi, Corbetta, and Jourdain.

Weak monotone rearrangement on the line

. The adapted weak topology is an extension of the weak topology for stochastic processes designed to adequately capture properties of underly-ing ﬁltrations. With the recent work of Bart–Beiglböck–P. [7] as starting point, the purpose of this note is to recover with topological arguments the intrigu-ing result by Backhoﬀ–Bartl–Beiglböck–Eder [5] that all adapted topologies in discrete time coincide. We also derive new characterizations of this topology including descriptions of its trace on the sets of Markov processes and processes equipped with their natural ﬁltration. To emphasize the generality of the argument, we also describe the classical weak topology for measures on R d by a weak Wasserstein metric based on the theory of weak optimal transport initiated by Gozlan–Roberto–Samson–Tetali [12].

A note on the adapted weak topology in discrete time

Motivated by applications to geometric inequalities, Gozlan, Roberto, Samson, and Tetali (J. Funct. Anal. 273 (2017) 3327–3405) introduced a transport problem for ‘weak’ cost functionals. Basic results of optimal transport theory can be extended to this setup in remarkable generality. In this article, we collect several problems from different areas that can be recast in the framework of weak transport theory, namely: the Schrodinger problem, the Brenier–Strassen theorem, optimal mechanism design, linear transfers, semimartingale transport. Our viewpoint yields a unified approach and often allows to strengthen the original results.

Applications of weak transport theory

The numerical solution of stochastic partial differential equations and numerical Bayesian estimation is computationally demanding. If the coefficients in a stochastic partial differential equation...

Cubature Formulas for Multisymmetric Functions and Applications to Stochastic Partial Differential Equations

The stochastic nonlinear Poisson--Boltzmann equation describes the electrostatic potential in a random environment in the presence of free charges and has applications in many fields. We show the e...

Gudmund Pammer

Papers

Weak monotone rearrangement on the line

A note on the adapted weak topology in discrete time

Applications of weak transport theory

Cubature Formulas for Multisymmetric Functions and Applications to Stochastic Partial Differential Equations

Existence, Uniqueness, and a Comparison of Nonintrusive Methods for the Stochastic Nonlinear Poisson--Boltzmann Equation