Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2603.28602

Lindbladian Simulation with Commutator Bounds

Xinzhao Wang, Shuo Zhou, Xiaoyang Wang, Yi-Cong Zheng, Shengyu Zhang, Tongyang Li

quant-phcs.DSSubmitted Mar 30, 2026• Announced Mar 31, 2026• 0 citations• Revised Aug 13, 2026 (v2)
v1 · Mar 30, 2026v2 · Aug 13, 2026
arXivHTMLPDF

Abstract

Trotter decomposition provides a simple approach to simulating open quantum systems by decomposing the Lindbladian into a sum of individual terms. While it is established that Trotter errors in Hamiltonian simulation depend on nested commutators of the summands, such a relationship remains poorly understood for Lindbladian dynamics. In this Letter, we derive commutator-based Trotter error bounds for Lindbladian simulation, yielding an $O(\sqrt{N})$ scaling in the number of Trotter steps for locally interacting systems on $N$ sites. When estimating observable averages, we apply Richardson extrapolation to achieve polylogarithmic precision while maintaining the commutator scaling. To bound the extrapolation remainder, we develop a general truncation bound for the Baker-Campbell-Hausdorff expansion that bypasses common convergence issues in physically relevant systems. For local Lindbladians, our results demonstrate that the Trotter-based methods outperform prior simulation techniques in system-size scaling while requiring only $O(1)$ ancillas. Numerical simulations further validate the predicted system-size and precision scaling.

Comments
28 pages, 5 figures

Related papers

Practical Estimation of Trotter Error for Hamiltonian Simulation

William Maxwell, Pablo A. M. Casares, Robert A. Lang +7 more

Trotter product formulas are a leading approach for Hamiltonian simulation on quantum computers, yet their practical performance has remained difficult to assess due to the challenge of accurately estimating the Trotter error. In this work, we develop new theoretical results, algorithms, and software tools that advance the state-of-the-art in Trotter error estimation by orders of magnitude in both scale and accuracy. On the theoretical side, we prove that in the asymptotic limit the error of a product formula depends on the diagonal elements of the Baker-Campbell-Hausdorff (BCH) error operator in the eigenbasis of the Hamiltonian, rather than its full spectral norm -- yielding an improved scaling for Hamiltonian simulation using product formulas. On the algorithmic side, we introduce a compact representation of the BCH expansion that reduces the number of commutators from $\mathcal{O}(n^3)$ to $\mathcal{O}(n)$ for second-order, and from $\mathcal{O}(n^5)$ to $\mathcal{O}(n^2)$ for fourth-order formulas on $n$ fragments, complemented by an importance sampling scheme to further reduce the computational cost. We provide implementations of these techniques in software and demonstrate their power on two applications: (i) X-ray absorption spectroscopy of an electronic Hamiltonian (Li$_4$Mn$_2$O) at up to 56 qubits using tensor networks; and (ii) vibronic dynamics of naphthalene at over 100 qubits using ML-MCTDH, where we find that naive analytical bounds overestimate the required number of Trotter steps by nearly five orders of magnitude. Our framework enables, for the first time, the accurate estimation of Trotter error at practically relevant system sizes, providing a foundation for fair algorithmic comparisons and rational design of product formulas.

quant-pharXiv:2606.30738•Jul 1, 2026•0 citations

Trotter error compensation with polylogarithmic precision and nested-commutator scaling without ancillas

Xinzhao Wang, Shuo Zhou, Ziruo Wang +5 more

Product formulas are among the most practical approaches to Hamiltonian simulation, requiring no ancillary qubits and exhibiting error bounds governed by nested commutators rather than only by Hamiltonian norms. Their circuit size, however, scales polynomially with the inverse precision. We develop a high-order nested-commutator compensation (HNCC) algorithm that preserves the main advantages of product formulas while achieving polylogarithmic precision dependence in the circuit size and the standard $\mathcal{O}(\varepsilon^{-2})$ sampling cost. HNCC uses a truncated Baker--Campbell--Hausdorff expansion to represent high-order Trotter errors by products of nested commutators and compensates these errors at the channel level through randomly sampled Pauli-rotation channels, avoiding Hadamard tests and ancillary qubits. For a fixed $K$-th order product formula applied to a $k$-local Hamiltonian on $N$ qubits with $Γ$ Pauli terms and local interaction strength $g_0$, HNCC estimates $\operatorname{tr}[Oe^{-i tH}ρe^{i tH}]$ to additive precision $\varepsilon\|O\|$ using $\mathcal{O}(\varepsilon^{-2})$ repetitions. Its maximum gate count per circuit is $\mathcal{O}\bigl( kN^{\frac{1}{2K+1}} Γ^{1-\frac{1}{2K+1}} \max\{Γ,N\log(1/\varepsilon)\}^{\frac{1}{2K+1}} (kg_0t\log(1/\varepsilon))^{1+\frac{1}{2K+1}} \bigr)$. Finite-size resource estimates for the periodic Heisenberg chain indicate that HNCC has the lowest estimated $T$-gate count per circuit among the product-formula-based methods considered.

