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

Quantum phase transitions in the K -layer Ising toric code

Lukas Schamriß, Lea Lenke, Matthias Mühlhauser, Kai Phillip Schmidt

Published May 26, 2022• 4 citations
arXivHTMLPDF

Abstract

We investigate the quantum phase diagram of the $K$-layer Ising toric code corresponding to $K$ layers of two-dimensional toric codes coupled by Ising interactions. While for small Ising interactions the system displays ${\mathbb{Z}}_{2}^{K}$ topological order originating from the toric codes in each layer, the system shows ${\mathbb{Z}}_{2}$ topological order in the high-Ising limit. The latter is demonstrated for general $K$ by deriving an effective low-energy model in $K\mathrm{th}$-order degenerate perturbation theory, which is given as an effective anisotropic single-layer toric code in terms of collective pseudospins 1/2 referring to the two ground states of isolated Ising chain segments. For the specific cases $K=3$ and $K=4$ we apply high-order series expansions to determine the gap series in the low- and high-Ising limit. Extrapolation of the elementary energy gaps gives convincing evidence that the ground-state phase diagram consists of a single quantum critical point and our findings suggest a quantum phase transition in the 3D Ising* universality class for both $K$ separating both types of topological order, which is consistent with former findings for the bilayer Ising toric code.

DOI
10.1103/physrevb.105.184425

Related papers

Quantum critical phase transition between two topologically ordered phases in the Ising toric code bilayer

Raymond Wiedmann, Lea Lenke, Matthias R. Walther +2 more

We demonstrate that two toric code layers on the square lattice coupled by an Ising interaction display two distinct phases with intrinsic topological order. The second-order quantum phase transition between the weakly coupled ${\mathbb{Z}}_{2}\ifmmode\times\else\texttimes\fi{}{\mathbb{Z}}_{2}$ and the strongly coupled ${\mathbb{Z}}_{2}$ topological order can be described by the condensation of bosonic quasiparticles from both sides and belongs to the three-dimensional ${\mathrm{Ising}}^{*}$ universality class. This can be shown by an exact duality transformation to the transverse-field Ising model on the square lattice, which builds on the existence of an extensive number of local ${\mathbb{Z}}_{2}$ conserved parities. These conserved quantities correspond to the product of two adjacent star operators on different layers. Notably, we show that the low-energy effective model derived about the limit of large Ising coupling is given by an effective single-layer toric code in terms of the conserved quantities of the Ising toric code bilayer. The two topological phases are further characterized by the topological entanglement entropy which serves as a nonlocal order parameter.

arXiv:2010.05982•Physical review. B./Physical review. B•Dec 17, 2020•19 citations

High-order series expansion of non-Hermitian quantum spin models

Lea Lenke, Matthias Mühlhauser, Kai Phillip Schmidt

We investigate the low-energy physics of non-Hermitian quantum spin models with $PT$ symmetry. To this end, we consider the one-dimensional Ising chain and the two-dimensional toric code in a non-Hermitian staggered field. For both systems, dual descriptions in terms of non-Hermitian staggered Ising interactions in a conventional transverse field exist. We perform high-order series expansions about the high- and low-field limit for both systems to determine the ground-state energy per site and the one-particle gap. The one-dimensional non-Hermitian Ising chain is known to be exactly solvable. Its ground-state phase diagram consists of second-order quantum phase transitions, which can be characterized by logarithmic singularities of the second derivative of the ground-state energy and, in the symmetry-broken phase, the gap closing of the low-field gap. In contrast, the gap closing from the high-field phase is not accessible perturbatively due to the complex energy and the occurrence of exceptional lines in the high-field gap expression. For the two-dimensional toric code in a non-Hermitian staggered field, we study the quantum robustness of the topologically ordered phase by the gap closing of the low-field gap. We find that the well-known second-order quantum phase transition of the toric code in a uniform field extends into a large portion of the non-Hermitian parameter space. However, the series expansions become unreliable for a dominant anti-Hermitian field. Interestingly, the analysis of the high-field gap reveals the potential presence of an intermediate region.

arXiv:2108.05760•Physical review. B./Physical review. B•Nov 19, 2021•21 citations

Topological and magnetic phase transitions in the bilayer Kitaev-Ising model

Aayush Vijayvargia, Urban F. P. Seifert, Onur Erten

