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 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.
The robustness of the topological color code, which is a class of error-correcting quantum codes, is investigated under the influence of a uniform magnetic field on the honeycomb lattice. Our study relies on two high-order series expansions using perturbative continuous unitary transformations in the limit of low and high fields, exact diagonalization, and a classical approximation. We show that the topological color code in a single parallel field is isospectral to the Baxter-Wu model in a transverse field on the triangular lattice. It is found that the topological phase is stable up to a critical field beyond which it breaks down to the polarized phase by a first-order phase transition. The results also suggest that the topological color code is more robust than the toric code in the parallel magnetic field.
We investigate the stability of the topological phase of the toric code model in the presence of a uniform magnetic field by means of variational and high-order series expansion approaches. We find that when this perturbation is strong enough, the system undergoes a topological phase transition whose first- or second-order nature depends on the field orientation. When this transition is of second order, it is in the Ising universality class except for a special line on which the critical exponent driving the closure of the gap varies continuously, unveiling a new topological universality class.
A quantum Ising model in a transverse field is of the simplest quantum many-body systems used for studying universal properties of quantum phase transitions. Interestingly, it is well known that such phase transitions can be mapped to topological phase transitions in toric code models. Therefore, one can expect that well-known properties of the transverse Ising model are used for characterizing topological phase transitions in toric code models. In this paper, we consider the magnetization of the Ising model and show, while it is a local order parameter, it is mapped to a nonlocal order parameter on the topological side. Consequently, we introduce a string order parameter for topological phase transitions in toric code models defined on different lattices with different dimensions in the presence of a uniform magnetic field. Since such a string order parameter is dual of the Ising order parameter with well-known properties, it leads to a topological (nonlocal) characterization of phase transition in the toric code model in uniform field. Our results show that although topological phase transitions do not follow a symmetry-breaking mechanism, it might still be possible to use concepts of symmetry-breaking phase transitions for topological ones by finding suitable mappings.
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 Z2K topological order originating from the toric codes in each layer, the system shows Z2 topological order in the high-Ising limit. The latter is demonstrated for general K by deriving an effective low-energy model in 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=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.
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.
Kitaev’s toric code has become one of the most studied models in physics and is highly relevant to the fields of both quantum error correction and condensed matter physics. Most notably, it is the simplest known model hosting an extended, deconfined topological bulk phase. To this day, it remains challenging to reliably and robustly probe topological phases, as many state-of-the-art order parameters are sensitive to specific models and even specific parameter regimes. With the emergence of powerful quantum simulators which are approaching the regimes of topological bulk phases, there is a timely need for experimentally accessible order parameters. Here we study the ground state physics of the parallel field toric code on the honeycomb, triangular, and cubic lattices using continuous-time quantum Monte Carlo. By extending the concept of experimentally accessible percolation-inspired order parameters (POPs) we show that electric and magnetic anyons are independently confined on the honeycomb and triangular lattices, unlike on the square lattice. Our work manifestly demonstrates that, even in the ground state, we must make a distinction between topological order and (de-)confinement. Moreover, we report multicritical points in the aforementioned confinement phase diagrams. Finally, we map out the topological phase diagrams on the honeycomb, triangular, and cubic lattices and compare the performance of the POPs with other topological order parameters. Our work paves the way for studies of confinement involving dynamical matter and the associated multicritical points in contemporary quantum simulation platforms for Z 2 lattice gauge theories.