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:1504.01566 · Physical Review E

Topology-driven phase transitions in the classical monomer-dimer-loop model

Sazi Li, Wei Li, Ziyu Chen

Published Jun 11, 2015• 1 citation
arXivHTMLPDF

Abstract

In this work, we investigate the classical loop models doped with monomers and dimers on a square lattice, whose partition function can be expressed as a tensor network (TN). In the thermodynamic limit, we use the boundary matrix product state technique to contract the partition function TN, and determine the thermodynamic properties with high accuracy. In this monomer-dimer-loop model, we find a second-order phase transition between a trivial monomer-condensation and a loop-condensation (LC) phase, which cannot be distinguished by any local order parameter, while nevertheless the two phases have distinct topological properties. In the LC phase, we find two degenerate dominating eigenvalues in the transfer-matrix spectrum, as well as a nonvanishing (nonlocal) string order parameter, both of which identify the topological ergodicity breaking in the LC phase and can serve as the order parameter for detecting the phase transitions.

DOI
10.1103/physreve.91.060104

Related papers

Kosterlitz-Thouless transitions and phase diagrams of the interacting monomer-dimer model on a checkerboard lattice

Sazi Li, Wei Li, Ziyu Chen

Using the tensor network approach, we investigate the monomer-dimer models on a checkerboard lattice, in which there are interactions (with strength $v$) between the parallel dimers on half of the plaquettes. For the fully packed interacting dimer model, we observe a Kosterlitz-Thouless (KT) transition between the low-temperature symmetry breaking and the high-temperature critical phases; for the doped monomer-dimer case with finite chemical potential $\ensuremath{\mu}$, we also find an order-disorder phase transition which is of second order instead. We use the boundary matrix product state approach to detect the KT and second-order phase transitions and obtain the phase diagrams $v\text{\ensuremath{-}}T$ and $\ensuremath{\mu}\text{\ensuremath{-}}T$. Moreover, for the noninteracting monomer-dimer model (setting $\ensuremath{\mu}=\ensuremath{\nu}=0$), we get an extraordinarily accurate determination of the free energy per site (negative of the monomer-dimer constant ${h}_{2}$) as $f=\ensuremath{-}0.662\phantom{\rule{0.16em}{0ex}}798\phantom{\rule{0.16em}{0ex}}972\phantom{\rule{0.16em}{0ex}}833\phantom{\rule{0.16em}{0ex}}746$ with the dimer density $n=0.638\phantom{\rule{0.16em}{0ex}}123\phantom{\rule{0.16em}{0ex}}109\phantom{\rule{0.16em}{0ex}}228\phantom{\rule{0.16em}{0ex}}547$, both of 15 correct digits.

arXiv:1504.01565•Physical Review E•Nov 5, 2014•4 citations

Topological sectors, dimer correlations, and monomers from the transfer-matrix solution of the dimer model

Neil Wilkins, Stephen Powell

We solve the classical square-lattice dimer model with periodic boundaries and in the presence of a field t that couples to the (vector) flux, by diagonalizing a modified version of Lieb's transfer matrix. After deriving the torus partition function in the thermodynamic limit, we show how the configuration space divides into topological sectors corresponding to distinct values of the flux. Additionally, we demonstrate in general that expectation values are t independent at leading order, and obtain explicit expressions for dimer occupation numbers, dimer-dimer correlation functions, and the monomer distribution function. The last of these is expressed as a Toeplitz determinant, whose asymptotic behavior for large monomer separation is tractable using the Fisher-Hartwig conjecture. Our results reproduce those previously obtained using Pfaffian techniques.

arXiv:2103.03260•Physical review. E•Jul 29, 2021•7 citations

Tensor network state approach to quantum topological phase transitions and their criticalities of Z2 topologically ordered states

Wen-Tao Xu, Guang-Ming Zhang

