Search papers, labs, and topics across Lattice.
The authors develop a perfect Born sampling framework for thermal tensor networks that resolves the long-standing trade-off between sampling efficiency and symmetry acceleration by sampling Abelian and non-Abelian quantum numbers directly from the purified density matrix. This enables unbiased, low-temperature exploration of strongly correlated quantum lattice models in both 1D and 2D geometries without suffering from symmetry-induced sampling bottlenecks. Using stochastic matrix product states and projected entangled pair states (stoMPS and stoPEPS), the method accurately tracks real-time dynamics, simulates width-10 square Hubbard cylinders, and reaches ultralow temperatures ($T/t = 1/64$) on triangular lattices to uncover scalar chiral order and kinetic ferromagnetism.
Thermal tensor networks have finally broken the fundamental trade-off between sampling efficiency and symmetry acceleration, scaling quantum lattice simulations to widths of $W=10$ and temperatures down to $T/t = 1/64$.
Accurate calculations of quantum lattice models at low temperatures constitute a major challenge in many-body physics. Stochastic sampling of tensor-network states offers a promising route to tackle this problem; however, existing schemes have long faced a fundamental dilemma---sampling efficiency and symmetry acceleration \textit{cannot} be achieved simultaneously. Here we propose a perfect Born sampling approach for thermal tensor networks, which performs importance sampling directly from the purified density matrix via the Born rule and incorporates Abelian and non-Abelian symmetries by sampling symmetry quantum numbers. We benchmark the method, realized as both stochastic matrix product states (stoMPS) and stochastic projected entangled pair states (stoPEPS), on large-scale quantum lattice models. Using stoMPS, we accurately simulate the square-lattice Hubbard model on cylinders up to width $W=10$, and study the triangular-lattice Hubbard model down to $T/t = 1/64$, revealing scalar chiral order at half filling and kinetic ferromagnetism upon electron doping. We further extend the stoMPS method to compute finite-temperature quantum dynamics, as demonstrated by the optical conductivity of the Hubbard model, and generalize it to stoPEPS, as showcased on the $20\times20$ square-lattice quantum Ising model at its quantum critical point. Our method combines high sampling efficiency with full symmetry acceleration, and can be used as a state-of-the-art framework for studying both equilibrium and dynamical properties down to ultralow temperatures.