Twice Hidden String Order and Competing Phases in the Spin-1/2 Kitaev-Gamma Ladder

,


INTRODUCTION
Over the course of investigating spin S=1/2 twodimensional (2D) honeycomb Kitaev materials [1][2][3][4][5], candidates of long-sought quantum spin liquids, an additional bond-dependent spin exchange term named Gamma (Γ) interaction [6] was found along with the bond-dependent Kitaev (K) interaction [7,8].Unlike the standard Heisenberg interaction, both the Kitaev and the Γ interactions are highly nontrivial and extremely frustrated.While the same sign of K and Γ cancels the frustration leading to magnetically ordered phases, the regions with different signs of these interactions are highly frustrated.[5] Since these two interactions are known to dominate over other symmetryallowed interactions in the emerging candidate material α-RuCl 3 , the combined Kitaev-Gamma (KΓ) model has attracted considerable attention in the theoretical community and it is widely accepted [9][10][11] that a realistic description of α-RuCl 3 should be sought in the regime with antiferromagnetic (AFM) Gamma (Γ>0) and ferromagnetic (FM) Kitaev (K<0) interactions.Understanding the phases that arise in this frustrated regime of the 2D honeycomb KΓ-model is therefore of crucial importance.
This particular region of the KΓ-model has been extensively studied by a range of numerical methods [5].Despite detailed studies, the nature of the phase next to the FM Kitaev spin liquid arising due to AFM Γ interaction remains controversial.
Motivated by such discrepancy, another approach to investigate the 2D limit of the KΓ model was taken by starting from low-dimensional models with the hope of furthering the understanding of the honeycomb model by determining the phases of n-leg (brick-wall) models.Despite the obvious challenge in connecting the two limits, it is reasonable to expect that potential spin liquid phases arising in the honeycomb model should correspond to regions where the n-leg models display disordered phases.Such an approach was employed earlier for the pure Kitaev model [22].Disordered phases in the anisotropic Kitaev 1-leg chain were found, and they were characterized by non-local string order parameters (SOPs) [22].It has also been shown that the isotropic Kitaev 2-leg ladder model exhibits a disordered phase, characterized by an unconventional SOP different from that of the anisotropic chain Kitaev phase [23].
The one-dimensional (1D) chain and ladder version of the KΓ model were investigated numerically with very high precision, as the reduced dimensionality allows accessing bigger system sizes [24][25][26][27][28].In the 1D chain model, it was found that the pure Gamma model belongs to a Luttinger liquid phase governed by the gapless hidden SU(2) Heisenberg chain, a fact revealed after a 6-site transformation, i.e, a duality mapping [24,29].A study of the same KΓ model on a quasi 1D ladder using DMRG and iDMRG techniques [18,30] found a magnetically disordered phase, possessing a small gap near the AFM pure Gamma limit.This phase surrounding the pure AFM Gamma point next to the FM Kitaev phase (denoted by FK) was referred to as the AΓ phase.Even though the AΓ phase in the ladder occurs in the same part of the phase diagram as the proposed KΓSL in the 2D limit, a distinct name was introduced, as < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 P i M = " > A A A B 6 H i c b V A 9 S w N B E J 3 z M 8 a v q K X N Y h C s w p 2 I W g Z s L C w S M B + Q H G F v M 5 e s 2 d s 7 d v f E c O Q X 2 F g o Y u t P s v P f u E m u 0 M Q H A 4 / 3 Z p i Z F y S C a + O 6 3 8 7 K 6 t r 6 x m Z h q 7 i 9 s 7 u 3 X z o 4 b O o 4 V Q w b L B a x a g d U o + A S G 4 Y b g e 1 E I Y 0 C g a 1 g d D P 1 W 4 + o N I / l v R k n 6 E d 0 I H n I G T V W q j / 1 S m W 3 4 s 5 A l o m X k z L k q P V K X 9 1 + z N I I p W G C a t 3 x 3 M T 4 G V W G M 4 G T Y j f V m F A 2 o g P s W C p p h N r P Z o d O y K l V + i S M l S 1 p y E z 9 P Z H R S O t x F N j O i J q h X v S m 4 n 9 e J z X h t Z 9 x m a Q G J Z s v C l N B T E y m X 5 M + V 8 i M G F t C m e L 2 V s K G V F F m b D Z F G 4 K 3 + P I y a Z 5 X v M u K V 7 8 o V + / y O A p w D C d w B h 5 c Q R V u o Q Y N Y I D w D K / w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + s z j Q w = < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y K U K 1 r G + N P o K 6 x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y x < l a t e x i t s h a 1 _ b a s e 6 4 = " L Z W 6 t 2 q 9 S i 0 y l 0 X p x 3 5 2 P R W n D y m W P 4 A + f z B + y 3 j Q 0 = < / l a t e x i t > y < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 w 5 j w 4 L 8 6 7 8 z F v X X H y m S P 4 A + f z B + 4 7 j Q 4 = < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " y 3 d 4 8 w p s z d l 6 c d + d j 0 V p w 8 p l j + A P n 8 w f 5 j I 0 H < / l a t e x i t > 18 < l a t e x i t s h a 1 _ b a s e 6 4 = " q 4 1 q g + q Y k 1 i q j Y J 3 0

