All you need is Trotter

Gumaro Rendon

2024-08-27

Abstract:The work here enables linear cost-scaling with evolution time $t$ while keeping ${\rm polylog} (1/\epsilon)$ scaling and no extra block-encoding qubits, where $\epsilon$ is the algorithmic error. This is achieved through product formulas, stable interpolation (Chebyshev), and to calculate the needed fractional queries, cardinal sine interpolation is used.

Quantum Physics

What problem does this paper attempt to address?

Trotter error time scaling separation via commutant decomposition

Yi-Hsiang Chen

2024-09-25

Abstract: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 $\mathcal{O}(t^{p+1}/r^p)$ and $\mathcal{O}(t^p/r^p)$ for a $p$th-order product formula evolving to time $t$ using $r$ partitions. Under a fixed step size $t/r$, it implies one would scale linearly with time $t$ and the other would be 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.

Quantum Physics
Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

James D. Watson,Jacob Watkins

2024-10-08

Abstract:Product formulae are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by use of these formulae. This work provides an improved, rigorous analysis of these techniques for the task of calculating time-evolved expectation values. We demonstrate that, to achieve error $\epsilon$ in a simulation of time $T$ using a $p^\text{th}$-order product formula with extrapolation, circuits depths of $O\left(T^{1+1/p} \textrm{polylog}(1/\epsilon)\right)$ are sufficient -- an exponential improvement in the precision over product formulae alone. Furthermore, we achieve commutator scaling, improve the complexity with $T$, and do not require fractional implementations of Trotter steps. Our results provide a more accurate characterisation of the algorithmic error mitigation techniques currently proposed to reduce Trotter error.

Quantum Physics,Strongly Correlated Electrons,Data Structures and Algorithms,Numerical Analysis
Trotterization in Quantum Theory

Grant Kluber

2023-10-23

Abstract:Trotterization in quantum mechanics is an important theoretical concept in handling the exponential of noncommutative operators. In this communication, we give a mathematical formulation of the Trotter Product Formula, and apply it to basic examples in which the utility of Trotterization is evident. Originally, this article was completed in December 2020 as a report under the mentorship of Esteban Cárdenas for the University of Texas at Austin Mathematics Directed Reading Program (DRP). However, the relevance of Trotterization in reducing quantum circuit complexity has warranted the release of a revised and more formal version of the original. Thus, we present a mathematical perspective on Trotterization, including a detailed sketch of a formal proof of the Trotter Product Formula.

Quantum Physics,Mathematical Physics
A Theory of Trotter Error

Andrew M. Childs,Yuan Su,Minh C. Tran,Nathan Wiebe,Shuchen Zhu

DOI: https://doi.org/10.1103/PhysRevX.11.011020

2021-02-04

Abstract:The Lie-Trotter formula, together with its higher-order generalizations, provides a direct approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of prior approaches based on truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, including simulations of second-quantized plane-wave electronic structure, $k$-local Hamiltonians, rapidly decaying power-law interactions, clustered Hamiltonians, the transverse field Ising model, and quantum ferromagnets, nearly matching or even outperforming the best previous results. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound as a byproduct. Our analysis reproduces known tight bounds for first- and second-order formulas. Our higher-order bound overestimates the complexity of simulating a one-dimensional Heisenberg model with an even-odd ordering of terms by only a factor of $5$, and is close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both asymptotic scaling and constant prefactor.

Quantum Physics,Strongly Correlated Electrons,Data Structures and Algorithms,Numerical Analysis,Chemical Physics
Improved Accuracy for Trotter Simulations Using Chebyshev Interpolation

Gumaro Rendon,Jacob Watkins,Nathan Wiebe

DOI: https://doi.org/10.22331/q-2024-02-26-1266

2024-02-23

Abstract:Quantum metrology allows for measuring properties of a quantum system at the optimal Heisenberg limit. However, when the relevant quantum states are prepared using digital Hamiltonian simulation, the accrued algorithmic errors will cause deviations from this fundamental limit. In this work, we show how algorithmic errors due to Trotterized time evolution can be mitigated through the use of standard polynomial interpolation techniques. Our approach is to extrapolate to zero Trotter step size, akin to zero-noise extrapolation techniques for mitigating hardware errors. We perform a rigorous error analysis of the interpolation approach for estimating eigenvalues and time-evolved expectation values, and show that the Heisenberg limit is achieved up to polylogarithmic factors in the error. Our work suggests that accuracies approaching those of state-of-the-art simulation algorithms may be achieved using Trotter and classical resources alone for a number of relevant algorithmic tasks.

