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

Classical Variational Optimization of Gate Sequences for Time Evolution of Large Quantum Systems

Refik Mansuroglu, Samuel A. Wilkinson, Ludwig Nützel, Michael J. Hartmann

Published Jun 7, 2021• 2 citations
OpenAlex

Abstract

The simulation of time evolution of large quantum systems is a classically challenging and often intractable task, making it a promising application for quantum computation. A Trotter-Suzuki approximation yields an implementation thereof, where a certain desired accuracy can be achieved by raising the gate count adequately. In this work, we introduce a variational algorithm which uses solutions of classical optimizations to predict efficient quantum circuits for time evolution of large quantum systems. Our strategy improves on the Trotter-Suzuki ansatz in accuracy and gate count by several orders of magnitude. In a NISQ-friendly setting, we can either significantly reduce the approximation error while using the same amount of gates as the Trotter-Suzuki ansatz or achieve a comparable accuracy with significantly fewer gates. We find, that our strategy outperforms the Trotter-Suzuki benchmark in translation invariant systems as well as for open boundaries.

Related papers

Variational Hamiltonian simulation for translational invariant systems via classical pre-processing

Refik Mansuroglu, Timo Eckstein, Ludwig Nützel +2 more

Abstract The simulation of time evolution of large quantum systems is a classically challenging and in general intractable task, making it a promising application for quantum computation. A Trotter–Suzuki approximation yields an implementation thereof, where a higher approximation accuracy can be traded for an increased gate count. In this work, we introduce a variational algorithm which uses solutions of classical optimizations to predict efficient quantum circuits for time evolution of translationally invariant quantum systems. Our strategy can improve upon the Trotter–Suzuki accuracy by several orders of magnitude. It translates into a reduction in gate count and hence gain in overall fidelity at the same algorithmic accuracy. This is important in noisy intermediate scale quantum-applications where the fidelity of the output state decays exponentially with the number of gates. The performance advantage of our classical assisted strategy can be extended to open boundaries with translational symmetry in the bulk. We can extrapolate our method to beyond classically simulatable system sizes, maintaining its total fidelity advantage over a Trotter–Suzuki approximation making it an interesting candidate for beyond classical time evolution.

arXiv:2106.03680•Quantum Science and Technology•Jan 10, 2023•37 citations

Beyond trotterization: Variational product formulas for quantum simulation

I. A. Assi, Michael Rudolf Vogl, Meenu Kumari +1 more

We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global action principle, where Euler-Lagrange equations govern their optimal dynamics. Unlike conventional wavefunction-based variational methods, our approach specifically targets the time evolution operation, and this allows a single set of optimized parameters to be applied to any initial state for a fixed Hamiltonian avoiding costly optimization procedures. Our method outperforms the standard Trotter-Suzuki formulas, typically achieving higher accuracy than higher-order Suzuki schemes. This translates directly to quantum computing applications, where it enables the design of quantum circuits with fewer gates, which reduces noise and improves precision. Although we focus on quantum dynamics, the method is broadly applicable to problems involving general time-evolution operators. Applied to various model Hamiltonians, our approach reduces errors by factors of 2 to 5 compared to Trotter-Suzuki decompositions, demonstrating its promise for accurate quantum simulation with improved efficiency. In certain cases, the variational ansatz achieves higher accuracy than more complex higher-order Suzuki formulas while reducing the gate count by nearly half within a single circuit layer. Furthermore, we derive approximate analytical expressions for the variational parameters up to cubic order in time, valid for generic Hamiltonians. These approximations enable long-time quantum simulations with improved accuracy over equivalent Suzuki decompositions, providing ready-to-use evolution formulas that match Suzuki's gate complexity while delivering better performance.

Physical review. B./Physical review. B•May 21, 2026•0 citations

Beyond Trotterization: Variational Product Formulas for Quantum Simulation

Ibsal Assi, Michael Vogl, Meenu Kumari +1 more

