The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L=12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.
The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L=12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.
We investigate how to simulate a quantum many body system using near term device known as noisy intermediate scale quantum (NISQ) hardware. We explore the two different types of hardware, analog and digital device, and discuss their potential for quantum many body simulation. For instance, one of the applications is to use the analog device to simulate the slow dynamics under a strong staggered field. However, noise is a major challenge for quantum hardware. Hence, we discuss the quantum error mitigation methods which are suitable for NISQ hardware. We apply the most commonly-used quantum error mitigation techniques to our many body simulation. We simulate quantum body scar dynamics using zero noise extrapolation techniques, Pauli twirling, and dynamical decoupling. We also calculate the corresponding correlation function using direct measurement without resorting to ancilla qubits. Furthermore, we use the same error mitigation techniques to explore dynamics of energy correlators. As additional error mitigation technique, we divide the noisy correlation results by their sum, which is ideally constant due to energy conservation. By doing so, we significantly improve our result. Finally, we consider another NISQ friendly quantum algorithm for imaginary time evolution. With enough samples from minimally entangled typical thermal state algorithms, we can obtain the thermal expectation value of observable and explore their thermodynamics. In chapter 2, we provide the background knowledge for the following 4 chapters. We introduce the rydberg atom system as the analog system and their special parameter regime for interesting physics. Moreover, we introduce the digital quantum hardware and trotterization that enable one to simulate many body dynamics on the quantum hardware. We also introduce the quantum error mitigation techniques to reduce the the impact of noise from the results. Finally, we introduce the variation quantum dynamics method for the dynamics simulation. In chapter 3, we explore the Rydberg blockade regime with a strong staggered field. With the strong staggered field, any computational basis states are thermalized. Here, we investigate the slow dynamics using a Schrieffer-Wolff (SW) transformation, which is similar to perturbation theory in calculating the different-order effective Hamiltonian in terms of energy scale. As numerical result shown, the naive Schrieffer-Wolff transformation is not enough to explain the some dynamics' relaxation time. Hence, here, we also consider another conservative quantity to perform the SW transformation. The corresponding analytical result agrees with the numerical result. In chapter 4, we simulate quantum many body scarred dynamics in IBM quantum superconducting hardware. We also use custom pulse gates, quantum error mitigation techniques, and quantum error suppression methods to help us handle the noise. The error mitigated result outperforms the raw result. With the all techniques, we also explore non-trivial correlation function which also manifest the long lived oscillation in the rydberg blockade regime. In chapter 5, we use the same quantum error mitigation techniques, which we used in ch. 5, to study the energy transport properties of the mixed Ising model at infinite temperature. Here, instead of using some random state to evaluate the necessary two point energy density correlations, we use Y basis product state sampling to evaluate them. Surprisingly, few samples are enough to estimate the correlation. In addition to the quantum error mitigation techniques, we also use the fact that sum of all correlation functions is constant to renormalize the noisy results. With this trick, we are able to extract an effective dynamical exponent, which enables us to classify the transport properties. In chapter 6, we combine a variational quantum algorithm and the minimally entangled typical thermal states approach to investigate the thermodynamics of 1+1-dimensional Z2 lattice gauge theory. In particular, with a confining field, the 1+1-dimensional Z2 lattice gauge theory cannot be solved exactly. We estimate the corresponding thermal particle number density and energy density to distinguish the confined and deconfined regime. Chapter 7 provides a short review of the main results described in the thesis and a discussion of future research directions that they motivate. We mainly give insight into the recent experimental breakthrough in quantum error correction. We also clarify that quantum error mitigation methods do not become totally useless when a fault-tolerant quantum computer is realized. Finally, we also comment on variational quantum algorithms and their potential applications.
Quantum mechanical problems are among the hardest to simulate and, in some cases, remain intractable even for the most powerful computers. Quantum computing has emerged as a new technological platform to address such challenges, with rapid advances in recent years. Yet, current quantum devices remain noisy and limited in scale. Hence, it is essential to identify classically hard problems of practical interest and tractable with existing quantum devices. Among potential applications, the real-time simulation of quantum systems is one of the most promising to deliver an early, practical quantum advantage. This doctoral thesis is therefore centered around simulating quantum dynamics on quantum devices. We first present an overview of the most relevant quantum algorithms for quantum dynamics, highlighting respective advantages and limitations. Further, we identify relevant problems within quantum dynamics that could benefit from quantum simulation in the near future. Second, we propose a method for benchmarking hardware and error mitigation algorithms that is based on well-understood theoretical results and suffers from no scaling issues. The resulting quality metric is intuitive and transferable to other applications. We successfully implement the scheme on up to 133 qubits, demonstrating coherent evolution up to a two-qubit gate depth of 28, featuring 1396 two-qubit gates, before noise becomes prevalent. Third, we propose a novel variant of probabilistic error amplification to implement open quantum dynamics, relying on characterizing and altering hardware noise to mimic the system-environment interaction under study. Lastly, we present two studies on state preparation and phase classification: first, a hybrid algorithm to prepare ground states of electron-phonon systems; second, a quantum machine learning-based approach to distinguish phases in previously prepared quantum states.
Abstract Hamiltonian simulation is one of the main proposed applications of near-term quantum computers, but its practical value depends on whether real hardware can preserve intended Hamiltonian dynamics once sampling limits, gate noise, decoherence, and circuit depth are taken into account. This study investigates whether the IBM quantum backend \texttt{ibm\_marrakesh} can reproduce the many-body dynamics shown by ideal simulation, using a 15-qubit one-dimensional tilted-field Ising chain as the benchmark system. The system was initialized with all qubits spin-up except for a localized spin-down excitation at the central qubit, q8, and evolved using a second-order Trotterized circuit with timestep Δt=0.1. The reference case, hz=0, was compared with the tilted-field case, hz=hx=0.35625, across three interaction regimes: J=0.90hx, J=hx, and J=1.40hx. Local magnetization, total magnetization, and connected spin correlations were used to track magnetization memory and finite-speed correlation spreading from the central excitation. The ideal simulation shows clear local and total magnetization evolution, as well as connected-correlation spreading outward from q8, giving a baseline for the intended Hamiltonian dynamics. The IBM hardware results reproduce the broad early-time magnetization trend but show compressed local-magnetization values, reduced separation between qubits, and weaker long-range connected-correlation structure. These differences indicate that the open-access quantum hardware preserves some early-time features of the intended dynamics, while longer-time correlation spreading is increasingly obscured by finite sampling, two-qubit gate errors, readout error, decoherence, and circuit depth. Overall, the study uses this spin-chain system as a direct benchmark for testing how much physically meaningful many-body behavior can still be extracted from noisy open-access quantum hardware.
Abstract Hamiltonian simulation is one of the main proposed applications of near-term quantum computers, but its practical value depends on whether real hardware can preserve intended Hamiltonian dynamics once sampling limits, gate noise, decoherence, and circuit depth are taken into account. This study investigates whether the IBM quantum backend \texttt{ibm\_marrakesh} can reproduce the many-body dynamics shown by ideal simulation, using a 15-qubit one-dimensional tilted-field Ising chain as the benchmark system. The system was initialized with all qubits spin-up except for a localized spin-down excitation at the central qubit, q8, and evolved using a second-order Trotterized circuit with timestep Δt=0.1. The reference case, hz=0, was compared with the tilted-field case, hz=hx=0.35625, across three interaction regimes: J=0.90hx, J=hx, and J=1.40hx. Local magnetization, total magnetization, and connected spin correlations were used to track magnetization memory and finite-speed correlation spreading from the central excitation. The ideal simulation shows clear local and total magnetization evolution, as well as connected-correlation spreading outward from q8, giving a baseline for the intended Hamiltonian dynamics. The IBM hardware results reproduce the broad early-time magnetization trend but show compressed local-magnetization values, reduced separation between qubits, and weaker long-range connected-correlation structure. These differences indicate that the open-access quantum hardware preserves some early-time features of the intended dynamics, while longer-time correlation spreading is increasingly obscured by finite sampling, two-qubit gate errors, readout error, decoherence, and circuit depth. Overall, the study uses this spin-chain system as a direct benchmark for testing how much physically meaningful many-body behavior can still be extracted from noisy open-access quantum hardware.
Quantum many-body scars are an intriguing dynamical regime in which quantum systems exhibit coherent dynamics and long-range correlations when prepared in certain initial states. We use this combination of coherence and many-body correlations to benchmark the performance of present-day quantum computing devices by using them to simulate the dynamics of an antiferromagnetic initial state in mixed-field Ising chains of up to 19 sites. In addition to calculating the dynamics of local observables, we also calculate the Loschmidt echo and a nontrivial unequal-time connected correlation function that witnesses long-range many-body correlations in the scarred dynamics. We find coherent dynamics to persist over up to 40 Trotter steps even in the presence of various sources of error. To obtain these results, we leverage a variety of error-mitigation techniques including noise tailoring, zero-noise extrapolation, dynamical decoupling, and physically motivated postselection of measurement results. Crucially, we also find that using pulse-level control to implement the Ising interaction yields a substantial improvement over the standard controlled-not-based compilation of this interaction. Our results demonstrate the power of error-mitigation techniques and pulse-level control to probe many-body coherence and correlation effects on present-day quantum hardware.
In a recent milestone experiment, Google's processor Sycamore heralded the era of "quantum supremacy" by sampling from the output of (pseudo-)random circuits. We show that such random circuits provide tailor-made building blocks for simulating quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. Specifically, we propose an algorithm consisting of a random circuit followed by a trotterized Hamiltonian time evolution to study hydrodynamics and to extract transport coefficients in the linear response regime. We numerically demonstrate the algorithm by simulating the buildup of spatiotemporal correlation functions in one- and two-dimensional quantum spin systems, where we particularly scrutinize the inevitable impact of errors present in any realistic implementation. Importantly, we find that the hydrodynamic scaling of the correlations is highly robust with respect to the size of the Trotter step, which opens the door to reach nontrivial time scales with a small number of gates. While errors within the random circuit are shown to be irrelevant, we furthermore unveil that meaningful results can be obtained for noisy time evolutions with error rates achievable on near-term hardware. Our work emphasizes the practical relevance of random circuits on NISQ devices beyond the abstract sampling task.