Due to the absence of local order parameters, it is a challenging task to characterize the quantum topological phase transitions between topologically ordered phases in two dimensions. In this paper, we construct a topologically ordered tensor network wave function with one parameter $\ensuremath{\lambda}$, describing both the toric code state $(\ensuremath{\lambda}=1)$ and double semion state $(\ensuremath{\lambda}=\ensuremath{-}1)$. Via calculating the correlation length defined from the one-dimensional quantum transfer operator of the wave function norm, we can map out the complete phase diagram in terms of the parameter $\ensuremath{\lambda}$, and three different quantum critical points (QCPs) at $\ensuremath{\lambda}=0,\ifmmode\pm\else\textpm\fi{}1.73$ are identified. The first one separates the toric code phase and double semion phase, while latter two describe the topological phase transitions from the toric code phase or double semion phase to the symmetry-breaking phase, respectively. When mapping the quantum tensor network wave function to the exactly solved statistical model, the norm of the wave function is identified as the partition function of the classical eight-vertex model, and both QCPs at $\ensuremath{\lambda}=\ifmmode\pm\else\textpm\fi{}1.73$ correspond to the eight-vertex model at the critical point $\ensuremath{\lambda}=\sqrt{3}$, while the QCP at $\ensuremath{\lambda}=0$ corresponds to the critical six-vertex model. Actually such a quantum-classical mapping cannot yield the complete low-energy excitations at these three QCPs. We further demonstrate that the full eigenvalue spectra of the transfer operators without/with the flux insertions can give rise to the complete quantum criticalities, which are described by the (2+0)-dimensional free boson conformal field theories (CFTs) compactified on a circle with the radius $R=\sqrt{6}$ at $\ensuremath{\lambda}=\ifmmode\pm\else\textpm\fi{}\sqrt{3}$ and $R=\sqrt{8/3}$ at $\ensuremath{\lambda}=0$. From the complete transfer operator spectra, the finite-size spectra of the CFTs for the critical eight-vertex model are obtained, and the topological sectors of anyonic excitations are yielded as well. Furthermore, for the QCP at $\ensuremath{\lambda}=0$, no anyon condensation occurs, but the emerged symmetries of the matrix product operators significantly enrich the topological sectors of the CFT spectra. Finally, we provide our understanding of the (2+0)-dimensional conformal quantum criticalities and their possible connection with the generic (2+1)-dimensional CFTs for quantum topological phase transitions.

arXiv:1807.08490•Physical review. B./Physical review. B•Oct 10, 2018•21 citations

Ashkin-Teller phase transition and multicritical behavior in a classical monomer-dimer model

Satoshi Morita, Hyun‐Yong Lee, Kedar Damle +1 more

We use Monte Carlo simulations and tensor network methods to study a classical monomer-dimer model on the square lattice with a hole (monomer) fugacity $z$, an aligning dimer-dimer interaction $u$ that favors columnar order, and an attractive dimer-dimer interaction $v$ between two adjacent dimers that lie on the same principal axis of the lattice. The Monte Carlo simulations of finite size systems rely on our grand-canonical generalization of the dimer worm algorithm, while the tensor network computations are based on a uniform matrix product ansatz for the eigenvector of the row-to-row transfer matrix that work directly in the thermodynamic limit. The phase diagram has nematic, columnar order and fluid phases, and a nonzero temperature multicritical point at which all three meet. For any fixed $v/u<\ensuremath{\infty}$, we argue that this multicritical point continues to be located at a nonzero hole fugacity ${z}_{\mathrm{mc}}(v/u)>0$; our numerical results confirm this theoretical expectation but find that ${z}_{\mathrm{mc}}(v/u)\ensuremath{\rightarrow}0$ very rapidly as $v/u\ensuremath{\rightarrow}\ensuremath{\infty}$. Our numerical results also confirm the theoretical expectation that the corresponding multicritical behavior is in the universality class of the four-state Potts multicritical point on critical line of the two-dimensional Ashkin-Teller model.

arXiv:2306.02578•Physical Review Research•Oct 19, 2023•4 citations

Transition of a Z3 topologically ordered phase to trivial and critical phases

Ching-Yu Huang, Tzu-Chieh Wei