We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global action principle, where Euler-Lagrange equations govern their optimal dynamics. Unlike conventional wavefunction-based variational methods, our approach specifically targets the time evolution operation and this allows a single set of optimized parameters to be applied to any initial state for a fixed Hamiltonian avoiding costly optimization procedures. Our method outperforms the standard Trotter-Suzuki formulas, typically achieving higher accuracy than higher-order Suzuki schemes. This translates directly to quantum computing applications, where it enables the design of quantum circuits with fewer gates which reduces noise and improves precision. Although we focus on quantum dynamics, the method is broadly applicable to problems involving general time-evolution operators. Applied to various model Hamiltonians, our approach reduces errors by factors of 2 to 5 compared to Trotter-Suzuki decompositions, demonstrating its promise for accurate quantum simulation with improved efficiency. In certain cases, the variational ansatz achieves higher accuracy than more complex higher-order Suzuki formulas while reducing the gate count by nearly half within a single circuit layer. Furthermore, we derive approximate analytical expressions for the variational parameters up to cubic order in time, valid for generic Hamiltonians. These approximations enable long-time quantum simulations with improved accuracy over equivalent Suzuki decompositions, providing ready-to-use evolution formulas that match Suzuki's gate complexity while delivering better performance.

quant-pharXiv:2511.15124•Nov 20, 2025•0 citations

Variational Optimization of the Time Evolution Operator in Correlated Systems

Kaelyn Ferris, Sergio E. Ulloa

Simulating the time evolution of quantum systems is a classically challenging and often intractable task due to the exponential scaling of its constituent operators. This scaling makes it a promising application for quantum computation, especially for correlated systems. Constructing an approximate time evolution operator can be done using the Trotter-Suzuki decomposition with an error that scales at most polynomially in the target simulation time and can be decreased at the cost of longer circuit depth. Variational approaches to classically optimize this approximation have shown significant improvement in accuracy for large quantum systems in both gate count and accuracy by several orders of magnitude. However, it is unclear if a similar improvement can be obtained for strongly correlated systems. We extend this variational approach to two example systems with strong correlations: a one-dimensional Fermi-Hubbard chain and an extended Fermi-Hubbard chain. The encoding between qubits and fermionic modes in these systems is generally nonlocal and so we propose an encoding scheme for these Hamiltonians with minimal overall operator distance between nearest neighbor connected qubits. We demonstrate that this variational method grants a similarly strong improvement in both gate count and accuracy and examine the performance of these optimized gate sequences using realistic noise models to estimate their performance on noisy quantum hardware. Finally, we discuss an approach to combine the optimized nonlocal gates of small chains to larger systems without the need for additional computational resources.

IEEE International Conference on Quantum Computing and Engineering (QCE)•Sep 1, 2022•1 citation

Problem-specific classical optimization of Hamiltonian simulation

Refik Mansuroglu, Felix Fischer, Michael J. Hartmann

Nonequilibrium time evolution of large quantum systems is a strong candidate for quantum advantage. Variational quantum algorithms have been put forward for this task, but their quantum optimization routines suffer from trainability and sampling problems. Here, we present a classical preprocessing routine for variational Hamiltonian simulation that circumvents the need for a quantum optimization by expanding rigorous error bounds in a perturbative regime for suitable time steps. The resulting cost function is efficiently computable on a classical computer. We show that there always exists potential for optimization with respect to a Trotter sequence of the same order and that the cost value has the same scaling as Trotter in simulation time and system size. Unlike previous work on classical preprocessing, the method is applicable to any Hamiltonian system independent of locality and interaction lengths. Via numerical experiments for spin-lattice models, we find that our approach significantly improves digital quantum simulation capabilities with respect to Trotter sequences for the same resources. For short times, we find accuracy improvements of more than three orders of magnitude for our method as compared to Trotter sequences of the same gate number. Moreover, for a given gate number and accuracy target, we find that the preoptimization we introduce enables simulation times that are consistently more than ten times longer for a target accuracy of 0.1%.