Quantum Physics
Making Trotterization adaptive and energy-self-correcting for NISQ devices and beyond

Hongzheng Zhao,Marin Bukov,Markus Heyl,Roderich Moessner

2023-08-15

Abstract:Simulation of continuous time evolution requires time discretization on both classical and quantum computers. A finer time step improves simulation precision, but it inevitably leads to increased computational efforts. This is particularly costly for today's noisy intermediate scale quantum computers, where notable gate imperfections limit the circuit depth that can be executed at a given accuracy. Classical adaptive solvers are well-developed to save numerical computation times. However, it remains an outstanding challenge to make optimal usage of the available quantum resources by means of adaptive time steps. Here, we introduce a quantum algorithm to solve this problem, providing a controlled solution of the quantum many-body dynamics of local observables. The key conceptual element of our algorithm is a feedback loop which self-corrects the simulation errors by adapting time steps, thereby significantly outperforming conventional Trotter schemes on a fundamental level and reducing the circuit depth. It even allows for a controlled asymptotic long-time error, where usual Trotterized dynamics is facing difficulties. Another key advantage of our quantum algorithm is that any desired conservation law can be included in the self-correcting feedback loop, which has potentially a wide range of applicability. We demonstrate the capabilities by enforcing gauge invariance which is crucial for a faithful and long-sought quantum simulation of lattice gauge theories. Our algorithm can be potentially useful on a more general level whenever time discretization is involved concerning, for instance, also numerical approaches based on time-evolving block decimation methods.

Quantum Physics,Statistical Mechanics,Strongly Correlated Electrons
Measuring Trotter error and its application to precision-guaranteed Hamiltonian simulations

Tatsuhiko N. Ikeda,Hideki Kono,Keisuke Fujii

2024-07-04

Abstract:Trotterization is the most common and convenient approximation method for Hamiltonian simulations on digital quantum computers, but estimating its error accurately is computationally difficult for large quantum systems. Here, we develop a method for measuring the Trotter error without ancillary qubits on quantum circuits by combining the $m$th- and $n$th-order ($m<n$) Trotterizations rather than consulting with mathematical error bounds. Using this method, we make Trotterization precision-guaranteed, developing an algorithm named Trotter$(m,n)$, in which the Trotter error at each time step is within an error tolerance $\epsilon$ preset for our purpose. Trotter$(m,n)$ is applicable to both time- independent and dependent Hamiltonians, and it adaptively chooses almost the largest stepsize $\mathrm{d}t$, which keeps quantum circuits shallowest within the error tolerance. Benchmarking it in a quantum spin chain, we find the adaptively chosen $\mathrm{d}t$ to be about ten times larger than that inferred from known upper bounds of Trotter errors.

Quantum Physics,Materials Science,Strongly Correlated Electrons,High Energy Physics - Lattice,Computational Physics
Lower Bounds for the Trotter Error

Alexander Hahn,Paul Hartung,Daniel Burgarth,Paolo Facchi,Kazuya Yuasa

2024-10-04

Abstract:In analog and digital simulations of practically relevant quantum systems, the target dynamics can only be implemented approximately. The Trotter product formula is the most common approximation scheme as it is a generic method which allows tuning accuracy. The Trotter simulation precision will always be inexact for non-commuting operators, but it is currently unknown what the minimum possible error is. This is an important quantity because upper bounds for the Trotter error are known to often be vast overestimates. Here, we present explicit lower bounds on the error, in norm and on states, allowing to derive minimum resource requirements. Numerical comparison with the true error shows that our bounds offer accurate and tight estimates.

Quantum Physics
Adaptive Trotterization for time-dependent Hamiltonian quantum dynamics using piecewise conservation laws

Hongzheng Zhao,Marin Bukov,Markus Heyl,Roderich Moessner

2024-06-20

