Pith. sign in

REVIEW 4 major objections 4 minor 14 references

Empirical Models of the Time Evolution of SPX Option Prices

T0 review · 4 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read A tiny neural net beats Black-Scholes on SPX puts.

desk verdict Solid empirical benchmark undermined by an overclaimed arbitrage-free result that its own tests refute. read the letter →

arxiv 2506.17511 v1 pith:LL5VZTZP submitted 2025-06-20 q-fin.PR q-fin.CP

classification q-fin.PRq-fin.CP
keywords neuralnetworkoptionpricingS&P500putoptionsBlack-ScholesbenchmarkGARCH(11)volatilityout-of-the-moneyno-arbitrageconditionsrandomforestSHAPexplainability
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper sets out to show that a deliberately small neural network—two hidden layers with four neurons each—can price S&P 500 put options better than the closed-form Black-Scholes-Merton formula and can do so without needing no-arbitrage conditions imposed by hand. The models are trained on roughly 7.4 million daily put observations from 1996 to 2022, with inputs such as moneyness, time to maturity, the index level, dividend yield, the risk-free rate, and a GARCH volatility forecast. The neural network is evaluated against a random forest and a linear regression under expanding and rolling training windows, and it is the only one that stays competitive when the Black-Scholes price is not supplied as an input. If the claim holds, the payoff is a simulation-ready pricing function that generates internally consistent option-price paths along simulated SPX paths for stress testing and tail-risk hedging.

What carries the argument

The load-bearing object is the feedforward pricing function $f_\theta(X)$: a network with two hidden layers of four ReLU neurons that maps option characteristics (strike, moneyness, time to maturity) and market state variables (SPX level, dividend yield, risk-free rate, GARCH(1,1) forecast volatility) to a put price. The GARCH(1,1) volatility estimate, re-estimated daily on a rolling 252-day window, supplies time-varying volatility information without relying on an implied-volatility surface. Training minimizes the Huber loss, which keeps the network from being dominated by a few high-priced options and preserves accuracy on the cheap deep-OTM contracts; the design that includes or excludes the Black-Scholes price as an input isolates how much of each model's performance is borrowed from the theoretical benchmark. SHAP values and a principal-component comparison are then used to argue that the network relies on strike, maturity, and index level in the way option theory would predict.

What would settle it

On a fixed test date, compute the network's prices on a dense strike-maturity grid and search for a static portfolio of puts, such as a butterfly or a calendar spread, with nonnegative payoffs and negative cost; if such a portfolio exists on more than a negligible fraction of dates, the claim that the network delivers arbitrage-free prices without constraints is false.

Watch

Extended reading notes

Core claim

On its own terms, the paper claims that an off-the-shelf two-layer feedforward network with four ReLU neurons per layer, trained with Adam, weight decay, early stopping, and the Huber loss, approximates the European put pricing function for SPX more accurately than Black-Scholes-Merton out of sample across 1996–2022. The advantage is largest for deep out-of-the-money puts and during heightened-volatility periods, and it persists whether or not the Black-Scholes price is included as a feature. The paper also claims that the trained network yields prices consistent with no-arbitrage conditions in the large majority of perturbed cases—about 93.5 percent monotone in strike, 95.1 percent convex in strike, and 82.9 percent monotone in time to maturity—without these conditions being enforced during training.

Load-bearing premise

The load-bearing premise is that European put prices must increase with time to maturity, which the paper treats as a no-arbitrage condition; with nonzero interest rates and dividend yields this monotonicity does not generally hold, so the violation statistics do not by themselves establish economic consistency.

Editorial extensions