arXiv:2306.07208•Physical Review Research•Oct 11, 2023•13 citations

Reducing the Gate Count with Efficient Trotter-Suzuki Schemes

Marko Maležič, Johann Ostmeyer

Hamiltonian formulations of lattice field theories provide access to real-time dynamics, but their simulation is difficult to implement efficiently. Trotter-Suzuki decompositions are at the center of time evolution computation, either on quantum hardware or classically, for instance with the use of tensor networks. While low-order Trotterizations remain the standard choice due to their simplicity, higher-order schemes offer the potential for improved efficiency. In this work we outline a short guide to Trotter-Suzuki schemes and their implementations in general. To help with this, we highlight new efficient schemes found by our optimization framework, and demonstrate their performance on the Heisenberg model.

hep-latCross-listedarXiv:2602.21145•Feb 25, 2026•0 citations
The simulation of time evolution of large quantum systems is a classically challenging and often intractable task, making it a promising application for quantum computation. A Trotter-Suzuki approximation yields an implementation thereof, where a certain desired accuracy can be achieved by raising the gate count adequately. In this work, we introduce a variational algorithm which uses solutions of classical optimizations to predict efficient quantum circuits for time evolution of large quantum systems. Our strategy improves on the Trotter-Suzuki ansatz in accuracy and gate count by several orders of magnitude. In a NISQ-friendly setting, we can either significantly reduce the approximation error while using the same amount of gates as the Trotter-Suzuki ansatz or achieve a comparable accuracy with significantly fewer gates. We find, that our strategy outperforms the Trotter-Suzuki benchmark in translation invariant systems as well as for open boundaries.
Abstract The simulation of time evolution of large quantum systems is a classically challenging and in general intractable task, making it a promising application for quantum computation. A Trotter–Suzuki approximation yields an implementation thereof, where a higher approximation accuracy can be traded for an increased gate count. In this work, we introduce a variational algorithm which uses solutions of classical optimizations to predict efficient quantum circuits for time evolution of translationally invariant quantum systems. Our strategy can improve upon the Trotter–Suzuki accuracy by several orders of magnitude. It translates into a reduction in gate count and hence gain in overall fidelity at the same algorithmic accuracy. This is important in noisy intermediate scale quantum-applications where the fidelity of the output state decays exponentially with the number of gates. The performance advantage of our classical assisted strategy can be extended to open boundaries with translational symmetry in the bulk. We can extrapolate our method to beyond classically simulatable system sizes, maintaining its total fidelity advantage over a Trotter–Suzuki approximation making it an interesting candidate for beyond classical time evolution.
We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global action principle, where Euler-Lagrange equations govern their optimal dynamics. Unlike conventional wavefunction-based variational methods, our approach specifically targets the time evolution operation, and this allows a single set of optimized parameters to be applied to any initial state for a fixed Hamiltonian avoiding costly optimization procedures. Our method outperforms the standard Trotter-Suzuki formulas, typically achieving higher accuracy than higher-order Suzuki schemes. This translates directly to quantum computing applications, where it enables the design of quantum circuits with fewer gates, which reduces noise and improves precision. Although we focus on quantum dynamics, the method is broadly applicable to problems involving general time-evolution operators. Applied to various model Hamiltonians, our approach reduces errors by factors of 2 to 5 compared to Trotter-Suzuki decompositions, demonstrating its promise for accurate quantum simulation with improved efficiency. In certain cases, the variational ansatz achieves higher accuracy than more complex higher-order Suzuki formulas while reducing the gate count by nearly half within a single circuit layer. Furthermore, we derive approximate analytical expressions for the variational parameters up to cubic order in time, valid for generic Hamiltonians. These approximations enable long-time quantum simulations with improved accuracy over equivalent Suzuki decompositions, providing ready-to-use evolution formulas that match Suzuki's gate complexity while delivering better performance.
We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global action principle, where Euler-Lagrange equations govern their optimal dynamics. Unlike conventional wavefunction-based variational methods, our approach specifically targets the time evolution operation and this allows a single set of optimized parameters to be applied to any initial state for a fixed Hamiltonian avoiding costly optimization procedures. Our method outperforms the standard Trotter-Suzuki formulas, typically achieving higher accuracy than higher-order Suzuki schemes. This translates directly to quantum computing applications, where it enables the design of quantum circuits with fewer gates which reduces noise and improves precision. Although we focus on quantum dynamics, the method is broadly applicable to problems involving general time-evolution operators. Applied to various model Hamiltonians, our approach reduces errors by factors of 2 to 5 compared to Trotter-Suzuki decompositions, demonstrating its promise for accurate quantum simulation with improved efficiency. In certain cases, the variational ansatz achieves higher accuracy than more complex higher-order Suzuki formulas while reducing the gate count by nearly half within a single circuit layer. Furthermore, we derive approximate analytical expressions for the variational parameters up to cubic order in time, valid for generic Hamiltonians. These approximations enable long-time quantum simulations with improved accuracy over equivalent Suzuki decompositions, providing ready-to-use evolution formulas that match Suzuki's gate complexity while delivering better performance.
Simulating the time evolution of quantum systems is a classically challenging and often intractable task due to the exponential scaling of its constituent operators. This scaling makes it a promising application for quantum computation, especially for correlated systems. Constructing an approximate time evolution operator can be done using the Trotter-Suzuki decomposition with an error that scales at most polynomially in the target simulation time and can be decreased at the cost of longer circuit depth. Variational approaches to classically optimize this approximation have shown significant improvement in accuracy for large quantum systems in both gate count and accuracy by several orders of magnitude. However, it is unclear if a similar improvement can be obtained for strongly correlated systems. We extend this variational approach to two example systems with strong correlations: a one-dimensional Fermi-Hubbard chain and an extended Fermi-Hubbard chain. The encoding between qubits and fermionic modes in these systems is generally nonlocal and so we propose an encoding scheme for these Hamiltonians with minimal overall operator distance between nearest neighbor connected qubits. We demonstrate that this variational method grants a similarly strong improvement in both gate count and accuracy and examine the performance of these optimized gate sequences using realistic noise models to estimate their performance on noisy quantum hardware. Finally, we discuss an approach to combine the optimized nonlocal gates of small chains to larger systems without the need for additional computational resources.
Nonequilibrium time evolution of large quantum systems is a strong candidate for quantum advantage. Variational quantum algorithms have been put forward for this task, but their quantum optimization routines suffer from trainability and sampling problems. Here, we present a classical preprocessing routine for variational Hamiltonian simulation that circumvents the need for a quantum optimization by expanding rigorous error bounds in a perturbative regime for suitable time steps. The resulting cost function is efficiently computable on a classical computer. We show that there always exists potential for optimization with respect to a Trotter sequence of the same order and that the cost value has the same scaling as Trotter in simulation time and system size. Unlike previous work on classical preprocessing, the method is applicable to any Hamiltonian system independent of locality and interaction lengths. Via numerical experiments for spin-lattice models, we find that our approach significantly improves digital quantum simulation capabilities with respect to Trotter sequences for the same resources. For short times, we find accuracy improvements of more than three orders of magnitude for our method as compared to Trotter sequences of the same gate number. Moreover, for a given gate number and accuracy target, we find that the preoptimization we introduce enables simulation times that are consistently more than ten times longer for a target accuracy of 0.1%.
Hamiltonian formulations of lattice field theories provide access to real-time dynamics, but their simulation is difficult to implement efficiently. Trotter-Suzuki decompositions are at the center of time evolution computation, either on quantum hardware or classically, for instance with the use of tensor networks. While low-order Trotterizations remain the standard choice due to their simplicity, higher-order schemes offer the potential for improved efficiency. In this work we outline a short guide to Trotter-Suzuki schemes and their implementations in general. To help with this, we highlight new efficient schemes found by our optimization framework, and demonstrate their performance on the Heisenberg model.