Abstract:Digital quantum simulation relies on Trotterization to discretize time evolution into elementary quantum gates. On current quantum processors with notable gate imperfections, there is a critical tradeoff between improved accuracy for finer timesteps, and increased error rate on account of the larger circuit depth. We present an adaptive Trotterization algorithm to cope with time-dependent Hamiltonians, where we propose a concept of piecewise "conserved" quantities to estimate errors in the time evolution between two (nearby) points in time; these allow us to bound the errors accumulated over the full simulation period. They reduce to standard conservation laws in the case of time-independent Hamiltonians, for which we first developed an adaptive Trotterization scheme [PRX Quantum 4, 030319]. We validate the algorithm for a time-dependent quantum spin chain, demonstrating that it can outperform the conventional Trotter algorithm with a fixed step size at a controlled error.

Quantum Physics,Quantum Gases,Statistical Mechanics
Efficient and practical Hamiltonian simulation from time-dependent product formulas

Jan Lukas Bosse,Andrew M. Childs,Charles Derby,Filippo Maria Gambetta,Ashley Montanaro,Raul A. Santos

2024-06-25

Abstract:In this work we propose an approach for implementing time-evolution of a quantum system using product formulas. The quantum algorithms we develop have provably better scaling (in terms of gate complexity and circuit depth) than a naive application of well-known Trotter formulas, for systems where the evolution is determined by a Hamiltonian with different energy scales (i.e., one part is "large" and another part is "small"). Our algorithms generate a decomposition of the evolution operator into a product of simple unitaries that are directly implementable on a quantum computer. Although the theoretical scaling is suboptimal compared with state-of-the-art algorithms (e.g., quantum signal processing), the performance of the algorithms we propose is highly competitive in practice. We illustrate this via extensive numerical simulations for several models. For instance, in the strong-field regime of the 1D transverse-field Ising model, our algorithms achieve an improvement of one order of magnitude in both the system size and evolution time that can be simulated with a fixed budget of 1000 arbitrary 2-qubit gates, compared with standard Trotter formulas.

Quantum Physics
Trotter error bounds and dynamic multi-product formulas for Hamiltonian simulation

Sergiy Zhuk,Niall Robertson,Sergey Bravyi

2024-02-09

Abstract:Multi-product formulas (MPF) are linear combinations of Trotter circuits offering high-quality simulation of Hamiltonian time evolution with fewer Trotter steps. Here we report two contributions aimed at making multi-product formulas more viable for near-term quantum simulations. First, we extend the theory of Trotter error with commutator scaling developed by Childs, Su, Tran et al. to multi-product formulas. Our result implies that multi-product formulas can achieve a quadratic reduction of Trotter error in 1-norm (nuclear norm) on arbitrary time intervals compared with the regular product formulas without increasing the required circuit depth or qubit connectivity. The number of circuit repetitions grows only by a constant factor. Second, we introduce dynamic multi-product formulas with time-dependent coefficients chosen to minimize a certain efficiently computable proxy for the Trotter error. We use a minimax estimation method to make dynamic multi-product formulas robust to uncertainty from algorithmic errors, sampling and hardware noise. We call this method Minimax MPF and we provide a rigorous bound on its error.

Quantum Physics,Mathematical Physics,Numerical Analysis,Optimization and Control
Simple Ways to improve Discrete Time Evolution

Johann Ostmeyer

DOI: https://doi.org/10.22323/1.453.0025

2023-09-07

Abstract:Suzuki-Trotter decompositions of exponential operators like $\exp(Ht)$ are required in almost every branch of numerical physics. Often the exponent under consideration has to be split into more than two operators, for instance as local gates on quantum computers. In this work, we demonstrate how highly optimised schemes originally derived for exactly two operators can be applied to such generic Suzuki-Trotter decompositions. After this first trick, we explain what makes an efficient decomposition and how to choose from the large variety available. Furthermore we demonstrate that many problems for which a Suzuki-Trotter decomposition might appear to be the canonical ansatz, are better approached with different methods like Taylor or Chebyshev expansions. In particular, we derive an efficient and numerically stable method to implement truncated polynomial expansions based on a linear factorisation using their complex zeros.

Quantum Physics,Strongly Correlated Electrons,High Energy Physics - Lattice,Computational Physics
Optimised Trotter Decompositions for Classical and Quantum Computing

Johann Ostmeyer

DOI: https://doi.org/10.1088/1751-8121/acde7a

2023-06-16