22
< l a t e x i t s h a 1 _ b a s e 6 4 = " q t l U a 2 P y j + y 8 O e / z 0 Z J T 7 O y j X 3 A + v g F q t J F b < / l a t e x i t > AK < l a t e x i t s h a 1 _ b a s e 6 4 = " P l E V 8 N e e w p i c L Y L x M J R W 4 N 4 A W r o = " > A A A B + n i c b V B P S 8 M w H E 3 n v z n / d X r 0 E h y C p 9 G K T o 9 D Q b w I E + w 2 2 E p J s 2 w L S 9 K S p M q o / S h e P C j i 1 U / i z W 9 j u v W g 0 w e B x 3 u / H 7 + X F 8 a M K u 0 4 X 1 Z p a X l l d a 2 8 X t n Y 3 N r e s a u 7 b R U l E h M P R y y S 3 R A p w q g g n q a a k W 4 s C e I h I 5 1 w c p n 7 n A n e 4 p e X S e u s 5 l 3 W L h 7 P q / X b o o 4 y O k L H 6 B R 5 6 A r V 0 Q N q o C b C a I y e 0 S t 6 s z L r x X q 3 P u a j J a v Y O U R / Y H 3 + A K u a k 7 I = < / l a t e x i t > A < l a t e x i t s h a 1 _ b a s e 6 4 = " r o N v j F a f x t 9 p f P w 6 S 1 a a t i S n A X v / y X t I / r b q N + e n 1 S a 1 4 U d Z T B P j g A R 8 A F Z 6 A J r k A L e A C D B / A E X s C r 9 W g 9 W 2 / W + 3 y 0 Z B U 7 e + A X r I 9 v c j i U I g = = < / l a t e x i t >

RS U6
< l a t e x i t s h a _ b a s e = " e i m K w + e A J y n Z P P e n H F g v M n a y u 8 W c p 6 s d 6 t j 2 V r z s p m T u E P r M 8 f d b + O 1 g = = < / l a t e x i t > 0.88271 (5) < l a t e x i t s h a 1 _ b a s e 6 4 = " J / r E q Z B Y M 9 N n m 8 w A m 8 g F c r t 5 6 t N + t 9 M V q x y p 0 D 8 A v W x z c 7 Q Z S U < / l a t e x i t > SPT < l a t e x i t s h a 1 _ b a s e 6 4 = " o m e e c 0 E Z a A strip of the KΓ honeycomb lattice corresponding to a two-leg KΓ ladder with alternating x and y bonds along the leg and z-bond between the chains, with the numbering of the sites used throughout the paper.The dotted z-bonds arise from imposing periodic boundary conditions in the direction perpendicular to the ladder.However, all z-bonds are taken to have equal strength.The dashed red line indicates the partition used for ρB while the dashed blue line indicates the partition used for ρA.b Schematic phase diagram of the KΓ ladder, Eq. ( 1), for Γ>0.
it is yet to be determined how the AΓ phase is connected to the 2D limit KΓSL.
While it is clear that there is no magnetic order in this phase, the precise nature of the AΓ phase has not yet been settled due to its complex nature.The presence of a gap indicates that the AΓ-phase is likely a symmetry protected topological (SPT) phase [31][32][33].If so, it is of importance to identify a corresponding SOP, edge states, and the symmetry that protects this phase, which are characteristic of the SPT, and to determine how this phase respond to an external magnetic field.If strong evidence for a non-trivial SPT nature of the AΓ phase can be established, it would establish a next step to the proposed KΓSL in the 2D limit.In the following sections, we will systematically examine these inquiries and provide comprehensive responses.To perform a thorough analysis, we start by reviewing the full phase diagram prior to focusing on the nature of the AΓ phase.As shown below, we also report two other disordered phases.

Model
The KΓ Hamiltonian is given by where (α, β) takes on the values (y, z)/(x, z)/(x, y) for γ = x/y/z, and ⟨i, j⟩ refers to the nearest neighbor sites.
An alternative representation of the honeycomb lattice is as a brick-wall lattice [22], and its two-leg limit with periodic boundary conditions is simply a ladder shown in Fig. 1a.The dotted bonds indicate Kitaev z-bonds arising from periodic boundary conditions, and we shall always take such bonds to be identical to the regular (solid) z-bonds, in which case the honeycomb strip can be viewed as a regular rectangular ladder.We parameterize the model by taking K = cos ϕ, and Γ= sin ϕ, and interpolate between the Kitaev and Γ interactions by varying ϕ.Our main interest is in the region with ϕ/π ∈ [0, π] where Γ>0 and the Kitaev term, K, changes from AFM to FM at ϕ=π/2, as this region is relevant to most two-dimensional (2D) Kitaev candidate materials.The total number of sites in the ladder (including both legs) is denoted by N .

Phase Diagram
A full phase diagram is shown in Fig. 1b.It is obtained by various quantities presented in the next subsection.Moving from ϕ=0 to π, the AFM Kitaev phase (denoted by AK), a FM phase denoted by FM U6 , AΓ, and FK phases are found consistent with the earlier works [18,30].At the special point ϕ = π/4, the FM U6 phase can be mapped to the ferromagnetic Heisenberg ladder by applying a local unitary U 6 transformation (see Supplementary Note 1), and is therefore gapless.However, for ϕ̸ =π/4 a small gap appear, as we show in Supplementary Note 3. Surprisingly, two additional phases denoted by SPT α and SPT β can be identified between the FM U6 and AΓ phases.As we will show below, they are magnetically disordered and display the characteristics of SPT phases, i.e., doubled entanglement spectrum.Beyond ϕ = π, the rung singlet phase denoted by RS U6 phase delineates the FK phase.
In addition to the expected Kitaev phases AK and FK the appearance of the FM U6 and RS U6 phases are well established in the KΓ honeycomb and n-leg models.After the local U 6 transformation [29], corresponding to local spin rotations, at ϕ=π/4 and 5π/4, i.e., K=Γ, the KΓ 2D honeycomb model is equivalent to the FM and AFM Heisenberg model, respectively.The application of the U 6 transformation is specified in Supplementary Note 1.To understand the nature of the other three phases, AΓ, SPT α and SPT β , we first performed a detailed analysis of the entanglement spectrum.