If this is right

  • A single small network can replace a daily recalibrated implied-volatility surface for pricing the SPX put cross-section and remains accurate out of sample.
  • Because options are priced directly from state variables, the model can produce simulated option prices along GARCH-simulated SPX paths without interpolating or extrapolating a volatility surface.
  • The neural network does not need the Black-Scholes price as an input, so it can be used in regimes where the theoretical model is misspecified.
  • The random forest degrades sharply without the Black-Scholes feature and linear regression generally stays below the Black-Scholes baseline, so the network is the only tested model that is both nonlinear and self-sufficient.
  • Residual no-arbitrage violations appear when inputs are perturbed far from the original point, so deployment would still benefit from regularization or post-hoc correction.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's monotonicity-in-maturity test is not a valid arbitrage bound for European puts when interest rates and dividend yields are nonzero, so the 17.1 percent violation rate is better interpreted as a property of the test than as evidence that the model admits arbitrage.
  • A proper economic-consistency check would test calendar-spread and butterfly arbitrage on a dense strike-maturity grid; the paper's local perturbation procedure only probes neighborhoods of observed points.
  • A natural extension is to train the same architecture on SPX calls or on a multi-asset index to see whether the accuracy and no-arbitrage findings generalize beyond puts.
  • The small network size suggests an implicit smoothness prior; pricing on a dense grid and inspecting the implied volatility smile would reveal whether the learned surface is itself arbitrage-free.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper develops statistical models for pricing S&P 500 index put options, using a GARCH(1,1) volatility estimate and a 1996–2022 OptionMetrics dataset. It compares a small feedforward neural network, a random forest, and linear regression under expanding and rolling training windows, with and without the Black-Scholes price as an input feature, reporting MAPE over test periods and moneyness segments. It also presents SHAP-based interpretability and a Section 7 evaluation of no-arbitrage conditions. The headline claim is that a two-hidden-layer, four-neuron-per-layer neural network trained with minimal tuning performs well against Black-Scholes-Merton and 'delivers arbitrage-free option prices without requiring these conditions to be imposed.'

Significance. If the central claims held, the paper would offer a practically useful result: a small, explainable network that both improves on Black-Scholes pricing and automatically satisfies no-arbitrage constraints. The empirical scope is a strength: roughly seven million filtered options over 27 years, two training-window schemes, three model classes, and extensive appendix tables. The SHAP and PCA analyses are also a genuine effort to connect the fitted models to economic intuition. However, the headline arbitrage-free claim is contradicted by the paper's own Section 7 results, and the performance comparisons are reported without uncertainty quantification. The paper's value as an empirical benchmark is diminished by these issues, and the central abstract claim is not supported as stated.

major comments (4)
  1. [Abstract; Section 7; Section 9] The abstract states that the neural network 'delivers arbitrage-free option prices without requiring these conditions to be imposed,' but Section 7's own tests report that, on 123,587 options, only 93.51% of predictions respect strike monotonicity, 95.09% respect strike convexity, and 82.92% respect time-to-maturity monotonicity. Strike monotonicity and convexity are model-free no-arbitrage bounds for European put prices, so the reported 6.49% and 4.91% violation rates directly contradict an unconditional arbitrage-free claim. Even if the TTM condition is set aside, the two valid conditions are still violated. The conclusion's phrase 'respect key no-arbitrage requirements' is therefore internally inconsistent with the evidence in Section 7, and the claim must either be removed or explicitly qualified as approximate/tolerance-based.
  2. [Section 5.3; Tables 6, 12, 18] The claim that the neural network 'consistently outperforms both linear regression and random forest benchmarks' is not supported by the appendix tables for in-the-money options. For example, in Table 6 the expanding-window ITM test MAPE for NN+ in the first window is 34.5% versus 5.0% for BS and 4.2% for RF+ in Table 18, and several later ITM windows show NN+ above both the BS baseline and RF+/LR+. Section 5.3's statement that the ITM analysis 'largely mirrors' the OTM findings is contradicted by these numbers. The comparisons are reported as point estimates with no confidence intervals or paired significance tests, so even the OTM outperformance claims are not statistically established.
  3. [Section 1; Section 9] The stated objective is a framework that 'enables the simulation of joint paths for asset prices and corresponding option prices,' and the conclusion describes the paper as providing a simulation framework. However, no simulation experiment appears anywhere in the manuscript: there are no generated paths, no distributional or tail-risk diagnostics, no stress-test application, and no code. The paper evaluates static out-of-sample pricing accuracy only, so the simulation capability is asserted rather than demonstrated.
  4. [Section 7] The time-to-maturity monotonicity test is not a valid no-arbitrage test for European put options when interest rates and dividend yields are nonzero, since put prices need not be increasing in maturity under those conditions. Counting the 17.1% TTM violation rate as evidence against economic consistency is therefore misleading. The Section 7 test design also perturbs only one input at a time and uses a $0.05 tolerance, so it is neither a sufficient nor a fully reported test of static arbitrage; a complete check would need simultaneous constraints across strikes and maturities, such as calendar-spread and butterfly conditions.