Abstract:Suzuki-Trotter decompositions of exponential operators like $\exp(Ht)$ are required in almost every branch of numerical physics. Often the exponent under consideration has to be split into more than two operators $H=\sum_k A_k$, for instance as local gates on quantum computers. We demonstrate how highly optimised schemes originally derived for exactly two operators $A_{1,2}$ can be applied to such generic Suzuki-Trotter decompositions, providing a formal proof of correctness as well as numerical evidence of efficiency. A comprehensive review of existing symmetric decomposition schemes up to order $n\le4$ is presented and complemented by a number of novel schemes, including both real and complex coefficients. We derive the theoretically most efficient unitary and non-unitary 4th order decompositions. The list is augmented by several exceptionally efficient schemes of higher order $n\le8$. Furthermore we show how Taylor expansions can be used on classical devices to reach machine precision at a computational effort at which state of the art Trotterization schemes do not surpass a relative precision of $10^{-4}$. Finally, a short and easily understandable summary explains how to choose the optimal decomposition in any given scenario.

Quantum Physics,Statistical Mechanics,Strongly Correlated Electrons,High Energy Physics - Lattice,Computational Physics
Optimal Trotterization in universal quantum simulators under faulty control

George C. Knee,William J. Munro

DOI: https://doi.org/10.1103/physreva.91.052327

IF: 2.971

2015-05-26

Physical Review A

Abstract:Universal quantum simulation may provide insights into those many-body systems that cannot be described classically and that cannot be efficiently simulated with current technology. The Trotter formula, which decomposes a desired unitary time evolution of the simulator into a stroboscopic sequence of repeated elementary evolutions, is a key algorithmic component which makes quantum simulation of dynamics tractable. The Trotter number n sets the time scale on which a computer running this algorithm is switched from one elementary evolution to another. In the ideal case, the precision of the simulation can be arbitrarily controlled by increasing n. We study a more realistic scenario where each gate is applied imperfectly. The resultant trade-off in errors leads to an ultimate limit on the precision of the simulation. We calculate the optimum Trotter number n* that achieves this limit, which is the minimum statistical distance from the actual simulation to the ideal one.

optics,physics, atomic, molecular & chemical
State-dependent Trotter Limits and their approximations

Daniel Burgarth,Niklas Galke,Alexander Hahn,Lauritz van Luijk

DOI: https://doi.org/10.1103/PhysRevA.107.L040201

2024-02-28

Abstract:The Trotter product formula is a key instrument in numerical simulations of quantum systems. However, computers cannot deal with continuous degrees of freedom, such as the position of particles in molecules, or the amplitude of electromagnetic fields. It is therefore necessary to discretize these variables to make them amenable to digital simulations. Here, we give sufficient conditions to conclude the validity of this approximate discretized physics. Essentially, it depends on the state-dependent Trotter error, for which we establish explicit bounds that are also of independent interest.

Quantum Physics,Mathematical Physics
Minimum Trotterization Formulas for a Time-Dependent Hamiltonian

Tatsuhiko N. Ikeda,Asir Abrar,Isaac L. Chuang,Sho Sugiura

DOI: https://doi.org/10.22331/q-2023-11-06-1168

2023-11-02

Abstract:When a time propagator $e^{\delta t A}$ for duration $\delta t$ consists of two noncommuting parts $A=X+Y$, Trotterization approximately decomposes the propagator into a product of exponentials of $X$ and $Y$. Various Trotterization formulas have been utilized in quantum and classical computers, but much less is known for the Trotterization with the time-dependent generator $A(t)$. Here, for $A(t)$ given by the sum of two operators $X$ and $Y$ with time-dependent coefficients $A(t) = x(t) X + y(t) Y$, we develop a systematic approach to derive high-order Trotterization formulas with minimum possible exponentials. In particular, we obtain fourth-order and sixth-order Trotterization formulas involving seven and fifteen exponentials, respectively, which are no more than those for time-independent generators. We also construct another fourth-order formula consisting of nine exponentials having a smaller error coefficient. Finally, we numerically benchmark the fourth-order formulas in a Hamiltonian simulation for a quantum Ising chain, showing that the 9-exponential formula accompanies smaller errors per local quantum gate than the well-known Suzuki formula.

Quantum Physics,Mathematical Physics
Average-Case Speedup for Product Formulas

Chi-Fang Chen,Fernando G. S. L. Brandão