Entanglement Spectrum
Our results for the entanglement spectrum, as well as for the susceptibility, χ e ϕ are shown in Figs. 2 and 3. Due to the complexity of the phase diagram, we split To the left of the AΓ-phase, two other phases SPT α and SPT β are clearly separated from the FM U6 and AΓ, as noted from their entanglement spectrum for ρ A .Furthermore, the double degeneracy of the entanglement spectrum shown in Fig. 2 (b) from ρ B is a clear signal of a SPT phase [33][34][35][36][37].Note that there is only a weak signature in χ e ϕ of the transition between the two phases, visible at ϕ=0.3956π in Fig. 2(c), corresponding to a small discontinuity in χ e ϕ .While the FM U6 phase appears abruptly at ϕ=0.3823π, the transitions between the AΓ phase and SPT β phase at ϕ ∼ 0.435π is not discernible in χ e ϕ .The blurriness of the transition is likely due to a field induced phase that pinches off to a single point at zero field at the AΓ-SPT β transition, thereby obscuring it.We note that, the quantum critical points (QCPs) are immediately noticeable in the entanglement spectra.
Moving to the right within the AΓ phase which encompass the point Γ =1, the transition to the FK phase from the AΓ phase occurs at ϕ=0.88271(5)π, which is clearly visible in the entanglement spectra as well as χ e ϕ as shown in in Fig. 3.The transition between the RS U6 and FK phases at ϕ=π is less apparent in the structure of the entanglement spectrum in Fig. 3 (a) and (b).However, from the results for χ e ϕ in Fig. 3(c) this transition is immediately visible, and it was also noted that the second derivative of λ 1 clearly detects the transition [30].
Similar to the AK phase, none of the phases SPT α , SPT β , AΓ and FK has any long-range magnetic ordering, nor is there any indication of nematic (quadropolar) or chiral ordering.As discussed in Supplementary Note 2, all four phases are gapped with a finite sizeable correlation length.The difference between them is captured in the entanglement spectrum.In the SPT α , SPT β and FK (and AK) phases, the entanglement spectrum have all entries doubled when considering ρ B .For the AΓ-phase, the same applies to the spectrum for ρ A .Since an entanglement spectrum where all eigenvalues have degeneracy larger than one is a signature of a topological non-trivial phase, in the next section, we investigate the projective symmetry analysis to confirm their non-trivial topology prior to presenting associated non-local SOPs.
< l a t e x i t s h a 1 _ b a s e 6 4 = " w K e t 7 r 9 s / 2 u j + y h w H D R t J S F 9 t y b k a 8 C U h 8 w E P x f B 9 5 y I 9 H d k x J w k 3 z 9 2 s r 7 7 3 / w Y c f r X 5 c + + T T z 6 6 v 3 a j f f J
For the KΓ ladder, it is a significant simplification to consider the transformed model obtained after applying the U 6 transformation.This transformation maps the original Hamiltonian H KΓ to the transformed ladder, denoted by H U6 KΓ .Under the U 6 transformation the x, y and z-bonds of the KΓ-ladder are transformed into anisotropic Heisenberg bonds, x ′ , y ′ and z ′ in the following manner: In the transformed ladder H U6 KΓ model, the definition of x ′ , y ′ and z ′ -bonds, can pictorially be represented as shown in Fig. 4(a).See Supplementary Note 1.We will consider two different open boundary conditions (OBC) as shown in Fig. 4. The unit cells depicted in panel (a) are referred to as regular unit cells with regular OBC, while the ones in panel (b) are referred to as slanted unit cells with slanted OBC.
In order to understand the symmetries of the matrix product state (MPS) wave-function, it is useful to write it in the canonical form [49][50][51][52]: where the M [n] jn are complex matrices and the M [n] , real, positive, square diagonal matrices.In the iDMRG formulation, the set of matrices on any unit cell becomes the same M [n] j =M j , M [n] = M for all n, although they may vary within the unit cell.For the translationally invariant state, it can be shown [34,53] that for any (site) symmetry operation g, represented in the spin basis by the unitary matrix, Σ jj ′ (g), the M j matrices must transform as [34,48]: where the unitary matrix U (g) commutes with the M matrices and e iθ is a phase factor.With D denoting the bond dimension, the U matrices form a D-dimensional projective representation of the symmetry group of the wave-function, and they can be determined from the unique eigenvector of the generalized transfer matrix [34,48] with eigenvalue |λ|=1, where the generalized transfer matrix is defined as If the largest eigenvalue is |λ| <1, the symmetry is not a property of the state being considered.The projective representation is reflected in the fact that if Σ(g)Σ(h)=Σ(gh), then where the phases ϕ(g, h) are characteristic of the topological phase.
Let us now consider the site symmetries, R x and R y defined as: The Hamiltonian H U6 KΓ is invariant under the operators l R γ l , with γ=x, y, z with distinct quantum numbers for the low-lying states.If these symmetries are respected, their representations can differ by a phase, ϕ(R x , R y ) that must obey e iϕ(Rx,Ry) =±1: Furthermore, the non-trivial value e iϕ(Rx,Ry) =-1 can only occur if all eigenvalues of the entanglement spectrum are at least twice degenerate [34].The phase factor can then be isolated by defining [48]: with D the bond dimension, with similar definitions for other pairs of operators in Eq. ( 8).In the above definitions it is understood that the transformations are applied throughout the lattice and in order to obtain the matrices U , generalized transfer matrices representing the relevant unit cell has to be considered.
The site symmetries, R x , R y and R z , forming the dihedral group, D 2 , are respected by both of the unit cells in Fig. 4(a) and 4(b).Using generalized transfer matrices obtained from unit cells of the shape shown in Fig 4(a) when studying the AΓ-phase and of the slanted shape shown in Fig. 4(b) when studying the FK, SPT α , and SPT β phases, we obtain for the SPT α , SPT β , AΓ, and FK phases.Similar analysis can be made for the time-reversal (TR) symmetry, defined by M j → j ′ e iπS y jj ′ M ⋆ j ′ , with ⋆ denoting complex conjugation.In this case, it can be established that [34] U TR U ⋆ TR =e iϕ(TR,TR) 1 where the phase ϕ(TR, TR) cannot be absorbed into the definition of U TR .One should note that for most other symmetries, with the notable exception of inversion, similar considerations will lead to U 2 =e iϕ 1 in which case the phase ϕ in fact can be absorbed into the definition of U .For instance, this is the case for U (R α ) discussed above.However, for time reversal, the phase factor e iϕ(TR,TR) can directly be extracted by defining [48]: and again one finds that ϕ(TR, TR)=0,π, so that O TR =±1.As an example, for the S=1 Heisenberg spin chain in the Haldane phase it is known that [34,35].Using generalized transfer matrices obtained from unit cells of the shape shown in Fig 4(a) for the AΓ-phase and of the slanted shape for the FK, SPT α , and SPT β phases, we obtain consistent with the presence of a doubled entanglement spectrum in all cases.As was the case for O Z2×Z2 , if the unit cells are interchanged, one finds instead O TR =1.

