SANOS Smooth strictly Arbitrage-free Non-parametric Option Surfaces

Anastasis Kratsios; Blanka Horvath; Hans Buehler; Raeid Saqur; Yannick Limmer

arxiv: 2601.11209 · v4 · pith:5SW6IXA5new · submitted 2026-01-16 · 💱 q-fin.CP · q-fin.MF

SANOS Smooth strictly Arbitrage-free Non-parametric Option Surfaces

Hans Buehler , Blanka Horvath , Anastasis Kratsios , Yannick Limmer , Raeid Saqur This is my paper

Pith reviewed 2026-05-25 07:14 UTC · model grok-4.3

classification 💱 q-fin.CP q-fin.MF

keywords option price surfacesarbitrage-freenon-parametriclinear programmingdiscrete local volatilitysmooth interpolationS&P 500 options

0 comments

The pith

A smooth generalization of linear interpolation produces strictly arbitrage-free option price surfaces with only positivity constraints.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a non-parametric method for option price surfaces that are smooth and strictly arbitrage-free across time and strike dimensions. This approach extends linear interpolation smoothly while preserving its simplicity and transparency. Calibration to market quotes is cast as a linear program that incorporates bid-ask spreads directly through penalties or inequalities. An equivalent form is derived using strictly positive discrete local volatility variables, which the authors present as the first construction of such surfaces needing only trivial positivity constraints. A sympathetic reader would care because the method delivers materially lower computational cost than most implied-volatility surface routines while remaining numerically efficient and flexible.

Core claim

The paper establishes that a smooth generalization of linear interpolation, calibrated by linear program to observed quotes, yields option price surfaces that are smooth and strictly arbitrage-free across time and strike. This surface admits an equivalent parameterization in terms of strictly positive discrete local volatility variables, allowing the arbitrage-free property to hold under only the constraint of positivity. The approach is illustrated using S&P 500 index options.

What carries the argument

The smooth generalization of linear interpolation calibrated by linear program, together with its equivalent parameterization via strictly positive discrete local volatility variables.

If this is right

Calibration becomes a linear program that directly handles bid-ask spreads.
Surfaces remain strictly arbitrage-free by construction in both time and strike.
Only positivity is required on the discrete local volatility parameters.
The method is numerically efficient and was demonstrated on S&P 500 index options.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The positivity-only parameterization may allow direct embedding into existing linear or convex optimization pipelines used in portfolio risk systems.
The same linear-program structure could be applied to construct surfaces for other underlyings where quote data arrive in real time.
Because the method starts from linear interpolation, existing codebases for that baseline could be upgraded with minimal changes to enforce the no-arbitrage property.

Load-bearing premise

The smooth generalization of linear interpolation remains strictly arbitrage-free across time and strike when calibrated by linear program to observed quotes.

What would settle it

Fitting the surface to real S&P 500 options quotes and then checking whether any calendar spread or butterfly spread violates no-arbitrage conditions would directly test the central claim.

read the original abstract

We present a simple, numerically efficient but highly flexible non-parametric method to construct representations of option price surfaces which are both smooth and strictly arbitrage-free across time and strike. The method can be viewed as a smooth generalization of the widely-known linear interpolation scheme, and retains the simplicity and transparency of that baseline. Calibration of the model to observed market quotes is formulated as a linear program, allowing bid-ask spreads to be incorporated directly via linear penalties or inequalities, and delivering materially lower computational cost than most of the currently available implied-volatility surface fitting routines. As a further contribution, we derive an equivalent parameterization of the proposed surface in terms of strictly positive "discrete local volatility" variables. This yields, to our knowledge, the first construction of smooth, strictly arbitrage-free option price surfaces while requiring only trivial parameter constraints (positivity). We illustrate the approach using S&P 500 index options

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The discrete local vol reparameterization is the actual advance, turning strict no-arbitrage into a positivity-only LP while keeping the surface smooth.

read the letter

The main takeaway is the reparameterization of the option surface through strictly positive discrete local volatility variables. This reduces the usual no-arbitrage inequalities to simple positivity constraints and lets them claim the first smooth, strictly arbitrage-free construction under only those trivial bounds. Calibration becomes a linear program that directly incorporates bid-ask spreads via penalties or inequalities, which is cheaper than most implied-vol fitting routines and keeps the method transparent, like a smoothed version of linear interpolation. The S&P 500 illustration shows it runs on real data without exotic machinery. That combination of smoothness, strict arbitrage-free property, and LP tractability is new relative to the cited literature. The approach is practical for anyone who needs reliable surfaces inside pricing or risk engines. The soft spot is that the abstract states the equivalence and the LP formulation but does not supply the explicit derivation steps or error bounds that would confirm the smoothed surface stays strictly arbitrage-free for every strike and maturity. Without those checks visible, it is hard to judge whether the positivity variables fully carry the no-arbitrage load after the smoothing step or whether hidden conditions appear in the full construction. The numerical example would also need to demonstrate that the LP solutions remain stable across different market regimes. This paper is for quant teams that build or maintain option surfaces day-to-day and want something lighter than current non-parametric fits. A reader who values implementable methods over heavy theory will find the LP and the new variables useful even if the proofs require tightening. It deserves a serious referee because the parameterization is a genuine simplification and the computational claim is testable; the technical details are the part that needs external verification rather than a reason to desk-reject.

Referee Report

2 major / 2 minor

