Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2607.28435

Classical Tensor Network and Quantum Fourier Transform Approaches for Large-Scale Carr-Madan Option Pricing

Sascha Hauck, Ivica Turkalj

quant-phSubmitted Jul 30, 2026• Announced Jul 31, 2026• 0 citations• Revised Aug 12, 2026 (v2)
v1 · Jul 30, 2026v2 · Aug 12, 2026
arXivHTMLPDF

Abstract

Fourier-based methods are among the most widely used techniques for pricing European options when the characteristic function of the underlying asset process is available. Their applicability to increasingly fine discretizations, however, is limited by the rapidly growing memory requirements of classical Fourier transforms, which become a computational bottleneck for large-scale pricing problems. In this work, we overcome this limitation by reformulating the Carr-Madan pricing framework using tensor networks. Specifically, we employ the Superfast Fourier Transform (SFFT), a compressed Tensor Train representation of the Quantum Fourier Transform (QFT), and apply it directly to tensorized option pricing without ever explicitly constructing exponentially large vectors or Fourier operators. This formulation also enables a direct comparison between the classical tensor network algorithm and its quantum counterpart through QFT-based option pricing on quantum simulators and quantum hardware. Numerical experiments for European call options demonstrate that the proposed SFFT method maintains pricing accuracy while substantially reducing memory requirements and achieving subexponential computational scaling compared with conventional FFT-based pricing. The accompanying quantum simulations and hardware executions enable a direct comparison between the classical tensor network formulation and its QFT-based quantum counterpart, showing that both approaches avoid the exponential scaling of conventional Fourier implementations and provide complementary perspectives on large-scale option pricing. Together, these results establish a unified framework connecting classical Fourier pricing, tensor network algorithms, and quantum computing approaches, demonstrating how tensorized Fourier methods can provide scalable alternatives for high-dimensional financial computations.

Related papers

Pricing of European Calls with the Quantum Fourier Transform

Ewen, Tom

The accurate valuation of financial derivatives plays a pivotal role in the finance industry. Although closed formulas for pricing are available for certain models and option types, exemplified by the European Call and Put options in the Black-Scholes Model, the use of either more complex models or more sophisticated options precludes the existence of such formulas, thereby requiring alternative approaches. The Monte Carlo simulation, an alternative approach effective in nearly all scenarios, has already been challenged by quantum computing techniques that leverage Amplitude Estimation. Despite its theoretical promise, this approach currently faces limitations due to the constraints of hardware in the Noisy Intermediate-Scale Quantum (NISQ) era. In this study, we introduce and analyze a quantum algorithm for pricing European call options across a broad spectrum of asset models. This method transforms a classical approach, which utilizes the Fast Fourier Transform (FFT), into a quantum algorithm, leveraging the efficiency of the Quantum Fourier Transform (QFT). Furthermore, we compare this novel algorithm with existing quantum algorithms for option pricing.

arXiv:2404.14115•Apr 22, 2024•0 citations

Quantum Computing in Option Pricing

Ewen, Tom

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.

•Jan 1, 2025•0 citations

A highly efficient tensor network algorithm for multi-asset Fourier options pricing

Michael J. Kastoryano, Nicola Pancotti

Risk assessment and in particular derivatives pricing is one of the core areas in computational finance and accounts for a sizeable fraction of the global computing resources of the financial industry. We outline a quantum-inspired algorithm for multi-asset options pricing. The algorithm is based on tensor networks, which have allowed for major conceptual and numerical breakthroughs in quantum many body physics and quantum computation. In the proof-of-concept example explored, the tensor network approach yields several orders of magnitude speedup over vanilla Monte Carlo simulations. We take this as good evidence that the use of tensor network methods holds great promise for alleviating the computation burden of risk evaluation in the financial and other industries, thus potentially lowering the carbon footprint these simulations incur today.

arXiv:2203.02804•Mar 5, 2022•2 citations

STN-GPR: A Singularity Tensor Network Framework for Efficient Option Pricing

Dominic Gribben, Carolina Allende, Alba Villarino +5 more

