Dimension Reduction in Martingale Optimal Transport: Geometry and Robust Option Pricing
Joshua Zoen-Git Hiew, Tongseok Lim, Brendan Pass, Marcelo Cruz de Souza
TL;DR
The paper analyzes robust option pricing via Vectorial Martingale Optimal Transport (VMOT) in multi-period, two-asset settings, showing that for $d=2$ and supermodular costs, the first-period joint distribution of the assets is monotone, enabling a dimension reduction that keeps the martingale structure intact. It develops a dual formulation and proves dual attainment, then uses convex analysis to establish the main monotonicity result, plus a perturbation argument to handle strict supermodularity; counterexamples demonstrate sharpness in higher dimensions ($d\ge 3$). The authors introduce a reduced-dimension Sinkhorn algorithm that exploits the monotone first-period structure, achieving up to roughly $99\%$ reductions in compute time while maintaining or improving accuracy, and they present several numeric experiments across two- and three-period problems. The work clarifies when dimension reduction is possible in VMOT, contrasts VMOT with MMOT, and provides practical guidance for robust pricing with scalable computation in the two-asset regime.
Abstract
This paper addresses the problem of robust option pricing within the framework of Vectorial Martingale Optimal Transport (VMOT). We investigate the geometry of VMOT solutions for $N$-period market models and demonstrate that, when the number of underlying assets is $d=2$ and the payoff is sub- or supermodular, the extremal model reduces to a single-factor structure in the first period. This structural result allows for a significant dimension reduction, transforming the problem into a more tractable format. We prove that this reduction is specific to the two-asset case and provide counterexamples showing it generally fails for $d \geq 3$. Finally, we exploit this monotonicity to develop a reduced-dimension Sinkhorn algorithm. Numerical experiments demonstrate that this structure-preserving approach reduces computational time by approximately 99\% compared to standard methods while improving accuracy.
