Variational quantum Monte Carlo (VMC) combined with neural-network quantum states offers a novel angle of attack on the curse-of-dimensionality encountered in a particular class of partial differential equations (PDEs); namely, the real- and imaginary time-dependent Schrödinger equation. In this paper, we present a simple generalization of VMC applicable to arbitrary time-dependent PDEs, showcasing the technique in the multi-asset Black-Scholes PDE for pricing European options contingent on many correlated underlying assets.
Variational quantum Monte Carlo (VMC) combined with neural-network quantum states offers a novel angle of attack on the curse-of-dimensionality encountered in a particular class of partial differential equations (PDEs); namely, the real- and imaginary time-dependent Schrödinger equation. In this paper, we present a simple generalization of VMC applicable to arbitrary time-dependent PDEs, showcasing the technique in the multi-asset Black-Scholes PDE for pricing European options contingent on many correlated underlying assets.
The variational quantum Monte Carlo (VQMC) method has received significant attention because of its ability to overcome the curse of dimensionality inherent in many-body quantum systems, by representing the exponentially complex quantum states variationally with machine learning models. We develop novel training strategies to improve the scalability of VQMC, and build parallelization frameworks for solving large-scale problems. The application of our method is extended to quantum chemistry and financial derivative pricing. For quantum chemistry, we build a pre-processing pipeline serving as an interface connecting molecular information and VQMC, and achieve remarkable performance in comparison with the classical approximate methods. On the other hand, we present a simple generalization of VQMC applicable to arbitrary linear PDEs, showcasing the technique in the Black-Scholes equation for pricing European contingent claims dependent on many underlying assets. We also introduce meta-learning and multi-fidelity active learning as exotic components to VQMC, which, under some reasonable assumptions on the problem formulation, can further improve the convergence and the sampling efficiency of our method.
Following the recent great advance of quantum computing technology, there are growing interests in its applications to industries, including finance. In this paper, we focus on derivative pricing based on solving the Black-Scholes partial differential equation by finite difference method (FDM), which is a suitable approach for some types of derivatives but suffers from the {\it curse of dimensionality}, that is, exponential growth of complexity in the case of multiple underlying assets. We propose a quantum algorithm for FDM-based pricing of multi-asset derivative with exponential speedup with respect to dimensionality compared with classical algorithms. The proposed algorithm utilizes the quantum algorithm for solving differential equations, which is based on quantum linear system algorithms. Addressing the specific issue in derivative pricing, that is, extracting the derivative price for the present underlying asset prices from the output state of the quantum algorithm, we present the whole of the calculation process and estimate its complexity. We believe that the proposed method opens the new possibility of accurate and high-speed derivative pricing by quantum computers.
Following the recent great advance of quantum computing technology, there are growing interests in its applications to industries, including finance. In this article, we focus on derivative pricing based on solving the Black–Scholes partial differential equation by the finite-difference method (FDM), which is a suitable approach for some types of derivatives but suffers from thecurse of dimensionality, that is, exponential growth of complexity in the case of multiple underlying assets. We propose a quantum algorithm for FDM-based pricing of multi-asset derivative with exponential speedup with respect to dimensionality compared with classical algorithms. The proposed algorithm utilizes the quantum algorithm for solving differential equations, which is based on quantum linear system algorithms. Addressing the specific issue in derivative pricing, that is, extracting the derivative price for the present underlying asset prices from the output state of the quantum algorithm, we present the whole of the calculation process and estimate its complexity. We believe that the proposed method opens the new possibility of accurate and high-speed derivative pricing by quantum computers.
A variational quantum algorithm for numerically solving partial differential equations (PDEs) on a quantum computer was proposed by Lubasch et al. [Phys. Rev. A 101, 010301(R) (2020)]. In this paper we generalize the method introduced by Lubasch et al. to cover a broader class of nonlinear PDEs as well as multidimensional PDEs and study the performance of the variational quantum algorithm on several example equations. Specifically, we show via numerical simulations that the algorithm can solve instances of the single-asset Black-Scholes equation with a nontrivial nonlinear volatility model, the double-asset Black-Scholes equation, the Buckmaster equation, and the deterministic Kardar-Parisi-Zhang equation. Our simulations use up to n=12 Ansatz qubits, computing PDE solutions with 2n grid points. We also perform proof-of-concept demonstrations with a trapped-ion quantum processor from IonQ [Nat. Commun. 10, 5464 (2019)], showing accurate computation of two representative expectation values needed for the calculation of a single time step of the nonlinear Black-Scholes equation. Through our classical simulations and demonstrations on quantum hardware, we identify and discuss several open challenges for using quantum variational methods to solve PDEs in a regime with a large number (much greater than 220) of grid points, but also a practical number of gates per circuit and circuit shots.
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\%.
Inspired by recent progress in quantum algorithms for ordinary and partial differential equations, we study quantum algorithms for stochastic differential equations (SDEs). Firstly we provide a quantum algorithm that gives a quadratic speed-up for multilevel Monte Carlo methods in a general setting. As applications, we apply it to compute expectation values determined by classical solutions of SDEs, with improved dependence on precision. We demonstrate the use of this algorithm in a variety of applications arising in mathematical finance, such as the Black-Scholes and Local Volatility models, and Greeks. We also provide a quantum algorithm based on sublinear binomial sampling for the binomial option pricing model with the same improvement.