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:2511.02809 · Physical review. B./Physical review. B

Majorana string simulation of nonequilibrium dynamics in two-dimensional lattice fermion systems

Matteo D'Anna, Jannes Nys, Juan Carrasquilla

cond-mat.quant-gascond-mat.str-elquant-phSubmitted Nov 4, 2025• Announced Nov 5, 2025• 0 citations
v1 · Nov 4, 2025v2 · Jul 14, 2026
arXivHTMLPDF

Abstract

The study of real-time dynamics of fermions remains one of the last frontiers beyond the reach of classical simulations and is key to our understanding of quantum behavior in chemistry and materials, with implications for quantum technology. Here we introduce a Heisenberg-picture algorithm that propagates observables expressed in a Majorana-string basis using a truncation scheme that preserves Trotter accuracy and aims at maintaining computational efficiency. The framework is exact for quadratic Hamiltonians--remaining restricted to a fixed low-weight sector determined by the physical observable--admits variational initial states, and can be extended to interacting regimes via systematically controlled truncations. We benchmark our approach on one- and two-dimensional Fermi-Hubbard quenches, comparing against tensor network methods (MPS and fPEPS) and recent experimental data. The method achieves high accuracy on timescales comparable to state-of-the-art variational techniques and experiments, demonstrating that controlled Majorana-string truncation is a practical tool for simulating two-dimensional fermionic dynamics.

DOI
10.1103/f6mn-jc3r

Related papers

Fast convergence of Majorana Propagation for weakly interacting fermions

Giorgio Facelli, Hamza Fawzi, Omar Fawzi

Simulating the time dynamics of an observable under Hamiltonian evolution is one of the most promising candidates for quantum advantage as we do not expect efficient classical algorithms for this problem except in restricted settings. Here, we introduce such a setting by showing that Majorana Propagation, a simple algorithm combining Trotter steps and truncations, efficiently finds a low-degree approximation of the time-evolved observable as soon as such an approximation exists. This provides the first provable guarantee about Majorana Propagation for Hamiltonian evolution. As an application of this result, we prove that Majorana Propagation can efficiently simulate the time dynamics of any sparse quartic Hamiltonian up to time $t_{\text{max}}(u)$ depending on the interaction strength $u$. For a time horizon $t \leq t_{\text{max}}(u)$, the runtime of the algorithm is $N^{O(\log(t/\varepsilon))}$ where $N$ is the number of Majorana modes and $\varepsilon$ is the error measured in the normalized Frobenius norm. Importantly, in the limit of small $u$, $t_{\text{max}}(u)$ goes to $+\infty$, formalizing the intuition that the algorithm is accurate at all times when the Hamiltonian is quadratic.

quant-pharXiv:2601.05226•Jan 9, 2026•0 citations

Efficient Classical Simulation of Weakly Interacting Fermion Dynamics

Chu Zhao, Iman Marvian, Yu Tong

We consider the task of simulating the real-time dynamics of weakly interacting fermionic systems. In particular, we focus on computing the expectation value of a local observable $A$ at time $t$. By analyzing the convergence of the perturbative expansion in the interaction strength $λ$ for the Heisenberg-picture observable, we propose a polynomial-time algorithm for estimating this expectation value in the weakly interacting regime $λ|t|^{2D+1}=\mathcal{O}(1)$, when the Hamiltonian is geometrically local on a $D$-dimensional lattice. Importantly, this condition is independent of the system size. If the goal is instead to approximate the time-evolved observable in normalized Frobenius norm, we extend the convergence regime to $λ|t|=\mathcal{O}(1)$ with quasi-polynomial runtime. When the non-interacting part exhibits Anderson localization, our polynomial-time algorithm can be extended up to $λ|t|=\mathcal{O}(1)$, modulo polylogarithmic factors. Our algorithm brings together ideas from continuous-time QMC, diagrammatic QMC, and Majorana Propagation, but with a new Heisenberg-picture operator-growth analysis that makes the sampling complexity rigorously controllable. This leads to provably efficient classical algorithms in regimes where the interaction is weak enough that the sampling variance remains bounded independently of system size. Together, these results identify broad regimes in which weak interactions, locality, and localization can be leveraged to make real-time fermionic dynamics classically tractable.

quant-pharXiv:2608.19448•Aug 21, 2026•0 citations

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

Lowering the temperature of two-dimensional fermionic tensor networks with cluster expansions

Anonymous, Atsushi Ueda, Yuchi He +2 more

