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:2412.10145 · Physical Review Letters

Roughening Dynamics of Interfaces in the Two-Dimensional Quantum Ising Model

Wladislaw Krinitsin, Niklas Tausendpfund, Matteo Rizzi, Markus Heyl, Markus Schmitt

Published May 21, 2025• 8 citations
arXivHTMLPDF

Abstract

The properties of interfaces are key to understanding the physics of matter. However, the study of quantum interface dynamics has remained an outstanding challenge. Here, we use large-scale tree tensor network simulations to identify the dynamical signature of an interface roughening transition within the ferromagnetic phase of the 2D quantum Ising model. For initial domain wall profiles we find extended prethermal plateaus for smooth interfaces, whereas above the roughening transition the domain wall decays quickly. Our results can be readily explored experimentally in Rydberg atomic systems.

DOI
10.1103/9bsk-x9rw

Related papers

Interface roughening in the 3-D Ising model with tensor networks

Atsushi Ueda, Lander Burgelman, Luca Tagliacozzo +1 more

Interfaces in three-dimensional many-body systems can exhibit rich phenomena beyond the corresponding bulk properties. In particular, they can fluctuate and give rise to massless low energy degrees of freedom even in the presence of a gapped bulk. In this work, we present the first tensor-network study of the paradigmatic interface roughening transition of the 3-D Ising model using highly asymmetric lattices that are infinite in the $(xy)$ direction and finite in $z$. By reducing the problem to an effective 2-D tensor network, we study how truncating the $z$ direction reshapes the physics of the interface. For a truncation based on open boundary conditions, we demonstrate that varying the interface width gives rise to either a $\mathbb{Z}_2$ symmetry breaking transition (for odd $L_z$) or a smooth crossover(for even $L_z$). For antiperiodic boundary conditions, we obtain an effective $\mathbb{Z}_q$ clock model description with $q=2L_z$ that exhibits an intermediate Luttinger liquid phase with an emergent $\U(1)$ symmetry.

arXiv:2601.07829•Jan 12, 2026•0 citations

Tree Tensor Network Simulation of Dynamical Quantum Phase Transitions in the 2D Transverse-Field Ising Model

Xiangyue Zhang, Dizhou Xie, Yongqiang Li

The discovery of dynamical quantum phase transitions (DQPTs) has fundamentally challenged the traditional view that phase transitions only occur in thermal equilibrium. Experimental platforms and 1D numerical methods, like matrix product states (MPS), have made great progress. However, exploring true 2D DQPTs remains difficult due to finite-size limitations and the geometric biases of quasi-1D cylinder mappings. Here, we bypass these limitations by deploying a tree tensor network (TTN) approach. This allows us to directly compute the quench dynamics of the transverse-field Ising model (TFIM) on an open 2D square lattice. Because the TTN architecture naturally mirrors 2D lattice connectivity, we can extract the global Loschmidt echo. Our simulations reveal that while deep quenches yield standard DQPTs, quenching within the ferromagnetic phase produces an anomalous dynamical response. In this regime, the rate function exhibits sharp non-analytic peaks even as the macroscopic order parameter maintains its initial sign. This decoupled behavior strongly indicates that local spin excitations drive 2D DQPTs, rather than the macroscopic domain-wall motions seen in 1D chains. These results provide a quantitative numerical baseline for understanding non-equilibrium quantum matter in higher dimensions.

Entropy•Apr 26, 2026•0 citations

Exploring the Relaxation Landscape of a 2D Quantum Magnet on a 256-Qubit Processor

Tiago Mendes-Santos, Joseph Vovrosh, Sergi Julià-Farré +33 more

How quantum matter relaxes far from equilibrium is a central open problem in many-body physics, and one for which analog quantum simulators are well positioned to move from confirming theory to discovering new physics. Here, we use a two-dimensional Rydberg atom array of 256 qubits to map the relaxation landscape of the two-dimensional transverse-field Ising model across its phase diagram. Beyond the expected rapid thermalization, we identify two further regimes. The first is a prethermal regime whose dynamics are governed by an effective XY model. The second, and most unexpected, is a crossover regime characterized by a slowdown in relaxation. This slowdown occurs precisely where state-of-the-art classical tensor-network methods lose control at late times, whereas the quantum simulation remains consistent across system sizes. These results establish Rydberg atom arrays as a platform for scientific discovery in nonequilibrium quantum many-body dynamics.

