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

Symmetrizing the constraints: Density matrix renormalization group for constrained lattice models

Ting-Tung Wang, Xiaoxue Ran, Zi Yang Meng

Published Jul 23, 2025• 1 citation
arXivHTMLPDF

Abstract

We develop a density matrix renormalization group (DMRG) algorithm for constrained quantum lattice models that successfully {\it{implements the local constraints as symmetries in the contraction of the matrix product states and matrix product operators}}. Such an implementation allows us to investigate a quantum dimer model in DMRG for any lattice geometry wrapped around a cylinder with substantial circumference. We have thence computed the ground state phase diagram of the quantum dimer model on triangular lattice, with the symmetry-breaking characteristics of the columnar solid phase and $\sqrt{12}\times\sqrt{12}$ valence bond solid phase fully captured, as well as the topological entanglement entropy of the $\mathbb{Z}_2$ quantum spin liquid phase that extends to the RK point on non-bipartite lattice accurately revealed. Our DMRG algorithm on constrained quantum lattice models opens new opportunities for matrix and tensor-based algorithms for these systems that have immediate relevance towards the frustrated quantum magnets and synthetic quantum simulators.

DOI
10.1103/z4z8-zps4

Related papers

Real-space parallel density matrix renormalization group

Edwin Miles Stoudenmire, Steven R. White

We demonstrate how to parallelize the density matrix renormalization group (DMRG) algorithm in real space through a straightforward modification of serial DMRG. This makes it possible to apply at least an order of magnitude more computational power to challenging simulations, greatly accelerating investigations of two-dimensional systems and large parameter spaces. We discuss details of the algorithm and present benchmark results including a study of valence-bond-solid order within the square-lattice ${Q}_{2}$ model and N\'eel order within the triangular lattice Heisenberg model. The parallel DMRG algorithm also motivates an alternative canonical form for matrix product states.

arXiv:1301.3494•Physical Review B•Apr 22, 2013•90 citations

Entanglement and the density matrix renormalization group in the generalized Landau paradigm

Laurens Lootens, Clement Delcamp, Frank Verstraete

The fields of entanglement theory and tensor networks have recently emerged as central tools for characterizing quantum phases of matter. Here we determine the entanglement structure of ground states of gapped symmetric quantum lattice models and use this to obtain the most efficient tensor network representation of those ground states. We do this by showing that degeneracies in the entanglement spectrum arise through a duality transformation of the original model to the unique dual model where the entire dual symmetry is spontaneously broken. Physically, this duality transformation amounts to a-potentially twisted-gauging of the unbroken symmetry in the original ground state. In general, the dual symmetries of the resulting models are generalized non-invertible symmetries that cannot be described by groups. This result has strong implications for the complexity of simulating many-body systems using variational tensor network methods. For every phase in the phase diagram, the dual representation of the ground state that completely breaks the symmetry minimizes both the entanglement entropy and the required number of variational parameters. We demonstrate the applicability of this idea by developing a generalized density matrix renormalization group algorithm that works on constrained Hilbert spaces and quantify the computational gains obtained over traditional tensor network methods in a perturbed Heisenberg model. Our work testifies to the usefulness of generalized non-invertible symmetries and their formal category theoretic description for the practical simulation of strongly correlated systems.

arXiv:2408.06334•Nature Physics•Aug 15, 2025•12 citations

Kitaev–Ising-J1-J2 Model: A Density Matrix Renormalization Group Study

A. V. Kapranov, Ramil S. Akzyanov

We numerically study the Kitaev honeycomb model with the additional XX Ising interaction between the nearest and the next nearest neighbors (Kitaev–Ising- $${{J}_{1}}$$ - $${{J}_{2}}$$ model), by using the density matrix renormalization group method. Such additional interaction correspond to the nearest and diagonal interactions on the square lattice. Phase diagram of the bare Kitaev model consist of low entangled commensurate magnetic phases and entangled Kitaev spin liquid. Anisotropic Ising interaction allows the entangled quantum paramagnetic phases in the phase diagram, which in the absence of the magnetic field previously was predicted for more complex type of interaction. We study the scaling law of the entanglement entropy and the bond dimension of the matrix product state with the size of the system. In addition, we propose an optimization algorithm to prevent density matrix renormalization group (DMRG) from getting stuck in the low-entangled phases.