DOI: https://doi.org/10.1007/s00220-023-04912-5

IF: 2.361

2024-02-06

Communications in Mathematical Physics

Abstract:Quantum simulation is a promising application of future quantum computers. Product formulas, or Trotterization, are the oldest and still remain an appealing method to simulate quantum systems. For an accurate product formula approximation, the state-of-the-art gate complexity depends on the number of terms in the Hamiltonian and a local energy estimate. In this work, we give evidence that product formulas, in practice, may work much better than expected. We prove that the Trotter error exhibits a qualitatively better scaling for the vast majority of input states, while the existing estimate is for the worst states. For general k -local Hamiltonians and higher-order product formulas, we obtain gate count estimates for input states drawn from any orthogonal basis. The gate complexity significantly improves over the worst case for systems with large connectivity. Our typical-case results generalize to Hamiltonians with Fermionic terms, with input states drawn from a fixed-particle number subspace, and with Gaussian coefficients (e.g., the SYK models). Technically, we employ a family of simple but versatile inequalities from non-commutative martingales called uniform smoothness , which leads to Hypercontractivity , namely p -norm estimates for k -local operators. This delivers concentration bounds via Markov's inequality. For optimality, we give analytic and numerical examples that simultaneously match our typical-case estimates and the existing worst-case estimates. Therefore, our improvement is due to asking a qualitatively different question, and our results open doors to the study of quantum algorithms in the average case.

physics, mathematical
Compilation of Trotter-Based Time Evolution for Partially Fault-Tolerant Quantum Computing Architecture

Yutaro Akahoshi,Riki Toshio,Jun Fujisaki,Hirotaka Oshima,Shintaro Sato,Keisuke Fujii

2024-10-23

Abstract:Achieving practical quantum speedup with limited resources is a crucial challenge in both academic and industrial communities. To address this, a partially fault-tolerant quantum computing architecture called ``space-time efficient analog rotation quantum computing architecture (STAR architecture)'' has been recently proposed. This architecture focuses on minimizing resource requirements while maximizing the precision of non-Clifford gates, essential for universal quantum computation. However, non-deterministic processes such as the repeat-until-success (RUS) protocol and state injection can introduce significant computational overhead. Therefore, optimizing the logical circuit to minimize this overhead by using efficient fault-tolerant operations is essential. This paper presents an efficient method for simulating the time evolution of the 2D Hubbard model Hamiltonian, a promising application of the STAR architecture. We present two techniques, parallel injection protocol and adaptive injection region updating, to reduce unnecessary time overhead specific to our architecture. By integrating these with the existing fSWAP technique, we develop an efficient Trotter-based time evolution operation for the 2D Hubbard model. Our analysis reveals an acceleration of over 10 times compared to naive serial compilation. This optimized compilation enables us to estimate the computational resources required for quantum phase estimation of the 2D Hubbard model. For devices with a physical error rate of $p_{\rm phys} = 10^{-4}$, we estimate that approximately $6.5 \times 10^4$ physical qubits are required to achieve faster ground state energy estimation of the $8\times8$ Hubbard model compared to classical computation.

Quantum Physics,Strongly Correlated Electrons

All you need is Trotter

Trotter error time scaling separation via commutant decomposition

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Trotterization in Quantum Theory

A Theory of Trotter Error

Improved Accuracy for Trotter Simulations Using Chebyshev Interpolation

Making Trotterization adaptive and energy-self-correcting for NISQ devices and beyond

Measuring Trotter error and its application to precision-guaranteed Hamiltonian simulations

Lower Bounds for the Trotter Error

Adaptive Trotterization for time-dependent Hamiltonian quantum dynamics using piecewise conservation laws

Efficient and practical Hamiltonian simulation from time-dependent product formulas

Trotter error bounds and dynamic multi-product formulas for Hamiltonian simulation

Simple Ways to improve Discrete Time Evolution

Optimised Trotter Decompositions for Classical and Quantum Computing

Optimal Trotterization in universal quantum simulators under faulty control

State-dependent Trotter Limits and their approximations

Minimum Trotterization Formulas for a Time-Dependent Hamiltonian

Average-Case Speedup for Product Formulas

Compilation of Trotter-Based Time Evolution for Partially Fault-Tolerant Quantum Computing Architecture