qDRIFT and SqDRIFT: Randomized Compilation

• qDRIFT and SqDRIFT are randomized compilation techniques that reduce circuit depth and gate count by stochastically sampling Hamiltonian terms.
• They employ power-series error analysis and Richardson extrapolation to achieve high-order accuracy and exponential improvements over traditional methods.
• SqDRIFT decouples gate cost from the number of Hamiltonian terms, making it effective for dense chemistry Hamiltonians and scalable on near-term quantum devices.

The qDRIFT Randomized Compilation protocol and its higher-order generalizations (“SqDRIFT”, including qFLO) represent a class of randomized circuit compilers for quantum simulation, delivering gate count and circuit depth reductions that are sharply favorable, especially for high-order accuracy and near-term hardware. These algorithms leverage stochastic sampling of Hamiltonian terms, power-series expansions of expectation values, and classical post-processing (notably, Richardson extrapolation), achieving asymptotically exponential improvements over traditional product formulas in simulation error scaling versus circuit depth. They stand out by decoupling gate cost from the number of Hamiltonian summands—a key advantage in chemistry and materials contexts—and are robust, ancilla-free, and compatible with high-noise devices (Watson, 2024).

1. Foundations: qDRIFT Randomized Compilation Protocol

The qDRIFT protocol targets unitary simulation $U=e^{-iHT}$ for time-independent Hamiltonians decomposed as $H = \sum_k h_k H_k$, where $\|H_k\|=1$ and $h_k>0$ (Campbell, 2018, Watson, 2024). The aggregate strength is $\lambda = \sum_k h_k$. Rather than deterministically implementing Trotter steps over all $L$ terms, qDRIFT executes $N$ rounds, each involving stochastic sampling:

• At each round, select index $j$ with probability $p_j = h_j / \lambda$.
• Apply unitary $V_j = \exp(-i\lambda H_j t)$, corresponding to time-step $H = \sum_k h_k H_k$0.

This yields the channel $H = \sum_k h_k H_k$1 per step, with error accumulating linearly over $H = \sum_k h_k H_k$2 steps. The circuit depth required to achieve diamond-norm error $H = \sum_k h_k H_k$3 is $H = \sum_k h_k H_k$4—crucially independent of $H = \sum_k h_k H_k$5 (Chen et al., 2020).

2. Higher-Order SqDRIFT (qFLO): Series Expansion and Richardson Extrapolation

SqDRIFT (also called qFLO) builds upon qDRIFT’s power-series error structure by explicitly extrapolating to higher orders via Richardson techniques (Watson, 2024):

• For step-size $H = \sum_k h_k H_k$6, the $H = \sum_k h_k H_k$7-fold channel is $H = \sum_k h_k H_k$8 with $H = \sum_k h_k H_k$9.
• Observable expectation $\|H_k\|=1$0 expands as $\|H_k\|=1$1.

Richardson extrapolation evaluates $\|H_k\|=1$2 at $\|H_k\|=1$3 tailored points and combines them $\|H_k\|=1$4 to cancel all terms up to $\|H_k\|=1$5, ensuring $\|H_k\|=1$6. Chebyshev node choices for $\|H_k\|=1$7 and corresponding weights $\|H_k\|=1$8 yield controlled norms $\|H_k\|=1$9, and at order $h_k>0$0, circuit depth per run is reduced to $h_k>0$1, exponentially better than the $h_k>0$2 depth for qDRIFT alone.

3. Error Analysis and Circuit Depth Scaling

The extrapolated observable error in SqDRIFT satisfies

$h_k>0$3

Selecting $h_k>0$4 and $h_k>0$5 so the geometric tail falls below $h_k>0$6 provides the depth bound

$h_k>0$7

and total gate count

$h_k>0$8

since $h_k>0$9 repeated runs per point are needed to estimate the sample mean to within $\lambda = \sum_k h_k$0 (Watson, 2024).

4. Comparison to qDRIFT and Deterministic Product Formulas

Method	Circuit Depth (per run)	Total Gate Count	Scaling With Terms $\lambda = \sum_k h_k$1
qDRIFT	$\lambda = \sum_k h_k$2	$\lambda = \sum_k h_k$3	Independent
SqDRIFT/qFLO (qFLO)	$\lambda = \sum_k h_k$4	$\lambda = \sum_k h_k$5	Independent
Trotter/Suzuki	$\lambda = \sum_k h_k$6	$\lambda = \sum_k h_k$7	Explicit dependence

SqDRIFT trades repeated short-depth experiments plus classical post-processing (no ancillas or control gates beyond native $\lambda = \sum_k h_k$8 unitaries) for dramatically reduced coherent depth (Watson, 2024). This reduction is especially pronounced when simulating dense chemistry Hamiltonians with $\lambda = \sum_k h_k$9, where qDRIFT/SqDRIFT’s circuit cost remains solely a function of $L$0 and not $L$1 (Campbell, 2018).

5. Extensions and Hybrid Randomized Compilers

Several generalizations fall within the SqDRIFT paradigm:

• Markov Chain Random Compilation: Extends qDRIFT to compile time-dependent $L$2 using continuous-time Markov chains to stochastically select Hamiltonian terms, allowing dwell times and jump rates to be tuned for optimal error scaling, with overall gate count $L$3 (Dubus et al., 2024).
• Stochastic Hamiltonian Sparsification: Interpolates between pure qDRIFT and randomized Trotter by optimally sparsifying Hamiltonian terms; convex optimization of term inclusion probabilities quadratically suppresses simulation error for given gate budget (Ouyang et al., 2019).
• Importance Sampling and Composite Channels: Sampling from cost-optimized distributions $L$4 (e.g., minimizing expected CNOT or $L$5-count) further suppresses overall gate cost while maintaining rigorous error bounds; composite protocols partition the Hamiltonian into deterministic (Trotter) and stochastic (qDRIFT) blocks (Kiss et al., 2022).
• Adaptive Random Sampling: Fluctuation-guided compilers update term-sampling probabilities in response to real-time state fluctuations, yielding improved error scaling when $L$6 (Wu et al., 12 Sep 2025).

6. Concentration Results and Cost Analysis

Rigorous martingale and concentration-inequality analyses establish that a single realization of the random product formula concentrates sharply around ideal evolution, with probability $L$7 requiring only

$L$8

gates for diamond-norm error $L$9 (Chen et al., 2020). For typical chemistry Hamiltonians, hundreds-to-thousands fold speedups are demonstrated compared to Trotter-Suzuki for relevant time regimes and precision targets, including phase estimation to chemical accuracy (Campbell, 2018). Cost-aware importance sampling further compresses gate counts near per-term hardware cost minima (Kiss et al., 2022).

7. Practical Considerations for Near-Term Quantum Devices

SqDRIFT methods require no ancilla qubits or advanced control operations. Circuit depth at each point is determined by $N$0, $N$1, and $N$2, not $N$3, enabling application to electronic structure problems with massive Hamiltonian decompositions (Watson, 2024). Robustness to imperfect state preparation and measurement noise follows from the well-conditioned structure and incoherent error averaging inherent to randomized compilation. Classical post-processing, including Richardson extrapolation, is computationally light; many short-depth circuits may be executed in parallel to exploit device throughput (Kiss et al., 2022, Watson, 2024). These features collectively render SqDRIFT highly suitable for early fault-tolerant and NISQ quantum hardware, particularly in contexts demanding scalable, chemistry-relevant quantum simulation.