Estimation and uncertainty quantification for the output from quantum simulators

Abstract

The problem of estimating certain distributions over $\{0,1{\}}^d$ is considered here. The distribution represents a quantum system of $d$ qubits, where there are non-trivial dependencies between the qubits. A maximum entropy approach is adopted to reconstruct the distribution from exact moments or observed empirical moments. The Robbins Monro algorithm is used to solve the intractable maximum entropy problem, by constructing an unbiased estimator of the un-normalized target with a sequential Monte Carlo sampler at each iteration. In the case of empirical moments, this coincides with a maximum likelihood estimator. A Bayesian formulation is also considered in order to quantify posterior uncertainty. Several approaches are proposed in order to tackle this challenging problem, based on recently developed methodologies. In particular, unbiased estimators of the gradient of the log posterior are constructed and used within a provably convergent Langevin-based Markov chain Monte Carlo method. The methods are illustrated on classically simulated output from quantum simulators.

References in corpus (16)

• Quantum Computing in the NISQ era and beyond
• A nonequilibrium equality for free energy differences
• Quantum state tomography via compressed sensing
• Experimental Quantum State Tomography of Optical Fields and Ultrafast Statistical Sampling
• Unknown Quantum States and Operations, a Bayesian View
• The pseudo-marginal approach for efficient Monte Carlo computations
• Quantum state estimation
• Optimal, reliable estimation of quantum states
• MCMC for doubly-intractable distributions
• Adaptive Bayesian Quantum Tomography
• Rejection-free Monte-Carlo sampling for general potentials
• Practical Bayesian Tomography
• On Russian Roulette Estimates for Bayesian Inference with Doubly-Intractable Likelihoods
• Exact sampling from non-attractive distributions using summary states
• Markov Chain Truncation for Doubly-Intractable Inference
• Multilevel Monte Carlo for Scalable Bayesian Computations