A complex but important challenge in understanding quantum mechanical phenomena is the simulation of quantum many-body dynamics. Although quantum computers offer significant potential to accelerate these simulations, their practical application is currently limited by noise and restricted scalability. In this work, we address these problems by introducing a hybrid ansatz combining the strengths of quantum and classical computational methods. Using Trotterization, we propose to evolve an initial state on the quantum computer according to a simplified Hamiltonian, focusing on terms that are difficult to simulate classically. A classical model then corrects the simulation by including the terms omitted in the quantum circuit. While the classical ansatz is optimized during the time evolution, the quantum circuit has no variational parameters. Derivatives can thus be calculated purely classically, avoiding challenges arising in the optimization of parametrized quantum circuits. We demonstrate three applications of this hybrid method. First, our approach allows us to avoid gates in the quantum circuit by restricting the quantum part of the ansatz to hardware-efficient terms of the Hamiltonian. Second, we can mitigate errors arising from the Trotterization of the time-evolution unitary. Finally, we can extend the system size while keeping the number of qubits in the quantum circuit constant by including additional degrees of freedom in the classical ansatz.
A complex but important challenge in understanding quantum mechanical phenomena is the simulation of quantum many-body dynamics. Although quantum computers offer significant potential to accelerate these simulations, their practical application is currently limited by noise and restricted scalability. In this work, we address these problems by introducing a hybrid ansatz combining the strengths of quantum and classical computational methods. Using Trotterization, we propose to evolve an initial state on the quantum computer according to a simplified Hamiltonian, focusing on terms that are difficult to simulate classically. A classical model then corrects the simulation by including the terms omitted in the quantum circuit. While the classical ansatz is optimized during the time evolution, the quantum circuit has no variational parameters. Derivatives can thus be calculated purely classically, avoiding challenges arising in the optimization of parametrized quantum circuits. We demonstrate three applications of this hybrid method. First, our approach allows us to avoid gates in the quantum circuit by restricting the quantum part of the ansatz to hardware-efficient terms of the Hamiltonian. Second, we can mitigate errors arising from the Trotterization of the time-evolution unitary. Finally, we can extend the system size while keeping the number of qubits in the quantum circuit constant by including additional degrees of freedom in the classical ansatz.
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.
The study of quantum many-body systems is crucial to our understanding of nature at the microscopic scale. As analytical solutions quickly become intractable, the development of efficient computational simulation techniques is indispensable to predict the behaviour of quantum systems. This task is particularly challenging as the dimension of the underlying Hilbert space grows exponentially with the system size. While classical computational methods have achieved remarkable success in specific regimes, their applicability is fundamentally limited by this exponential scaling. Quantum computers offer a natural platform for representing and simulating quantum systems, but current devices are constrained by noise and limited qubit counts, drastically restricting the algorithms that can be run. In this context, hybrid quantum-classical algorithms provide a promising approach to combine the strengths of both paradigms. This thesis develops hybrid quantum-classical algorithms for the simulation of quantum many-body systems. The key idea is to perform classically difficult tasks on near-term quantum devices, while leveraging classical computational power to control, correct, and extend the obtained results. Within this framework, we introduce two complementary approaches for the simulation of quantum dynamics. The first method employs circuit knitting techniques to distribute variational quantum simulations across multiple quantum devices. The entanglement between weakly coupled clusters of highly correlated spins is recovered in classical post-processing, where the overhead is controlled through a constrained optimisation problem. The second method combines a simplified, Trotterised quantum dynamics simulation with a parametrised classical ansatz trained to correct for the approximations made on the quantum device. We show that this method can be employed to correct Trotter errors, implement long-range gates with hardware-efficient circuits, or simulate larger systems by incorporating additional degrees of freedom purely classically. In addition, we present a hybrid approach to computing thermal properties of quantum many-body systems via the Quantum Finite Temperature Lanczos Method. A Krylov subspace is generated on a quantum computer by combining the efficient preparation of typical states with real-time evolution. A classical computer then diagonalises the system in the low-dimensional subspace and estimates the finite-temperature expectation values through trace estimation. The results presented in this thesis demonstrate that hybrid quantum-classical algorithms can expand the range of systems that can be studied on near-term quantum devices. Conversely, we show that incorporating quantum subroutines into classical simulation techniques can improve their efficiency. These results suggest that, rather than treating classical and quantum simulation as competing paradigms, their careful integration can provide a more effective route toward practical quantum many-body simulation.
The advent of near-term digital quantum computers could offer us an exciting opportunity to investigate quantum many-body phenomena beyond that of classical computing. To make the best use of the hardware available, it is paramount that we have methods that accurately simulate Hamiltonian dynamics for limited circuit depths. In this paper, we propose a method to classically optimize unitary brickwall circuits to approximate quantum time evolution operators. Our method is scalable in system size through the use of tensor networks. We demonstrate that, for various three-body Hamiltonians, our approach produces quantum circuits that can outperform trotterization in both their accuracy and the quantum circuit depth needed to implement the dynamics, with the exact details being dependent on the Hamiltonian. We also explain how to choose an optimal time step that minimizes the combined errors of the quantum device and the brickwall circuit approximation. Published by the American Physical Society 2024
Computational methods are central to modern physics, yet many problems, especially those involving Quantum Many-Body Systems (QMBSs), remain challenging at scale. The core difficulty of simulating QMBSs is the exponential growth of the solution space with system size. This renders exact treatments intractable and pushes approximate methods to their limits when tackling two- and higher-dimensional systems and long-time dynamics. Quantum computers offer a quantum-native platform for representing and manipulating such systems. Recent algorithmic and technological advancements have enabled experiments on quantum processors in the quantum utility regime, a scale at which exact classical reproduction is no longer feasible. This represents the first step towards realizing quantum advantage, where quantum computers solve specific, practically relevant problems more efficiently than any classical approach. While advancing towards the goal of quantum advantage, useful computations on quantum processors are expected to be realised through hybrid quantum-classical algorithms, which combine quantum experiments and classical computational resources. A prominent example are variational quantum algorithms, where a Parametrized Quantum Circuit (PQC) prepares a trial quantum state whose parameters are iteratively updated by a classical optimisation routine to minimize a target cost function (e.g., the energy of a target Hamiltonian). However, the optimisation of such PQCs is affected by trainability issues, such as barren plateaus in the parameter landscape, which can render the optimisation increasingly difficult for larger system sizes. This motivates the development of novel hybrid approaches overcoming the limits of PQCs. In this thesis, we introduce a new class of approaches based on hybrid Tensor Networks (TNs). A hybrid TN is obtained by integrating quantum states prepared on a quantum processor within TNs -- one of the state-of-the-art approaches for simulating quantum systems. We first show how hybrid TNs can be used for ground-state searches within a variational setting. Then, we further extend the framework to non-variational hybrid TN ansätze for which the quantum tensors are prepared via quantum annealing and subsequently post-processed by classical TNs. In the context of quantum machine learning, which also relies on PQCs, we further investigate heuristic best practices that mitigate the same exponential concentration effects in kernel-based models for anomaly detection. Another class of hybrid algorithms are quantum error-mitigation workflows, where classical post-processing is used to reconstruct expectation values from noise-affected circuit executions. Such methods generally feature a trade-off between the bias and the variance: the price to pay for mitigating noise in a bias-free fashion is that the number of measurements usually grows exponentially with the depth of the circuit in order to achieve a fixed accuracy. In this thesis, we develop a low-resource error-mitigation strategy that reduces the measurement overhead in practical digital Hamiltonian-simulation workflows. We apply this method to Trotterized dynamics in a lattice-gauge-theory model, enabling the extraction of expectation values with improved accuracy at time-scales that allow the observation of post-collision scattering dynamics.
Understanding the dynamics of quantum systems is crucial in many areas of physics, but simulating many-body systems presents significant challenges due to the large Hilbert space to navigate and the exponential growth of computational overhead. Quantum computers offer a promising platform to overcome these challenges, particularly for simulating the time evolution with Hamiltonians. Trotterization is a widely used approach among available algorithms in this regard, and well suited for near-term quantum devices. However, it introduces algorithmic Trotter errors due to the non-commutativity of Hamiltonian components. Several techniques such as multi-product formulas have been developed to mitigate Trotter errors, but often require deep quantum circuits, which can introduce additional physical errors. In this work, we propose a resource-efficient scheme to reduce the algorithmic Trotter error with relatively shallow circuit depth. We develop a profiling method by introducing an auxiliary parameter to estimate the error effects in expectation values, enabling significant error suppression with a fixed number of Trotter steps. Our approach offers an efficient way of quantum simulation on near-term quantum processors with shallow circuits.