We develop a tensor-network surrogate for option pricing, targeting large-scale portfolio revaluation problems arising in market risk management (e.g., VaR and Expected Shortfall computations). The method involves representing high-dimensional price surfaces in tensor-train (TT) form using TT-cross approximation, constructing the surrogate directly from black-box price evaluations without materializing the full training tensor. For inference, we use a Laplacian kernel and derive TT representations of the kernel matrix and its closed-form inverse in the noise-free setting, enabling TT-based Gaussian process regression without dense matrix factorization or iterative linear solves. We found that hyperparameter optimization consistently favors a large kernel length-scale and show that in this regime the GPR predictor reduces to multilinear interpolation for off-grid inputs; we also derive a low-rank TT representation for this limit. We evaluate the approach on five-asset basket options over an eight dimensional parameter space (asset spot levels, strike, interest rate, and time to maturity). For European geometric basket puts, the tensor surrogate achieves lower test error at shorter training times than standard GPR by scaling to substantially larger effective training sets. For American arithmetic basket puts trained on LSMC data, the surrogate exhibits more favorable scaling with training-set size while providing millisecond-level evaluation per query, with overall runtime dominated by data generation.

q-fin.PRCross-listedarXiv:2603.26318•Mar 30, 2026

Monte-Carlo Option Pricing in Quantum Parallel

Robert Scriba, Yuying Li, Jingbo B. Wang

Abstract Financial derivative pricing is a significant challenge in finance, involving the valuation of instruments like options based on underlying assets. While some cases have simple solutions, many require complex classical computational methods like Monte Carlo simulations and numerical techniques. However, as derivative complexities increase, these methods face limitations in computational power. Cases involving Non-Vanilla Basket pricing, American Options, and derivative portfolio risk analysis need extensive computations in higher-dimensional spaces, posing challenges for classical computers. Quantum computing presents a promising avenue by harnessing quantum superposition and entanglement, allowing the handling of high-dimensional spaces effectively. In this paper, we introduce a self-contained and all-encompassing quantum algorithm that operates without reliance on oracles or presumptions. More specifically, we develop an effective stochastic method for simulating exponentially many potential asset paths in quantum parallel, leading to a highly accurate final distribution of stock prices. Furthermore, we demonstrate how this algorithm can be extended to price more complex options and analyze risk within derivative portfolios.

arXiv:2505.09459•Sep 9, 2025•1 citation

End-to-End PDE-Based Quantum Algorithms for Multi-Asset Option Pricing under Local and Stochastic Volatility

Nikita Guseynov, Nana Liu, Chi Seng Pun +1 more

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=2^n$ 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 $\widetilde{O}(d^2 N^{2+d/2})$ for local-volatility Black--Scholes and $\widetilde{O}(d^2 N^{d+2})$ for Heston. Relative to grid-based finite-difference baselines, these scalings correspond to polynomial improvement factors $N^{d/2}$ and $N^d$, 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.