quant-pharXiv:2607.11856•Jul 14, 2026•0 citations

Entanglement-Dependent Error Bounds for Hamiltonian Simulation

Prateek P. Kulkarni

We establish tight connections between entanglement entropy and the approximation error in Trotter-Suzuki product formulas for Hamiltonian simulation. Product formulas remain the workhorse of quantum simulation on near-term devices, yet standard error analyses yield worst-case bounds that can vastly overestimate the resources required for structured problems. For systems governed by geometrically local Hamiltonians with maximum entanglement entropy $S_\text{max}$ across all bipartitions, we prove that the first-order Trotter error scales as $\mathcal{O}(t^2 S_\text{max} \operatorname{polylog}(n)/r)$ rather than the worst-case $\mathcal{O}(t^2 n/r)$, where $n$ is the system size and $r$ is the number of Trotter steps. This yields improvements of $\tildeΩ(n^2)$ for one-dimensional area-law systems and $\tildeΩ(n^{3/2})$ for two-dimensional systems. We extend these bounds to higher-order Suzuki formulas, where the improvement factor involves $2^{pS^*/2}$ for the $p$-th order formula. We further establish a separation result demonstrating that volume-law entangled systems fundamentally require $\tildeΩ(n)$ more Trotter steps than area-law systems to achieve the same precision. This separation is tight up to logarithmic factors. Our analysis combines Lieb-Robinson bounds for locality, tensor network representations for entanglement structure, and novel commutator-entropy inequalities that bound the expectation value of nested commutators by the Schmidt rank of the state. These results have immediate applications to quantum chemistry, condensed matter simulation, and resource estimation for fault-tolerant quantum computing.

quant-pharXiv:2602.00555•Feb 3, 2026•0 citations

Trotter error timescaling separation via commutant decomposition

Yi‐Hsiang Chen

Suppressing the Trotter error in dynamical quantum simulation typically requires running deeper circuits, posing a great challenge for noisy near-term quantum devices. Studies have shown that the empirical error is usually much smaller than the one suggested by existing bounds, implying the actual circuit cost required is much less than the ones based on those bounds. Here, we improve the estimate of the Trotter error over existing bounds, by introducing a general framework of commutant decomposition that separates disjoint error components that have fundamentally different scaling with time. In particular we identify two error components that each scale as $O({\ensuremath{\tau}}^{p}t)$ and $O({\ensuremath{\tau}}^{p})$ for a $p\mathrm{th}$-order product formula evolving to time $t$ using a fixed step size $\ensuremath{\tau}$, implying one would scale linearly with time $t$ and the other would be a constant of $t$. We show that this formalism not only straightforwardly reproduces previous results but also provides a better error estimate for higher-order product formulas. We demonstrate the improvement both analytically and numerically. We also apply the analysis to observable error relating to the heating in Floquet dynamics and thermalization, which is of independent interest.

arXiv:2409.16634•Physical Review A•Feb 11, 2025•1 citation

Joint Mitigation of Algorithmic and Physical Errors in Noisy Hamiltonian Simulation

Shuo Zhou, Xinzhao Wang, Ruiqi Zhang +4 more