Topologically ordered quantum systems have robust physical properties, such as quasiparticle statistics, and ground-state degeneracy, which do not depend on the microscopic details of the Hamiltonian. We consider topological phase transitions under a deformation such as an effective string tension on a ${\mathbb{Z}}_{3}$ topological state. This is studied in terms of the gauge-symmetry-preserved quantum state renormalization group, first proposed by He, Moradi, and Wen [Phys. Rev. B 90, 205114 (2014)]. In this approach modular matrices $S$ and $T$ can be obtained and used as order parameters to characterize the topological properties of the phase and determine phase transitions. From a mapping to a classical two-dimensional Potts model on the square lattice, the critical string tension, at which the transition to a topologically trivial phase takes place, can be obtained analytically and agrees with the numerically determined value. Such a transition can be generalized to a ${\mathbb{Z}}_{N}$ topological model under a string tension and determined in the same way. With different deformations, the ${\mathbb{Z}}_{3}$ topological phase can also be driven to a critical phase, which contains, in the large deformation limit, the wave function analogous to the Rokhsar-Kivelson point in the quantum dimer model in one case and the fully packed loop model in another case.

arXiv:1503.03419•Physical Review B•Aug 6, 2015•8 citations

Quantum phase transition between symmetry enriched topological phases in tensor-network states

Lukas Haller, Wen-Tao Xu, Yu-Jie Liu +1 more

Quantum phase transitions between different topologically ordered phases exhibit rich structures and are generically challenging to study in microscopic lattice models. In this paper, we propose a tensor-network solvable model that allows us to tune between different symmetry enriched topological (SET) phases. Concretely, we consider a decorated two-dimensional toric code model for which the ground state can be expressed as a two-dimensional tensor-network state with bond dimension $D=3$ and two tunable parameters. We find that the time-reversal (TR) symmetric system exhibits three distinct phases: (i) an SET toric code phase in which anyons transform nontrivially under TR, (ii) a toric code phase in which TR does not fractionalize, and (iii) a topologically trivial phase that is adiabatically connected to a product state. We characterize the different phases using the topological entanglement entropy and a membrane order parameter that distinguishes the two SET phases. Along the phase boundary between the SET toric code phase and the toric code phase, the model has an enhanced $U(1)$ symmetry and the ground state is a quantum critical loop gas wavefunction whose squared norm is equivalent to the partition function of the classical $O(2)$ model. By duality transformations, this tensor-network solvable model can also be used to describe transitions between SET double-semion phases and between ${\mathbb{Z}}_{2}\ifmmode\times\else\texttimes\fi{}{\mathbb{Z}}_{2}^{T}$ symmetry protected topological phases in two dimensions.