Summary. The paper presents SANOS, a non-parametric method to construct smooth, strictly arbitrage-free option price surfaces across strikes and maturities. It is positioned as a smooth generalization of linear interpolation, with calibration to market quotes formulated as a linear program that directly incorporates bid-ask spreads, and an equivalent reparameterization in terms of strictly positive discrete local volatility variables that is claimed to enforce the no-arbitrage property by construction with only trivial positivity constraints. The approach is illustrated on S&P 500 index options data.

Significance. If the central construction holds, the method would supply a computationally lightweight, transparent alternative to existing implied-volatility surface routines while guaranteeing strict no-arbitrage without complex constraints; the LP calibration and positivity-only parameterization would be practically useful strengths for pricing and risk applications.

major comments (2)

[Abstract / parameterization section] Abstract and the section introducing the discrete-local-volatility parameterization: the central claim that positivity of the discrete local volatility variables alone ensures the resulting smooth surface is strictly arbitrage-free for all strikes and maturities rests on an unshown equivalence and verification; no explicit derivation steps, error bounds, or exhaustive check against the full set of static and dynamic arbitrage conditions are supplied.
[LP calibration section] The linear-program calibration section: while the formulation is stated, there is no analysis confirming that the optimal solution under positivity constraints preserves the strict no-arbitrage property of the underlying interpolation scheme when the surface is evaluated at arbitrary (non-grid) strikes and maturities.

minor comments (2)

Notation for the smooth interpolation weights and the discrete local volatility variables should be introduced with a single consistent table or equation block to improve readability.
The numerical example on S&P 500 options would benefit from an explicit statement of the number of quotes, the chosen grid, and a quantitative measure (e.g., maximum violation of butterfly or calendar-spread conditions) to demonstrate the claimed strict no-arbitrage.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. Below we respond point by point to the major comments and indicate the revisions we will make.

read point-by-point responses

Referee: [Abstract / parameterization section] Abstract and the section introducing the discrete-local-volatility parameterization: the central claim that positivity of the discrete local volatility variables alone ensures the resulting smooth surface is strictly arbitrage-free for all strikes and maturities rests on an unshown equivalence and verification; no explicit derivation steps, error bounds, or exhaustive check against the full set of static and dynamic arbitrage conditions are supplied.

Authors: Section 3 derives the equivalence by showing that the smooth interpolant can be re-expressed exactly in terms of the discrete local volatility variables and that positivity of those variables is necessary and sufficient for the surface to satisfy the static no-arbitrage inequalities at every strike and maturity. We agree that the derivation steps can be expanded and that an explicit verification table against the standard static arbitrage conditions would improve clarity. We will add these elements in the revised manuscript. Dynamic arbitrage conditions lie outside the scope of a static surface construction. revision: yes
Referee: [LP calibration section] The linear-program calibration section: while the formulation is stated, there is no analysis confirming that the optimal solution under positivity constraints preserves the strict no-arbitrage property of the underlying interpolation scheme when the surface is evaluated at arbitrary (non-grid) strikes and maturities.

Authors: Because every feasible point of the LP corresponds to a strictly positive discrete-local-volatility vector, and the parameterization maps any such vector to a surface that is arbitrage-free at all (not merely grid) strikes and maturities, the optimal solution automatically inherits the property. We will insert a short paragraph in the LP section that makes this preservation explicit and notes that the continuous nature of the parameterization extends the guarantee beyond the calibration grid. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via explicit parameterization

full rationale

The paper constructs a surface via smooth generalization of linear interpolation, then reparameterizes it explicitly in terms of strictly positive discrete local volatility variables. The no-arbitrage property is stated to follow directly from this parameterization under the positivity constraint alone, with calibration performed as an LP on external market quotes. No step reduces a claimed prediction to a fitted input by construction, no self-citation is invoked as load-bearing justification for uniqueness or ansatz, and the central equivalence (positivity implies strict arbitrage-freeness) is presented as a definitional property of the chosen coordinates rather than an independent derivation that loops back to itself. The method is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on standard no-arbitrage axioms from option pricing theory plus the new discrete-local-volatility parameterization; no free parameters are introduced beyond the market data itself.

axioms (1)

domain assumption Option prices must satisfy no-arbitrage conditions (positive butterfly spreads, positive calendar spreads) to be strictly arbitrage-free.
Invoked throughout the abstract to guarantee the strict arbitrage-free property of the constructed surface.

invented entities (1)

discrete local volatility variables no independent evidence
purpose: Provide an equivalent parameterization of the option surface that requires only positivity constraints.
Introduced in the abstract as the key technical device enabling trivial constraints while preserving smoothness and no-arbitrage.

pith-pipeline@v0.9.0 · 5697 in / 1352 out tokens · 63720 ms · 2026-05-25T07:14:27.223887+00:00 · methodology

discussion (0)

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Structure-Preserving Reconstruction of Convex Lipschitz Functionals on Hilbert Spaces from Finite Samples
math.FA 2026-05 unverdicted novelty 7.0

Any convex L-Lipschitz functional on a compact convex subset of a separable Hilbert space can be uniformly approximated to arbitrary accuracy by an explicit convex L-Lipschitz reconstruction from finitely many linear ...
Robust Volatility Index Calculation with OTM Option-implied Probability
q-fin.MF 2026-05 unverdicted novelty 6.0

A construction of continuous European option prices from OTM bid-ask data with fewer parameters enables robust model-free volatility index calculation under no-arbitrage.