← Blog/blog/circular-convolution-calendar-arbitrage

A wrapped convolution can manufacture option arbitrage

An implied-volatility surface is a map of option prices across strike and maturity. A visually smooth surface can still quote an impossible market: a butterfly can cost less than zero, or a longer-lived option can cost less than its shorter-lived twin. Luo and colleagues avoid policing those failures with a loss penalty. Their decoder emits a probability distribution whose prices are valid by construction.

01

Price a distribution, not a volatility pixel

An option strike is the price at which its holder may buy the asset. Forward normalization divides away carry so the future asset ratio M has expected value one. The normalized call price is then E[(M−z)⁺], where z is strike divided by the forward. Because this payoff is convex, its expectation cannot fall as a martingale gains another independent increment.

AddMart's learned Transformer reads nine sparse quotes and predicts non-negative mixtures for 30 log-price increments. The decoder convolves those increments, exponentiates and normalizes them to mean one, prices calls, then inverts Black's formula to implied volatility. The learning network is flexible; the decoder supplies the financial structure.

02

The tiny martingale we can audit by hand

Take one log increment on the grid −1, 0, 1 with probabilities 0.25, 0.50, 0.25. After two independent increments, exact convolution places mass on five sums from −2 through 2. The TypeScript core below runs the page; Python and C++ are line-for-line translations of the same index rule.

def convolve(a, b, circular=False):
    size = len(a) if circular else len(a) + len(b) - 1
    out = [0.0] * size
    for i, ai in enumerate(a):
        for j, bj in enumerate(b):
            index = (i + j) % size if circular else i + j
            out[index] += ai * bj
    total = sum(out)
    return [mass / total for mass in out]
Illustrative deterministic toy. Exact convolution expands to all five attainable log sums. Circular convolution keeps only three bins, so sums outside the grid wrap onto unrelated states.

The wrap is not harmless compression. It changes which states carry tail probability, then each maturity is renormalized to have mean one. That preserves a valid marginal distribution at each date, but it no longer preserves the coupling that made the dates a martingale.

maturity 1maturity 2: exactmaturity 2: circular
Illustrative core-computed call curves for strikes 0.5–1.8. Exact maturity-two prices stay above maturity one. In the right wing, the wrapped decoder crosses below it.
exact calendar residualcircular calendar residual
Calendar residual means later price minus earlier price. Negative is an arbitrage violation. At strike 1.5, circular convolution creates a residual of -1.74 price points; exact convolution remains non-negative.
03

Why exact convolution preserves the theorem

If L₁,…,Lⱼ are independent log increments, AddMart sets Mⱼ equal to exp(ΣL)/E[exp(ΣL)]. The next multiplicative factor is non-negative, independent of the past, and has mean one. Therefore E[Mⱼ₊₁ | past] = Mⱼ. Jensen's inequality says every convex payoff—including a call— has non-decreasing expected value across that step.

full-support atoms
Core-computed support growth for the paper's 2,048-point increment grid. Exact support needs G(X−1)+1 atoms: 61,411 at 30 maturity buckets. Circular convolution is tempting precisely because it avoids this growth.
04

The paper's ablation separates fit from validity

The authors compare exact full-support convolution with two shortcuts: cropping the expanding support and using circular FFT convolution. Reliable-region IV-MAE barely distinguishes them. The arbitrage audit does.

DecoderSPY IV-MAENVDA IV-MAESPY calendar viol.NVDA calendar viol.
linear exact0.561.840.000%0.000%
linear crop0.561.823.855%0.023%
circular0.551.763.652%0.070%
Paper-reported means over five seeds. MAE is in volatility points; violation columns are percentages.
Paper-reported SPY reliable-region IV-MAE. The circular shortcut is microscopically better on this accuracy metric; lower is better.
Paper-reported SPY calendar-violation rate. The accuracy winner is structurally invalid on 3.652% of audited comparisons; lower is better.
05

What the guarantee does—and does not—cover

The formal result is global over the modeled strike domain and from time zero through the last maturity bucket. The market data are listed equity and ETF options, but the construction treats liquid out-of-the-money quotes as European-equivalent Black–76 IVs using parity-implied forwards. That limits early-exercise and dividend effects; it does not prove the original American contracts themselves obey a European pricing model.

The main experiment also excludes same-day-expiry options, removes maturities below one day, and defines a reliable accuracy region that excludes contracts with fewer than 14 days remaining. The nine reference quotes use a structured three-maturity by three-delta layout. Extra references help but not monotonically, and zero-shot transfer deteriorates when assets have different volatility levels or smile shapes.

06

What to probe next

First, audit violations by strike and maturity rather than averaging them: the toy shows why aliases can hide in a wing. Second, report the runtime and memory trade-off of exact support at the same batch size. Third, stress the quote layout, same-day expiries, and dividend-heavy names. Finally, compare tail functionals against later realized outcomes—not only a smoothed option-implied reference.

References

  1. Luo, H., Lin, S., Luo, X., Minh, A. B., Zhang, Y., and Liu, C. (2026). Arbitrage-Free Implied Volatility Surface Reconstruction via Additive Martingales. arXiv:2610.25671
  2. Breeden, D. T. and Litzenberger, R. H. (1978). Prices of State-Contingent Claims Implicit in Option Prices. Journal of Business, 51(4)
  3. Kellerer, H. G. (1972). Markov-Komposition und eine Anwendung auf Martingale. Mathematische Annalen, 198