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

arXiv:1809.09827 · Journal of the Physical Society of Japan

Quantum Annealing of Pure and Random Ising Chains Coupled to a Bosonic Environment

Sei Suzuki, Hiroki Oshiyama, Naokazu Shibata

Published Mar 1, 2019• 14 citations
arXivHTMLPDF

Abstract

We study quantum annealing of pure and random transverse Ising chains that are coupled to boson baths, using a novel numerical method based on the combination of the quasi-adiabatic propagation path integral (QUAPI) and the matrix product state (MPS) formalisms. We present numerical results on systems up to 64 spins and conclude that the baths disturb quantum annealing of pure and random Ising chains. The numerical method to compute the reduced density matrix is described in detail.

DOI
10.7566/jpsj.88.061003

Related papers

Numerically efficient quasi-adiabatic propagator path integral approach with two independent non-commuting baths

Roman Ovcharenko, Benjamin Philipp Fingerhut

Path integral methods, such as the quasi-adiabatic propagator path integral (QUAPI), are widely used in general-purpose and highly accurate numerical benchmark simulations of open quantum systems, particularly in regimes inaccessible to perturbative methods. Nevertheless, the applicability of the QUAPI method to realistic systems of interest is restricted by the exponentially growing computer memory requirements with respect to the size of the quantum system and the time range of non-Markovian correlation effects. This exponential "wall" becomes even more severe for multiple non-commuting fluctuating environments. In the present work, we address the numerical efficiency and accuracy of approximations that have been introduced for the QUAPI method with a single general environment, for the case of two independent non-commuting environments where one of them is considered as a pure dephasing environment. In particular, we consider a sharply defined cutoff of the memory time, path filtering, and mask assisted coarse graining of influence functional coefficients as approximations. We demonstrate that commonly applied numerical techniques, such as path filtering, cannot be straightforwardly transferred to the two-bath case even in the weak-coupling and quasi-Markovian limits. On the other hand, the sharply defined memory cutoff can be accurately handled with the mask assisted coarse graining approach. Our findings demonstrate that if system coupling operators to different baths do not commute, the additive nature of the statistically independent environments may be misleading. In particular, the quasi-Markovian nature of a pure dephasing bath is lost once there simultaneously exists another non-commuting source of fluctuations.

arXiv:2503.21693•The Journal of Chemical Physics•Jul 7, 2025•3 citations

Simulated quantum annealing as a simulator of nonequilibrium quantum dynamics

Yuki Bando, Hidetoshi Nishimori

Simulated quantum annealing based on the path-integral Monte Carlo is one of the most common tools to simulate quantum annealing on classical hardware. Nevertheless, it is in principle highly nontrivial whether or not this classical algorithm can correctly reproduce the quantum dynamics of quantum annealing, particularly in the diabatic regime. We study this problem numerically through the generalized Kibble-Zurek mechanism of defect distribution in the simplest ferromagnetic one-dimensional transverse-field Ising model with and without coupling to the environment. We find that, in the absence of coupling to the environment, simulated quantum annealing correctly describes the annealing-time dependence of the average number of defects, but a detailed analysis of the defect distribution shows clear deviations from the theoretical prediction. When the system is open (coupled to the environment), the average number of defects does not follow the theoretical prediction but is qualitatively compatible with the numerical result by the infinite-time-evolving block decimation combined with the quasiadiabatic propagator path integral, which is valid in a very-short-time region. The distribution of defects in the open system turns out to be not far from the theoretical prediction. It is surprising that the classical stochastic dynamics of simulated quantum annealing ostensibly reproduces some aspects of the quantum dynamics. However, a serious problem is that it is hard to predict for which physical quantities in which system it is reliable. Those results suggest the necessity to exert a good amount of caution in using simulated quantum annealing to study the detailed quantitative aspects of the dynamics of quantum annealing.

arXiv:2106.00928•Physical Review A•Aug 16, 2021•32 citations

Classical Simulation and Theory of Quantum Annealing in a Thermal Environment

Hiroki Oshiyama, Sei Suzuki, Naokazu Shibata

We study quantum annealing in the quantum Ising model coupled to a thermal environment. When the speed of quantum annealing is sufficiently slow, the system evolves following the instantaneous thermal equilibrium. This quasistatic and isothermal evolution, however, fails near the end of annealing because the relaxation time grows infinitely, therefore yielding excess energy from the thermal equilibrium. We develop a phenomenological theory based on this picture and derive a scaling relation of the excess energy after annealing. The theoretical results are numerically confirmed using a novel non-Markovian method that we recently proposed based on a path-integral representation of the reduced density matrix and the infinite time evolving block decimation. In addition, we discuss crossovers from weak to strong coupling as well as from the adiabatic to quasistatic regime, and propose experiments on the D-Wave quantum annealer.

arXiv:2102.02570•Physical Review Letters•Apr 27, 2022•11 citations

Optimal working point in dissipative quantum annealing

Luca Arceci, Simone Barbarino, Davide Rossini +1 more

We study the effect of a thermal environment on the quantum annealing dynamics of a transverse-field Ising chain. The environment is modeled as a single Ohmic bath of quantum harmonic oscillators weakly interacting with the total transverse magnetization of the chain in a translationally invariant manner. We show that the density of defects generated at the end of the annealing process displays a minimum as a function of the annealing time, the so-called optimal working point, only in rather special regions of the bath temperature and coupling strength plane. We discuss the relevance of our results for current and future experimental implementations with quantum annealing hardware.