Product-formula Hamiltonian simulation is naturally suited to near-term quantum processors, but its accuracy is set by two competing errors: finite-step Trotter bias and physical hardware noise. We introduce a joint extrapolation strategy that ties the tunable per-layer noise strength to the Trotter step size, $λ(s)=c(sT)^{p+1}$ for a $p$-th order product formula. Along this one-dimensional path, the leading physical-noise and Trotter corrections over the full evolution both enter at order $s^p$ and can be canceled by a single Richardson extrapolation. Building on a previously established finite-order Baker--Campbell--Hausdorff truncation bound, we derive a provably commutator-scaling resource guarantee for the joint extrapolation. For local Hamiltonians with local Lindbladian noise, the protocol mitigates physical noise together with Trotter error with only a constant asymptotic overhead relative to noiseless Trotter extrapolation, provided the required noise strengths lie above the intrinsic device-noise floor. We experimentally demonstrate the protocol for Ising dynamics on a superconducting quantum computer using learned noise amplification. Complementary 100-qubit Sparse Pauli Dynamics (SPD) simulations achieve comparable accuracy to two-dimensional Richardson extrapolation with fewer circuit settings.

quant-pharXiv:2609.11508•Sep 11, 2026•0 citations

On the commutator scaling in Hamiltonian simulation with multi-product formulas

Kaoru Mizuta

A multi-product formula (MPF) is a promising approach for Hamiltonian simulation efficiently both in the system size N and the inverse allowable error 1 / ε by combining Trotterization and the linear combination of unitaries (LCU). It achieves poly-logarithmic cost in 1 / ε like LCU [G. H. Low, V. Kliuchnikov, N. Wiebe, (2019)]. The efficiency in N is expected to come from the commutator scaling in Trotterization, and this appears to be confirmed by the error bound of MPF expressed by nested commutators [J. Aftab, D. An, K. Trivisa, (2024)]. However, we point out that the efficiency of MPF in the system size N is not exactly resolved yet in that the present error bound expressed by nested commutators is incompatible with the size-efficient complexity reflecting the commutator scaling. The problem is that q -fold nested commutators with arbitrarily large q are involved in their requirement and error bound. The benefit of commutator scaling by locality is absent, and the cost efficient in N becomes prohibited in general. In this paper, we show an alternative commutator-scaling error of MPF and derive its size-efficient cost properly inheriting the advantage in Trotterization. The requirement and the error bound in our analysis, derived by techniques from the Floquet-Magnus expansion, have a certain truncation order in the nested commutators and can fully exploit the locality. We prove that Hamiltonian simulation by MPF certainly achieves the cost whose system-size dependence is as large as Trotterization while keeping the p o l y l o g ( 1 / ε ) -scaling like the LCU. Our results will provide improved or accurate error and cost also for various algorithms using interpolation or extrapolation of Trotterization.