Representing the time-evolution operator as a tensor network constitutes a key ingredient in several algorithms for studying quantum lattice systems at finite temperature or in a nonequilibrium setting. For a Hamiltonian composed of strictly short-range interactions, the Suzuki-Trotter decomposition is the main technique for strategy-the cluster expansion. This approach naturally preserves internal and lattice symmetries and can more easily be extended to higher-order representations or longer-range interactions. We extend the cluster expansion to two-dimensional fermionic systems and employ it to construct projected entangled-pair operator (PEPO) approximations of Gibbs states. We also discuss and benchmark different truncation schemes for multiplying layers of PEPOs together. Applying the resulting framework to a two-dimensional spinless fermion model with attractive interactions, we resolve a clear phase boundary at finite temperature.

Physical review. B./Physical review. B•Jul 23, 2026•0 citations

Lowering the temperature of two-dimensional fermionic tensor networks with cluster expansions

Sander De Meyer, Atsushi Ueda, Yuchi He +2 more

Representing the time-evolution operator as a tensor network constitutes a key ingredient in several algorithms for studying quantum lattice systems at finite temperature or in a non-equilibrium setting. For a Hamiltonian composed of strictly short-ranged interactions, the Suzuki-Trotter decomposition is the main technique for obtaining such a representation. In [B.~Vanhecke, L.~Vanderstraeten and F.~Verstraete, Physical Review A, L020402 (2021)], an alternative strategy, the cluster expansion, was introduced. This approach naturally preserves internal and lattice symmetries and can more easily be extended to higher-order representations or longer-ranged interactions. We extend the cluster expansion to two-dimensional fermionic systems, and employ it to construct projected entangled-pair operator (PEPO) approximations of Gibbs states. We also discuss and benchmark different truncation schemes for multiplying layers of PEPOs together. Applying the resulting framework to a two-dimensional spinless fermion model with attractive interactions, we resolve a clear phase boundary at finite temperature.

cond-mat.str-elCross-listedarXiv:2602.22113•Feb 26, 2026•1 citation

Hamiltonian truncation tensor networks for quantum field theories

Philipp Schmoll, Jan Naumann, Alexander Nietner +2 more

Understanding the equilibrium properties and out of equilibrium dynamics of quantum field theories are key aspects of fundamental problems in theoretical particle physics and cosmology. However, their classical simulation is highly challenging. In this work, we introduce a tensor network method for the classical simulation of continuous quantum field theories that is suitable for the study of low-energy eigenstates and out-of-equilibrium time evolution. The method is built on Hamiltonian truncation and tensor network techniques, bridging the gap between two successful approaches. One of the key developments is the exact construction of matrix product state representations of global projectors, crucial for the implementation of interacting theories. Despite featuring a relatively high computational effort, our method dramatically improves predictive precision compared to exact diagonalisation-based Hamiltonian truncation, allowing the study of so far unexplored parameter regimes and dynamical effects. We corroborate trust in the accuracy of the method by comparing it with exact theoretical results for ground state properties of the sine-Gordon model. We then proceed with discussing $(1+1)$-dimensional quantum electrodynamics, the massive Schwinger model, for which we accurately locate its critical point and study the growth and saturation of momentum-space entanglement in sudden quenches.