quant-pharXiv:2605.26610•May 27, 2026•0 citations
Fourier-based methods are among the most widely used techniques for pricing European options when the characteristic function of the underlying asset process is available. Their applicability to increasingly fine discretizations, however, is limited by the rapidly growing memory requirements of classical Fourier transforms, which become a computational bottleneck for large-scale pricing problems.
In this work, we overcome this limitation by reformulating the Carr-Madan pricing framework using tensor networks. Specifically, we employ the Superfast Fourier Transform (SFFT), a compressed Tensor Train representation of the Quantum Fourier Transform (QFT), and apply it directly to tensorized option pricing without ever explicitly constructing exponentially large vectors or Fourier operators. This formulation also enables a direct comparison between the classical tensor network algorithm and its quantum counterpart through QFT-based option pricing on quantum simulators and quantum hardware. Numerical experiments for European call options demonstrate that the proposed SFFT method maintains pricing accuracy while substantially reducing memory requirements and achieving subexponential computational scaling compared with conventional FFT-based pricing. The accompanying quantum simulations and hardware executions enable a direct comparison between the classical tensor network formulation and its QFT-based quantum counterpart, showing that both approaches avoid the exponential scaling of conventional Fourier implementations and provide complementary perspectives on large-scale option pricing.
Together, these results establish a unified framework connecting classical Fourier pricing, tensor network algorithms, and quantum computing approaches, demonstrating how tensorized Fourier methods can provide scalable alternatives for high-dimensional financial computations.
The accurate valuation of financial derivatives plays a pivotal role in the finance industry. Although closed formulas for pricing are available for certain models and option types, exemplified by the European Call and Put options in the Black-Scholes Model, the use of either more complex models or more sophisticated options precludes the existence of such formulas, thereby requiring alternative approaches. The Monte Carlo simulation, an alternative approach effective in nearly all scenarios, has already been challenged by quantum computing techniques that leverage Amplitude Estimation. Despite its theoretical promise, this approach currently faces limitations due to the constraints of hardware in the Noisy Intermediate-Scale Quantum (NISQ) era. In this study, we introduce and analyze a quantum algorithm for pricing European call options across a broad spectrum of asset models. This method transforms a classical approach, which utilizes the Fast Fourier Transform (FFT), into a quantum algorithm, leveraging the efficiency of the Quantum Fourier Transform (QFT). Furthermore, we compare this novel algorithm with existing quantum algorithms for option pricing.
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.
Risk assessment and in particular derivatives pricing is one of the core areas in computational finance and accounts for a sizeable fraction of the global computing resources of the financial industry. We outline a quantum-inspired algorithm for multi-asset options pricing. The algorithm is based on tensor networks, which have allowed for major conceptual and numerical breakthroughs in quantum many body physics and quantum computation. In the proof-of-concept example explored, the tensor network approach yields several orders of magnitude speedup over vanilla Monte Carlo simulations. We take this as good evidence that the use of tensor network methods holds great promise for alleviating the computation burden of risk evaluation in the financial and other industries, thus potentially lowering the carbon footprint these simulations incur today.
We develop a tensor-network surrogate for option pricing, targeting large-scale portfolio revaluation problems arising in market risk management (e.g., VaR and Expected Shortfall computations). The method involves representing high-dimensional price surfaces in tensor-train (TT) form using TT-cross approximation, constructing the surrogate directly from black-box price evaluations without materializing the full training tensor. For inference, we use a Laplacian kernel and derive TT representations of the kernel matrix and its closed-form inverse in the noise-free setting, enabling TT-based Gaussian process regression without dense matrix factorization or iterative linear solves. We found that hyperparameter optimization consistently favors a large kernel length-scale and show that in this regime the GPR predictor reduces to multilinear interpolation for off-grid inputs; we also derive a low-rank TT representation for this limit. We evaluate the approach on five-asset basket options over an eight dimensional parameter space (asset spot levels, strike, interest rate, and time to maturity). For European geometric basket puts, the tensor surrogate achieves lower test error at shorter training times than standard GPR by scaling to substantially larger effective training sets. For American arithmetic basket puts trained on LSMC data, the surrogate exhibits more favorable scaling with training-set size while providing millisecond-level evaluation per query, with overall runtime dominated by data generation.
Abstract Financial derivative pricing is a significant challenge in finance, involving the valuation of instruments like options based on underlying assets. While some cases have simple solutions, many require complex classical computational methods like Monte Carlo simulations and numerical techniques. However, as derivative complexities increase, these methods face limitations in computational power. Cases involving Non-Vanilla Basket pricing, American Options, and derivative portfolio risk analysis need extensive computations in higher-dimensional spaces, posing challenges for classical computers. Quantum computing presents a promising avenue by harnessing quantum superposition and entanglement, allowing the handling of high-dimensional spaces effectively. In this paper, we introduce a self-contained and all-encompassing quantum algorithm that operates without reliance on oracles or presumptions. More specifically, we develop an effective stochastic method for simulating exponentially many potential asset paths in quantum parallel, leading to a highly accurate final distribution of stock prices. Furthermore, we demonstrate how this algorithm can be extended to price more complex options and analyze risk within derivative portfolios.
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=2nN=2^nN=2n grid points per spatial direction and ddd 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)\widetilde{O}(d^2 N^{2+d/2})O(d2N2+d/2) for local-volatility Black--Scholes and O~(d2Nd+2)\widetilde{O}(d^2 N^{d+2})O(d2Nd+2) for Heston. Relative to grid-based finite-difference baselines, these scalings correspond to polynomial improvement factors Nd/2N^{d/2}Nd/2 and NdN^dNd, 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.