Monte Carlo integration using quantum computers has been widely investigated, including applications to concrete problems. It is known that quantum algorithms based on quantum amplitude estimation (QAE) can compute an integral with a smaller number of iterative calls of the quantum circuit which calculates the integrand, than classical methods call the integrand subroutine. However, the issues about the iterative operations in the integrand circuit have not been discussed so much. That is, in the high-dimensional integration, many random numbers are used for calculation of the integrand and in some cases similar calculations are repeated to obtain one sample value of the integrand. In this paper, we point out that we can reduce the number of such repeated operations by a combination of the nested QAE and the use of pseudorandom numbers (PRNs), if the integrand has the separable form with respect to contributions from distinct random numbers. The use of PRNs, which the authors originally proposed in the context of the quantum algorithm for Monte Carlo, is the key factor also in this paper, since it enables parallel computation of the separable terms in the integrand. Furthermore, we pick up one use case of this method in finance, the credit portfolio risk measurement, and estimate to what extent the complexity is reduced.
Monte Carlo integration using quantum computers has been widely investigated, including applications to concrete problems. It is known that quantum algorithms based on quantum amplitude estimation (QAE) can compute an integral with a smaller number of iterative calls of the quantum circuit which calculates the integrand, than classical methods call the integrand subroutine. However, the issues about the iterative operations in the integrand circuit have not been discussed so much. That is, in the high-dimensional integration, many random numbers are used for calculation of the integrand and in some cases similar calculations are repeated to obtain one sample value of the integrand. In this paper, we point out that we can reduce the number of such repeated operations by a combination of the nested QAE and the use of pseudorandom numbers (PRNs), if the integrand has the separable form with respect to contributions from distinct random numbers. The use of PRNs, which the authors originally proposed in the context of the quantum algorithm for Monte Carlo, is the key factor also in this paper, since it enables parallel computation of the separable terms in the integrand. Furthermore, we pick up one use case of this method in finance, the credit portfolio risk measurement, and estimate to what extent the complexity is reduced.
Quantum algorithms for Monte Carlo integration (QMCI), which are based on quantum amplitude estimation (QAE), provide quadratic speed-up compared with classical counterparts, and are therefore widely investigated, along with applications to industries such as finance. One typical feature of problems in finance is high dimensionality, or, in other words, the necessity to create numerous random numbers (RNs) to calculate the integrand. In such a situation, the original implementation of QMCI, in which states encoding the probability distributions of the RNs are generated on different quantum registers, requires a large number of qubits. In this poster, we propose an implementation of QMCI, in which pseudo-random numbers (PRNs) are sequentially generated on one register, and therefore the qubit number is tremendously reduced. Moreover, we show that, when the integrand has the form such that contributions from different RNs are separable into different terms, we can also achieve time complexity reduction with respect to the number of dimensions by combining the nested QAE and use of PRNs.
It is known that quantum computers can speed up Monte Carlo simulation compared to classical counterparts. There are already some proposals of application of the quantum algorithm to practical problems, including quantitative finance. In many problems in finance to which Monte Carlo simulation is applied, many random numbers are required to obtain one sample value of the integrand, since those problems are extremely high-dimensional integrations, for example, risk measurement of credit portfolio. This leads to the situation that the required qubit number is too large in the naive implementation where a quantum register is allocated per random number. In this paper, we point out that we can reduce qubits keeping quantum speed up if we perform calculation similar to the classical one, that is, estimate the average of integrand values sampled by a pseudo-random-number generator (PRNG) implemented on a quantum circuit. We present not only the overview of the idea but also concrete implementation of PRNG and application to credit risk measurement. Actually, reduction of qubits is a trade-off against increase of circuit depth. Therefore, full reduction might be impractical, but such a trade-off between speed and memory space will be important in adjustment of calculation setting considering machine specs, if large-scale Monte Carlo simulation by quantum computer is in operation in the future.
The approximation of an integral with numerical integration methods is a relevant problem in many fields, from finance to chemistry. Monte Carlo methods are a family of methods widely used for numerical integration but which are not always applicable to industrial tasks. Instead, there are quantum computing methods such as Quantum Amplitude Estimation able to offer a speed-up compared to the Monte Carlo methods for numerical integration which promise a significant impact, especially in finance. The Quantum Amplitude Estimation algorithm is a fundamental quantum algorithm with the potential to achieve a quadratic speed-up. Unfortunately, being competitive over classical Monte Carlo methods would require hardware performances currently not available. Therefore, alternative methods have been developed in order to exploit fewer resources while maintaining an advantageous theoretical scaling. We present a comparison between the standard QAE algorithm and two NISQ-friendly versions of the QAE on a numerical integration task. The Monte Carlo technique of Metropolis-Hastings instead is used as a classical benchmark. The scaling of the different algorithms are evaluated in terms of estimation error as a function of the number of samples, computational time, and length of the quantum circuits required by the quantum solutions. The conclusion highlights that an alternative approach is preferable with respect to employng the phase estimation routine. Indeed, the maximum likelihood estimation guarantees the best trade-off between the length of the quantum circuits and the precision in the integral estimation.
Monte Carlo simulations underpin computational quantitative finance, where they are used for derivative pricing, risk quantification and scenario analysis whenever closed-form solutions are unavailable. Their main limitation is the slow O(N^-1/2) convergence of classical sampling. Quantum amplitude estimation promises a quadratic improvement in query complexity over classical Monte Carlo, reducing the sample requirement from O(eps^-2) to O(eps^-1) under idealised conditions. This thesis measures the extent to which the quadratic improvement in query complexity survives a concrete risk task under gate noise when compared to classical Monte Carlo. A four-obligor correlated portfolio, whose loss distribution can be computed exactly, is encoded into quantum circuits, and its value at risk and expected shortfall are estimated by iterative amplitude estimation under exactly simulated depolarising noise on every two-qubit gate. Six configurations of 11 to 15 qubits vary the portfolio, the factor grid and the arithmetic construction. At matched loss-model evaluation budgets and without noise, amplitude estimation reduces the RMSE of the tail probability against classical Monte Carlo by a factor of 13.6 at the largest budget tested. The error budget is dominated by the circuit encoding, not the estimator. As a result, the value at risk at the 95% level picks the wrong loss level in none of the noiseless repetitions, while the encoding understates expected shortfall by 0.177 loss units. Under noise, the signal decays with every amplification round. The number of rounds that still help falls from 64 at an error rate of 3x10^-6 to none at 3x10^-2. This work concludes that amplification benefits only while the product of gate count and error rate stays small against one. Because every obligor added enlarges the circuit, the error rate at which amplification remains useful falls as the portfolio grows, and advantage starts depending on circuit compression instead of estimator design.
Abstract Quantum algorithms for the pricing of financial derivatives have been discussed in recent papers. However, the pricing model discussed in those papers is too simple for practical purposes. It motivates us to consider how to implement more complex models used in financial institutions. In this paper, we consider the local volatility (LV) model, in which the volatility of the underlying asset price depends on the price and time. As in previous studies, we use the quantum amplitude estimation (QAE) as the main source of quantum speedup and discuss the state preparation step of the QAE, or equivalently, the implementation of the asset price evolution. We compare two types of state preparation: One is the amplitude encoding (AE) type, where the probability distribution of the derivative’s payoff is encoded to the probabilistic amplitude. The other is the pseudo-random number (PRN) type, where sequences of PRNs are used to simulate the asset price evolution as in classical Monte Carlo simulation. We present detailed circuit diagrams for implementing these preparation methods in fault-tolerant quantum computation and roughly estimate required resources such as the number of qubits and T-count.
Numerical integration with Monte Carlo methods is a central computational task in many scientific and industrial applications, including financial derivative pricing and risk management. Classical Monte Carlo algorithms are computationally demanding: achieving an accuracy ε typically requires a number of function evaluations scaling as O(1/ε2). Quantum-accelerated Monte Carlo methods based on quantum amplitude estimation can in principle quadratically improve this dependence. However, \textit{quasi}-Monte Carlo methods have not been explored in the quantum context. In this work, we introduce a quantum quasi-Monte Carlo algorithm that combines low-discrepancy nets with quantum amplitude estimation. The proposed method prepares the quasi-random point set coherently in superposition. The method does not yield an asymptotic improvement over classical quasi-Monte Carlo, since the total error separates into a discretization error, determined by the finite net, and a quantum estimation error. Instead, we explore a pre-asymptotic advantage window: for a target accuracy that would classically require 2q low discrepancy points, one can prepare a higher-resolution net of size 2Q, with Q>q, in superposition and reach the same accuracy using significantly fewer function queries. This window can be controlled by tuning the circuit resolution and amplitude-estimation parameters, making the approach relevant for practical regimes where the number of queries is finite rather than asymptotically large.