Simulation-Free Fidelity Estimation via Quantum Output Order Statistics

Basic Information

• Paper ID: 2510.13026
• Title: Simulation-Free Fidelity Estimation via Quantum Output Order Statistics
• Author: Tobias Micklitz (Centro Brasileiro de Pesquisas Físicas, Rio de Janeiro, Brazil)
• Classification: quant-ph cond-mat.dis-nn
• Publication Date: October 16, 2025
• Paper Link: https://arxiv.org/abs/2510.13026

Abstract

This paper proposes a simulation-free method for estimating the fidelity of large-scale quantum circuits based on order statistics of measurement output probabilities. The method requires only the highest probability output bitstring (the most frequently observed measurement result) and is built upon exact analytical results for order statistics of Haar-random quantum states derived in this work. By analyzing corrections under depolarizing noise, the authors propose a scalable fidelity estimator, validated on Google's 12-qubit Sycamore experiment and further supported by numerical simulations.

Research Background and Motivation

Core Problem

As universal quantum processors approach the regime where classical simulation becomes infeasible, marking the onset of quantum advantage, verifying the correctness of quantum computations or benchmarking their performance becomes a fundamental challenge, particularly when the ideal output distribution cannot be efficiently simulated.

Problem Significance

1. Quantum Advantage Verification: In quantum advantage experiments, verifying the correctness of quantum computation results is essential
2. Noise Assessment: Understanding and quantifying noise in quantum circuits is crucial for improving quantum computing performance
3. Scalability Requirements: Existing methods face computational complexity challenges in large-scale quantum systems

Limitations of Existing Methods

1. Quantum state tomography and direct fidelity estimation work well for small or structured systems but cannot scale beyond ~20 qubits
2. Cross-entropy benchmarking relies on classical simulation of ideal outputs, becoming computationally prohibitive for circuits larger than ~30 qubits
3. Requires complete output distribution information, which is difficult to obtain experimentally

Research Motivation

Develop simulation-free fidelity estimation techniques that remain practical in the ≳O(20) qubit regime, where classical simulation becomes expensive and direct fidelity estimation is typically infeasible.

Core Contributions

1. Theoretical Contribution: Derives exact analytical expressions for order statistics of Haar-random quantum states, extending previous results limited to k=1
2. Methodological Innovation: Proposes a simulation-free fidelity estimation method based on order statistics, requiring only the highest probability measurement result
3. Noise Modeling: Analyzes corrections to order statistics under depolarizing noise, providing simple rescaling relationships
4. Experimental Validation: Validates the method on Google Sycamore 12-qubit experimental data
5. Scalability: Develops count-based likelihood estimators that overcome finite sampling limitations, supporting larger-scale systems

Detailed Methodology

Task Definition

Given a chaotic quantum circuit on N qubits, estimate the circuit's fidelity f ∈ 0,1 under depolarizing noise by measuring the most frequently occurring output bitstring.

Theoretical Foundation

Order Statistics of Haar-Random States

For a fully chaotic N-qubit quantum state |ψ⟩ = ∑_^D z_k|k⟩ (D = 2^N), the output probabilities are p_k = |z_k|^2. The joint distribution of these probabilities is uniform on the D-dimensional simplex:

$P_{CUE}(p_1, p_2, ..., p_D) = (D-1)!\delta(p_1 + p_2 + ... + p_D - 1)$

The probability distribution for the k-th largest probability taking value x is:

$P_k(x,N) = \mathcal{N} \sum_{j=k}^{j_{max}} \binom{D-k}{j-k} (-1)^j (1-jx)^{D-2}$

where j_ = min(D, ⌊1/x⌋), which is the core technical result of this work.

Noise Model

Under depolarizing noise, each ideal output probability undergoes an affine shift toward the uniform value 1/D:

$p_k(f) = fp_k + (1-f)/D$

The effect of noise on order statistics is captured through simple rescaling:

$P_k(x;N,f) = P_k(x_f, N), \quad fx_f = x - (1-f)/D$

Likelihood Function Construction

Probability-Based Likelihood

For M circuit implementations, each with K maximum measurement probabilities {p^m_k}, the likelihood function is:

$\ln\Lambda(f;N,\{p^m_k\}) = \sum_{m=1}^M \sum_{k \in K^*} \ln P_k(p^m_k;N,f)$

Count-Based Likelihood

To overcome finite sampling limitations, a Poisson likelihood is constructed using raw measurement counts n_k:

$\ln\Lambda(f;N,\{n_k\}) = \sum_{k \in K^*} (n_k \ln p_k(f) - Sp_k(f))$

where S is the total number of measurements.

Technical Innovations

1. Order Statistics Decomposition: Decomposes the Porter-Thomas distribution into ordered statistical components, revealing transitions from extremal (Gumbel-like) to median (Gaussian) to tail (exponential) regimes
2. Noise Invariance: Discovers simple rescaling relationships under depolarizing noise that preserve ordering
3. Multi-Scale Estimation: Supports averaging strategies across fixed ranks over circuits and multiple ranks within single circuits
4. Adaptive Sampling: Count-based estimators remain robust even under severe sampling limitations

Experimental Setup

Datasets

1. Google Sycamore 12-Qubit Data: 20 circuit implementations, each with 5×10^5 measurements
2. Numerical Simulations: Chaotic circuits with N = 20-42 qubits under depolarizing noise with known fidelity

Evaluation Metrics

• Relative Error: |f_ - f_|/f_
• Likelihood Peak Position: Fidelity value at maximum likelihood estimate
• Required Sampling: Minimum number of measurements to achieve given precision threshold