arXiv:2305.02432•Physical Review Research•Oct 23, 2023•18 citations
In this work, we investigate the classical loop models doped with monomers and dimers on a square lattice, whose partition function can be expressed as a tensor network (TN). In the thermodynamic limit, we use the boundary matrix product state technique to contract the partition function TN, and determine the thermodynamic properties with high accuracy. In this monomer-dimer-loop model, we find a second-order phase transition between a trivial monomer-condensation and a loop-condensation (LC) phase, which cannot be distinguished by any local order parameter, while nevertheless the two phases have distinct topological properties. In the LC phase, we find two degenerate dominating eigenvalues in the transfer-matrix spectrum, as well as a nonvanishing (nonlocal) string order parameter, both of which identify the topological ergodicity breaking in the LC phase and can serve as the order parameter for detecting the phase transitions.
Using the tensor network approach, we investigate the monomer-dimer models on a checkerboard lattice, in which there are interactions (with strength vvv) between the parallel dimers on half of the plaquettes. For the fully packed interacting dimer model, we observe a Kosterlitz-Thouless (KT) transition between the low-temperature symmetry breaking and the high-temperature critical phases; for the doped monomer-dimer case with finite chemical potential $\ensuremath{\mu}$, we also find an order-disorder phase transition which is of second order instead. We use the boundary matrix product state approach to detect the KT and second-order phase transitions and obtain the phase diagrams $v\text{\ensuremath{-}}T$ and $\ensuremath{\mu}\text{\ensuremath{-}}T$. Moreover, for the noninteracting monomer-dimer model (setting $\ensuremath{\mu}=\ensuremath{\nu}=0$), we get an extraordinarily accurate determination of the free energy per site (negative of the monomer-dimer constant h2{h}_{2}h2​) as $f=\ensuremath{-}0.662\phantom{\rule{0.16em}{0ex}}798\phantom{\rule{0.16em}{0ex}}972\phantom{\rule{0.16em}{0ex}}833\phantom{\rule{0.16em}{0ex}}746$ with the dimer density n=0.638123109228547n=0.638\phantom{\rule{0.16em}{0ex}}123\phantom{\rule{0.16em}{0ex}}109\phantom{\rule{0.16em}{0ex}}228\phantom{\rule{0.16em}{0ex}}547n=0.638123109228547, both of 15 correct digits.
We solve the classical square-lattice dimer model with periodic boundaries and in the presence of a field t that couples to the (vector) flux, by diagonalizing a modified version of Lieb's transfer matrix. After deriving the torus partition function in the thermodynamic limit, we show how the configuration space divides into topological sectors corresponding to distinct values of the flux. Additionally, we demonstrate in general that expectation values are t independent at leading order, and obtain explicit expressions for dimer occupation numbers, dimer-dimer correlation functions, and the monomer distribution function. The last of these is expressed as a Toeplitz determinant, whose asymptotic behavior for large monomer separation is tractable using the Fisher-Hartwig conjecture. Our results reproduce those previously obtained using Pfaffian techniques.
Due to the absence of local order parameters, it is a challenging task to characterize the quantum topological phase transitions between topologically ordered phases in two dimensions. In this paper, we construct a topologically ordered tensor network wave function with one parameter $\ensuremath{\lambda}$, describing both the toric code state $(\ensuremath{\lambda}=1)$ and double semion state $(\ensuremath{\lambda}=\ensuremath{-}1)$. Via calculating the correlation length defined from the one-dimensional quantum transfer operator of the wave function norm, we can map out the complete phase diagram in terms of the parameter $\ensuremath{\lambda}$, and three different quantum critical points (QCPs) at $\ensuremath{\lambda}=0,\ifmmode\pm\else\textpm\fi{}1.73$ are identified. The first one separates the toric code phase and double semion phase, while latter two describe the topological phase transitions from the toric code phase or double semion phase to the symmetry-breaking phase, respectively. When mapping the quantum tensor network wave function to the exactly solved statistical model, the norm of the wave function is identified as the partition function of the classical eight-vertex model, and both QCPs at $\ensuremath{\lambda}=\ifmmode\pm\else\textpm\fi{}1.73$ correspond to the eight-vertex model at the critical point $\ensuremath{\lambda}=\sqrt{3}$, while the QCP at $\ensuremath{\lambda}=0$ corresponds to the critical six-vertex model. Actually such a quantum-classical mapping cannot yield the complete low-energy excitations at these three QCPs. We further demonstrate that the full eigenvalue spectra of the transfer operators without/with the flux insertions can give rise to the complete quantum criticalities, which are described by the (2+0)-dimensional free boson conformal field theories (CFTs) compactified on a circle with the radius R=6R=\sqrt{6}R=6​ at $\ensuremath{\lambda}=\ifmmode\pm\else\textpm\fi{}\sqrt{3}$ and R=8/3R=\sqrt{8/3}R=8/3​ at $\ensuremath{\lambda}=0$. From the complete transfer operator spectra, the finite-size spectra of the CFTs for the critical eight-vertex model are obtained, and the topological sectors of anyonic excitations are yielded as well. Furthermore, for the QCP at $\ensuremath{\lambda}=0$, no anyon condensation occurs, but the emerged symmetries of the matrix product operators significantly enrich the topological sectors of the CFT spectra. Finally, we provide our understanding of the (2+0)-dimensional conformal quantum criticalities and their possible connection with the generic (2+1)-dimensional CFTs for quantum topological phase transitions.
We use Monte Carlo simulations and tensor network methods to study a classical monomer-dimer model on the square lattice with a hole (monomer) fugacity zzz, an aligning dimer-dimer interaction uuu that favors columnar order, and an attractive dimer-dimer interaction vvv between two adjacent dimers that lie on the same principal axis of the lattice. The Monte Carlo simulations of finite size systems rely on our grand-canonical generalization of the dimer worm algorithm, while the tensor network computations are based on a uniform matrix product ansatz for the eigenvector of the row-to-row transfer matrix that work directly in the thermodynamic limit. The phase diagram has nematic, columnar order and fluid phases, and a nonzero temperature multicritical point at which all three meet. For any fixed $v/u<\ensuremath{\infty}$, we argue that this multicritical point continues to be located at a nonzero hole fugacity zmc(v/u)>0{z}_{\mathrm{mc}}(v/u)>0zmc​(v/u)>0; our numerical results confirm this theoretical expectation but find that ${z}_{\mathrm{mc}}(v/u)\ensuremath{\rightarrow}0$ very rapidly as $v/u\ensuremath{\rightarrow}\ensuremath{\infty}$. Our numerical results also confirm the theoretical expectation that the corresponding multicritical behavior is in the universality class of the four-state Potts multicritical point on critical line of the two-dimensional Ashkin-Teller model.
Topologically ordered quantum systems have robust physical properties, such as quasiparticle statistics, and ground-state degeneracy, which do not depend on the microscopic details of the Hamiltonian. We consider topological phase transitions under a deformation such as an effective string tension on a Z3{\mathbb{Z}}_{3}Z3​ topological state. This is studied in terms of the gauge-symmetry-preserved quantum state renormalization group, first proposed by He, Moradi, and Wen [Phys. Rev. B 90, 205114 (2014)]. In this approach modular matrices SSS and TTT can be obtained and used as order parameters to characterize the topological properties of the phase and determine phase transitions. From a mapping to a classical two-dimensional Potts model on the square lattice, the critical string tension, at which the transition to a topologically trivial phase takes place, can be obtained analytically and agrees with the numerically determined value. Such a transition can be generalized to a ZN{\mathbb{Z}}_{N}ZN​ topological model under a string tension and determined in the same way. With different deformations, the Z3{\mathbb{Z}}_{3}Z3​ topological phase can also be driven to a critical phase, which contains, in the large deformation limit, the wave function analogous to the Rokhsar-Kivelson point in the quantum dimer model in one case and the fully packed loop model in another case.
Quantum phase transitions between different topologically ordered phases exhibit rich structures and are generically challenging to study in microscopic lattice models. In this paper, we propose a tensor-network solvable model that allows us to tune between different symmetry enriched topological (SET) phases. Concretely, we consider a decorated two-dimensional toric code model for which the ground state can be expressed as a two-dimensional tensor-network state with bond dimension D=3D=3D=3 and two tunable parameters. We find that the time-reversal (TR) symmetric system exhibits three distinct phases: (i) an SET toric code phase in which anyons transform nontrivially under TR, (ii) a toric code phase in which TR does not fractionalize, and (iii) a topologically trivial phase that is adiabatically connected to a product state. We characterize the different phases using the topological entanglement entropy and a membrane order parameter that distinguishes the two SET phases. Along the phase boundary between the SET toric code phase and the toric code phase, the model has an enhanced U(1)U(1)U(1) symmetry and the ground state is a quantum critical loop gas wavefunction whose squared norm is equivalent to the partition function of the classical O(2)O(2)O(2) model. By duality transformations, this tensor-network solvable model can also be used to describe transitions between SET double-semion phases and between ${\mathbb{Z}}_{2}\ifmmode\times\else\texttimes\fi{}{\mathbb{Z}}_{2}^{T}$ symmetry protected topological phases in two dimensions.