TABLE I. Summary of projective analysis.
The superscript regular refers to the unit cell from Fig. 4a, while the superscript slanted refers to the unit cell from Fig. 4b.Negative values indicate that the state transforms non-trivially.
A summary of our results from the projective analysis are provided in Table I, negative values indicate that the state transforms non-trivially.For all 4 phases, it is seen that a unit cell can be chosen for which the state transforms non-trivially under both the TR and O Z2×Z2 symmetries.
Based on the above analysis of entanglement spectrum and projective symmetry, we conclude that AΓ phase is an SPT phase.It is then important to further identify its SOP that differentiates this phase from the other disordered phases.

Twice Hidden String Order
To establish a non-local string order parameter (SOP) characterizing the AΓ phase, we need to exploit a non-local unitary transformation that maps the original Hamiltonian with OBC to a new Hamiltonian that exhibits a local long-range order [54][55][56].We found that it is difficult to identify such non-local transformation starting from the original H KΓ , but it can be achieved by first applying the U 6 transformation to arrive at H U6 KΓ .It is then possible to define a non-local unitary transformation W , mapping H U6 KΓ to a new local Hamiltonian.We denote the resulting Hamiltonian, where four-spin terms appear, by H KQ .
For the parameters relevant for the AΓ-phase, H KQ exhibits long-range order in the spin-spin correlation functions, corresponding to a local order parameter.Due to the application of two separate unitary transformations, one might consider the resulting order to be twice hidden.
The non-local unitary operator W for a N -site ladder with OBC that maps H U6 KΓ to H KQ takes the following form.
With the individual w(j, k) given as follows: and W † =W .The OBC are here crucial for the mapping to be exact.Evidently, all w(j, k) are unitary and therefore also W , and [w(j, k), w(l, m)] = 0 ∀ j, k, l, m.Note that this is a different labelling of the unitary operator introduced in Ref. [23,30].It can also be shown that other combinations of the spin operators S α appearing in Eq. ( 15) lead to equivalent unitary operators, for instance, e iπ(S y j +S y j+1 )•(S x k +S x k+1 ) is another valid choice.However, the specific choice made in Eq. ( 15) will influence the type of ordering that is observed in H KQ , as well as the specific form of H KQ .Schematically, the transformations can be viewed as shown in Fig. 5.The detailed form of H KQ after the W transformation on H U6 KΓ is presented in Supplementary Note 4.
We first discuss the ordering in H KQ where we denote the spins by S ′′α , where the double prime represents the two transformations of spin from the original Hamiltonian.It is then convenient to study correlation functions of the form 4⟨S ′′α (0)S ′′α (r)⟩ along   the legs of the ladder with r measured in lattice spacings along the leg.To avoid boundary effects, r=0 is usually taken to correspond to a site in the bulk of the chain.In Fig. 6 we show results for 4⟨S ′′y (0)S ′′y (r)⟩ starting from site 47 with ϕ=0.85π.Long-range order is clearly present.Similar results can be obtained for the other leg of the ladder as well as for 4⟨S ′′z (0)S ′′z (r)⟩.However, due to the choice of spin operators in the definition of w in Eq. ( 15), there is no ordering in 4⟨S ′′x (0)S ′′x (r)⟩.
Using the inverse of the non-local unitary operator W from Eq. ( 14) the above results for 4⟨S ′′y (0)S ′′y (r)⟩ is reproduced as a non-local string order correlation function in H U6 KΓ where we denote the spin variables by S ′α =σ ′α /2.We find ⟨S ′′y in Fig. 6 with the string correlation function in ( 16) for H U6 KΓ , a relabelling of the sites needs to be done that we have skipped for clarity.
We can now apply the inverse U 6 transformation (see Supplementary Note 1) to the above expressions for O y (r) to determine the string order correlation functions that define the AΓ-phase in the original H KΓ : with the U −1 6 transformation detailed in Supplementary Note 1.The SOP in H KΓ Hamiltonian is then given by It is interesting to note that for a small r, for example r = 7 Ξ y (r = 7) corresponds to the plaquette operator found in the pure Gamma model in the honeycomb lattice.[57] In Fig. 7 we show iDMRG results for Ξ y (orange circles).The inset shows iDMRG results for Ξ y (r) versus r at ϕ=0.85π which at large r can be compared to the results for H KQ shown in Fig. 6.Note that the results in Fig. 6 show the correlations along a leg, whereas the inset in Fig. 7 show results from both legs without including the sign and the relabelling of the sites.Due to the different methods used, there is not an exact equivalence for very small values of r.
We emphasize that the appearance of a non-local SOP in the AΓ phase of H KΓ is equivalent to the presence of long-range order in H KQ .Hence, Ξ y is non-zero throughout the AΓ-phase and goes to zero at the critical points delineating this phase.It is absent in the other disordered phases, SPT α , SPT β , and FK, and thus uniquely defines the AΓ phase.
With the identification of the AΓ-phase with regular long-range ordering in the H KQ model, it is natural to ask if a regular (local) order parameter can also be identified for the SPT α -and SPT β -phases in the H KQ .However, all local order parameters that we have investigated have not shown any ordering in the SPT α -and SPT β -phases for H KQ .Extending the string-order correlation function defined for the S=1 Haldane chain, a heuristic stringorder correlation function has been proposed for S=1/2 ladders by pairing two S=1/2 [58,59].Following the numbering of Fig. 1a, if τ α i =S α 2i + S α 2i+1 are the sum of two diagonally situated spins, one defines [58,59]: The associated SOP is non-zero in the phase surrounding ϕ=0 in H KQ [30], corresponding to the AFM Kitaev (AK) phase in H KΓ .The magenta points in Fig. 7 show our results for O z even for the H KQ model, which clearly is non-zero in the SPT α -and SPT β -phases.This is consistent with the nonexistence of a local order parameter in these two phases for H KQ .We note that, due to the heuristic nature of O z even , it is not clear how to associate it with a local order in a related model.Since H KΓ and H U6 KΓ are related by a local unitary transformation, any ordering in the SPT α -and SPT βphases in either model would immediately be apparent in both, and we have not observed any for either model.
Building on the above results for O z even for H KQ we propose a closely related heuristic string order correlation function for H U6 KΓ in the following way: Define 2i+1 as the difference of two diagonally situated spins, following the numbering from Fig. 1a.We then have with τ α i =S α i,1 + S α i+1,2 as above.Results for P z even versus ϕ obtained from iDMRG calculations with H U6 KΓ are shown in Fig. 7 as the light blue points.The FK, SPT α and SPT β -phases are clearly defined by a nonzero P even .By applying U −1 6 to the definition of P even it is straightforward to perform similar calculations using H KΓ by evaluating U −1 6 (P α even (r)).Even though the definition P even is heuristic, we interpret this result as a verification of the SPT nature of the FK, SPT α -and SPT β -phases.
For the FK-phase, it is also instructive to consider an even simpler heuristic string order correlation function, Z defined as follows [22]: We consider this string correlator for H U6 KΓ or equivalently to H KΓ through application of U −1 6 .Results for Z are shown in Fig. 7 versus ϕ (green points) as obtained from iDMRG calculations.Throughout the FK-phase Z is almost identical to 1, dropping to zero at the transition to the AΓ-phase.However, we note that Z remain sizable throughout much of the RS U6 -phase, reflecting its heuristic nature.
Edge states and response to magnetic field Another signature of SPT phases is the presence of edge states under OBC related to a ground state degeneracy.For the SPT phases in the KΓ ladder, it is clear from the degeneracy of the entanglement spectrum that we need to consider different shapes of clusters (regular vs. slanted OBCs) for the different SPT phases.In this section, we exclusively study the original Hamiltonian H KΓ and do not consider H U6 KΓ nor H KQ .For AΓ phase, we use N =4n with the regular OBC and for the remaining SPT phases, we use the slanted OBC with N =4n + 2 in order to have an equal number of the different bond types.
We first demonstrate the presence of edge states in the SPT phases.For the AΓ phase, results for the 16 lowest states with the regular OBC at ϕ=0.85π are obtained using ED (see Supplementary Figure 7).Four low-lying states below the gap are clearly present.With increasing N , these 4 states quickly become degenerate while the gap stabilizes at a finite value.Similar results can be obtained for the AK and FK phase using the slanted cluster, and a degeneracy of 4 is also observed for this case (see Supplementary Figure 8).For the SPT α -and SPT β -phases, it has not been possible to produce reliable results in the same manner, likely due to the significantly larger correlation lengths.
To understand the nature of the edge-states, let us explore how the AK, FK, and AΓ phases respond to an external magnetic field.An external magnetic field introduces an additional term in the Hamiltonian of the form H ′ =g L µ B B • S tot , where S tot = i S i , g L is the Landé factor and µ B the Bohr magneton.Following Ref. [38] we denote the 4 states ψ 1 , ψ 2 , ψ 3 and ψ 4 and consider S tot,α in this four-fold degenerate space by defining Here, the components of the total spin S tot,α are usually taken to be identical to x, y, z but given the underlying honeycomb structure we shall find it useful to instead consider α=â, b, ĉ corresponding to the three directions , [1][2][3][4][5][6][7][8][9][10] and [111] that correspond to the perpendicular and parallel to the z-bond, and perpendicular to the plane of the honeycomb (or nleg brick-wall), respectively.One then finds that the eigenvalues of the matrices (S γβ tot,α ) for four corresponding degenerate states are simply given by (s α , −s α , 0, 0).
Our ED results for s α are shown in Fig. 8.For the AK phase at ϕ= −0.08π, shown in Fig. 8(a), we find s a ∼3/4, s b ∼1/2 and s c ≃0.85.Similarly, for the FKphase at ϕ=0.95π, shown in Fig. 8(b), we find s a =1, s b ∼2 and s c =3/4.In both cases, we expect some variation in the values of s α as ϕ is tuned.We note that for both the FK and AΓ phase, the values of s α quickly saturate at a small, finite value as N is increased.This is indicative of excitations localized at the edges as opposed to an actual magnetically ordered groundstate which should show s α continually growing with N .A calculation of s α for the SPT α and SPT β -phases p / 2 a v Q 1 7 3 S / Z + 3 7 I K F u y l c S S R c s y 0 E J G Y s j k O U e f P l I 8 J D / L v Z F j h x E A / 9 6 4 + c G t D z / 6 + P Y n e 5 9 + 9 v k X X 9 6 5 + 9 W L 0 B 8 H F n 9 u + Y 4 f v O q Z I X d s j z + P 7 M j h r 0 Y B N 9 2 e w 1 / 2 z n 6 I 6 1 + + 4 U F o + 9 5 v 0 c W I d 1 1 z 6 N k D 2 z I j U f T 6 7 i / / G T 0 + t L 3 J w B 6 O A z 7 d U 5 L P k X 3 2 d m R b U V x 4 E l q m w w 9 A h 3 T 3 F P E y + n w g z j c 7 + 6 T n j P l 0 r / a S 6 w c v U A f S D v 0 V 3 X 9 y m P S g 2 6 1 v W t + 2 9 l u w x V p P W j + 1 j l r P W 9 b R z 0 e / H 7 0 9 m h z / c f z n 8 V / H f 8 9 N b 9 5 I f L 5 u X X o d / / M / i P i q Q Q = = < / l a t e x i t > do not yield clear results for the range of N available with ED, as discussed in Supplementary Note 6.In a realistic experimental setting, the presence of impurities will always create finite open segments of ladders, with a resulting Curie-law behavior.The response to an applied magnetic field is in that case highly anisotropic and the low temperature Curie-law response should show a strong directional dependence [38,41,43] with χ a (T ), χ b (T ) and χ c (T ) clearly distinguishable.The results for the AΓ-phase, shown in Fig. 8(c), are even more intriguing.They are not only more anisotropic, but only s b is nonzero, approaching a value close to 2 at ϕ=0.85π.At a slightly different point in the AΓ-phase with ϕ=0.8π we instead find s b ≃4/3 but again only s b is non-zero.This implies that the phase does not respond at all to a field applied along the â and ĉ directions, effectively g a , g c ≃0.
To gain a clearer picture of how the ladder in the AΓ phase responds to a magnetic field applied along the b-direction, we have performed ED calculations in the presence of a small field in the b-direction.The resulting site dependent magnetization ⟨S a,b,c i ⟩ can then easily be obtained for the â, b and ĉ-directions.Our results are shown in Fig. 9 for N =24 as obtained from the ground-state with a small field in the b-direction of 0.002 at ϕ/π=0.85 in the AΓ-phase.The green, blue and cyan colors represent positive expectation values, whereas orange, red and pink colors indicate negative values, with the size of the points proportional to the expectation value.The response along the â, ĉ directions is completely symmetric in the positive and negative directions, yielding i S a,c i =0.However, along the bdirection the response is much larger, and we clearly find (See Fig. S9) (c).With the field along the b-direction, sizeable excitations are visible at both ends of the ladder.The response of the edge states to a magnetic field correlates with the lifting (or absence of lifting) of the degeneracy in the entanglement spectrum.As previously discussed, non-trivial indices in the projective analysis can only arise if the degeneracy of all states in the entanglement spectrum (ES) is larger than one.This implies that if a finite strength of the perturbation is needed to remove the degeneracy, then a phase transition does not occur until that strength is reached and the phase persists till that point.On the other hand, if the degeneracy is lifted for any non-zero strength of the perturbation, the symmetry protection is broken without an associated phase transition.We can then investigate the response of the degeneracy of the ES to a magnetic field in the â, b and ĉ directions.This is shown in Fig. 10 for the AK, FK, and AΓ phases.For simplicity, we focus exclusively on the difference in the two largest eigenvalues ∆=λ 1 −λ 2 , and in each case we employ the reduced density matrix that has a two-fold degeneracy at zero field.
As can be seen in Fig. 10(a) the degeneracy in the ES for the AK phase is immediately lifted by a field in any of the three directions which correlates with the response of the edge-states (See Fig. S9(a) in the SI).However, the response is rather weak, and a relatively large field has to be applied to see a significant splitting.For the FK phase, we have a similar effect as shown in Fig. 10(b), but in this case the response is much stronger.However, for the AΓ phase, where we show results in Fig. 10(c), it is clear that a field in the â direction of around h a =0.06 is needed to lift the degeneracy of the ES.For a field applied in the ĉ direction, a significantly larger field is needed.On the other hand, a field in the b direction immediately induces a large splitting in the entanglement spectrum even for infinitesimal field strengths.Since the degeneracy remains intact in the â and ĉ directions, we conclude that the SPT character of the AΓ phase persists with respect to a field applied in the â and ĉ directions.
The AΓ phase is then protected by the product of timereversal (TR) and π rotation around the b-axis (R b ), TR×R b , the only remaining symmetry [60] when the field is in the âĉ plane, but broken when it is along the b-axis.Hence, if a field is applied in the âĉ plane, a transition to the trivial polarized state can only occur at finite field strengths with potentially other phases intervening before the polarized state is encountered.Several such transitions were observed for the AΓ phase (denoted by KΓSL) for a field in the ĉ direction [18,30].On the other hand, the FK phase is not protected by the TR × R b symmetry, and if a field is applied in the ĉ direction the ES degeneracy is lost as shown in Fig. 10(b).However, the field induced FK phase can still be distinguished from the polarized state.