arXiv:2502.08389•Journal of Experimental and Theoretical Physics Letters•Apr 1, 2026•0 citations

Frustrated spin systems, an MPS approach

S. N. Saadatmand

In spin lattices, frustration happens when there exists a competition between different terms of the Hamiltonian, which cannot be simultaneously minimized in all local physical bases. Frustration can destabilize the bipartite Neel order, which is the conventional groundstate of antiferromagnetic Heisenberg-type models in two dimensions. Hence, the focus of two-dimensional quantum magnetism studies has been directed toward elucidating the zero-temperature properties of Heisenberg-type models on the Archimedean lattices that exhibit geometrical frustration. This type of frustration is of great interest since it can lead to exotic forms of the quantum matter, such as topological spin-liquid (SL) and many-sublattice long-range-ordered phases. The spin-1/2 triangular Heisenberg model (THM) is a prototypical model with geometrical frustration, which has the highest coordination number of all Archimedean lattices. However, there exist two major, distinct difficulties for numerical methods to efficiently determine the groundstate properties of the THM: the immense dimension of the Hilbert space and the geometrical frustration phenomenon. In this thesis, to deal with such difficulties, we employ the state-of-the- art non-Abelian matrix product states (MPS) and density matrix renormalization group (DMRG) methods. We exploit the symmetries of the Hamiltonian to reduce the working dimension of the Hilbert space and, to lessen the effects of the geometrical frustration, we depend on the local entanglement (area-law) nature of the MPS as it scales with the size of a lattice subspace with reduced spatial dimensions. In such a way, we thoroughly examine the phase diagram of the THM with nearest and next-nearest neighbor (NN and NNN) exchange couplings on finite- and infinite-length cylinders of widths up to 12 sites. In addition, we develop new numerical tools to analyze topological and symmetry-broken phases in the context of the SU(2)-symmetric translation-invariant MPS algorithm. We establish that on infinite cylinders, the anyonic sectors of a topological order can be classified through the fractionalization of global symmetries and degeneracy patterns of the entanglement spectrum (ES). On the other hand, we detect symmetry-breaking phase transitions by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters, even though symmetry implies that the order parameter itself is strictly zero. Furthermore, the appearance of Nambu-Goldstone modes in the excitation spectrum can be detected using “tower of states” (TOS) levels in the momentum-resolved ES. By applying the new tools to the THM, we discover a variety of exotic groundstates such as a Z2-gauge toric-code-type topological order (with vison and spinon excitations), a NNN Majumdar-Gosh, and a long-range 120◦-ordered state.

•Jun 28, 2017•2 citations

Detection and characterization of symmetry-broken long-range orders in the spin- 12 triangular Heisenberg model

S. N. Saadatmand, Ian P. McCulloch

We present new numerical tools to analyze symmetry-broken phases in the context of $\mathrm{SU}(2)$-symmetric translation-invariant matrix product states (MPS) and density-matrix renormalization-group (DMRG) methods for infinite cylinders, and determine the phase diagram of the geometrically frustrated triangular Heisenberg model with nearest- and next-nearest-neighbor (NN and NNN) interactions. The appearance of Nambu-Goldstone modes in the excitation spectrum is characterized by ``tower of states'' levels in the momentum-resolved entanglement spectrum. Symmetry-breaking phase transitions are detected by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters (which we call the Binder ratio), even though symmetry implies that the order parameter itself is strictly zero. Using this approach, we have identified a ${120}^{\ensuremath{\circ}}$ order, a columnar order, and an algebraic spin liquid (specific to width-6 systems), alongside the previously studied topological spin liquid phase. For the latter, we also demonstrate robustness against chiral perturbations.

arXiv:1704.03418•Physical review. B./Physical review. B•Aug 9, 2017•41 citations

Chiral Spin Liquid Phase of the Triangular Lattice Hubbard Model: A Density Matrix Renormalization Group Study

