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:2607.05178

Efficient classical simulation of two-dimensional long-range systems: Rydberg arrays and beyond

Jia-Lin Chan, Tao Xiang, Yantao Wu

quant-phcond-mat.quant-gascond-mat.str-elSubmitted Jul 6, 2026• Announced Jul 7, 2026• 0 citations
arXivHTMLPDF

Abstract

In variational Monte Carlo (VMC) calculations of $N$-site quantum systems with arbitrary all-to-all two-body interactions, evaluating the local energy generally costs $O(N^3)$. We introduce a new framework that reduces this cost to $O(N)$ for tensor network states, capable of scalable and accurate computation of real-time dynamics and ground states. As a result, we obtain accurate simulations of the adiabatic real-time protocol of a $10\times10$ dipolar XY model realized in a Rydberg simulator [C. Chen et al., Nature 616, 691 (2023)], which was previously beyond the reach of classical simulation. Going beyond quantum experiments, we also directly perform ground state VMC to compare with the adiabatic state preparation. Our work demonstrates tensor network VMC as a powerful classical simulator for long-range quantum platforms such as Rydberg and ion-trap simulators, which are currently in urgent need of scalable classical benchmarking tools. As a separate technical contribution, we resolve the pathology of evolving from product states within of tensor network VMC.

Related papers

Variational approach to open quantum systems with long-range competing interactions

Dawid A. Hryniuk, Marzena Hanna Szymańska

Competition between short- and long-range interactions underpins many emergent phenomena in nature. Despite rapid progress in their experimental control, computational methods capable of accurately simulating open quantum many-body systems with complex long-ranged interactions at scale remain scarce. Here, we address this limitation by introducing an efficient and scalable approach to dissipative quantum lattices in one and two dimensions, combining matrix product operators and time-dependent variational Monte Carlo. We showcase the versatility, effectiveness, and unique methodological advantages of our algorithm by simulating the non-equilibrium dynamics and steady states of spin-$$\frac{1}{2}$$ lattices with competing algebraically-decaying interactions for as many as N = 200 sites, revealing the emergence of spatially-modulated magnetic order far from equilibrium. This approach offers promising prospects for advancing our understanding of the complex non-equilibrium properties of a diverse variety of experimentally-realizable quantum systems with long-ranged interactions, including Rydberg atoms, ultracold dipolar molecules, and trapped ions. Far-from-equilibrium many-body quantum systems with complex non-local interactions are of significant fundamental and practical importance. Here the authors demonstrate how variational Monte Carlo in combination with compact tensor network parameterization of the density matrix enables accurate and scalable simulation of such open quantum systems.

arXiv:2510.01543•Communications Physics•Jan 10, 2026•1 citation

Real-Time Dynamics in Two Dimensions with Tensor Network States via Time-Dependent Variational Monte Carlo Method

Anonymous, Jannes Nys

Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned tVMC equation, yielding robust solutions using the efficient Cholesky decomposition. This enables long-time evolution in previously inaccessible regimes. We explain how the difficulties in the traditional approach, particularly those associated with gauge redundancies, are resolved within tVMC. We demonstrate the power and generality of the method through five representative real-time local quench dynamics in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) vison propagation in a pure Z 2 gauge theory; (III) vison confinement dynamics in a Z 2 lattice gauge theory coupled to Higgs field; (IV) fractionalized charge transport in a fractional Chern insulator; and (V) superfluidity and critical velocity in interacting bosons. All simulations are performed on ≥ 10 × 10 lattices with evolution times beyond T = 10 using modest computational resources. Where exact benchmarks exist (case I), PEPS-tVMC matches free-fermion dynamics with high accuracy up to T = 12 . In addition, we also simulate the paradigmatic dynamics of the Ising model following a global quench at the critical transverse field, and obtain, with modest bond dimension, agreement with previous results. These results demonstrate PEPS-tVMC as a practical and versatile tool for real-time quantum dynamics in two dimensions. The method significantly extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems in real time and provides a valuable computational counterpart to emerging quantum simulators. As a by-product in the development, we also present a new form of minSR, which is more stable and offers a new perspective on tVMC.

PRX Quantum•Jul 8, 2026•1 citation

Real-Time Dynamics in Two Dimensions with Tensor Network States via Time-Dependent Variational Monte Carlo

