Estimating the energy spectra of quantum many-body systems is a fundamental task in quantum physics, with applications ranging from chemistry to condensed matter. Algorithmic shadow spectroscopy is a recent method that leverages randomized measurements on time-evolved quantum states to extract spectral information. However, implementing accurate time evolution with low-depth circuits remains a key challenge for near-term quantum hardware. In this work, we propose a hybrid quantum-classical protocol that integrates Time Evolution via Probabilistic Angle Interpolation (TE-PAI) into the shadow spectroscopy framework. TE-PAI enables the simulation of time evolution using shallow stochastic circuits while preserving unbiased estimates through quasiprobability sampling. We construct the combined estimator and derive its theoretical properties. Through numerical simulations, we demonstrate that our method accurately resolves energy gaps and exhibits enhanced robustness to gate noise compared to standard Trotter-based shadow spectroscopy. We further validate the protocol experimentally on up to 20 qubits using IBM quantum hardware. This makes TE-PAI shadow spectroscopy a promising tool for spectral analysis on noisy intermediate-scale quantum (NISQ) devices.
Estimating the energy spectra of quantum many-body systems is a fundamental task in quantum physics, with applications ranging from chemistry to condensed matter. Algorithmic shadow spectroscopy is a recent method that leverages randomized measurements on time-evolved quantum states to extract spectral information. However, implementing accurate time evolution with low-depth circuits remains a key challenge for near-term quantum hardware. In this work, we propose a hybrid quantum-classical protocol that integrates Time Evolution via Probabilistic Angle Interpolation (TE-PAI) into the shadow spectroscopy framework. TE-PAI enables the simulation of time evolution using shallow stochastic circuits while preserving unbiased estimates through quasiprobability sampling. We construct the combined estimator and derive its theoretical properties. Through numerical simulations, we demonstrate that our method accurately resolves energy gaps and exhibits enhanced robustness to gate noise compared to standard Trotter-based shadow spectroscopy. We further validate the protocol experimentally on up to 20 qubits using IBM quantum hardware. This makes TE-PAI shadow spectroscopy a promising tool for spectral analysis on noisy intermediate-scale quantum (NISQ) devices.
Extracting energy spectra from quantum Hamiltonians is a fundamental task for quantum simulation, yet remains challenging on noisy intermediate-scale quantum (NISQ) devices. We propose Measured Quantum Time Evolution (MQTE), an ancilla-free algorithm that estimates energy gaps by applying real-time evolution to a reference state and measuring time-resolved probabilities via repeated projective measurements. Spectral analysis of these signals reveals oscillation frequencies corresponding to eigenvalue differences. Crucially, MQTE exhibits inherent robustness to quantum hardware noise and sampling errors: these disturbances manifest as a white-noise background, which does not distort the underlying spectral structure but rather obscures the frequency information. By increasing the number of measurement samples, the intensity of the background white noise can be suppressed, thereby recovering the original spectral content. We validate the algorithm's performance via numerical simulations on one- and two-dimensional Heisenberg models, demonstrating accurate extraction of energy gaps and resilience against both sampling and circuit-level noise. Experimental implementation on the superconducting quantum processor Tianyan-176-II further confirms the practical feasibility and noise tolerance of MQTE under real hardware conditions. This work provides a robust and scalable framework for quantum spectral estimation in the NISQ era.
Abstract Simulating time evolution under quantum Hamiltonians is one of the most natural applications of quantum computers. We introduce TE-PAI, which simulates time evolution exactly by sampling random quantum circuits for the purpose of estimating observable expectation values at the cost of an increased circuit repetition. The approach builds on the Probabilistic Angle Interpolation (PAI) technique and we prove that it simulates time evolution without discretisation or algorithmic error while achieving shallow circuit depths with optimal scaling that saturates the Lieb–Robinson bound. Another significant advantage of TE-PAI is that it only requires executing random circuits that consist of Pauli rotation gates of only two kinds of rotation angles ± Δ and π , along with measurements. While TE-PAI is highly beneficial for NISQ devices, we additionally develop an optimised early fault-tolerant implementation using catalyst circuits and repeat-until-success teleportation, concluding that the approach requires orders of magnitude fewer T-states than conventional techniques, such as Trotterization—we estimate 3 × 10 5 T states are sufficient for the fault-tolerant simulation of a 100-qubit Heisenberg spin Hamiltonian. Furthermore, TE-PAI allows for a highly configurable trade-off between circuit depth and measurement overhead by adjusting the rotation angle Δ arbitrarily. We expect that the approach will be a major enabler in the late NISQ and early fault-tolerant periods as it can compensate circuit-depth and qubit-number limitations through an increased circuit repetition.
We present shadow spectroscopy as a simulator-agnostic quantum algorithm for estimating energy gaps using very few circuit repetitions (shots) and no extra resources (ancilla qubits) beyond performing time evolution and measurements. The approach builds on the fundamental feature that every observable property of a quantum system must evolve according to the same harmonic components: we can reveal them by postprocessing classical shadows of time-evolved quantum states to extract a large number of time-periodic signals N o ∝ 10 8 , whose frequencies correspond to Hamiltonian energy differences with precision limited as ϵ ∝ 1 / T for simulation time T . We provide strong analytical guarantees that (a) quantum resources scale as O ( log N o ) , while the classical computational complexity is linear O ( N o ) , (b) the signal-to-noise ratio increases with the number of processed signals as ∝ N o , and (c) spectral peak positions are immune to reasonable levels of noise. We demonstrate our approach on model spin systems and the excited-state conical intersection of molecular CH 2 and verify that our method is indeed intuitively easy to use in practice, robust against gate noise, amiable to a new type of algorithmic-error mitigation technique, and uses relatively few shots given a reasonable initial state is supplied—we demonstrate that even 10 shots per time step can be sufficient. Finally, we measured a high-quality, experimental shadow spectrum of a spin chain on readily available IBM quantum computers, achieving the same precision as in noise-free simulations without using any advanced error mitigation, and verified scalability in tensor-network simulations of up to 100-qubit systems.
Accurately estimating observables on noisy quantum devices remains a central challenge for near-term quantum algorithms. While quantum error mitigation techniques can reduce noise-induced bias, they often rely on unverifiable assumptions about the circuit noise, and cannot guarantee the magnitude of residual bias error. Here, rather than using classical resources to mitigate a noisy quantum circuit execution, we propose a hybrid algorithm that uses quantum resources to improve the accuracy of approximate classical Pauli-path simulation. Our protocol, Quantum Enhanced Pauli Propagation (QuEPP), uses Clifford perturbation theory (CPT) to construct a classically simulable ensemble of Clifford circuits from the low-order terms in CPT, which directly provide the approximate classical Pauli-path simulation of the target circuit. Noisy quantum expectation values of this ensemble are then used to infer a global rescaling factor that corrects quantum execution of the target circuit, providing higher-order contributions absent from the truncated low-order classical simulation. This approach requires no noise characterization, applies to arbitrary circuits, and provides a provable route to asymptotically unbiased estimates. Using IBM Heron hardware, we demonstrate QuEPP on 2D random mirror circuits of up to 49 qubits and circuit depth 80, as well as Trotterized Hamiltonian evolution, showing consistent improvements beyond classical CPT and unmitigated quantum results. QuEPP offers a simple, scalable, and model-free framework for enabling accurate quantum computation in the pre-fault-tolerant era.
Extracting information efficiently from quantum systems is crucial for quantum information processing. Classical shadows enable predicting many properties of arbitrary quantum states using few measurements. While random single-qubit measurements are experimentally friendly and suitable for learning low-weight Pauli observables, they perform poorly for nonlocal observables. Introducing a shallow random quantum circuit before measurements improves sample efficiency for high-weight Pauli observables and low-rank properties. However, in practice, these circuits can be noisy and bias the measurement results. Here, we propose the robust shallow shadows, which employs Bayesian inference to learn and mitigate noise in postprocessing. We analyze noise effects on sample complexity and the optimal circuit depth. We provide theoretical guarantees for the success of error mitigation under a wide class of noise processes. Experimental validation on a superconducting quantum processor confirms the advantage of our method, even in the presence of realistic noise, over single-qubit measurements for predicting diverse state properties, such as fidelity and entanglement entropy. Our protocol thus offers a scalable, robust, and sample-efficient method for quantum state characterization on near-term quantum devices. Fast and reliable characterisation of quantum systems is a key part of quantum technologies development. Here, the authors propose and demonstrate a way to embed noise characterisation in classical shadow estimation of nonlocal properties, enabling an efficient way to extract information from noisy quantum systems.
Abstract Quantum phase estimation (QPE) is a fundamental tool in quantum computing, facilitating efficient simulations of complex problems in quantum chemistry and materials science. While most phase estimation algorithms are deterministic, recent advancements indicate that incorporating randomness can enhance performance. This study introduces a framework for randomized QPE that merges the benefits of randomized compilation with phase estimation algorithms based on quantum signal processing. Our proposed algorithms effectively reduce circuit depths by eliminating the need for precise Hamiltonian time evolution, making them advantageous for digital quantum computers estimating the eigenvalue and eigenvector properties of Hamiltonians. Notably, our findings show that the quantum stochastic drift protocol (qDRIFT)-based randomized algorithm surpasses the original phase estimation with qDRIFT, especially in scaling inverse failure probabilities. We also establish that a circuit depth of O ( log ( M ) ) suffices for estimating M distinct observables. The protocol is executed through multiple iterations of the randomized algorithms combined with classical shadow techniques. Overall, our framework retains many advantages of the randomized compilation technique, making it a compelling solution for challenges in quantum chemistry.