quant-pharXiv:2608.07178•Aug 10, 2026•0 citations

Superfluid-insulator and roughening transitions in domain walls

Şebnem Güneş Söyler, Barbara Capogrosso-Sansone, Nikolay V. Prokof'ev +1 more

We have performed quantum Monte Carlo simulations to investigate the superfluid behavior of one- and two-dimensional interfaces separating checkerboard solid domains. The system is described by the hard-core Bose-Hubbard Hamiltonian with nearest-neighbor interaction. In accordance with Burovski et al. [Phys. Rev. Lett. 94, 165301 (2005)] we find that (i) the interface remains superfluid in a wide range of interaction strength before it undergoes a superfluid-insulator transition; (ii) in one dimension, the transition is of the Kosterlitz-Thouless type and is accompanied by the roughening transition, driven by proliferation of charge-$1∕2$ quasiparticles; (iii) in two dimensions, the transition belongs to the three-dimensional $\mathrm{U}(1)$ universality class and the interface remains smooth. Similar phenomena are expected for domain walls in quantum antiferromagnets.

arXiv:0707.3930•Physical Review A•Oct 29, 2007•9 citations

Interface dynamics in the two-dimensional quantum Ising model

Federico Balducci, Andrea Gambassi, Alessio Lerose +2 more

In a recent paper [Phys. Rev. Lett. 129, 120601 (2022)], we have shown that the dynamics of interfaces, in the symmetry-broken phase of the two-dimensional ferromagnetic quantum Ising model, displays a robust form of ergodicity breaking. In this paper, we elaborate more on the issue. First, we discuss two classes of initial states on the square lattice, the dynamics of which is driven by complementary terms in the effective Hamiltonian and may be solved exactly: (a) Strips of consecutive neighboring spins aligned in the opposite direction of the surrounding spins and (b) a large class of initial states, characterized by the presence of a well-defined ``smooth'' interface separating two infinitely extended regions with oppositely aligned spins. The evolution of the latter states can be mapped onto that of an effective one-dimensional fermionic chain, which is integrable in the infinite-coupling limit. In this case, deep connections with noteworthy results in mathematics emerge, as well as with similar problems in classical statistical physics. We present a detailed analysis of the evolution of these interfaces both on the lattice and in a suitable continuum limit, including the interface fluctuations and the dynamics of entanglement entropy. Second, we provide analytical and numerical evidence supporting the conclusion that the observed nonergodicity---arising from Stark localization of the effective fermionic excitations---persists away from the infinite-Ising-coupling limit, and we highlight the presence of a timescale $T\ensuremath{\sim}{e}^{cLlnL}$ for the decay of a region of large linear size $L$. The implications of our work for the classic problem of the decay of a false vacuum are also discussed.