Yantao Wu, Jannes Nys

Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned tVMC equation. This enables long-time evolution in previously inaccessible regimes. We explain how the difficulties in the traditional approach, particularly those associated with gauge redundancies, are resolved within tVMC. We demonstrate the power and generality of the method through five representative real-time local quench dynamics in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) vison propagation in a pure Z2 gauge theory; (III) vison confinement dynamics in a Z2 lattice gauge theory coupled to Higgs field; (IV) fractionalized charge transport in a fractional Chern insulator; and (V) superfluidity and critical velocity in interacting bosons. All simulations are performed on >= 10 x 10 lattices with evolution times beyond T = 10 using modest computational resources. In addition, we also simulate the paradigmatic dynamics of the Ising model following a global quench at the critical transverse field, and obtain, with modest bond dimension, agreement with previous results. The method significantly extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems in real-time and provides a valuable computational counterpart to emerging quantum simulators. As a by-product in the development, we also present a new form of minSR, which is more stable and offers a new perspective on tVMC.

cond-mat.str-elCross-listedarXiv:2512.06768•Dec 9, 2025•0 citations

Variational Monte Carlo with large patched transformers

Kyle Sprague, Stefanie Czischek

Abstract Large language models, like transformers, have recently demonstrated immense powers in text and image generation. This success is driven by the ability to capture long-range correlations between elements in a sequence. The same feature makes the transformer a powerful wavefunction ansatz that addresses the challenge of describing correlations in simulations of qubit systems. Here we consider two-dimensional Rydberg atom arrays to demonstrate that transformers reach higher accuracies than conventional recurrent neural networks for variational ground state searches. We further introduce large, patched transformer models, which consider a sequence of large atom patches, and show that this architecture significantly accelerates the simulations. The proposed architectures reconstruct ground states with accuracies beyond state-of-the-art quantum Monte Carlo methods, allowing for the study of large Rydberg systems in different phases of matter and at phase transitions. Our high-accuracy ground state representations at reasonable computational costs promise new insights into general large-scale quantum many-body systems.

arXiv:2306.03921•Communications Physics•Mar 11, 2024•32 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

Data-enhanced variational Monte Carlo simulations for Rydberg atom arrays

Stefanie Czischek, M. Schuyler Moss, Matthew Radzihovsky +2 more

Rydberg atom arrays are programmable quantum simulators capable of preparing interacting qubit systems in a variety of quantum states. Due to long experimental preparation times, obtaining projective measurement data can be relatively slow for large arrays, which poses a challenge for state reconstruction methods such as tomography. Today, novel ground-state wave-function Ans\"atze like recurrent neural networks (RNNs) can be efficiently trained not only from projective measurement data, but also through Hamiltonian-guided variational Monte Carlo (VMC). In this paper, we demonstrate how pretraining modern RNNs on even small amounts of data significantly reduces the convergence time for a subsequent variational optimization of the wave function. This suggests that essentially any amount of measurements obtained from a state prepared in an experimental quantum simulator could provide significant values for neural-network-based VMC strategies.