DISCUSSION
Our initial inquiry in this paper pertains to the nature of the AΓ phase and whether there exists a defining quantity for its characterization.For example, the Kitaev phases (AK and FK) in the ladder display the character of SPT phases.It is likely that AΓ is another SPT phase.If so, we expect all the signatures of the SPT such as the degeneracy of the entanglement spectrum, ground state degeneracy under OBC, and the presence of a SOP.Using the iDMRG, DMRG, and ED techniques, we indeed found that the entanglement spectrum is degenerate and there exists four-fold ground state degeneracy under the regular OBC in the AΓ phase.It is interesting to note that the same results were obtained for the Kitaev phases, AK and FK, but under the slanted OBC.
Despite such clear signatures of the SPT, determining the corresponding SOP in the AΓ phase has been challenging.We found that the string order correlation function is related to ordinary local order in a regular correlation function in a model H KQ obtained only after two consecutive unitary transformations.Hence, we term this order as 'twice' hidden.
To understand the symmetry that protects the AΓ phase, we also investigated the effects of the external magnetic field.From the magnetic field response, we noted that the AΓ phase is completely inert to the magnetic field when the field is applied in the â and ĉ directions, which correspond to perpendicular to the zbond and the ladder plane, respectively.Accordingly, the entanglement spectrum degeneracy remains intact when the field is applied in the â and ĉ directions.This is in contrast to the effect of applying the magnetic field along the b direction, i.e., parallel to the z-bond, which immediately lifts the degeneracy of the entanglement spectrum.We note that the product of TR and R b , TR × R b symmetry, is preserved when the magnetic field is applied in the âĉ-plane, which is valid for the generic honeycomb Kitaev model beyond the ladder [60].Thus, we conclude that the AΓ is protected by the TR × R b symmetry.Another intriguing implication from the magnetic field study is that the edge state in the AΓ phase is not a isotropic free spin-1/2 unlike the standard S=1 Haldane SPT.They act like spinless modes under the field in â and ĉ directions.Further studies are needed to fully understand the nature of the zero-energy modes at the boundary of the system with the regular OBC.
In the context of Kitaev materials, let us revisit our motivations for investigating the AΓ phase in the ladder model.As previously mentioned in the introduction, the majority of d 5 Kitaev materials prominently feature FM Kitaev and AFM Γ interactions.However, ongoing debates persist regarding the specific phase that arises in this region.Several numerical studies have suggested the presence of magnetic disorder [12][13][14][15][16][17][18][19], while others have indicated a magnetically ordered phase, such as a zig-zag order [21].Should a zig-zag order indeed be manifest in the 2D limit, we would anticipate observing the same ordering pattern in the 2-leg ladder, as the magnetic unit cell of the zig-zag can be captured in the ladder geometry.This is indeed the case for the Kitaev-Heisenberg ladder model, where the zig-zag, stripy, and FM ordered phases reported in the 2D honeycomb clusters are found in the ladder geometry [23].
Our findings have substantiated the presence of disordered state in the AΓ phase of the 2-leg ladder, categorizing it as a SPT phase characterized by a SOP.The 2D limit can be constructed by stacking the ladders, and one possibility of the resulting 2D phase is a stacked SPTs with edge modes known as a weak-SPT [61].However, the coupling between the ladders may generate a new phase or critical point.It is interesting to note that the evolution from a stacked weak-SPT chains to a gapless critical point supporting edge modes that do not hybridize with bulk modes was reported in the extended anisotropic Kitaev model approaching from the dimer limit [62].Our findings hint at the possibility that as the 2D limit is approached, the AΓ phase may become a 2D spin liquid, denoted as the KΓ spin liquid.However, we cannot rule out a possibility of large unit cell [63] or incommensurate [20] magnetic orders whose magnetic unit cells are beyond the ladder geometry, and a definitive resolution to this question remains a subject for future investigation.