arXiv:2209.08992•Physical review. B./Physical review. B•Jan 23, 2023•16 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
The properties of interfaces are key to understanding the physics of matter. However, the study of quantum interface dynamics has remained an outstanding challenge. Here, we use large-scale tree tensor network simulations to identify the dynamical signature of an interface roughening transition within the ferromagnetic phase of the 2D quantum Ising model. For initial domain wall profiles we find extended prethermal plateaus for smooth interfaces, whereas above the roughening transition the domain wall decays quickly. Our results can be readily explored experimentally in Rydberg atomic systems.
Interfaces in three-dimensional many-body systems can exhibit rich phenomena beyond the corresponding bulk properties. In particular, they can fluctuate and give rise to massless low energy degrees of freedom even in the presence of a gapped bulk. In this work, we present the first tensor-network study of the paradigmatic interface roughening transition of the 3-D Ising model using highly asymmetric lattices that are infinite in the (xy)(xy)(xy) direction and finite in zzz. By reducing the problem to an effective 2-D tensor network, we study how truncating the zzz direction reshapes the physics of the interface. For a truncation based on open boundary conditions, we demonstrate that varying the interface width gives rise to either a Z2\mathbb{Z}_2Z2​ symmetry breaking transition (for odd LzL_zLz​) or a smooth crossover(for even LzL_zLz​). For antiperiodic boundary conditions, we obtain an effective Zq\mathbb{Z}_qZq​ clock model description with q=2Lzq=2L_zq=2Lz​ that exhibits an intermediate Luttinger liquid phase with an emergent $\U(1)$ symmetry.
The discovery of dynamical quantum phase transitions (DQPTs) has fundamentally challenged the traditional view that phase transitions only occur in thermal equilibrium. Experimental platforms and 1D numerical methods, like matrix product states (MPS), have made great progress. However, exploring true 2D DQPTs remains difficult due to finite-size limitations and the geometric biases of quasi-1D cylinder mappings. Here, we bypass these limitations by deploying a tree tensor network (TTN) approach. This allows us to directly compute the quench dynamics of the transverse-field Ising model (TFIM) on an open 2D square lattice. Because the TTN architecture naturally mirrors 2D lattice connectivity, we can extract the global Loschmidt echo. Our simulations reveal that while deep quenches yield standard DQPTs, quenching within the ferromagnetic phase produces an anomalous dynamical response. In this regime, the rate function exhibits sharp non-analytic peaks even as the macroscopic order parameter maintains its initial sign. This decoupled behavior strongly indicates that local spin excitations drive 2D DQPTs, rather than the macroscopic domain-wall motions seen in 1D chains. These results provide a quantitative numerical baseline for understanding non-equilibrium quantum matter in higher dimensions.
How quantum matter relaxes far from equilibrium is a central open problem in many-body physics, and one for which analog quantum simulators are well positioned to move from confirming theory to discovering new physics. Here, we use a two-dimensional Rydberg atom array of 256 qubits to map the relaxation landscape of the two-dimensional transverse-field Ising model across its phase diagram. Beyond the expected rapid thermalization, we identify two further regimes. The first is a prethermal regime whose dynamics are governed by an effective XY model. The second, and most unexpected, is a crossover regime characterized by a slowdown in relaxation. This slowdown occurs precisely where state-of-the-art classical tensor-network methods lose control at late times, whereas the quantum simulation remains consistent across system sizes. These results establish Rydberg atom arrays as a platform for scientific discovery in nonequilibrium quantum many-body dynamics.
We have performed quantum Monte Carlo simulations to investigate the superfluid behavior of one- and two-dimensional interfaces separating checkerboard solid domains. The system is described by the hard-core Bose-Hubbard Hamiltonian with nearest-neighbor interaction. In accordance with Burovski et al. [Phys. Rev. Lett. 94, 165301 (2005)] we find that (i) the interface remains superfluid in a wide range of interaction strength before it undergoes a superfluid-insulator transition; (ii) in one dimension, the transition is of the Kosterlitz-Thouless type and is accompanied by the roughening transition, driven by proliferation of charge-1∕21∕21∕2 quasiparticles; (iii) in two dimensions, the transition belongs to the three-dimensional U(1)\mathrm{U}(1)U(1) universality class and the interface remains smooth. Similar phenomena are expected for domain walls in quantum antiferromagnets.
In a recent paper [Phys. Rev. Lett. 129, 120601 (2022)], we have shown that the dynamics of interfaces, in the symmetry-broken phase of the two-dimensional ferromagnetic quantum Ising model, displays a robust form of ergodicity breaking. In this paper, we elaborate more on the issue. First, we discuss two classes of initial states on the square lattice, the dynamics of which is driven by complementary terms in the effective Hamiltonian and may be solved exactly: (a) Strips of consecutive neighboring spins aligned in the opposite direction of the surrounding spins and (b) a large class of initial states, characterized by the presence of a well-defined ``smooth'' interface separating two infinitely extended regions with oppositely aligned spins. The evolution of the latter states can be mapped onto that of an effective one-dimensional fermionic chain, which is integrable in the infinite-coupling limit. In this case, deep connections with noteworthy results in mathematics emerge, as well as with similar problems in classical statistical physics. We present a detailed analysis of the evolution of these interfaces both on the lattice and in a suitable continuum limit, including the interface fluctuations and the dynamics of entanglement entropy. Second, we provide analytical and numerical evidence supporting the conclusion that the observed nonergodicity---arising from Stark localization of the effective fermionic excitations---persists away from the infinite-Ising-coupling limit, and we highlight the presence of a timescale $T\ensuremath{\sim}{e}^{cLlnL}$ for the decay of a region of large linear size LLL. The implications of our work for the classic problem of the decay of a false vacuum are also discussed.
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.