We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped ${\mathbb{Z}}_{2}$ spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated ${120}^{\ensuremath{\circ}}$ compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{z}$ axis or lying in the $xy$ plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.

arXiv:2311.06245•Physical review. B./Physical review. B•Jan 30, 2024•6 citations

Quantum robustness of the toric code in a parallel field on the honeycomb and triangular lattice

Viktor Kott, Matthias Mühlhauser, Jan Alexander Koziol +1 more

We investigate the quantum robustness of the topological order in the toric code on the honeycomb lattice in the presence of a uniform parallel field. For a field in z z -direction, the low-energy physics is in the flux-free sector and can be mapped to the transverse-field Ising model on the honeycomb lattice. One finds a second-order quantum phase transition in the 3D Ising ^\star ⋆ universality class for both signs of the field. The same is true for a postive field in x x -direction where an analogue mapping in the charge-free sector yields a ferromagnetic transverse-field Ising model on the triangular lattice and the phase transition is still 3D Ising ^\star ⋆ . In contrast, for negative x x -field, the charge-free sector is mapped to the highly frustrated antiferromagnetic transverse-field Ising model on the triangular lattice which is known to host a quantum phase transition in the 3D XY ^\star ⋆ universality class. Further, the charge-free sector does not always contain the low-energy physics for negative x x -fields and a first-order phase transition to the polarized phase in the charge-full sector takes place at larger negative field values. We quantify the location of this transition by comparing quantum Monte Carlo simulations and high-field series expansions. The full extension of the topological phase in the presence of x x - and z z -fields is determined by perturbative linked-cluster expansions using a full graph decomposition. Extrapolating the high-order series of the charge and the flux gap allows to estimate critical exponents of the gap closing. This analysis indicates that the topological order breaks down by critical lines of 3D Ising ^\star ⋆ and 3D XY ^\star ⋆ type with interesting potential multi-critical crossing points. All findings for the toric code on the honeycomb lattice can be transferred exactly to the toric code on the triangular lattice.

arXiv:2402.15389•SciPost Physics•Aug 15, 2024•4 citations

The toric code under antiferromagnetic isotropic Heisenberg interactions

Won Jang, Robert Peters, Thore Posske

We investigate the impact of an isotropic antiferromagnetic Heisenberg perturbation on the toric code, focusing on the resulting quantum phase transition and the nature of the phase that emerges beyond topological order. Using neural-network quantum states (NQS), we compute ground states over a wide range of Heisenberg couplings while fully respecting the exact symmetries of the model. In the weak-coupling regime, the numerical results are in excellent agreement with an effective low-energy description derived from a Schrieffer-Wolff (SW) transformation, providing analytic control over the perturbative breakdown of topological order. We show that the Heisenberg perturbation only renormalizes local operators at low orders, whereas mixing between topological sectors occurs only at a perturbative order proportional to the system size. At intermediate values of the Heisenberg interaction, the topological phase breaks down. We estimate the critical point through a combination of the fidelity susceptibility and the logarithmic susceptibility of non-contractible Wilson loops for various system sizes. Furthermore, we utilize the topological entanglement entropy to provide a comprehensive characterization of the phase transition. Beyond the transition, an antiferromagnetic $\pm X/\pm Z$ Néel phase emerges, characterized by a fourfold-degenerate symmetry-broken manifold, which is explicitly probed using staggered-magnetization-based diagnostics. Our results show how local two-spin interactions, which naturally arise in realistic implementations of the toric code, drive the breakdown of topological order. Moreover, we establish the SW approach as a systematic framework for analyzing such perturbations in combination with variational many-body methods.

cond-mat.str-elCross-listedarXiv:2603.05707•Physical Review Research•Mar 9, 2026•0 citations

Quantum phase transition from Z2×Z2 to Z2 topological order

Mohammad Hossein Zarei