minor comments (4)
  1. [Figure 10 caption] The caption for Figure 10 describes the panels as 'NN trained with Black-Scholes information' and 'NN trained without,' but the figure panels are labeled RF+ and RF-; the caption should refer to the random forest model.
  2. [Equation (10)] The SHAP formula is written with an awkward summation notation and a garbled subscript; it should be rewritten as a sum over all subsets S not containing i, with the standard Shapley weights.
  3. [Section 8.1] The sentence 'MAPE increases as options move deeper OTM, MAPE increases significantly across all models' contains a duplicated subject and should be edited for clarity.
  4. [General] The manuscript does not provide a data availability statement or code repository; given the volume of numerical results, making the preprocessing and training pipelines available would materially improve reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the paper's pricing models are fitted to observed market prices and evaluated out-of-sample; the no-arbitrage claim is an empirical test, not a derivation from its own inputs.

full rationale

The models are trained on market option prices and tested on held-out periods, so the predicted prices are not equal to any fitted input by construction. The GARCH parameters are estimated from SPX index returns, which is an external input, and the Black-Scholes price is used only as a feature or baseline; including it does not force the neural network output to equal the Black-Scholes value because the target remains the observed market price. No uniqueness theorem, ansatz-smuggling self-citation, or renamed empirical pattern carries the central claim. The abstract's 'arbitrage-free' statement is contradicted in part by the paper's own Section 7 violation counts, but that is a correctness or internal-consistency concern, not circularity: the claim is testable and, if the numbers are accurate, false as stated, rather than true by definition. The paper also acknowledges in Section 8.2 that the neural network 'can still produce such violations under certain training conditions,' which further confirms the claim is empirical rather than definitionally forced. Self-citations such as Brini and Lenz 2024 appear only in related-work context and are not load-bearing. Overall, the derivation chain is self-contained with respect to the data and benchmark comparisons, so the appropriate circularity score is 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new theoretical entities. It relies on standard GARCH estimation, off-the-shelf ML models, and a Black-Scholes benchmark. The free parameters are the GARCH coefficients and manually chosen hyperparameters; the key assumptions are the GARCH data-generating process and the validity of the no-arbitrage tests used in Section 7.

free parameters (3)
  • GARCH(1,1) parameters (mu, a0, a1, b1) = Re-estimated daily on rolling 252-day window
    Used to compute the GARCH volatility input and the Black-Scholes benchmark price. If misspecified, the BS baseline is handicapped and the NN advantage may be inflated.
  • NN hyperparameters (2 hidden layers, 4 neurons per layer, Huber delta=1, learning rate 1e-4, weight decay 1e-3, batch… = Chosen by hand, not tuned
    The claim of strong performance with minimal tuning depends on these fixed choices; other architectures could perform differently.
  • RF hyperparameters (100 trees, max depth 10) = Chosen by hand
    Random forest performance and stability depend on these settings; they are not derived from first principles.
assumptions (4)
  • domain assumption SPX log returns follow a GARCH(1,1) process with standard normal innovations.
    Invoked in Section 3, Equations (1) and (2), to produce volatility forecasts used as model inputs and in the Black-Scholes benchmark.
  • domain assumption No-arbitrage bounds are characterized by monotonicity in strike and maturity and convexity in strike, as tested in Section 7.
    The paper treats TTM monotonicity as a universal bound for European puts, which is not generally valid with dividends.
  • domain assumption OptionMetrics quotes from 1996-2022 with the stated filters are representative of SPX put market prices.
    Section 3 relies on this for all training and evaluation.
  • domain assumption MAPE is an appropriate error metric for comparing option pricing performance across orders of magnitude.
    Section 5 motivates MAPE, but it is known to be dominated by small-price options and is not a risk-neutral or economic loss measure.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Empirical Models of the Time Evolution of SPX Option Prices." pith.science (2026). https://pith.science/paper/LL5VZTZP

@misc{pith2026250617511,
  author       = {Pith},
  title        = {Pith review of: Empirical Models of the Time Evolution of SPX Option Prices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LL5VZTZP}},
  note         = {Machine review of arXiv:2506.17511}
}
read the original abstract

The key objective of this paper is to develop an empirical model for pricing SPX options that can be simulated over future paths of the SPX. To accomplish this, we formulate and rigorously evaluate several statistical models, including neural network, random forest, and linear regression. These models use the observed characteristics of the options as inputs -- their price, moneyness and time-to-maturity, as well as a small set of external inputs, such as the SPX and its past history, dividend yield, and the risk-free rate. Model evaluation is performed on historical options data, spanning 30 years of daily observations. Significant effort is given to understanding the data and ensuring explainability for the neural network. A neural network model with two hidden layers and four neurons per layer, trained with minimal hyperparameter tuning, performs well against the theoretical Black-Scholes-Merton model for European options, as well as two other empirical models based on the random forest and the linear regression. It delivers arbitrage-free option prices without requiring these conditions to be imposed.

Figures