arXiv:2203.04988•Physical review. B./Physical review. B•May 9, 2022•22 citations
In variational Monte Carlo (VMC) calculations of NNN-site quantum systems with arbitrary all-to-all two-body interactions, evaluating the local energy generally costs O(N3)O(N^3)O(N3). We introduce a new framework that reduces this cost to O(N)O(N)O(N) for tensor network states, capable of scalable and accurate computation of real-time dynamics and ground states. As a result, we obtain accurate simulations of the adiabatic real-time protocol of a 10×1010\times1010×10 dipolar XY model realized in a Rydberg simulator [C. Chen et al., Nature 616, 691 (2023)], which was previously beyond the reach of classical simulation. Going beyond quantum experiments, we also directly perform ground state VMC to compare with the adiabatic state preparation. Our work demonstrates tensor network VMC as a powerful classical simulator for long-range quantum platforms such as Rydberg and ion-trap simulators, which are currently in urgent need of scalable classical benchmarking tools. As a separate technical contribution, we resolve the pathology of evolving from product states within of tensor network VMC.
Competition between short- and long-range interactions underpins many emergent phenomena in nature. Despite rapid progress in their experimental control, computational methods capable of accurately simulating open quantum many-body systems with complex long-ranged interactions at scale remain scarce. Here, we address this limitation by introducing an efficient and scalable approach to dissipative quantum lattices in one and two dimensions, combining matrix product operators and time-dependent variational Monte Carlo. We showcase the versatility, effectiveness, and unique methodological advantages of our algorithm by simulating the non-equilibrium dynamics and steady states of spin-12\frac{1}{2}21​ lattices with competing algebraically-decaying interactions for as many as N = 200 sites, revealing the emergence of spatially-modulated magnetic order far from equilibrium. This approach offers promising prospects for advancing our understanding of the complex non-equilibrium properties of a diverse variety of experimentally-realizable quantum systems with long-ranged interactions, including Rydberg atoms, ultracold dipolar molecules, and trapped ions. Far-from-equilibrium many-body quantum systems with complex non-local interactions are of significant fundamental and practical importance. Here the authors demonstrate how variational Monte Carlo in combination with compact tensor network parameterization of the density matrix enables accurate and scalable simulation of such open quantum systems.
Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned tVMC equation, yielding robust solutions using the efficient Cholesky decomposition. This enables long-time evolution in previously inaccessible regimes. We explain how the difficulties in the traditional approach, particularly those associated with gauge redundancies, are resolved within tVMC. We demonstrate the power and generality of the method through five representative real-time local quench dynamics in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) vison propagation in a pure Z 2 gauge theory; (III) vison confinement dynamics in a Z 2 lattice gauge theory coupled to Higgs field; (IV) fractionalized charge transport in a fractional Chern insulator; and (V) superfluidity and critical velocity in interacting bosons. All simulations are performed on ≥ 10 × 10 lattices with evolution times beyond T = 10 using modest computational resources. Where exact benchmarks exist (case I), PEPS-tVMC matches free-fermion dynamics with high accuracy up to T = 12 . In addition, we also simulate the paradigmatic dynamics of the Ising model following a global quench at the critical transverse field, and obtain, with modest bond dimension, agreement with previous results. These results demonstrate PEPS-tVMC as a practical and versatile tool for real-time quantum dynamics in two dimensions. The method significantly extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems in real time and provides a valuable computational counterpart to emerging quantum simulators. As a by-product in the development, we also present a new form of minSR, which is more stable and offers a new perspective on tVMC.
Reliably simulating two-dimensional many-body quantum dynamics with projected entangled pair states (PEPS) has long been a difficult challenge. In this work, we overcome this barrier for low-energy quantum dynamics by developing a stable and efficient time-dependent variational Monte Carlo (tVMC) framework for PEPS. By analytically removing all gauge redundancies of the PEPS manifold and exploiting tensor locality, we obtain a numerically well-conditioned tVMC equation. This enables long-time evolution in previously inaccessible regimes. We explain how the difficulties in the traditional approach, particularly those associated with gauge redundancies, are resolved within tVMC. We demonstrate the power and generality of the method through five representative real-time local quench dynamics in two dimensions: (I) chiral edge propagation in a free-fermion Chern insulator; (II) vison propagation in a pure Z2 gauge theory; (III) vison confinement dynamics in a Z2 lattice gauge theory coupled to Higgs field; (IV) fractionalized charge transport in a fractional Chern insulator; and (V) superfluidity and critical velocity in interacting bosons. All simulations are performed on >= 10 x 10 lattices with evolution times beyond T = 10 using modest computational resources. In addition, we also simulate the paradigmatic dynamics of the Ising model following a global quench at the critical transverse field, and obtain, with modest bond dimension, agreement with previous results. The method significantly extends the reach of classical tensor-network simulations for studying elementary excitations in quantum many-body systems in real-time and provides a valuable computational counterpart to emerging quantum simulators. As a by-product in the development, we also present a new form of minSR, which is more stable and offers a new perspective on tVMC.
Abstract Large language models, like transformers, have recently demonstrated immense powers in text and image generation. This success is driven by the ability to capture long-range correlations between elements in a sequence. The same feature makes the transformer a powerful wavefunction ansatz that addresses the challenge of describing correlations in simulations of qubit systems. Here we consider two-dimensional Rydberg atom arrays to demonstrate that transformers reach higher accuracies than conventional recurrent neural networks for variational ground state searches. We further introduce large, patched transformer models, which consider a sequence of large atom patches, and show that this architecture significantly accelerates the simulations. The proposed architectures reconstruct ground states with accuracies beyond state-of-the-art quantum Monte Carlo methods, allowing for the study of large Rydberg systems in different phases of matter and at phase transitions. Our high-accuracy ground state representations at reasonable computational costs promise new insights into general large-scale quantum many-body systems.
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.
Rydberg atom arrays are programmable quantum simulators capable of preparing interacting qubit systems in a variety of quantum states. Due to long experimental preparation times, obtaining projective measurement data can be relatively slow for large arrays, which poses a challenge for state reconstruction methods such as tomography. Today, novel ground-state wave-function Ans\"atze like recurrent neural networks (RNNs) can be efficiently trained not only from projective measurement data, but also through Hamiltonian-guided variational Monte Carlo (VMC). In this paper, we demonstrate how pretraining modern RNNs on even small amounts of data significantly reduces the convergence time for a subsequent variational optimization of the wave function. This suggests that essentially any amount of measurements obtained from a state prepared in an experimental quantum simulator could provide significant values for neural-network-based VMC strategies.