A. Numerical Methods
We use a fully parallelized implementation [64] of the Lanczos algorithm to perform the exact diagonalizations (ED) of ladders with up to N =30 using both open and periodic boundary conditions.

B. Energy susceptibility
To determine the phase diagram we study the susceptibility derived from the ground-state energy per spin, e 0 : At a quantum critical point (QCP) it is known [72] that, for a finite system of size N , the energy susceptibility diverges as Here ν and z are the correlation and dynamical critical exponents and d is the dimension.Hence, χ e only diverges at the phase transition if the critical exponent ν is smaller than 2/(d + z).For the present case we have d=1 and if we assume z=1, then a divergence is observed only if ν < 1.

C. Entanglement Entropy and Spectrum, Schmidt gap
When studying the ladder shown in Fig. 1a it is important to realize that there are different ways of partitioning the system in two partitions of size x and N − x.This is crucial when considering the bipartite von Neumann entanglement entropy, EE, as well as for the entanglement spectrum [73] of central importance for understanding topological properties [34,[73][74][75].Both are obtained from the spectrum of the reduced density matrix, ρ x , of either one of the two partitions.Here we focus on two specific partitions shown in Fig. 1a as the red and blue dashed lines.With the numbering in Fig. 1a, they correspond to either an odd (N/2-1, red) or even (N/2, blue) number of sites in the partitions.We refer to the density matrix derived from the former case with N/2-1 as ρ B and to the latter case with N/2 as ρ A .We mainly focus on the case where the number of sites in the partition is close to the mid-point, either N/2 − 1 (ρ B ) or N/2 (ρ A ). but when considering the bipartite entanglement entropy, we let the number in the partition vary but only consider an even number of sites in the sub system corresponding to moving the blue dashed line in Fig. 1a along the ladder.For a subsystem, A, of size x the entanglement entropy is defined by: EE(x) = −Trρ x ln ρ x . ( Our results for EE(x) can be found in Supplementary Note 2. The eigenvalues, l α , of the reduced density matrix, ρ x , correspond to the Schmidt decomposition, l α =λ 2 α and thereby the entanglement spectrum [73], which then will depend on whether ρ A or ρ B from Fig. 1a is used.The Schmidt gap is then defined as λ 1 − λ 2 .lengths hinder a clear interpretation of the results.Even though we find non-zero values for all s α the trend with system size N is not clear, nor systematic, and we cannot determine if any of them remain non-zero in the thermodynamic limit or attain a finite value.In both cases, we also expect some variation in the values of s α as ϕ is tuned, further complicating the analysis.

11 <
l a t e x i t s h a 1 _ b a s e 6 4 = " 2 9 5 N o 6 e Y 3 5 u n 0 2 1 7 9 d 0 c 5 P 6 p o w 8 u I I 6 3 E I D m s A g h G d 4 h T d n 7 L w 4 7 8 7 H o r X g 5 D P H 8 A f O 5 w / x + Y 0 C < / l a t e x i t >

3 a
< l a t e x i t s h a 1 _ b a s e 6 4 = " 2 5 a V N e e 4 b x 3 o s E D n R v w 7 O n d y s 6 c = " > A A A B 7 3 i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e 5 K f B y D

1 <
5 9 3 5 m L e u O P n M E f y B 8 / k D f d 2 M v w = = < / l a t e x i t > l a t e x i t s h a 1 _ b a s e 6 4 = " 9 o g P 3 O A p f 5 Q + Y C o Q o z Y 2 2 6 a B P W 0 = " > A A A B 8 n i c b V D L S g M x F M 3 U V 6 2 v q k s 3 w S K 4 K j P i a 1 l 1 I 7 i p Y B 8 w H U o m z b S h y W R I 7 g h l 6 G e 4 c a G I W 7 / G n X 9 j p p 2 F V g 8 E D u f c S 8 4 9 Y S K 4 A d f 9 c k p L y y u r a + X 1 y s b m 1 H N U c 6 L 8 + 5 8 L E o L T t 5 z D H / g f P 4 A R y W N s g = = < / l a t e x i t > 0.086< l a t e x i t s h a 1 _ b a s e 6 4 = " J 2 F 3 4 V v 8 y C 8 U u E D N I Z p L z B G 5 2 E I = " > A A A B 8 n i c b V D L S g M x F M 3 U V 6 2 v q k s 3 w S K 4 K j P i a 1 k U R H B T w T 5 g O p R M m m l D k 8 m Q 3 B H K 0 M 9 w 4 0 I R t 3 6 N O / / G T D s L rR 4 I H M 6 5 l 5 x 7 w k R w A 6 7 7 5 Z S W l l d W 1 8 r r l Y 3 N r e 2 d 6 u 5 e 2 6 h U U 9 a i S i j d D Y l h g s e s B R w E 6 y a a E R k K 1 g n H 1 7 n f e W T a c B U / w C R h g S T D m E e c E r C S 3 5 M E R l p m N 3 f T f r X m 1 t 0 Z 8 F / i F a S G C j T 7 1 c / e Q N F U s h i o I M b 4 n p t A k B E N n A o 2 r f R S w x J C x 2 T I f E t j I p k J s l n k K T 6 y y g B H S t s X A 5 6 p P z c y I o 2 Z y N B O 5 h H N o p e L / 3 l + C t F l k P E 4 S Y H F d P 5 R l A o M C u f 3 4 w H X j I K Y W E K o 5 j Y r p i O i C Q X b U s W W 4 C 2 e / J e 0 T + r e e f 3 s / r T W u C r q K K M D d I i O k Y c u U A P d o i Z q I Y o U e k I v 6 N U B 5 9 l 5 c 9 7 n o y W n 2 N l H v + B 8 f A N y U p F g < / l a t e x i t > FK < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 6 e 8 b e m I y G T a b e H O k 8 + K 8 O x / z 1 h U n n z m C P 3 A + f w B 5 L o 3 E < / l a t e x i t > 0.435 < l a t e x i t s h a 1 _ b a s e 6 4 = " v t W c d y 5 J b e j X / D n

4 bFIG. 1 .
FIG. 1.The KΓ ladder and corresponding phase diagram.aA strip of the KΓ honeycomb lattice corresponding to a two-leg KΓ ladder with alternating x and y bonds along the leg and z-bond between the chains, with the numbering of the sites used throughout the paper.The dotted z-bonds arise from imposing periodic boundary conditions in the direction perpendicular to the ladder.However, all z-bonds are taken to have equal strength.The dashed red line indicates the partition used for ρB while the dashed blue line indicates the partition used for ρA.b Schematic phase diagram of the KΓ ladder, Eq. (1), for Γ>0.

FIG. 2 .
FIG.2.Entanglement spectrum and susceptibility for ϕ/π<0.5.iDMRG results for: a Entanglement spectrum from ρA, red dashed line in Fig.1a.b Entanglement spectrum from ρB, blue dashed line in Fig.1a.The numbers refer to the degeneracy of the eigenvalue.c χ e ϕ .Two distinct phases, SPTα and SPT β , are visible between the FMU 6 and AΓ phase.The red shading between ϕ/π=0.428-0.442denotes a transitional region of limited convergence due to a field instability.

FIG. 3 .
FIG.3.Entanglement spectrum and susceptibility for ϕ/π>0.5.iDMRG results for: a Entanglement spectrum from ρA, red dashed line in Fig.1a.b Entanglement spectrum from ρB, blue dashed line in Fig.1a.The numbers refer to the degeneracy of the eigenvalue.c χ e ϕ .The AΓ FK and RSU 6 phases are clearly delineated.

FIG. 4 .
FIG. 4. Unit cells of the U6 transformed ladder.Two unit cells of H U 6 KΓ , the KΓ ladder after the U6 local transformation with N =6n sites.a Two regular unit cells with regular open boundary conditions (OBC).b Two slanted unit cells with slanted OBC.
< l a t e x i t s h a 1 _ b a s e 6 4 = " d + 7 I 8 u e p 9 k y / H P v v g C 4 e Q I y + B B 4 = " > A A A B 6 3 i c b 4 h l d 4 c 8 b O i / P u f M x b C 0 4 + c w h / 4 H z + A O v Y j P k = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " 5 p L D s n I j B 5 u 5 m P Y X 3 h 5 X S D b s h x

6 KFIG. 5 .
FIG. 5. Pictorial view of KG transformations.The two transformations from the original KΓ model to H U 6 KΓ , and the subsequent transformation to HKQ are sketched.The type of order parameter for each Hamiltonian is indicated.

FIG. 6 . 6 |FIG. 7 .
FIG.6.Spin correlations in the HKQ model at ϕ=0.85π.DMRG results with N =384 for the correlation function 4⟨S ′′y (0)S ′′y (r)⟩ versus distance, r along one leg of the ladder.r=0 corresponds to site 47.Results are for HKQ with ϕ=0.85π and r is measured in lattice spacings along the leg.

<
l a t e x i t s h a 1 _ b a s e 6 4 = " A d 9 w c m C w B o 4 S R N n c c x + C P a n w w w I = " > A A B T A X i c 7 V x b b 9 s 2 F H a 7 b u u y W 7 s 9 7 k

FIG. 9 .
FIG. 9. Ground state magnetization in an infinitesimal field in the AΓ phase.ED results with N =24 for ⟨S a,b,c i ⟩ in the lowest state in the AG phase at ϕ/π=0.85 with a field in the b-direction of strength 0.002.The green, blue and cyan colors indicate positive values while orange, red, and pink indicate negative values.The size of the circles are proportional to the value of ⟨S a,b,c i ⟩. a ⟨S a i ⟩, b ⟨S b i ⟩, c ⟨S c i ⟩.