Figures reproduced from arXiv: 2506.17511 by the authors.

Figure 1
Figure 1. Empirical distributions of options data observations across different dimensions. The left panel shows the [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Evolution over time of the MAPE values for OTM options trained using an expanding window schema [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Evolution over time of the MAPE values for OTM options trained using an expanding window schema [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (27 more)
Figure 4
Figure 4. Figure 4: Evolution over time of the MAPE values for OTM options trained using an expanding window schema [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: Evolution over time of the MAPE values for OTM options trained using a rolling window schema with [PITH_FULL_IMAGE:figures/full_fig_p014_5.png]
Figure 6
Figure 6. Figure 6: Evolution over time of the MAPE values for OTM options trained using a rolling window schema with [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]
Figure 7
Figure 7. Figure 7: Evolution over time of the MAPE values for OTM options trained using a rolling window schema with [PITH_FULL_IMAGE:figures/full_fig_p015_7.png]
Figure 8
Figure 8. Figure 8: Boxplots of the mean absolute percentage error (MAPE) for OTM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: Boxplots of the mean absolute percentage error (MAPE) for OTM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p017_9.png]
Figure 10
Figure 10. Figure 10: Boxplots of the mean absolute percentage error (MAPE) for OTM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p018_10.png]
Figure 11
Figure 11. Figure 11: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using [PITH_FULL_IMAGE:figures/full_fig_p020_11.png]
Figure 12
Figure 12. Figure 12: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using [PITH_FULL_IMAGE:figures/full_fig_p020_12.png]
Figure 13
Figure 13. Figure 13: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using [PITH_FULL_IMAGE:figures/full_fig_p021_13.png]
Figure 14
Figure 14. Figure 14: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using a [PITH_FULL_IMAGE:figures/full_fig_p022_14.png]
Figure 15
Figure 15. Figure 15: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using a [PITH_FULL_IMAGE:figures/full_fig_p023_15.png]
Figure 16
Figure 16. Figure 16: Evolution over time of the MAPE values for OTM options grouped by moneyness ranges trained using a [PITH_FULL_IMAGE:figures/full_fig_p024_16.png]
Figure 17
Figure 17. Figure 17: Boxplots of the mean absolute percentage error (MAPE) for OTM options grouped by moneyness ranges [PITH_FULL_IMAGE:figures/full_fig_p025_17.png]
Figure 18
Figure 18. Figure 18: Boxplots of the mean absolute percentage error (MAPE) for OTM options grouped by moneyness ranges [PITH_FULL_IMAGE:figures/full_fig_p026_18.png]
Figure 19
Figure 19. Figure 19: Boxplots of the mean absolute percentage error (MAPE) for OTM options grouped by moneyness ranges [PITH_FULL_IMAGE:figures/full_fig_p027_19.png]
Figure 20
Figure 20. Figure 20: Evolution over time of the MAPE values for ITM options trained using an expanding window schema. [PITH_FULL_IMAGE:figures/full_fig_p029_20.png]
Figure 21
Figure 21. Figure 21: Evolution over time of the MAPE values for ITM options grouped by moneyness ranges trained using a [PITH_FULL_IMAGE:figures/full_fig_p030_21.png]
Figure 22
Figure 22. Figure 22: Boxplots of the mean absolute percentage error (MAPE) for ITM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p031_22.png]
Figure 23
Figure 23. Figure 23: Boxplots of the mean absolute percentage error (MAPE) for ITM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p032_23.png]
Figure 24
Figure 24. Figure 24: Boxplots of the mean absolute percentage error (MAPE) for ITM options under two data-splitting [PITH_FULL_IMAGE:figures/full_fig_p033_24.png]
Figure 25
Figure 25. Figure 25: SHAP Analysis for OTM Models. Each row in the grid represents a model tested on a specific [PITH_FULL_IMAGE:figures/full_fig_p035_25.png]
Figure 26
Figure 26. Figure 26: SHAP Analysis for ITM Models. Each row in the grid represents a model tested on a specific [PITH_FULL_IMAGE:figures/full_fig_p036_26.png]
Figure 27
Figure 27. Figure 27: Principal Component and SHAP Analysis. particular, similar behavior of the second and third principal components is observed, termed the slope and curvature components in Litterman and Scheinkman [1991]. This PCA plot reveals parabolic and cubic shapes in the second a…
Figure 28
Figure 28. Figure 28: Option price variations for different samples. [PITH_FULL_IMAGE:figures/full_fig_p038_28.png]
Figure 29
Figure 29. Figure 29: Distribution of monotonicity violation distances. Left: strike price violations measured in steps of $5. [PITH_FULL_IMAGE:figures/full_fig_p039_29.png]
Figure 30
Figure 30. Figure 30: 3D distribution of violation distance and magnitude. Each surface shows how often violations occur at [PITH_FULL_IMAGE:figures/full_fig_p039_30.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 6 canonical work pages

  1. [3]

    doi: https://doi.org/10.1016/j.jfineco.2021.12.007

    ISSN 0304-405X. doi: https://doi.org/10.1016/j.jfineco.2021.12.007. URL https://www.sciencedirect.com/science/article/pii/S0304405X21005249. Jay Cao, Jacky Chen, and John Hull. A neural network approach to understanding implied volatility movements.Quantitative Finance, 20(9):1405–1413,

  2. [6]

    The us treasury yield curve: 1961 to the present

    Refet S G¨ urkaynak, Brian Sack, and Jonathan H Wright. The us treasury yield curve: 1961 to the present. Journal of monetary Economics, 54(8):2291–2304,

  3. [8]

    Deep learning calibration of option pricing models: some pitfalls and solutions.arXiv preprint arXiv:1906.03507,

    Andrey Itkin. Deep learning calibration of option pricing models: some pitfalls and solutions.arXiv preprint arXiv:1906.03507,

  4. [9]

    Scaling laws for neural language models.arXiv preprint arXiv:2001.08361,

    Jared Kaplan, Sam McCandlish, Tom Henighan, Tom B Brown, Benjamin Chess, Rewon Child, Scott Gray, Alec Radford, Jeffrey Wu, and Dario Amodei. Scaling laws for neural language models.arXiv preprint arXiv:2001.08361,

  5. [13]

    Machine learning for option pricing: an empirical investigation of network architectures.arXiv preprint arXiv:2307.07657,

    Laurens Van Mieghem, Antonis Papapantoleon, and Jonas Papazoglou-Hennig. Machine learning for option pricing: an empirical investigation of network architectures.arXiv preprint arXiv:2307.07657,

  6. [14]

    The last six months always serve as the test set

    46 A Tables with Performance Results In each table, the Date column indicates the period covered, spanning from the beginning of the training set to the end of the test set. The last six months always serve as the test set. The format YY/MM - YY/MM is used, where the first date marks the start of the training set and the second denotes the end of the six-...

  7. [1986]

    Tim Bollerslev, Michael Gibson, and Hao Zhou

    doi: 10.1016/0304-4076(86)90063-1. Tim Bollerslev, Michael Gibson, and Hao Zhou. Dynamic estimation of volatility risk premia and investor risk aversion from option-implied and realized volatilities.Journal of econometrics, 160(1):235–245,

  8. [1994]

    Machine learning methods for pricing financial derivatives.arXiv preprint arXiv:2406.00459,

    Lei Fan and Justin Sirignano. Machine learning methods for pricing financial derivatives.arXiv preprint arXiv:2406.00459,

Show all 14 references
  1. [2002]

    Application of machine learning in option pricing: A review

    Wenda Li. Application of machine learning in option pricing: A review. In2022 7th International Conference on Social Sciences and Economic Development (ICSSED 2022), pages 209–214. Atlantis Press,

  2. [2016]

    Black–scholes option pricing using machine learning

    Shreyan Sood, Tanmay Jain, Nishant Batra, and HC Taneja. Black–scholes option pricing using machine learning. InProceedings of International Conference on Data Science and Applications: ICDSA 2022, Volume 1, pages 481–493. Springer,

  3. [2017]

    Dynamics of the implied volatility surface: Theory and empirical evidence

    45 Jacinto Marabel Romo. Dynamics of the implied volatility surface: Theory and empirical evidence. Quantitative Finance DOI, 10(14697688.2012):686668,

  4. [2020]

    Deep learning for exotic option valuation.arXiv preprint arXiv:2103.12551,

    Jay Cao, Jacky Chen, John Hull, and Zissis Poulos. Deep learning for exotic option valuation.arXiv preprint arXiv:2103.12551,

  5. [2021]

    Differential machine learning.arXiv preprint arXiv:2005.02347,

    Brian Huge and Antoine Savine. Differential machine learning.arXiv preprint arXiv:2005.02347,

  6. [2022]

    Empirical performance of alternative option pricing models.The Journal of finance, 52(5):2003–2049,

    Gurdip Bakshi, Charles Cao, and Zhiwu Chen. Empirical performance of alternative option pricing models.The Journal of finance, 52(5):2003–2049,

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.