arXiv:1804.04251•Physical review. B./Physical review. B•Aug 27, 2018•27 citations

Numerical and analytical approach to the quantum dynamics of two coupled spins in bosonic baths

Alessandro Sergi, Ilya Sinayskiy, Francesco Petruccione

The quantum dynamics of a spin chain interacting with multiple bosonic baths is described in a mixed Wigner-Heisenberg representation. The formalism is illustrated by simulating the time evolution of the reduced density matrix of two coupled spins, where each spin is also coupled to its own bath of harmonic oscillators. In order to prove the validity of the approach, an analytical solution in the Born-Markov approximation is found. The agreement between the two methods is shown.

arXiv:0903.5181•Physical Review A•Jul 15, 2009•27 citations

Simulating dynamics of the two-dimensional transverse-field Ising model: a comparative study of large-scale classical numerics

Joseph Vovrosh, Sergi Julià-Farré, Wladislaw Krinitsin +11 more

The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.

quant-pharXiv:2511.19340•Physical Review Research•Nov 25, 2025•1 citation
We study quantum annealing of pure and random transverse Ising chains that are coupled to boson baths, using a novel numerical method based on the combination of the quasi-adiabatic propagation path integral (QUAPI) and the matrix product state (MPS) formalisms. We present numerical results on systems up to 64 spins and conclude that the baths disturb quantum annealing of pure and random Ising chains. The numerical method to compute the reduced density matrix is described in detail.
Path integral methods, such as the quasi-adiabatic propagator path integral (QUAPI), are widely used in general-purpose and highly accurate numerical benchmark simulations of open quantum systems, particularly in regimes inaccessible to perturbative methods. Nevertheless, the applicability of the QUAPI method to realistic systems of interest is restricted by the exponentially growing computer memory requirements with respect to the size of the quantum system and the time range of non-Markovian correlation effects. This exponential "wall" becomes even more severe for multiple non-commuting fluctuating environments. In the present work, we address the numerical efficiency and accuracy of approximations that have been introduced for the QUAPI method with a single general environment, for the case of two independent non-commuting environments where one of them is considered as a pure dephasing environment. In particular, we consider a sharply defined cutoff of the memory time, path filtering, and mask assisted coarse graining of influence functional coefficients as approximations. We demonstrate that commonly applied numerical techniques, such as path filtering, cannot be straightforwardly transferred to the two-bath case even in the weak-coupling and quasi-Markovian limits. On the other hand, the sharply defined memory cutoff can be accurately handled with the mask assisted coarse graining approach. Our findings demonstrate that if system coupling operators to different baths do not commute, the additive nature of the statistically independent environments may be misleading. In particular, the quasi-Markovian nature of a pure dephasing bath is lost once there simultaneously exists another non-commuting source of fluctuations.
Simulated quantum annealing based on the path-integral Monte Carlo is one of the most common tools to simulate quantum annealing on classical hardware. Nevertheless, it is in principle highly nontrivial whether or not this classical algorithm can correctly reproduce the quantum dynamics of quantum annealing, particularly in the diabatic regime. We study this problem numerically through the generalized Kibble-Zurek mechanism of defect distribution in the simplest ferromagnetic one-dimensional transverse-field Ising model with and without coupling to the environment. We find that, in the absence of coupling to the environment, simulated quantum annealing correctly describes the annealing-time dependence of the average number of defects, but a detailed analysis of the defect distribution shows clear deviations from the theoretical prediction. When the system is open (coupled to the environment), the average number of defects does not follow the theoretical prediction but is qualitatively compatible with the numerical result by the infinite-time-evolving block decimation combined with the quasiadiabatic propagator path integral, which is valid in a very-short-time region. The distribution of defects in the open system turns out to be not far from the theoretical prediction. It is surprising that the classical stochastic dynamics of simulated quantum annealing ostensibly reproduces some aspects of the quantum dynamics. However, a serious problem is that it is hard to predict for which physical quantities in which system it is reliable. Those results suggest the necessity to exert a good amount of caution in using simulated quantum annealing to study the detailed quantitative aspects of the dynamics of quantum annealing.
We study quantum annealing in the quantum Ising model coupled to a thermal environment. When the speed of quantum annealing is sufficiently slow, the system evolves following the instantaneous thermal equilibrium. This quasistatic and isothermal evolution, however, fails near the end of annealing because the relaxation time grows infinitely, therefore yielding excess energy from the thermal equilibrium. We develop a phenomenological theory based on this picture and derive a scaling relation of the excess energy after annealing. The theoretical results are numerically confirmed using a novel non-Markovian method that we recently proposed based on a path-integral representation of the reduced density matrix and the infinite time evolving block decimation. In addition, we discuss crossovers from weak to strong coupling as well as from the adiabatic to quasistatic regime, and propose experiments on the D-Wave quantum annealer.
We study the effect of a thermal environment on the quantum annealing dynamics of a transverse-field Ising chain. The environment is modeled as a single Ohmic bath of quantum harmonic oscillators weakly interacting with the total transverse magnetization of the chain in a translationally invariant manner. We show that the density of defects generated at the end of the annealing process displays a minimum as a function of the annealing time, the so-called optimal working point, only in rather special regions of the bath temperature and coupling strength plane. We discuss the relevance of our results for current and future experimental implementations with quantum annealing hardware.
The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.