arXiv:2507.06557•Quantum•Jan 19, 2026•2 citations
Trotter decomposition provides a simple approach to simulating open quantum systems by decomposing the Lindbladian into a sum of individual terms. While it is established that Trotter errors in Hamiltonian simulation depend on nested commutators of the summands, such a relationship remains poorly understood for Lindbladian dynamics. In this Letter, we derive commutator-based Trotter error bounds for Lindbladian simulation, yielding an O(N)O(\sqrt{N})O(N​) scaling in the number of Trotter steps for locally interacting systems on NNN sites. When estimating observable averages, we apply Richardson extrapolation to achieve polylogarithmic precision while maintaining the commutator scaling. To bound the extrapolation remainder, we develop a general truncation bound for the Baker-Campbell-Hausdorff expansion that bypasses common convergence issues in physically relevant systems. For local Lindbladians, our results demonstrate that the Trotter-based methods outperform prior simulation techniques in system-size scaling while requiring only O(1)O(1)O(1) ancillas. Numerical simulations further validate the predicted system-size and precision scaling.
Trotter product formulas are a leading approach for Hamiltonian simulation on quantum computers, yet their practical performance has remained difficult to assess due to the challenge of accurately estimating the Trotter error. In this work, we develop new theoretical results, algorithms, and software tools that advance the state-of-the-art in Trotter error estimation by orders of magnitude in both scale and accuracy. On the theoretical side, we prove that in the asymptotic limit the error of a product formula depends on the diagonal elements of the Baker-Campbell-Hausdorff (BCH) error operator in the eigenbasis of the Hamiltonian, rather than its full spectral norm -- yielding an improved scaling for Hamiltonian simulation using product formulas. On the algorithmic side, we introduce a compact representation of the BCH expansion that reduces the number of commutators from O(n3)\mathcal{O}(n^3)O(n3) to O(n)\mathcal{O}(n)O(n) for second-order, and from O(n5)\mathcal{O}(n^5)O(n5) to O(n2)\mathcal{O}(n^2)O(n2) for fourth-order formulas on nnn fragments, complemented by an importance sampling scheme to further reduce the computational cost. We provide implementations of these techniques in software and demonstrate their power on two applications: (i) X-ray absorption spectroscopy of an electronic Hamiltonian (Li4_44​Mn2_22​O) at up to 56 qubits using tensor networks; and (ii) vibronic dynamics of naphthalene at over 100 qubits using ML-MCTDH, where we find that naive analytical bounds overestimate the required number of Trotter steps by nearly five orders of magnitude. Our framework enables, for the first time, the accurate estimation of Trotter error at practically relevant system sizes, providing a foundation for fair algorithmic comparisons and rational design of product formulas.
Product formulas are among the most practical approaches to Hamiltonian simulation, requiring no ancillary qubits and exhibiting error bounds governed by nested commutators rather than only by Hamiltonian norms. Their circuit size, however, scales polynomially with the inverse precision. We develop a high-order nested-commutator compensation (HNCC) algorithm that preserves the main advantages of product formulas while achieving polylogarithmic precision dependence in the circuit size and the standard O(ε−2)\mathcal{O}(\varepsilon^{-2})O(ε−2) sampling cost. HNCC uses a truncated Baker--Campbell--Hausdorff expansion to represent high-order Trotter errors by products of nested commutators and compensates these errors at the channel level through randomly sampled Pauli-rotation channels, avoiding Hadamard tests and ancillary qubits. For a fixed KKK-th order product formula applied to a kkk-local Hamiltonian on NNN qubits with ΓΓΓ Pauli terms and local interaction strength g0g_0g0​, HNCC estimates tr⁡[Oe−itHρeitH]\operatorname{tr}[Oe^{-i tH}ρe^{i tH}]tr[Oe−itHρeitH] to additive precision ε∥O∥\varepsilon\|O\|ε∥O∥ using O(ε−2)\mathcal{O}(\varepsilon^{-2})O(ε−2) repetitions. Its maximum gate count per circuit is O(kN12K+1Γ1−12K+1max⁡{Γ,Nlog⁡(1/ε)}12K+1(kg0tlog⁡(1/ε))1+12K+1)\mathcal{O}\bigl( kN^{\frac{1}{2K+1}} Γ^{1-\frac{1}{2K+1}} \max\{Γ,N\log(1/\varepsilon)\}^{\frac{1}{2K+1}} (kg_0t\log(1/\varepsilon))^{1+\frac{1}{2K+1}} \bigr)O(kN2K+11​Γ1−2K+11​max{Γ,Nlog(1/ε)}2K+11​(kg0​tlog(1/ε))1+2K+11​). Finite-size resource estimates for the periodic Heisenberg chain indicate that HNCC has the lowest estimated TTT-gate count per circuit among the product-formula-based methods considered.
We establish tight connections between entanglement entropy and the approximation error in Trotter-Suzuki product formulas for Hamiltonian simulation. Product formulas remain the workhorse of quantum simulation on near-term devices, yet standard error analyses yield worst-case bounds that can vastly overestimate the resources required for structured problems. For systems governed by geometrically local Hamiltonians with maximum entanglement entropy SmaxS_\text{max}Smax​ across all bipartitions, we prove that the first-order Trotter error scales as O(t2Smaxpolylog⁡(n)/r)\mathcal{O}(t^2 S_\text{max} \operatorname{polylog}(n)/r)O(t2Smax​polylog(n)/r) rather than the worst-case O(t2n/r)\mathcal{O}(t^2 n/r)O(t2n/r), where nnn is the system size and rrr is the number of Trotter steps. This yields improvements of Ω~(n2)\tildeΩ(n^2)Ω~(n2) for one-dimensional area-law systems and Ω~(n3/2)\tildeΩ(n^{3/2})Ω~(n3/2) for two-dimensional systems. We extend these bounds to higher-order Suzuki formulas, where the improvement factor involves 2pS∗/22^{pS^*/2}2pS∗/2 for the ppp-th order formula. We further establish a separation result demonstrating that volume-law entangled systems fundamentally require Ω~(n)\tildeΩ(n)Ω~(n) more Trotter steps than area-law systems to achieve the same precision. This separation is tight up to logarithmic factors. Our analysis combines Lieb-Robinson bounds for locality, tensor network representations for entanglement structure, and novel commutator-entropy inequalities that bound the expectation value of nested commutators by the Schmidt rank of the state. These results have immediate applications to quantum chemistry, condensed matter simulation, and resource estimation for fault-tolerant quantum computing.
Suppressing the Trotter error in dynamical quantum simulation typically requires running deeper circuits, posing a great challenge for noisy near-term quantum devices. Studies have shown that the empirical error is usually much smaller than the one suggested by existing bounds, implying the actual circuit cost required is much less than the ones based on those bounds. Here, we improve the estimate of the Trotter error over existing bounds, by introducing a general framework of commutant decomposition that separates disjoint error components that have fundamentally different scaling with time. In particular we identify two error components that each scale as $O({\ensuremath{\tau}}^{p}t)$ and $O({\ensuremath{\tau}}^{p})$ for a pthp\mathrm{th}pth-order product formula evolving to time ttt using a fixed step size $\ensuremath{\tau}$, implying one would scale linearly with time ttt and the other would be a constant of ttt. We show that this formalism not only straightforwardly reproduces previous results but also provides a better error estimate for higher-order product formulas. We demonstrate the improvement both analytically and numerically. We also apply the analysis to observable error relating to the heating in Floquet dynamics and thermalization, which is of independent interest.
Product-formula Hamiltonian simulation is naturally suited to near-term quantum processors, but its accuracy is set by two competing errors: finite-step Trotter bias and physical hardware noise. We introduce a joint extrapolation strategy that ties the tunable per-layer noise strength to the Trotter step size, λ(s)=c(sT)p+1λ(s)=c(sT)^{p+1}λ(s)=c(sT)p+1 for a ppp-th order product formula. Along this one-dimensional path, the leading physical-noise and Trotter corrections over the full evolution both enter at order sps^psp and can be canceled by a single Richardson extrapolation. Building on a previously established finite-order Baker--Campbell--Hausdorff truncation bound, we derive a provably commutator-scaling resource guarantee for the joint extrapolation. For local Hamiltonians with local Lindbladian noise, the protocol mitigates physical noise together with Trotter error with only a constant asymptotic overhead relative to noiseless Trotter extrapolation, provided the required noise strengths lie above the intrinsic device-noise floor. We experimentally demonstrate the protocol for Ising dynamics on a superconducting quantum computer using learned noise amplification. Complementary 100-qubit Sparse Pauli Dynamics (SPD) simulations achieve comparable accuracy to two-dimensional Richardson extrapolation with fewer circuit settings.
A multi-product formula (MPF) is a promising approach for Hamiltonian simulation efficiently both in the system size N and the inverse allowable error 1 / ε by combining Trotterization and the linear combination of unitaries (LCU). It achieves poly-logarithmic cost in 1 / ε like LCU [G. H. Low, V. Kliuchnikov, N. Wiebe, (2019)]. The efficiency in N is expected to come from the commutator scaling in Trotterization, and this appears to be confirmed by the error bound of MPF expressed by nested commutators [J. Aftab, D. An, K. Trivisa, (2024)]. However, we point out that the efficiency of MPF in the system size N is not exactly resolved yet in that the present error bound expressed by nested commutators is incompatible with the size-efficient complexity reflecting the commutator scaling. The problem is that q -fold nested commutators with arbitrarily large q are involved in their requirement and error bound. The benefit of commutator scaling by locality is absent, and the cost efficient in N becomes prohibited in general. In this paper, we show an alternative commutator-scaling error of MPF and derive its size-efficient cost properly inheriting the advantage in Trotterization. The requirement and the error bound in our analysis, derived by techniques from the Floquet-Magnus expansion, have a certain truncation order in the nested commutators and can fully exploit the locality. We prove that Hamiltonian simulation by MPF certainly achieves the cost whose system-size dependence is as large as Trotterization while keeping the p o l y l o g ( 1 / ε ) -scaling like the LCU. Our results will provide improved or accurate error and cost also for various algorithms using interpolation or extrapolation of Trotterization.