Although the topological order is known as a quantum order in quantum many-body systems, it seems that there is not a one-to-one correspondence between topological phases and quantum phases. As a well-known example, it has been shown that all one-dimensional (1D) quantum phases are topologically trivial [X. Chen et al., Phys. Rev. B 83, 035107 (2011)]. By such facts, it seems a challenging task to understand when a quantum phase transition between different topological models necessarily reveals different topological classes of them. In this paper, we make an attempt to consider this problem by studying a phase transition between two different quantum phases which have a universal topological phase. We define a Hamiltonian as interpolation of the toric code model with ${Z}_{2}$ topological order and the color code model with ${Z}_{2}\ifmmode\times\else\texttimes\fi{}{Z}_{2}$ topological order on a hexagonal lattice. We show such a model is exactly mapped to many copies of 1D quantum Ising model in transverse field by rewriting the Hamiltonian in a new complete basis. Consequently, we show that the universal topological phase of the color code model and the toric code model reflects in the 1D nature of the phase transition. We also consider the expectation value of Wilson loops by a perturbative calculation and show that behavior of the Wilson loop captures the nontopological nature of the quantum phase transition. The result on the point of phase transition also shows that the color code model is strongly robust against the toric code model.

arXiv:1601.06506•Physical Review A•Apr 4, 2016•8 citations
We investigate the quantum phase diagram of the KKK-layer Ising toric code corresponding to KKK layers of two-dimensional toric codes coupled by Ising interactions. While for small Ising interactions the system displays Z2K{\mathbb{Z}}_{2}^{K}Z2K​ topological order originating from the toric codes in each layer, the system shows Z2{\mathbb{Z}}_{2}Z2​ topological order in the high-Ising limit. The latter is demonstrated for general KKK by deriving an effective low-energy model in KthK\mathrm{th}Kth-order degenerate perturbation theory, which is given as an effective anisotropic single-layer toric code in terms of collective pseudospins 1/2 referring to the two ground states of isolated Ising chain segments. For the specific cases K=3K=3K=3 and K=4K=4K=4 we apply high-order series expansions to determine the gap series in the low- and high-Ising limit. Extrapolation of the elementary energy gaps gives convincing evidence that the ground-state phase diagram consists of a single quantum critical point and our findings suggest a quantum phase transition in the 3D Ising* universality class for both KKK separating both types of topological order, which is consistent with former findings for the bilayer Ising toric code.
We demonstrate that two toric code layers on the square lattice coupled by an Ising interaction display two distinct phases with intrinsic topological order. The second-order quantum phase transition between the weakly coupled ${\mathbb{Z}}_{2}\ifmmode\times\else\texttimes\fi{}{\mathbb{Z}}_{2}$ and the strongly coupled Z2{\mathbb{Z}}_{2}Z2​ topological order can be described by the condensation of bosonic quasiparticles from both sides and belongs to the three-dimensional Ising∗{\mathrm{Ising}}^{*}Ising∗ universality class. This can be shown by an exact duality transformation to the transverse-field Ising model on the square lattice, which builds on the existence of an extensive number of local Z2{\mathbb{Z}}_{2}Z2​ conserved parities. These conserved quantities correspond to the product of two adjacent star operators on different layers. Notably, we show that the low-energy effective model derived about the limit of large Ising coupling is given by an effective single-layer toric code in terms of the conserved quantities of the Ising toric code bilayer. The two topological phases are further characterized by the topological entanglement entropy which serves as a nonlocal order parameter.
We investigate the low-energy physics of non-Hermitian quantum spin models with PTPTPT symmetry. To this end, we consider the one-dimensional Ising chain and the two-dimensional toric code in a non-Hermitian staggered field. For both systems, dual descriptions in terms of non-Hermitian staggered Ising interactions in a conventional transverse field exist. We perform high-order series expansions about the high- and low-field limit for both systems to determine the ground-state energy per site and the one-particle gap. The one-dimensional non-Hermitian Ising chain is known to be exactly solvable. Its ground-state phase diagram consists of second-order quantum phase transitions, which can be characterized by logarithmic singularities of the second derivative of the ground-state energy and, in the symmetry-broken phase, the gap closing of the low-field gap. In contrast, the gap closing from the high-field phase is not accessible perturbatively due to the complex energy and the occurrence of exceptional lines in the high-field gap expression. For the two-dimensional toric code in a non-Hermitian staggered field, we study the quantum robustness of the topologically ordered phase by the gap closing of the low-field gap. We find that the well-known second-order quantum phase transition of the toric code in a uniform field extends into a large portion of the non-Hermitian parameter space. However, the series expansions become unreliable for a dominant anti-Hermitian field. Interestingly, the analysis of the high-field gap reveals the potential presence of an intermediate region.
We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped Z2{\mathbb{Z}}_{2}Z2​ spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated ${120}^{\ensuremath{\circ}}$ compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{z}$ axis or lying in the xyxyxy plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.
We investigate the quantum robustness of the topological order in the toric code on the honeycomb lattice in the presence of a uniform parallel field. For a field in z z -direction, the low-energy physics is in the flux-free sector and can be mapped to the transverse-field Ising model on the honeycomb lattice. One finds a second-order quantum phase transition in the 3D Ising ^\star ⋆ universality class for both signs of the field. The same is true for a postive field in x x -direction where an analogue mapping in the charge-free sector yields a ferromagnetic transverse-field Ising model on the triangular lattice and the phase transition is still 3D Ising ^\star ⋆ . In contrast, for negative x x -field, the charge-free sector is mapped to the highly frustrated antiferromagnetic transverse-field Ising model on the triangular lattice which is known to host a quantum phase transition in the 3D XY ^\star ⋆ universality class. Further, the charge-free sector does not always contain the low-energy physics for negative x x -fields and a first-order phase transition to the polarized phase in the charge-full sector takes place at larger negative field values. We quantify the location of this transition by comparing quantum Monte Carlo simulations and high-field series expansions. The full extension of the topological phase in the presence of x x - and z z -fields is determined by perturbative linked-cluster expansions using a full graph decomposition. Extrapolating the high-order series of the charge and the flux gap allows to estimate critical exponents of the gap closing. This analysis indicates that the topological order breaks down by critical lines of 3D Ising ^\star ⋆ and 3D XY ^\star ⋆ type with interesting potential multi-critical crossing points. All findings for the toric code on the honeycomb lattice can be transferred exactly to the toric code on the triangular lattice.
We investigate the impact of an isotropic antiferromagnetic Heisenberg perturbation on the toric code, focusing on the resulting quantum phase transition and the nature of the phase that emerges beyond topological order. Using neural-network quantum states (NQS), we compute ground states over a wide range of Heisenberg couplings while fully respecting the exact symmetries of the model. In the weak-coupling regime, the numerical results are in excellent agreement with an effective low-energy description derived from a Schrieffer-Wolff (SW) transformation, providing analytic control over the perturbative breakdown of topological order. We show that the Heisenberg perturbation only renormalizes local operators at low orders, whereas mixing between topological sectors occurs only at a perturbative order proportional to the system size. At intermediate values of the Heisenberg interaction, the topological phase breaks down. We estimate the critical point through a combination of the fidelity susceptibility and the logarithmic susceptibility of non-contractible Wilson loops for various system sizes. Furthermore, we utilize the topological entanglement entropy to provide a comprehensive characterization of the phase transition. Beyond the transition, an antiferromagnetic ±X/±Z\pm X/\pm Z±X/±Z Néel phase emerges, characterized by a fourfold-degenerate symmetry-broken manifold, which is explicitly probed using staggered-magnetization-based diagnostics. Our results show how local two-spin interactions, which naturally arise in realistic implementations of the toric code, drive the breakdown of topological order. Moreover, we establish the SW approach as a systematic framework for analyzing such perturbations in combination with variational many-body methods.
Although the topological order is known as a quantum order in quantum many-body systems, it seems that there is not a one-to-one correspondence between topological phases and quantum phases. As a well-known example, it has been shown that all one-dimensional (1D) quantum phases are topologically trivial [X. Chen et al., Phys. Rev. B 83, 035107 (2011)]. By such facts, it seems a challenging task to understand when a quantum phase transition between different topological models necessarily reveals different topological classes of them. In this paper, we make an attempt to consider this problem by studying a phase transition between two different quantum phases which have a universal topological phase. We define a Hamiltonian as interpolation of the toric code model with Z2{Z}_{2}Z2​ topological order and the color code model with ${Z}_{2}\ifmmode\times\else\texttimes\fi{}{Z}_{2}$ topological order on a hexagonal lattice. We show such a model is exactly mapped to many copies of 1D quantum Ising model in transverse field by rewriting the Hamiltonian in a new complete basis. Consequently, we show that the universal topological phase of the color code model and the toric code model reflects in the 1D nature of the phase transition. We also consider the expectation value of Wilson loops by a perturbative calculation and show that behavior of the Wilson loop captures the nontopological nature of the quantum phase transition. The result on the point of phase transition also shows that the color code model is strongly robust against the toric code model.