We propose a hybrid approach to simulate quantum many-body dynamics by combining Trotter-based quantum evolution with classical dynamic mode decomposition. In quantum simulation, the interest often lies in estimating observables rather than explicitly obtaining the wave function. Instead of characterizing a quantum state, our method predicts observables of the quantum state at later times by using data from a set of short-time measurements from a quantum computer. For a fixed set of measurements, the global error bound for our method scales at most as $\mathcal{O}(t^{3/2})$, with constants determined by the DMD approximation error and the conditioning of the learned modes. We apply our method to quench dynamics in the Hubbard model and nearest neighbor spin systems and show that the observable properties can be predicted up to a reasonable error by controlling the number of data points obtained from the quantum measurements.
We propose a hybrid approach to simulate quantum many-body dynamics by combining Trotter-based quantum evolution with classical dynamic mode decomposition. In quantum simulation, the interest often lies in estimating observables rather than explicitly obtaining the wave function. Instead of characterizing a quantum state, our method predicts observables of the quantum state at later times by using data from a set of short-time measurements from a quantum computer. For a fixed set of measurements, the global error bound for our method scales at most as O(t3/2), with constants determined by the DMD approximation error and the conditioning of the learned modes. We apply our method to quench dynamics in the Hubbard model and nearest neighbor spin systems and show that the observable properties can be predicted up to a reasonable error by controlling the number of data points obtained from the quantum measurements.
Quantum computers offer the potential to efficiently simulate the dynamics of quantum systems, a task whose difficulty scales exponentially with system size on classical devices. To assess the potential for near-term quantum computers to simulate many-body systems we develop a formalism to straightforwardly compute bounds on the number of Trotter steps needed to accurately simulate the time evolution of fermionic lattice models based on the first-order commutator scaling. We apply this formalism to two closely related many-body models prominent in condensed matter physics, the Hubbard and t-J models. We find that, while a naive comparison of the Trotter depth first seems to favor the Hubbard model, careful consideration of the model parameters and the allowable error for accurate simulation leads to a substantial advantage in favor of the t-J model. These results and formalism set the stage for significant improvements in quantum simulation costs.
Digital quantum simulation of many-body dynamics relies on Trotterization to decompose the target time evolution into elementary quantum gates operating at a fixed equidistant time discretization. Recent advances have outlined protocols enabling more efficient adaptive Trotter protocols, which have been shown to exhibit a controlled error in the dynamics of local observables and correlation functions. However, it has remained open to which extent the errors on the actual generator of the dynamics, i.e., the target many-body Hamiltonian, remain controlled. Here, we propose to use quantum Hamiltonian learning to numerically obtain the effective Hamiltonian and apply it on the recently introduced ADA-Trotter algorithm as a concrete demonstration. Our key observation is that deviations from the target generator remain bounded on all simulation times. This result suggests that the ADA-Trotter not only generates reliable digital quantum simulation of local dynamics, but also controllably approximates the global quantum state of the target system. Our proposal is sufficiently general and readily applicable to other adaptive time-evolution algorithms.
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.
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.
Quantum algorithms for quantum dynamics simulations are traditionally based on implementing a Trotter approximation of the time-evolution operator. This approach typically relies on deep circuits and is therefore hampered by the substantial limitations of available noisy and near-term quantum hardware. On the other hand, variational quantum algorithms (VQAs) have become an indispensable alternative, enabling small-scale simulations on present-day hardware. However, despite the recent development of VQAs for quantum dynamics, a detailed assessment of their efficiency and scalability is yet to be presented. To fill this gap, we applied a VQA based on McLachlan's principle to simulate the dynamics of a spin-boson model subject to varying levels of realistic hardware noise as well as in different physical regimes, and discuss the algorithm's accuracy and scaling behavior as a function of system size. We observe a good performance of the variational approach used in combination with a general, physically motivated wave function ansatz, and compare it to the conventional first-order Trotter evolution. Finally, based on this, we make scaling predictions for the simulation of a classically intractable system. We show that, despite providing a clear reduction of quantum gate cost, the variational method in its current implementation is unlikely to lead to a quantum advantage for the solution of time-dependent problems.
In this thesis, we present three innovations that can be used to simulate quantum dynamics on both quantum and classical computers. Our primary innovation is an efficient quantum algorithm for simulating time-dependent Hamiltonian evolution of general input states on a quantum computer. Given conditions on the smoothness of the Hamiltonian, the complexity of the algorithm is close to linear in the evolution time, and therefore is comparable to algorithms for time-independent Hamiltonians. In addition, we show how the complexity can be reduced by optimizing the time steps. The complexity of the algorithm is quantified by calls to an oracle, which yields information about the Hamiltonian. In contrast to previous work, which allowed an oracle query to yield an arbitrary number of bits or qubits, we assign a cost for each bit or qubit accessed, revealing hitherto unnoticed simulation costs. We also account for discretization errors in the time and the representation of the Hamiltonian. This work is enabled by our development of a decomposition scheme based on Lie-Trotter- Suzuki product formulae for ordered operator exponentials. We show, by counterexample, that Lie-Trotter-Suzuki approximations may be lower order than expected when applied to problems that have singularities or discontinuous derivatives of appropriate order. We address this problem by presenting a set of criteria that is sufficient for the validity of these approximations, prove convergence and provide upper bounds for the approximation error. We also provide a new analysis of the evolution of quantum adiabatic systems in this thesis. Our work contains the tightest upper bounds that are currently known, and the first tight lower bounds, for the error in the adiabatic approximation applied to fixed Hamiltonians in the limit of slow evolution. We use our results to address the Marzlin- Sanders counterexample Hamiltonian. Specifically, we provide sufficiency conditions for the validity of the adiabatic approximation and demonstrate that the Marzlin-Sanders counterexample fails to satisfy them. We also find a new class of time-dependent Hamiltonians that are related to the Marzlin0Sanders counterexample Hamiltonian, but nonetheless obey the adiabatic approximation for sufficiently slow evolutions.