Comparison Methods

• Cross-entropy benchmarking (as reference standard)
• Comparison of different rank selection strategies

Implementation Details

• K = 500 ranked results retained to maintain moderate computational cost
• Low-rank large-D approximations to avoid numerical instability
• Bounded scalar optimization to solve for f ∈ 0,1

Experimental Results

Main Results

Sycamore 12-Qubit Verification

• Fidelity Estimation Range: f ≃ 0.45-0.50, consistent with cross-entropy benchmarking
• Cross-Circuit Averaging: Fixed-rank fidelity estimates fall within narrow range f ≈ 0.47-0.49
• Single-Circuit Analysis: Individual circuit fidelities vary f ≈ 0.44-0.52, reflecting inter-circuit noise variations

Scalability Verification

• Error Scaling: Estimation error scales as 1/N with system size, matching theoretical predictions
• Likelihood Width: Scales as ~f²/(√kM), becoming sharper with system size
• Rank Robustness: Sparse rank sampling deviates <0.01 compared to full rank case

Finite Sampling Performance

Sampling Requirements

Traditional methods require sampling to scale as: $S \gtrsim \frac{2^N}{\varepsilon_{rel}^2 N}$

This method achieves reliable fidelity estimation well below this threshold.

Improvement Factor

Count-based estimators improve sampling requirements by factor f/K, where K is the number of retained ranks.

Ablation Studies

• Rank Set Comparison: k = {1,2,3,5,6} vs k = {3,4,5,6}, no significant difference in estimation accuracy
• Rank Number Impact: Doubling rank numbers reduces required sampling by approximately half
• Noise Level Dependence: Performance validation at different fidelities f = 0.1, 0.3, 0.5

Related Work

Quantum Fidelity Estimation

• Quantum State Tomography: Complete reconstruction but not scalable
• Direct Fidelity Estimation: Efficient but limited to small systems
• Cross-Entropy Benchmarking: Requires classical simulation, computationally expensive

Quantum Randomness Statistics

• Porter-Thomas Distribution: Universal characteristic of chaotic quantum systems
• Quantum Speckle Patterns: Unique fingerprints of interference structures
• Extreme Value Statistics: Applications in quantum systems

Noise Modeling

• Depolarizing Channels: Primary error source in large-scale circuits
• Random Compilation: Averaged error behavior of Pauli twirling
• Coherent vs Incoherent Errors: Effects of different noise types

Conclusions and Discussion

Main Conclusions

1. Theoretical Breakthrough: First complete derivation of analytical expressions for order statistics of Haar-random states
2. Practical Method: Provides scalable simulation-free fidelity estimation tool
3. Experimental Validation: Demonstrates method effectiveness on real quantum experimental data
4. Technical Advantages: Overcomes sampling and computational limitations of traditional methods

Limitations

1. Noise Model Assumptions: Limited to uniform depolarizing noise; real hardware may have coherent or correlated error components
2. Chaotic Circuit Requirement: Method specifically designed for highly entangled chaotic states, may not apply to structured circuits
3. Finite Precision: Final accuracy still limited by finite sampling, particularly as N increases
4. Rank Correlations: Correlations between ranks within single implementations may introduce likelihood bias

Future Directions

1. Higher-Order Statistics: Explore higher-order statistics to improve estimation accuracy
2. Bayesian Inference: More advanced inference techniques
3. Rank Optimization: Optimize rank selection strategies
4. Extended Noise Models: Handle more complex noise structures
5. Other Diagnostic Applications: Apply order statistics to other quantum diagnostic tasks

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Provides complete mathematical derivations and exact analytical results
2. Methodological Innovation: First systematic application of order statistics to quantum fidelity estimation
3. Practical Value: Addresses actual verification needs in quantum advantage regime
4. Comprehensive Validation: Supported by both real experimental data and large-scale numerical simulations
5. Strong Scalability: Method performs better at larger system scales

Weaknesses

1. Limited Applicability: Primarily applicable to chaotic quantum circuits; applicability to structured circuits unknown
2. Simplified Noise Model: Depolarizing noise assumption may be oversimplified
3. Rough Correlation Handling: Treatment of inter-rank correlations is relatively crude
4. Limited Experimental Scale: Sycamore validation restricted to 12 qubits; larger-scale experimental validation lacking

Impact

1. Academic Contribution: Provides new theoretical tools and methods for quantum computation verification
2. Practical Value: Direct application to quantum advantage experiments and NISQ device benchmarking
3. Reproducibility: Detailed method description and complete theoretical derivations facilitate reproduction and extension
4. Interdisciplinary Significance: Connects quantum physics, statistics, and computational science

Applicable Scenarios

1. Quantum Advantage Verification: Performance verification of large-scale quantum circuits
2. NISQ Device Benchmarking: Fidelity assessment of medium-scale noisy quantum devices
3. Quantum Algorithm Development: Performance monitoring of chaotic quantum algorithms
4. Quantum Hardware Optimization: Guiding hardware improvements through fidelity feedback

References

This paper cites 33 important references spanning quantum computing, random matrix theory, statistical physics and other fields, particularly:

• Original papers on Google's quantum advantage experiments
• Classical literature on Porter-Thomas distribution and quantum chaos
• State-of-the-art methods in quantum fidelity estimation
• Mathematical foundations of extreme value statistics

Overall Assessment: This is a high-quality paper combining theory and experiment, proposing an innovative quantum fidelity estimation method with significant theoretical value and practical importance. The method is rigorous, validation is comprehensive, and it makes important contributions to the quantum computation verification field.