Aaron Szasz, Johannes Motruk, Michael P. Zaletel +1 more

Motivated by experimental studies that have found signatures of a quantum spin liquid phase in organic crystals whose structure is well described by the two-dimensional triangular lattice, we study the Hubbard model on this lattice at half filling using the infinite-system density matrix renormalization group (iDMRG) method. On infinite cylinders with finite circumference, we identify an intermediate phase between observed metallic behavior at low interaction strength and Mott insulating spin-ordered behavior at strong interactions. Chiral ordering from spontaneous breaking of time-reversal symmetry, a fractionally quantized spin Hall response, and characteristic level statistics in the entanglement spectrum in the intermediate phase provide strong evidence for the existence of a chiral spin liquid in the full two-dimensional limit of the model.

•Apr 1, 2020•50 citations
We develop a density matrix renormalization group (DMRG) algorithm for constrained quantum lattice models that successfully {\it{implements the local constraints as symmetries in the contraction of the matrix product states and matrix product operators}}. Such an implementation allows us to investigate a quantum dimer model in DMRG for any lattice geometry wrapped around a cylinder with substantial circumference. We have thence computed the ground state phase diagram of the quantum dimer model on triangular lattice, with the symmetry-breaking characteristics of the columnar solid phase and 12×12\sqrt{12}\times\sqrt{12}12​×12​ valence bond solid phase fully captured, as well as the topological entanglement entropy of the Z2\mathbb{Z}_2Z2​ quantum spin liquid phase that extends to the RK point on non-bipartite lattice accurately revealed. Our DMRG algorithm on constrained quantum lattice models opens new opportunities for matrix and tensor-based algorithms for these systems that have immediate relevance towards the frustrated quantum magnets and synthetic quantum simulators.
We demonstrate how to parallelize the density matrix renormalization group (DMRG) algorithm in real space through a straightforward modification of serial DMRG. This makes it possible to apply at least an order of magnitude more computational power to challenging simulations, greatly accelerating investigations of two-dimensional systems and large parameter spaces. We discuss details of the algorithm and present benchmark results including a study of valence-bond-solid order within the square-lattice Q2{Q}_{2}Q2​ model and N\'eel order within the triangular lattice Heisenberg model. The parallel DMRG algorithm also motivates an alternative canonical form for matrix product states.
The fields of entanglement theory and tensor networks have recently emerged as central tools for characterizing quantum phases of matter. Here we determine the entanglement structure of ground states of gapped symmetric quantum lattice models and use this to obtain the most efficient tensor network representation of those ground states. We do this by showing that degeneracies in the entanglement spectrum arise through a duality transformation of the original model to the unique dual model where the entire dual symmetry is spontaneously broken. Physically, this duality transformation amounts to a-potentially twisted-gauging of the unbroken symmetry in the original ground state. In general, the dual symmetries of the resulting models are generalized non-invertible symmetries that cannot be described by groups. This result has strong implications for the complexity of simulating many-body systems using variational tensor network methods. For every phase in the phase diagram, the dual representation of the ground state that completely breaks the symmetry minimizes both the entanglement entropy and the required number of variational parameters. We demonstrate the applicability of this idea by developing a generalized density matrix renormalization group algorithm that works on constrained Hilbert spaces and quantify the computational gains obtained over traditional tensor network methods in a perturbed Heisenberg model. Our work testifies to the usefulness of generalized non-invertible symmetries and their formal category theoretic description for the practical simulation of strongly correlated systems.
We numerically study the Kitaev honeycomb model with the additional XX Ising interaction between the nearest and the next nearest neighbors (Kitaev–Ising- J1{{J}_{1}}J1​ - J2{{J}_{2}}J2​ model), by using the density matrix renormalization group method. Such additional interaction correspond to the nearest and diagonal interactions on the square lattice. Phase diagram of the bare Kitaev model consist of low entangled commensurate magnetic phases and entangled Kitaev spin liquid. Anisotropic Ising interaction allows the entangled quantum paramagnetic phases in the phase diagram, which in the absence of the magnetic field previously was predicted for more complex type of interaction. We study the scaling law of the entanglement entropy and the bond dimension of the matrix product state with the size of the system. In addition, we propose an optimization algorithm to prevent density matrix renormalization group (DMRG) from getting stuck in the low-entangled phases.
In spin lattices, frustration happens when there exists a competition between different terms of the Hamiltonian, which cannot be simultaneously minimized in all local physical bases. Frustration can destabilize the bipartite Neel order, which is the conventional groundstate of antiferromagnetic Heisenberg-type models in two dimensions. Hence, the focus of two-dimensional quantum magnetism studies has been directed toward elucidating the zero-temperature properties of Heisenberg-type models on the Archimedean lattices that exhibit geometrical frustration. This type of frustration is of great interest since it can lead to exotic forms of the quantum matter, such as topological spin-liquid (SL) and many-sublattice long-range-ordered phases. The spin-1/2 triangular Heisenberg model (THM) is a prototypical model with geometrical frustration, which has the highest coordination number of all Archimedean lattices. However, there exist two major, distinct difficulties for numerical methods to efficiently determine the groundstate properties of the THM: the immense dimension of the Hilbert space and the geometrical frustration phenomenon. In this thesis, to deal with such difficulties, we employ the state-of-the- art non-Abelian matrix product states (MPS) and density matrix renormalization group (DMRG) methods. We exploit the symmetries of the Hamiltonian to reduce the working dimension of the Hilbert space and, to lessen the effects of the geometrical frustration, we depend on the local entanglement (area-law) nature of the MPS as it scales with the size of a lattice subspace with reduced spatial dimensions. In such a way, we thoroughly examine the phase diagram of the THM with nearest and next-nearest neighbor (NN and NNN) exchange couplings on finite- and infinite-length cylinders of widths up to 12 sites. In addition, we develop new numerical tools to analyze topological and symmetry-broken phases in the context of the SU(2)-symmetric translation-invariant MPS algorithm. We establish that on infinite cylinders, the anyonic sectors of a topological order can be classified through the fractionalization of global symmetries and degeneracy patterns of the entanglement spectrum (ES). On the other hand, we detect symmetry-breaking phase transitions by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters, even though symmetry implies that the order parameter itself is strictly zero. Furthermore, the appearance of Nambu-Goldstone modes in the excitation spectrum can be detected using “tower of states” (TOS) levels in the momentum-resolved ES. By applying the new tools to the THM, we discover a variety of exotic groundstates such as a Z2-gauge toric-code-type topological order (with vison and spinon excitations), a NNN Majumdar-Gosh, and a long-range 120◦-ordered state.
We present new numerical tools to analyze symmetry-broken phases in the context of SU(2)\mathrm{SU}(2)SU(2)-symmetric translation-invariant matrix product states (MPS) and density-matrix renormalization-group (DMRG) methods for infinite cylinders, and determine the phase diagram of the geometrically frustrated triangular Heisenberg model with nearest- and next-nearest-neighbor (NN and NNN) interactions. The appearance of Nambu-Goldstone modes in the excitation spectrum is characterized by ``tower of states'' levels in the momentum-resolved entanglement spectrum. Symmetry-breaking phase transitions are detected by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters (which we call the Binder ratio), even though symmetry implies that the order parameter itself is strictly zero. Using this approach, we have identified a ${120}^{\ensuremath{\circ}}$ order, a columnar order, and an algebraic spin liquid (specific to width-6 systems), alongside the previously studied topological spin liquid phase. For the latter, we also demonstrate robustness against chiral perturbations.
Motivated by experimental studies that have found signatures of a quantum spin liquid phase in organic crystals whose structure is well described by the two-dimensional triangular lattice, we study the Hubbard model on this lattice at half filling using the infinite-system density matrix renormalization group (iDMRG) method. On infinite cylinders with finite circumference, we identify an intermediate phase between observed metallic behavior at low interaction strength and Mott insulating spin-ordered behavior at strong interactions. Chiral ordering from spontaneous breaking of time-reversal symmetry, a fractionally quantized spin Hall response, and characteristic level statistics in the entanglement spectrum in the intermediate phase provide strong evidence for the existence of a chiral spin liquid in the full two-dimensional limit of the model.