Skip to main content
Article
100%
HD

Imagine having to decide the exact moment to exercise a complex financial option where hundreds of variables are at play. Classical methods break down: too many calculations, too much memory. Enter the Deep Penalty Method (DPM), a new algorithm that fuses deep neural networks with partial differential equation mathematics.

The trick? Transform the "optimal stopping" problem β€” typical of American options β€” into a PDE with a penalty term. This term forces the solution to respect the early-exercise constraint without explicitly tracking the free boundary, notoriously hard in high dimensions.

Then the Deep BSDE framework kicks in: a neural network learns to represent the solution of the associated stochastic equation, leveraging Monte Carlo simulations. The result is a method that scales well with dimensions: tests on Bermudan max-call options up to 100 risk factors show surprising accuracy and speed.

There's a catch: the error depends on the penalty parameter Ξ» and the time step h. Theory bounds the total error by the loss function plus O(1/Ξ») + O(Ξ»h) + O(√h). Choosing Ξ» β‰ˆ 1/√h balances the terms and yields a convergence rate of order 1/2. In practice, doubling time resolution halves the error, provided Ξ» is adjusted accordingly.

A breakthrough for quantitative finance and beyond: wherever sequential decisions in high-dimensional spaces arise, from optimal control to robotics, DPM opens previously impractical paths.

Bonus
100%
HD

The Deep Penalty Method (DPM) tackles high-dimensional optimal stopping through a dual approximation. First, the penalty method transforms the variational inequality of the free-boundary problem into a semilinear PDE: βˆ‚β‚œu + Lu + Ξ» max(0, g - u) = 0, where g is the obstacle (exercise payoff), L the infinitesimal generator of the underlying diffusion, and Ξ» > 0 the penalty parameter. As Ξ» β†’ ∞, the penalized solution u^Ξ» converges to the true solution at rate O(1/Ξ»).

Second, the penalized PDE is solved via the Deep BSDE framework (Weinan et al., 2017). The solution u(t,Xβ‚œ) is represented by a deep neural network with parameters ΞΈ, and the loss function derived from the time-discretized BSDE scheme is minimized: L(ΞΈ) = 𝔼[|u^ΞΈ(T,X_T) - g(X_T)|² + Ξ£β‚™ |u^ΞΈ(tβ‚™,X_{tβ‚™}) - u^ΞΈ(tβ‚™β‚Šβ‚,X_{tβ‚™β‚Šβ‚}) - f(tβ‚™,X_{tβ‚™},u^ΞΈ,βˆ‡u^ΞΈ)Ξ”t|²]. The driver f incorporates the penalty Ξ» max(0, g - u).

Total error analysis combines three contributions: penalty error O(1/λ), time-discretization error O(√h) from the Euler-Maruyama scheme, and a coupling term O(λh) arising from the penalty nonlinearity in discretization. The optimum balances λ ∼ h^{-1/2}, yielding overall convergence O(√h).

In numerical tests on 100-asset Bermudan max-call options (multi-dimensional Black-Scholes model), DPM achieves relative errors < 1% versus Longstaff-Schwartz Monte Carlo regression and sparse-grid PDE benchmarks, with runtimes of minutes on a single GPU. The network uses a ResNet architecture with 4 residual blocks, 128 neurons per layer, ReLU activation, trained with Adam (lr=1e-3) for 20k iterations, batch size 256.

Ask AI

Explore key aspects of this story with 3 quick insights

Source

Authors:
Yunfei PengPengyu WeiWei WeiXiaole Xue
Date Published:
Original Article arxiv.org
Fact-check verified on 27/10/2026 Primary Study / Original Source
Read our Editorial Guidelines & AI Ethics Policy
Creative Commons License CC BY-NC-SA 4.0

Keep exploring

Discover the full archive of science bites or head back to home

Riconnessione al server... Aggiornamento della sessione in corso
Connessione interrotta Impossibile stabilire il collegamento
Sessione aggiornata Ricarica per visualizzare i contenuti