Bayesian quantum phase estimation beyond the eigenstate approximation: sample complexity, Cramér–Rao bounds, and adaptive measurement design for non-eigenstate and mixed-state inputs on NISQ devices

Background Quantum phase estimation (QPE) is a foundational subroutine in quantum algorithms that ranges from Shor’s factoring scheme to variational quantum eigensolvers and quantum metrology protocols. Classical post-processing replacements for the inverse quantum Fourier transform (iQFT) —specifically Kitaev-style iterative estimation and Bayesian update schemes—substantially reduce circuit depth, making them practical for noisy intermediate-scale quantum (NISQ) hardware. Despite this promise, the statistical performance of these methods under the conditions that actually arise in practice—non-eigenstate inputs characterised by eigenstate overlap p ≡ O = | ψ | ϕ | 2 < 1 and mixed-state purity γ ≡ P = Tr ( ρ 2 ) < 1 —remains uncharacterised in the literature. Methods We develop a rigorous information-theoretic framework for Bayesian QPE (BQPE) beyond the eigenstate assumption. Starting from the binary measurement likelihood attenuated by the joint factor O P , we derive exact classical Fisher information expressions and apply the Cramér–Rao inequality to obtain tight variance lower bounds. We then construct an adaptive measurement protocol based on von Mises posterior updates, prove its asymptotic efficiency in the frequentist sense, and establish matching upper and lower bounds on sample complexity. All theoretical results are validated against 2,000 -trial Monte Carlo simulations and cross-checked against published nitrogen-vacancy (NV) centre and photonic experimental data. Results We prove that the fundamental variance floor for any unbiased frequentist estimator is ( 1 − p γ ) / ( N p γ n ̄ 2 ) 1 / ( N ( p γ ) 2 n ̄ 2 ) , where N