No Arabic abstract
Quantum Monte Carlo and quantum simulation are both important tools for understanding quantum many-body systems. As a classical algorithm, quantum Monte Carlo suffers from the sign problem, preventing its applications to most fermion systems and real time dynamics. In this paper, we introduce a novel non-variational algorithm using quantum simulation as a subroutine to accelerate quantum Monte Carlo by easing the sign problem. The quantum subroutine can be implemented with shallow circuits and, by incorporating error mitigation, reduce the Monte Carlo variance by orders of magnitude even when the circuit noise is significant. As such, the proposed quantum algorithm is applicable to near-term noisy quantum hardware.
The simulation of strongly correlated many-electron systems is one of the most promising applications for near-term quantum devices. Here we use a class of eigenvalue solvers (presented in Phys. Rev. Lett. 126, 070504 (2021)) in which a contraction of the Schrodinger equation is solved for the two-electron reduced density matrix (2-RDM) to resolve the energy splittings of ortho-, meta-, and para-isomers of benzyne ${textrm C_6} {textrm H_4}$. In contrast to the traditional variational quantum eigensolver, the contracted quantum eigensolver solves an integration (or contraction) of the many-electron Schrodinger equation onto the two-electron space. The quantum solution of the anti-Hermitian part of the contracted Schrodinger equation (qACSE) provides a scalable approach with variational parameters that has its foundations in 2-RDM theory. Experimentally, a variety of error mitigation strategies enable the calculation, including a linear shift in the 2-RDM targeting the iterative nature of the algorithm as well as a projection of the 2-RDM onto the convex set of approximately $N$-representable 2-RDMs defined by the 2-positive (DQG) $N$-representability conditions. The relative energies exhibit single-digit millihartree errors, capturing a large part of the electron correlation energy, and the computed natural orbital occupations reflect the significant differences in the electron correlation of the isomers.
Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic environment on a quantum simulator may be achieved by coupling independent ancilla qubits that facilitate energy exchange in an appropriate manner with the system and mimic an environment. This approach requires a large, and possibly exponential number of ancillary degrees of freedom which is impractical. In contrast, we develop a digital quantum algorithm that simulates interaction with an environment using a small number of ancilla qubits. By combining periodic modulation of the ancilla energies, or spectral combing, with periodic reset operations, we are able to mimic interaction with a large environment and generate thermal states of interacting many-body systems. We evaluate the algorithm by simulating preparation of thermal states of the transverse Ising model. Our algorithm can also be viewed as a quantum Markov chain Monte Carlo (QMCMC) process that allows sampling of the Gibbs distribution of a multivariate model. To demonstrate this we evaluate the accuracy of sampling Gibbs distributions of simple probabilistic graphical models using the algorithm.
Recently we find several candidates of quantum algorithms that may be implementable in near-term devices for estimating the amplitude of a given quantum state, which is a core subroutine in various computing tasks such as the Monte Carlo methods. One of those algorithms is based on the maximum likelihood estimate with parallelized quantum circuits; in this paper, we extend this method so that it can deal with the realistic noise effect. The validity of the proposed noise model is supported by an experimental demonstration on an IBM Q device, which accordingly enables us to predict the basic requirement on the hardware components (particularly the gate error) in quantum computers to realize the quantum speedup in the amplitude estimation task.
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.
Various quantum applications can be reduced to estimating expectation values, which are inevitably deviated by operational and environmental errors. Although errors can be tackled by quantum error correction, the overheads are far from being affordable for near-term technologies. To alleviate the detrimental effects of errors, quantum error mitigation techniques have been proposed, which require no additional qubit resources. Here, we benchmark the performance of a quantum error mitigation technique based on probabilistic error cancellation in a trapped-ion system. Our results clearly show that effective gate fidelities exceed physical fidelities, i.e. we surpass the break-even point of eliminating gate errors, by programming quantum circuits. The error rates are effectively reduced from $(1.10pm 0.12)times10^{-3}$ to $(1.44pm 5.28)times10^{-5}$ and from $(0.99pm 0.06)times10^{-2}$ to $(0.96pm 0.10)times10^{-3}$ for single- and two-qubit gates, respectively. Our demonstration opens up the possibility of implementing high-fidelity computations on a near-term noisy quantum device.