Optimized Lie-Trotter-Suzuki decompositions for two and three non-commuting terms

Abstract

Lie-Trotter-Suzuki decompositions are an efficient way to approximate operator exponentials $\exp (tH)$ when $H$ is a sum of $n$ (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. We provide optimized decompositions up to order $t^6$. The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For $n=2$ terms, several of the optima we find are close to those in McLachlan, SlAM J. Sci. Comput. 16, 151 (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

30 pages, 8 figures, 8 tables; added results, figures, and references, extended discussion

References in corpus (13)

• The density-matrix renormalization group in the age of matrix product states
• Real time evolution using the density matrix renormalization group
• An Open-System Quantum Simulator with Trapped Ions
• Simulating Hamiltonian dynamics with a truncated Taylor series
• Time-evolution methods for matrix-product states
• Criticality, the area law, and the computational power of PEPS
• Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions
• The iTEBD algorithm beyond unitary evolution
• Applications of quantum Monte Carlo methods in condensed systems
• Time-step targetting methods for real-time dynamics using DMRG
• Time evolution algorithms for Matrix Product States and DMRG
• Lanczos algorithm with Matrix Product States for dynamical correlation functions
• One-dimensional quantum systems at finite temperatures can be simulated efficiently on classical computers

Cited by in corpus (29)

• Stability of superdiffusion in nearly integrable spin chains
• The Folded Spin-1/2 XXZ Model: II. Thermodynamics and Hydrodynamics with a Minimal Set of Charges
• Time-dependent study of disordered models with infinite projected entangled pair states
• Variational Hamiltonian simulation for translational invariant systems via classical pre-processing
• Many-body localization in tilted and harmonic potentials
• What is a quantum shock wave?
• Symmetric cluster expansions with tensor networks
• Quantum-classical eigensolver using multiscale entanglement renormalization
• Comparing numerical methods for hydrodynamics in a one-dimensional lattice spin model
• Problem specific classical optimization of Hamiltonian simulation
• Measurement catastrophe and ballistic spread of charge density with vanishing current
• Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases
• Molecular Structure Optimization based on Electrons-Nuclei Quantum Dynamics Computation
• Accelerating Quantum Optimal Control of Multi-Qubit Systems with Symmetry-Based Hamiltonian Transformations
• Novel Trotter formulas for digital quantum simulation
• Non-Hermiticity in quantum nonlinear optics through symplectic transformations
• Many-body effects on high-harmonic generation in Hubbard ladders
• Breakdown of superdiffusion in perturbed quantum integrable spin chains and ladders
• Missing Rung Problem in Vibrational Ladder Climbing
• Recipes for the Digital Quantum Simulation of Lattice Spin Systems
• Optimize quantum simulation using a force-gradient integrator
• Phase-space stochastic quantum hydrodynamics for interacting Bose gases
• Unifying framework for quantum simulation algorithms for time-dependent Hamiltonian dynamics
• A Quantum-Inspired Algorithm for Wave Simulation Using Tensor Networks
• Spinsim: a GPU optimized python package for simulating spin-half and spin-one quantum systems
• Scaling of contraction costs for entanglement renormalization algorithms including tensor Trotterization and variational Monte Carlo
• Crystalline topological defects within response theory
• Hybrid quantum gap estimation algorithm using a filtered time series
• Pauli weight requirement of the matrix elements in time-evolved local operators: dependence beyond the equilibration temperature