arXiv:2312.12506•Dec 19, 2023•1 citation
The study of real-time dynamics of fermions remains one of the last frontiers beyond the reach of classical simulations and is key to our understanding of quantum behavior in chemistry and materials, with implications for quantum technology. Here we introduce a Heisenberg-picture algorithm that propagates observables expressed in a Majorana-string basis using a truncation scheme that preserves Trotter accuracy and aims at maintaining computational efficiency. The framework is exact for quadratic Hamiltonians--remaining restricted to a fixed low-weight sector determined by the physical observable--admits variational initial states, and can be extended to interacting regimes via systematically controlled truncations. We benchmark our approach on one- and two-dimensional Fermi-Hubbard quenches, comparing against tensor network methods (MPS and fPEPS) and recent experimental data. The method achieves high accuracy on timescales comparable to state-of-the-art variational techniques and experiments, demonstrating that controlled Majorana-string truncation is a practical tool for simulating two-dimensional fermionic dynamics.
Simulating the time dynamics of an observable under Hamiltonian evolution is one of the most promising candidates for quantum advantage as we do not expect efficient classical algorithms for this problem except in restricted settings. Here, we introduce such a setting by showing that Majorana Propagation, a simple algorithm combining Trotter steps and truncations, efficiently finds a low-degree approximation of the time-evolved observable as soon as such an approximation exists. This provides the first provable guarantee about Majorana Propagation for Hamiltonian evolution. As an application of this result, we prove that Majorana Propagation can efficiently simulate the time dynamics of any sparse quartic Hamiltonian up to time tmax(u)t_{\text{max}}(u)tmax​(u) depending on the interaction strength uuu. For a time horizon t≤tmax(u)t \leq t_{\text{max}}(u)t≤tmax​(u), the runtime of the algorithm is NO(log⁡(t/ε))N^{O(\log(t/\varepsilon))}NO(log(t/ε)) where NNN is the number of Majorana modes and ε\varepsilonε is the error measured in the normalized Frobenius norm. Importantly, in the limit of small uuu, tmax(u)t_{\text{max}}(u)tmax​(u) goes to +∞+\infty+∞, formalizing the intuition that the algorithm is accurate at all times when the Hamiltonian is quadratic.
We consider the task of simulating the real-time dynamics of weakly interacting fermionic systems. In particular, we focus on computing the expectation value of a local observable AAA at time ttt. By analyzing the convergence of the perturbative expansion in the interaction strength λλλ for the Heisenberg-picture observable, we propose a polynomial-time algorithm for estimating this expectation value in the weakly interacting regime λ∣t∣2D+1=O(1)λ|t|^{2D+1}=\mathcal{O}(1)λ∣t∣2D+1=O(1), when the Hamiltonian is geometrically local on a DDD-dimensional lattice. Importantly, this condition is independent of the system size. If the goal is instead to approximate the time-evolved observable in normalized Frobenius norm, we extend the convergence regime to λ∣t∣=O(1)λ|t|=\mathcal{O}(1)λ∣t∣=O(1) with quasi-polynomial runtime. When the non-interacting part exhibits Anderson localization, our polynomial-time algorithm can be extended up to λ∣t∣=O(1)λ|t|=\mathcal{O}(1)λ∣t∣=O(1), modulo polylogarithmic factors. Our algorithm brings together ideas from continuous-time QMC, diagrammatic QMC, and Majorana Propagation, but with a new Heisenberg-picture operator-growth analysis that makes the sampling complexity rigorously controllable. This leads to provably efficient classical algorithms in regimes where the interaction is weak enough that the sampling variance remains bounded independently of system size. Together, these results identify broad regimes in which weak interactions, locality, and localization can be leveraged to make real-time fermionic dynamics classically tractable.
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.
Representing the time-evolution operator as a tensor network constitutes a key ingredient in several algorithms for studying quantum lattice systems at finite temperature or in a nonequilibrium setting. For a Hamiltonian composed of strictly short-range interactions, the Suzuki-Trotter decomposition is the main technique for strategy-the cluster expansion. This approach naturally preserves internal and lattice symmetries and can more easily be extended to higher-order representations or longer-range interactions. We extend the cluster expansion to two-dimensional fermionic systems and employ it to construct projected entangled-pair operator (PEPO) approximations of Gibbs states. We also discuss and benchmark different truncation schemes for multiplying layers of PEPOs together. Applying the resulting framework to a two-dimensional spinless fermion model with attractive interactions, we resolve a clear phase boundary at finite temperature.
Representing the time-evolution operator as a tensor network constitutes a key ingredient in several algorithms for studying quantum lattice systems at finite temperature or in a non-equilibrium setting. For a Hamiltonian composed of strictly short-ranged interactions, the Suzuki-Trotter decomposition is the main technique for obtaining such a representation. In [B.~Vanhecke, L.~Vanderstraeten and F.~Verstraete, Physical Review A, L020402 (2021)], an alternative strategy, the cluster expansion, was introduced. This approach naturally preserves internal and lattice symmetries and can more easily be extended to higher-order representations or longer-ranged interactions. We extend the cluster expansion to two-dimensional fermionic systems, and employ it to construct projected entangled-pair operator (PEPO) approximations of Gibbs states. We also discuss and benchmark different truncation schemes for multiplying layers of PEPOs together. Applying the resulting framework to a two-dimensional spinless fermion model with attractive interactions, we resolve a clear phase boundary at finite temperature.
Understanding the equilibrium properties and out of equilibrium dynamics of quantum field theories are key aspects of fundamental problems in theoretical particle physics and cosmology. However, their classical simulation is highly challenging. In this work, we introduce a tensor network method for the classical simulation of continuous quantum field theories that is suitable for the study of low-energy eigenstates and out-of-equilibrium time evolution. The method is built on Hamiltonian truncation and tensor network techniques, bridging the gap between two successful approaches. One of the key developments is the exact construction of matrix product state representations of global projectors, crucial for the implementation of interacting theories. Despite featuring a relatively high computational effort, our method dramatically improves predictive precision compared to exact diagonalisation-based Hamiltonian truncation, allowing the study of so far unexplored parameter regimes and dynamical effects. We corroborate trust in the accuracy of the method by comparing it with exact theoretical results for ground state properties of the sine-Gordon model. We then proceed with discussing (1+1)(1+1)(1+1)-dimensional quantum electrodynamics, the massive Schwinger model, for which we accurately locate its critical point and study the growth and saturation of momentum-space entanglement in sudden quenches.