The solution of option-pricing problems may turn out to be computationally demanding due to non-linear and path-dependent payoffs, the high dimensionality arising from multiple underlying assets, and sophisticated models of price dynamics. In this context, quantum computing has been proposed as a means to address these challenges efficiently. Prevailing approaches either simulate the stochastic differential equations governing the forward dynamics of underlying asset prices or directly solve the backward pricing partial differential equation. Here, we present an end-to-end quantum algorithmic framework that solves the Kolmogorov forward (Fokker-Planck) partial differential equation for local-volatility models by mapping it to a Hamiltonian-simulation problem via the Schrödingerisation technique. The algorithm specifies how to prepare the initial quantum state, perform Hamiltonian simulation, and how to efficiently recover the option price via a swap test. In particular, the efficiency of the final solution recovery is an important advantage of solving the forward versus the backward partial differential equation. Thus, our end-to-end framework offers a potential route toward quantum advantage for challenging option-pricing tasks. In particular, we obtain a polynomial advantage in grid size for the discretization of a single dimension. Nevertheless, the true power of our methodology lies in pricing high-dimensional systems, such as baskets of options, because the quantum framework admits an exponential speedup with respect to dimension, overcoming the classical curse of dimensionality.
The solution of option-pricing problems may turn out to be computationally demanding due to non-linear and path-dependent payoffs, the high dimensionality arising from multiple underlying assets, and sophisticated models of price dynamics. In this context, quantum computing has been proposed as a means to address these challenges efficiently. Prevailing approaches either simulate the stochastic differential equations governing the forward dynamics of underlying asset prices or directly solve the backward pricing partial differential equation. Here, we present an end-to-end quantum algorithmic framework that solves the Kolmogorov forward (Fokker-Planck) partial differential equation for local-volatility models by mapping it to a Hamiltonian-simulation problem via the Schrödingerisation technique. The algorithm specifies how to prepare the initial quantum state, perform Hamiltonian simulation, and how to efficiently recover the option price via a swap test. In particular, the efficiency of the final solution recovery is an important advantage of solving the forward versus the backward partial differential equation. Thus, our end-to-end framework offers a potential route toward quantum advantage for challenging option-pricing tasks. In particular, we obtain a polynomial advantage in grid size for the discretization of a single dimension. Nevertheless, the true power of our methodology lies in pricing high-dimensional systems, such as baskets of options, because the quantum framework admits an exponential speedup with respect to dimension, overcoming the classical curse of dimensionality.
Multi-asset option pricing under local- and stochastic-volatility models leads naturally to high-dimensional parabolic PDEs. We develop an end-to-end quantum PDE framework for European option pricing under local-volatility Black--Scholes and Heston models. The framework takes classical contract and model data as input and returns classical estimates of selected option values. We solve the pricing PDEs after finite-difference discretization on spatial grids. For N=2n grid points per spatial direction and d assets, the end-to-end gate complexity for single-point recovery, counted in elementary CNOT gates and one-qubit Pauli-axis rotations, has leading grid-size dependence O(d2N2+d/2) for local-volatility Black--Scholes and O(d2Nd+2) for Heston. Relative to grid-based finite-difference baselines, these scalings correspond to polynomial improvement factors Nd/2 and Nd, respectively. These estimates translate to Clifford+T resources via standard compilation. We complement the complexity analysis with numerical benchmarks against standard classical methods. In the Heston setting, the framework recovers option prices across strikes together with the associated implied-volatility smile/skew. Overall, this work provides a complete end-to-end quantum pricing pipeline with explicit resource accounting and theoretical performance guarantees.
Pricing financial derivatives, in particular European-style options at different time-maturities and strikes, means a relevant problem in finance. The dynamics describing the price of vanilla options when constant volatilities and interest rates are assumed is governed by the Black-Scholes model, a linear parabolic partial differential equation with terminal value given by the payoff of the option contract and no additional boundary conditions. Here, we present a digital quantum algorithm to solve the Black-Scholes equation on a quantum computer by mapping it to the Schr\"odinger equation. The non-Hermitian nature of the resulting Hamiltonian is solved by embedding its propagator into an enlarged Hilbert space by using only one additional ancillary qubit. Moreover, due to the choice of periodic boundary conditions, given by the definition of the discretized momentum operator, we duplicate the initial condition, which substantially improves the stability and performance of the protocol. The algorithm shows a feasible approach for using efficient Hamiltonian simulation techniques as quantum signal processing to solve the price dynamics of financial derivatives on a digital quantum computer. Our approach differs from those based on Monte Carlo integration, exclusively focused on sampling the solution assuming the dynamics is known. We report expected accuracy levels comparable to classical numerical algorithms by using nine qubits to simulate its dynamics on a fault-tolerant quantum computer, and an expected success probability of the post-selection procedure due to the embedding protocol above 60%.
Pricing financial derivatives, in particular European-style options at different time-maturities and strikes, is a relevant financial problem. The dynamics describing the price of vanilla options when constant volatilities and interest rates are assumed, is governed by the Black-Scholes model, a linear parabolic partial differential equation with terminal value given by the pay-off of the option contract and no additional boundary conditions. Here, we present a digital quantum algorithm to solve Black-Scholes equation on a quantum computer for a wide range of relevant financial parameters by mapping it to the Schrodinger equation. The non-Hermitian nature of the resulting Hamiltonian is solved by embedding the dynamics into an enlarged Hilbert space, which makes use of only one additional ancillary qubit. Moreover, we employ a second ancillary qubit to transform initial condition into periodic boundary conditions, which substantially improves the stability and performance of the protocol. This algorithm shows a feasible approach for pricing financial derivatives on a digital quantum computer based on Hamiltonian simulation, technique which differs from those based on Monte Carlo simulations to solve the stochastic counterpart of the Black Scholes equation. Our algorithm remarkably provides an exponential speedup since the terms in the Hamiltonian can be truncated by a polynomial number of interactions while keeping the error bounded. We report expected accuracy levels comparable to classical numerical algorithms by using 10 qubits and 94 entangling gates on a fault-tolerant quantum computer, and an expected success probability of the post-selection procedure due to the embedding protocol above 60\%.
Stochastic differential equations (SDEs) model noisy dynamical systems across quantitative finance, chemistry, physics, and climate science. Classical solvers-Euler-Maruyama schemes, Monte Carlo methods, and finite-difference PDE solvers for Fokker-Planck equations-suffer from exponential growth in computational cost with dimension and accuracy requirements. We present a family of quantum algorithms for simulating and solving SDEs with improved scaling. Our approach combines Hamiltonian simulation for drift terms, quantum random number generation for noise increments, and quantum amplitude estimation for computing expectations. For linear and semilinear SDEs with Lipschitz coefficients, we prove that relevant quantities (option prices, expected hitting times, moments) can be estimated to accuracy ε in time Õ(poly(d)polylog(1/ε)), versus classical O(poly(d)/ε 2). For nonlinear SDEs, we introduce variational quantum schemes leveraging parameterized circuits, with a priori error bounds and trainability analysis avoiding barren plateaus. We analyze multi-asset Black-Scholes models, Heston stochastic volatility, and Langevin dynamics for molecular coarse-graining, quantifying resource requirements (qubits, gate depth, oracle calls). Numerical simulations on small instances demonstrate polynomial-to-super-polynomial empirical speedups in accuracy for fixed runtime, establishing SDE solving as a promising quantum-accelerated task with practical applications.
Options are financial instruments whose payoffs are determined by the value of an underlying asset. These underlying assets are typically stocks, interest rates, or exchange rates, but can also include commodities such as energy or agricultural products. Options serve diverse purposes, including hedging against price volatility and speculating on market movements. Additionally, they play a critical role in assessing certain positions on the balance sheets of insurance companies, particularly in the context of solvency capital requirement calculations. As an option’s payoff lies in the future and is stochastic in nature, finding a fair price today is a non-trivial task. The focus of this work is finding quantum algorithms for calculating these prices and compare them to algorithms designed for classical computers. This work follows two different directions on how to improve on classical methods for option pricing. The first direction is improving on Monte-Carlo (MC) methods and the second one is adapting Fourier based pricing algorithms. The Amplitude Estimation (AE) algorithm promises a quadratic speedup compared to classical MC. This fact makes this algorithm an interesting candidate for improvements by quantum algorithms. It was shown in the literature, that simple options with very crude discretization of the underlying stochastic models, can be implemented on a quantum device and evaluated with AE. One big challenge is the efficient implementation of the option’s payoff profile. To solve this problem colleagues and I have introduced a multi-objective genetic algorithm to automatically find well performing and efficient circuits to implement the payoff of non-path-dependent European options on one or more underlyings. Another strategy for pricing options is based on the fact, that in many models the characteristic function of the underlying is numerically more tractable than its distribution function. There are methods that make use of this fact and the performance of the Fast Fourier Transform (FFT) to efficiently calculate prices for European call options for many strikes in one go. The quantum version of the FFT, the Quantum Fourier Transform (QFT), is in some sense, even faster than its classical counterpart. We have adapted the existing classical method to make use of QFT, to benefit from its performance. Finally, the different strategies are benchmarked and compared regarding their accuracy, their potential for scaling and their adaptability to different styles of derivatives and models. The result of this comparison are requirements on the quantum hardware for the methods to be viable as well as use cases where each method has its benefits over the others.
This paper explores advancements in quantum algorithms for derivative pricing of exotics, a computational pipeline of fundamental importance in quantitative finance. For such cases, the classical Monte Carlo integration procedure provides the state-of-the-art provable, asymptotic performance: polynomial in problem dimension and quadratic in inverse-precision. While quantum algorithms are known to offer quadratic speedups over classical Monte Carlo methods, end-to-end speedups have been proven only in the simplified setting over the Black-Scholes geometric Brownian motion (GBM) model. This paper extends existing frameworks to demonstrate novel quadratic speedups for more practical models, such as the Cox-Ingersoll-Ross (CIR) model and a variant of Heston's stochastic volatility model, utilizing a characteristic of the underlying SDEs which we term fast-forwardability. Additionally, for general models that do not possess the fast-forwardable property, we introduce a quantum Milstein sampler, based on a novel quantum algorithm for sampling Lévy areas, which enables quantum multi-level Monte Carlo to achieve quadratic speedups for multi-dimensional stochastic processes exhibiting certain correlation types.
We also present an improved analysis of numerical integration for derivative pricing, leading to substantial reductions in the resource requirements for pricing GBM and CIR models. Furthermore, we investigate the potential for additional reductions using arithmetic-free quantum procedures. Finally, we critique quantum partial differential equation (PDE) solvers as a method for derivative pricing based on amplitude estimation, identifying theoretical barriers that obstruct achieving a quantum speedup through this approach. Our findings significantly advance the understanding of quantum algorithms in derivative pricing, addressing key challenges and open questions in the field.