the study of the disruption of the ordering tendency in low - dimensional antiferromagnets at zero temperature and the subsequent emergence of unconventional phases is an active topic of research , on both theoretical or experimental fronts .
frustrating interactions , quantum fluctuations or imposition of a magnetic field are among the main ingredients to destabilize the antiferromagnetic long - range order  @xcite .
considerable efforts have been devoted to studying and characterizing the resulting phases of matter , which may assume another type of non - magnetic ordering ( such as crystalline states of singlet valence bonds ) , or may not exhibit any kind of local order at all ( _ quantum spin liquids _
@xcite ) .
two - dimensional non - bipartite lattice with antiferromagnetic interactions between spin-1/2 moments are amongst the most studied , as antiferromagnetic order is most severely challenged by quantum fluctuations there .
the recent numerical findings of topological spin liquid phases in the realistic su(2 ) kagom - lattice antiferromagnet  @xcite , or its xxz version  @xcite provide examples of stabilization of such exotic phases of matter .
previously , @xmath3 topological spin liquids were found in toy models , such as constrained quantum dimer models ( qdm )  @xcite on the frustrated triangular  @xcite or kagom  @xcite lattices .
these simplified models advantageously allow for a better analytical handle and understanding of the physics of spin liquids , but their connection to su(2 ) microscopic hamiltonians is not always direct .
another route to realizing exotic phases of matter is to start from a model with highly degenerate manifold of ground states ( _ ice manifold _ ) obeying local constraints and then derive an effective hamiltonian describing the emergent degrees of freedom in this manifold @xcite . following this strategy ,
balents , fisher and girvin ( bfg ) introduced a spin-1/2 xxz model with extended interactions on the kagom lattice  @xcite . using a generalized qdm derived for low energies
, these authors argued that this system hosted a topological gapped @xmath3 spin liquid phase  @xcite , which was further confirmed numerically in subsequent works  @xcite .
one of the hallmarks of a @xmath3 liquid phase is the existence of a topological correction to entanglement entropy , which was shown to be present in the bfg model  @xcite .    in considering the nature of the ( gapped ) phases resulting from the destruction of antiferromagnetic order ,
one important guiding principle is provided by a general statement about the relation between the presence / absence of an excitation gap and commensurability .
for instance , it can be rigorously shown @xcite that a given spin system ( in the absence of magnetic field ) is gapless when the value of total spin per unit cell is an half - odd - integer ( this is the case for , _
e.g. _ , the spin-1/2 heisenberg model on the kagom lattice ) and the ground state is non - degenerate .
that is , if we have a featureless liquid - like ground state with an excitation gap , the state necessarily has a _ hidden _ degeneracy ( probably of topological origin ) .
a simpler version of this argument @xcite applies to any spin-@xmath4 systems with u(1 ) symmetry and is of direct relevance to magnetization plateaux ( see @xcite for a review ) : a unique ground state with a gap is possible only when the number of spins @xmath5 within a lattice unit cell and the magnetization @xmath6 per site satisfy @xmath7 .
recently , a field - theoretical meaning was given to this relation  @xcite and it was predicted that when the @xmath8 is a simple ( non - integral ) rational number , we may have gapped featureless liquid phases ( _ i.e. _ , spin - liquid plateaus ) as well as more conventional plateaus with magnetic superstructures .
the search for microscopic models that potentially host these spin - liquid plateaus formed in high magnetic fields is quite intriguing in the light of both the field control of spin liquids @xcite and a recent theoretical report @xcite .    in this paper
, we will investigate the ground - state phases of the bfg model in the presence of a magnetic field , which , as we will see , exhibits large magnetization plateaus at magnetization per site @xmath0 and @xmath1 as well as the one at @xmath9 ( the gapped spin liquid reported in ref .
@xcite ) . on these plateaus ,
the low - energy physics is also well captured by effective constrained models on a triangular lattice : a generalized qdm for the @xmath9 plateau @xcite , the usual qdm for @xmath10 , and a quantum loop model ( qlm ) for @xmath0 .
while much is known already about the nature of the @xmath10 plateau thanks to numerous extensive studies  @xcite done for the qdm on the triangular lattice , the phase diagram of the qlm , which is of direct relevance to the @xmath0 plateau , is not known to our best knowledge .
we will investigate it in details using simple analytical considerations and quantum monte carlo ( qmc ) simulations .
the paper is organized as follows . in sec .
[ sec : model ] , we present the spin 1/2 model considered in this study and provide the low - energy mapping to constrained models valid for the magnetization plateaus at @xmath0 ( qlm ) and @xmath1 ( qdm ) .
we also review some key results obtained previously when no magnetic field is present .
[ sec : loop_model ] is devoted to a complete mapping of the phase diagram of the qlm on the triangular lattice .
for the qlm with only the kinetic term , that is relevant to the original spin model , we find that the ground state of the qlm displays crystalline behavior with a spontaneous breaking of six - fold rotation symmetry , but not of translation symmetry .
we will then turn , in sec .
[ sec : spin_model ] , to large - scale numerical simulations of the original microscopic spin model at @xmath0 using direct qmc simulations .
computation of the diagonal spin structure factor confirms the existence of a stripe ordering of down spins with a moderately large correlation length .
the presence of this crystalline phase is further corroborated using rnyi entanglement entropy .
conclusions are given in sec .
[ sec : conclusion ] .
typical examples of the third - nearest - neighbor interaction between different hexagons [ that are not included in the hamiltonian ] are shown by green dashed lines .
others are obtained by space - group operations . ]    in the search for simple microscopic models ( e.g. with two - spin interactions ) hosting quantum spin liquid phases , bfg introduced a spin-1/2 xxz model on the kagom lattice  @xcite with the aim to reproduce , in a certain regime of parameters , the physics of a ( generalized ) quantum dimer model , known for hosting a @xmath3 topological spin liquid phase . we reproduce here their construction for completeness and obtain the effective hamiltonians for the plateaus .
the hamiltonian is defined on the hexagonal plaquettes of the kagom lattice ( see fig .  [
fig : kagome_lattice ] ) and reads : @xmath11{hex } } } } \left [ j_{z}\sum_{\langle i , j \rangle\in { \raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } s^z_i s^z_j - h\sum_{i \in { \raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } s^z_i \right ] , \\ & h_{xy } = \sum_{{\raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } j_{xy } \sum_{\langle i , j \rangle\in { \raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } ( s^{+}_{i}s^{-}_{j } + s^{-}_{i}s^{+}_{j } ) , \end{split } \label{eq : bfg_model}\ ] ] where @xmath12 is the spin-1/2 operator .
@xmath13 is the diagonal part of the hamiltonian made up of the ising interactions between first , second and third neighbors _ within _ each plaquette ( see fig .  [
fig : kagome_lattice ] for all the corresponding links ) as well as the zeeman term for a magnetic field @xmath14 along the @xmath15-axis .
note that the third - neighbor interactions between sites belonging to different hexagons , shown by the green dashed lines in fig .
[ fig : kagome_lattice ] , are not included ( see sec .  [
sec : mapping - to - qdm ] for the effect of these neglected interactions ) .
the second part @xmath16 contains the spin - flip terms between sites within the same hexagon .
originally , bfg considered all such terms on each plaquette , but , as will be clear below , setting @xmath17 for the second and third neighbors does not alter the physical behavior [ see eq .  ] .
compared to the original hamiltonian defined in ref .  , we added the magnetic field term and the coefficient of the @xmath18 terms is modified by a factor of @xmath19 .
the model can also be reformulated in terms of hardcore bosons , with the boson density playing the role of the magnetization . in this language ,
the plateaus at @xmath0 and @xmath10 that we will investigate correspond to @xmath1 and @xmath20 filling , respectively ( and their respective particle - hole values ) .     and
different values of @xmath21 , in the ground state ( @xmath22 is taken as the unit of energy ) .
inset : magnetization curve for @xmath23 ( ising limit ) .
the three plateaus at magnetizations @xmath24 are clearly visible for the lowest value of the off - diagonal coupling constant @xmath25 , and gradually disappear as @xmath26 gets larger . ]    in the following we will be interested in the strongly anisotropic limit @xmath27 where we will see that the dimer physics emerges . since the diagonal part in can be written as : @xmath28{hex } } } } ( s^z_{{\raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex}}}})^2 - h\sum_{{\raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } s^z_{{\raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } - \frac{3}{4 } n_{{\raisebox{-0.3ex}{\includegraphics[scale=0.1]{hex } } } } j_z , \label{eq : ice_limit}\ ] ] with @xmath29{hex}}}}$ ] the number of hexagons , every state satisfying the constraint @xmath30{hex}}}}=0,1,2,3 $ ] ( depending on the value of @xmath14 ) on all the plaquettes is a ground state of the classical part @xmath13 .
first excited states are separated by a gap of magnitude @xmath18 . as a consequence ,
the magnetization curve @xmath31 ( with normalization @xmath32 ) displays plateaus of width @xmath18 at the values @xmath33 ( see the inset of fig .  [
fig : mofh ] ) and saturation at @xmath19 ( @xmath34 where @xmath35 is the number of down spins per hexagon ) which are expected to survive for a finite value of @xmath21 . as an illustration , we present in fig .
[ fig : mofh ] the magnetization curve in the ground state for a moderate system size @xmath36 and finite small values of @xmath21 , as obtained from qmc simulations ( see sec .  [ sec : spin_model ] for details and a finite - size scaling analysis of the size of the @xmath0 plateau ) .
the above construction leads to magnetization plateaus where the ground states degeneracies scale exponentially with the system size  @xcite .
dynamics induced by the off - diagonal terms will lift this degeneracy , as presented below .
a simple way to visualize the local constraints @xmath30{hex}}}}=0,1,2,3 $ ] ( for all hexagons ) in the ising limit is to faithfully map the ground state configurations to those of dimers : draw a dimer on a link of the triangular lattice formed by the centers of hexagons when the middle of the link is occupied by a down spin , and leave the same link empty for an up spin .
the magnetization plateaus @xmath34 are then reproduced by having a constraint such that every triangular lattice site must be touched by @xmath37 dimers .
this procedure is illustrated in fig .
[ fig : dimer_mapping ] for the @xmath0 ( @xmath38 ) plateau with two dimers per site , which allows us to visualize the ground state spin configurations of the diagonal part @xmath39 as fully - packed self - avoiding loops on the triangular lattice .     and
@xmath0 in terms of dimers ( red ) on the underlying triangular lattice ( blue ) .
loop configurations are formed due to the constraint of two dimers per triangular lattice site . ]
the generalized quantum dimer models appear when treating @xmath16 in degenerate perturbation theory .
clearly , first - order processes do not contribute , since applying the spin - flip term only once produces configurations with two excited hexagons , thus outside of the ground state manifold . in second - order , two types of processes are allowed .
the first process simply flips the same pair of spins twice , merely shifting the total energy by @xmath40{hex}}}}j_{xy}^2/j_z$ ] .
the second process involves a flippable bow - tie plaquette , as pictured in green in fig .
[ fig : dimer_hopping ] .
the two pairs of antiparallel spins are flipped , generating another state fulfilling the ground state constraint .
the second - order effective hamiltonian is thus given by the ring - exchange hamiltonian  @xcite : @xmath41 where the sum runs over all flippable bow - tie plaquettes ( labels @xmath42 and @xmath43 denote the exterior sites of the bow - ties ) and the ring - exchange amplitude is given by @xmath44 .
quite importantly , the sign of @xmath21 is immaterial which allows numerical simulations of the original microscopic hamiltonian ( [ eq : bfg_model ] ) with qmc techniques , by considering the model with a ferromagnetic @xmath45 coupling .
also in order to obtain the ring - exchange hamiltonian , only the xy - interactions on the nearest - neighbor ( @xmath46 ) and the next - nearest - neighbor ( @xmath47 ) bonds are crucial : @xmath48 from this , one can see , as stated previously , that suppressing the further - neighbor xy interactions in @xmath16 [ eq .  ]
does not alter physics qualitatively and simply rescales @xmath49 by a factor 1/2 .
we exploit this property in sec .
[ sec : spin_model ] .     on a flippable bow - tie plaquette ( green ) for @xmath50 .
the antiparallel spins of two flippable pairs are exchanged , corresponding to a flip of two parallel dimers on the underlying triangular lattice ( blue ) . ]    in the dimer language , the ring move on a bow - tie of the original lattice corresponds to the flip of two parallel dimers on a diamond of the dual triangular lattice ( fig .
[ fig : dimer_hopping ] ) .
this is nothing but the off - diagonal term of the quantum dimer model  @xcite on a triangular lattice  @xcite , which , as usual , is conveniently extended with a diagonal potential ( @xmath51 ) term which counts the number of flippable plaquettes ( diamonds ) : @xmath52{plaq } } } \left ( | \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer1 } } \rangle   \langle \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer2 } } | + | \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer2 } } \rangle   \langle \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer1 } } | \right ) \\ & + v \sum_{{\includegraphics[scale=0.25]{plaq } } } \left ( | \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer1 } } \rangle   \langle \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer1 } } | + | \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer2 } } \rangle   \langle \raisebox{-0.65ex}{\includegraphics[scale=0.45]{dimer2 } } | \right ) .
\end{split } \label{eq : qdm_triangular}\ ] ] in the ising limit @xmath53 , low - energy physics of the bfg model is expected to be described by this hamiltonian for @xmath54 .
the following four - spin interaction @xcite @xmath55 ( @xmath56 label the four spins on a bow - tie in a anti - clockwise way ) reduces , in the same limit , to the @xmath51-term ( @xmath57 ) as it counts the number of flippable plaquettes . for the @xmath10 plateau corresponding to one dimer ( or , one @xmath58 spin ) per hexagon , one may use instead the _ inter - hexagon _ ( two - spin ) third - neighbor interaction @xmath59 that is not contained in model to realize the dimer interaction @xmath60 .
thus , one sees that _ the effective model is identical for the three magnetization plateaux , the only difference coming from the constraint on the hilbert space_. for the constraint of one dimer per site ( corresponding to the plateau @xmath10 ) , we have the usual qdm  @xcite , while in the case with two dimers per site ( the plateau at @xmath0 ) , a qlm is obtained as announced in the introduction .
the original model derived by bfg @xcite describing the @xmath9 plateau corresponds to three dimers per site .
table  [ tab : effective - model - summary ] summarizes the nature of the effective models for the different plateaus .
the mapping to dimer models [ including the diagonal @xmath51-term in eq .  ] allows us to benefit from the accumulated knowledge on this family of models  @xcite .
qdms generically admit a rokhsar - kivelson  @xcite ( rk ) point @xmath61 , where the ground - state wave - function is exactly known . for the qdm on the triangular lattice ,
the ground state is exactly shown to be a topological gapped @xmath3 liquid phase  @xcite , which furthermore extends in the region @xmath62 .        in the three - dimer model ,
bfg confirmed the presence of a topological phase at the rk point  @xcite , showing in particular that visons are gapped , implying the presence of deconfined fractionalized spinons  @xcite .
the authors of ref .   also speculated that the @xmath3 liquid phase could extend beyond the @xmath54 point corresponding to the original microscopic spin model , in contrast to the one - dimer case  @xcite for which the liquid phase ends at @xmath63  @xcite .
an exact diagonalization study  @xcite of the hamiltonian ( [ eq : ring_hamiltonian ] ) shows the presence of a vison gap in a regime including the bfg point ( @xmath54 ) , thereby supporting the above suggestion .
however , the presence of the spin liquid was only firmly established using qmc simulations on large systems .
first , a complete numerical phase diagram of the model ( [ eq : bfg_model ] ) with ferromagnetic @xmath21 and @xmath64 was obtained  @xcite : at strongly negative @xmath65 , a planar ( _ i.e. _ , superfluid in the bosonic language ) phase accompanied by the breaking of the u(1 ) symmetry is found , as expected . when the magnitude of the spin - flip term was increased , a continuous phase transition to an apparently featureless insulating phase was observed .
the phase diagram was also extended to finite temperature  @xcite .
finally , a direct evidence of the topological nature of the insulator is obtained by computing the topological entanglement ( rnyi ) entropy  @xcite , which was evaluated as @xmath66 as expected for a @xmath3 spin liquid phase . as a consequence of the condensation of fractionalized excitations ,
the transition between the planar phase and the featureless insulator is an exotic 3d @xmath67 quantum critical point  @xcite .
while the nature of the phases encountered at the @xmath9 and @xmath10 plateaus is well understood from this previous set of study , the ground state physics on the @xmath0 plateau has never been investigated to our best knowledge . in the following ,
we focus on this situation by performing numerical simulations first on the effective qlm ( sec .
[ sec : loop_model ] ) , and then on the microscopic spin model ( sec .
[ sec : spin_model ] ) .
.[tab : effective - model - summary ] three plateaus of @xmath68 bfg model and the corresponding effective models . [ cols="<,^,^,^",options="header " , ]
we consider here the effective model ( [ eq : qdm_triangular ] ) in the case of two dimers per site , where allowed states are represented by configurations of self - avoiding loops such as represented in fig .
[ fig : dimer_mapping ] . motivated by the original spin model eq .   at the @xmath0 plateau , we are primarily interested in the nature of the ground state of the effective qlm at @xmath69 , however we will also investigate the nature of the surrounding phases in the phase diagram .
we first give , in sec .
[ sec : loop_general ] , general arguments on the structure of the phase diagram as well as on the topological properties of the loop configurations , before complementing this analysis in sec .
[ sec : loop_qmc ] with qmc simulations of the loop model . anticipating the results obtained , we present in fig .  [
fig : loop_phase_diag ] the ground - state phase diagram of the qlm on the triangular lattice to illustrate the following discussion .      when @xmath70 , the ground - state energy is minimized by configurations where loops form straight lines along one of the three lattice directions , as shown in the left part of the phase diagram presented in fig .
[ fig : loop_phase_diag ] .
as aligned up and down spins form alternating stripes in the corresponding configurations of the original spin models ( see the right panel of fig .
[ fig : qmccorr ] ) , we call this a _ stripe phase_. in this gapped ordered phase , some of the rotations and reflections of the triangular lattice are spontaneously broken ( _ i.e. _ , the point group changes from @xmath71 to @xmath72 ) , leading to the three - fold degenerate ground states ( in the qdm , the corresponding columnar phase further breaks translation symmetry ) .    at the rk point
@xmath73 , the ground state is given by the equal - weighted sum of all fully - packed loop configurations on the triangular lattice .
to our best knowledge , the equivalent classical problem has not been studied before . as the triangular lattice is non - bipartite
, we expect that the loop segment correlations are short - ranged , indicating a _ liquid _ phase presumably with a gap , as also found for the rk point of the qdm  @xcite .     and @xmath74 .
the high symmetry points @xmath75 , @xmath76 and @xmath77 required for the non - flippable states are represented . ]
when @xmath78 , the ground states are readily found to have zero energy and correspond to non - flippable loop configurations ( _ i.e. _ for which the kinetic term vanishes ) , which is again similar to the staggered phase of the qdm . while the construction of those states is quite straightforward in the qdm , we have found that for @xmath79 systems , the qlm accepts two different families of non - flippable configurations , both of which are shown in the right part of the phase diagram of fig .
[ fig : loop_phase_diag ] .
the first family is obtained by starting from the staggered states of the qdm @xcite , which are 12-fold degenerate .
on the left , one of those 12 states is represented in red , with only two of the three types of links ( corresponding to the three lattice directions 0 , 1 , and 2 in fig .
[ fig : brillouin_zone ] ) occupied .
a non - flippable loop configuration can then be generated by adding parallel ( blue ) dimers on half of the remaining type of links , for which there are two possibilities ( a one - step translation in the direction of these parallel dimers also creates a different non - flippable state ) .
the same prescription can be applied starting from the other qdm staggered states , leading to a total degeneracy of @xmath80 .
the second family is obtained by tiling the lattice with the shortest possible triangular loops , as depicted on the right part of fig .
[ fig : loop_phase_diag ] .
the degeneracy of this family is readily found to be 6 .
considering the brillouin zone drawn in the right panel of fig .
[ fig : brillouin_zone ] , the first family of non - flippable states requires the cluster to have both @xmath81 and @xmath82 points ( i.e. @xmath83 multiple of 4 ) , while the second one needs the @xmath84 point ( i.e. @xmath83 multiple of 3 ) .
we finally remark that the nature of non - flippable states ( or absence thereof ) can be different for other types of clusters ( rectangular clusters for instance ) .    from these considerations
, the point @xmath54 lies between the two limiting cases @xmath85 ( stripe ) and @xmath86 ( liquid ) .
as both the stripe and liquid phases are gapped , they should extend in a finite region of phase space close to these limits , and it is not clear _ a priori _ in which phase the point @xmath54 is located . of course ,
it is also possible that one or several intervening crystalline phases arise : natural candidates are equivalent of plaquette , or @xmath87 phases observed in qdms  @xcite .
this can only be determined by exact numerical calculations , which will be presented below in sec .
[ sec : loop_qmc ] .
invariants @xmath88 associated to the parity of the number of loop segments crossed by the line . in this example , the number of loop segments pointing in the vertical @xmath89 ( horizontal @xmath90 ) direction which cross the line is odd ( even ) , corresponding to @xmath91 ( @xmath92 ) . this number does not depend on the precise location of the line ( parallel lines shifted by one or several lattice units result in the same parity ) .
the flipping terms of the hamiltonian , acting for instance on the shaded plaquette ( other flips are possible in this configuration ) , do not change these two parities . note that a loop segment located on a diagonal bond of a plaquette contributes to both @xmath93 and @xmath94 . ]
it is interesting to note at this stage that on a manifold with non - zero genus ( cylinder , torus , etc ) , the qlm model possesses topological sectors which can be defined in the same way as in the qdm case . as illustrated in fig .
[ fig : loop_winding ] for a torus , the parity of the number of loop segments crossing a dashed line in either the @xmath90 or @xmath89 directions of the torus is conserved by the dynamics of the hamiltonian . on the torus ,
this defines a pair of @xmath3 invariants @xmath95 and @xmath96 , and thus divides the full hilbert space into four distinct topological sectors , which we label by @xmath88 .
the existence of these four topological sectors is the key to identify the @xmath3 spin liquid phase close to @xmath73 , as _
e.g. _ in the qdm on the triangular lattice  @xcite .
one may wonder if the @xmath3 spin liquid found here is the same as that in the usual triangular - lattice qdm since even ( odd ) number of dimers exiting from each site in the qlm ( qdm ) suggests a mapping to an even ( odd ) ising gauge theory  @xcite . however , this difference is not important in considering the underlying nature of the spin liquids .
in fact , by explicitly constructing the ground state of the toric code  @xcite on a triangular lattice  , one can show that both spin liquids share the same quasiparticle ( anyon ) contents and exhibit the same ground - state degeneracy originating from the @xmath2 flux , which we will use in the following .    in the topological @xmath2 spin liquid phase
, the four topological sectors ( on a torus ) corresponds to the four different ways of threading the @xmath2 flux through the two periods of the system that can not be eliminated by finite - order perturbations @xcite .
therefore , the four sectors , in the topological spin liquid phase , are separated from each other by exponentially small gaps ( _ topological gaps _ ) and get degenerate in the thermodynamic limit . as was previously observed for the qdm  @xcite
, we expect the liquid phase to be detected by the closing of the above topological gaps .    on the other hand , in the ordered stripe phase , topological sectors which do not contain the pure stripe configurations ( shown in the left part of fig .
[ fig : loop_phase_diag ] ) can not optimize the attractive @xmath51-term and should have a higher energy : correspondingly , the topological gap should be finite even in the thermodynamic limit , providing us with a practical method to distinguish the two phases .
one should however be careful in identifying the sectors to which the stripe configurations belong .
indeed on a torus with @xmath97 , the qlm can be studied for even and odd @xmath83 ( in contrast to the qdm that can only host even-@xmath83 samples ) .
when @xmath83 is odd , the three stripe configurations shown in fig .  [ fig : loop_phase_diag ] are located in the three different sectors @xmath98 , @xmath99 and @xmath100 which are strictly identical ( the number of configurations is the same and the corresponding blocks of the hamiltonian are identical ) : we define the topological gap for odd-@xmath83 samples as @xmath101 , where @xmath102 is the ground state energy in the @xmath88 sector . for even @xmath83 ,
the three stripe configurations belong to the @xmath103 sector , and all other sectors are identical , and we use @xmath104 .
we now turn to qmc simulations of the qlm on a torus ( for clusters with the geometry represented in fig .
[ fig : brillouin_zone ] , with @xmath97 and @xmath105 sites ) , as defined in eq .  .
we use the reptation qmc  @xcite , which projects a given initial configuration belonging to each topological sector onto the ground state in the same sector .
we use a formulation of reptation qmc similar to the one presented in ref .  : to accelerate the convergence , we work in a continuous - time formalism as well as with a guiding wave function that provides a finite fugacity for each flippable plaquette of the triangular lattice ( an optimized value for the fugacity is obtained from a short preliminary run ) .
we simulated even-@xmath83 samples up to @xmath106 ( and one data point at @xmath107 ) , and odd-@xmath83 samples up to @xmath108 .
as a function of @xmath109 , for different sample sizes @xmath83 .
the finite - size effects are different for even ( @xmath110 ; left panel ) and odd ( @xmath111 ; right panel ) values of @xmath83 . in the thermodynamic limit
, we expect that in both cases , @xmath112 if finite for the stripe phase @xmath113 , and vanishes for the @xmath2 spin liquid phase @xmath114 .
left inset : finite - size scaling as a function of @xmath115 for the even-@xmath83 samples provides an estimate @xmath116 of the value where the gap first vanishes .
right inset : zoom close to the rk point for odd-@xmath83 samples . ]
we first present the results for the topological gap @xmath112 , as obtained by computing the ground state energy in each non - equivalent sector using the standard mixed estimator  @xcite .
[ fig : loop_gap ] shows @xmath112 as a function of @xmath109 for even ( left panel ) and odd ( right panel ) @xmath83 , for different system sizes . for even samples ,
the topological gap @xmath110 has a monotonous behavior : it is clearly finite for low values of @xmath109 and vanishes with system size close enough to the rk point @xmath73 .
finite - size scaling as a function of @xmath115 ( left inset of fig .
[ fig : loop_gap ] ) indicates that the gap vanishes around @xmath116 , separating the stripe phase for @xmath117 from the topological @xmath2 liquid phase for @xmath118 . for the largest samples ,
our data at @xmath119 and @xmath120 show that the ground state energies in different topological sectors are identical within error bars ( we can not resolve the exponentially small splitting in the liquid phase ) .
the gap for the odd-@xmath83 samples displays an apparently different behavior , with the topological sector @xmath103 having a higher energy for low enough @xmath109 , resulting in @xmath121 as the @xmath51-term favors the sectors @xmath98 , @xmath99 and @xmath100 . on the other hand , when the system is close enough to the rk point , the sector @xmath103 has lower energy ( resulting in @xmath122 in our definition of @xmath111 ) for finite - size samples .
if we assume that the ordering of the four levels in the spin liquid phase is the same as that for the even-@xmath83 samples , the sign change marks the transition .
the value of @xmath109 at which @xmath123 varies with system size and appears to extrapolate to a consistent value for the critical point @xmath116 ( see right inset of fig .  [
fig : loop_gap ] ) . here
again , for the largest odd-@xmath83 samples , the energies in all topological sectors are indistinguishable within error bars for @xmath119 and @xmath120 , indicating the topological degeneracy .
the topological gap for both odd and even sample is clearly finite for @xmath69 ( with a value @xmath124 , as estimated from the even-@xmath83 data , for which there is almost no size dependence ) , which indicates that this point is located in an ordered phase ( presumably , in the stripe phase ) . in order to confirm this
, we computed the loop segment occupations ( using middle - slices in the reptile representation  @xcite ) in the topological sector @xmath103 ( which hosts all three stripe configurations ) for the even-@xmath83 samples .
this allows to determine the stripe structure factor , which we define as : @xmath125 where @xmath126 if there is a loop segment at site @xmath127 occupying a link in the direction @xmath128(@xmath129 ; the three values taken by @xmath128 correspond to the three different directions as represented in the left bottom part of fig .  [
fig : brillouin_zone ] ) and @xmath130 otherwise . in the second line of eq .
, we have introduced the notation @xmath131 for the total number of bonds in the direction @xmath128 . the structure factor divided by the sample size
is expected to converge to the square of the stripe order parameter @xmath132 in the thermodynamic limit : @xmath133 , with a non - zero @xmath134 if the ground state exhibits the long - range stripe order .     as a function of inverse system size , for different values of @xmath109 . ]
the stripe structure factor @xmath135 is shown in fig .  [
fig : loop_sf ] for different values of the coupling @xmath109 as a function of inverse size @xmath136 ( for @xmath137 ) .
we expect that @xmath135 should saturate in an exponential way ( due to the finite correlation length ) to @xmath138 in the stripe phase , and to @xmath139 in a phase with no stripe order ( at the critical point , it should decay to @xmath139 with a non - trivial power ) . due to the moderate system sizes accessible , we can not expect to see the exponential decay : the extrapolation to the thermodynamic limit is indeed delicate for large values of @xmath109 .
considering the data for @xmath140 ( see fig .
[ fig : loop_sf ] ) , we observe that the structure factor extrapolates to a non - zero value . from this data
set , we obtained an estimate for the critical point @xmath141 , which is less precise , though consistent , than the one determined above by the topological gap .     of the quantifier of @xmath142 anisotropy of histograms of loop occupations , as a function of @xmath109 for three even lattices sizes @xmath137 .
insets : color - coded normalized histograms of loop occupations numbers ( using a color interpolation scheme ) , deep in the stripe phase ( left inset , @xmath36 , @xmath143 ) and closer to the critical point ( right inset , @xmath36 , @xmath144 ) .
vertical axis of histograms is @xmath145 ( ranging from @xmath146 to @xmath147 ) , horizontal axis is @xmath148 ( from @xmath149 to @xmath150 ) .
the histograms are enclosed in an equilateral triangle pictured close to the right inset . ]
note that , in analogy with the situation in the qdm , one may suspect that other types of symmetry - breaking order ( such as plaquette order ) could also exhibit a non - zero value of @xmath132 . to clarify this , it is useful to consider the histograms of the occupations of loop segments appearing in the monte carlo simulations , as was done in similar simulations for valence bond crystals  @xcite . to have a two - dimensional representation of the histograms for the loop segments occupying one of the three lattice directions , we introduce @xmath151 and @xmath152 ( with @xmath153 being the total number of bonds in the direction @xmath128 ) .
the three line configurations shown in the left part of fig .  [ fig : loop_phase_diag ] ( having respectively @xmath154 or @xmath155 equal to @xmath156 ) occupy , in the @xmath157 representation , the three vertices @xmath158 of an equilateral triangle in which the histogram is enclosed , as represented in fig .  [
fig : loop_cos ] .
the left inset of fig .
[ fig : loop_cos ] presents the histogram obtained for a moderate system size @xmath36 and @xmath143 , which is already located deep inside the stripe phase : we clearly observe a @xmath159-symmetric feature with three peaks ( corresponding to high occupations ) close to the corners of the triangle corresponding to the three line configurations .
the distance of these peaks from the origin @xmath160 gives a measure of the order parameter @xmath132 .    on the other hand ,
the _ angles _
@xmath161 associated to the three peaks in this two - dimensional representation also carry important information on the precise nature of the phase .
for instance , the angles associated with the three stripe patterns are @xmath139 , @xmath162 and @xmath163 .
other candidate phases ( such as the equivalent of plaquette phases ) would show up in the histograms as sharp peaks at different angles .
we find that , as @xmath109 gets closer to the critical point @xmath164 , the typical extents ( radius ) of the histograms get smaller ( as expected from the decrease of the order parameter ) , and that at the same time the shape of the histograms change from a @xmath159-symmetric one to a circular - symmetric one as is exemplified in the right inset of fig .
[ fig : loop_cos ] .
this is reminiscent of the @xmath165 valence bond histograms observed in the vicinity of deconfined quantum critical points  @xcite or even at moderate proximity of the rk point of the qdm on the square lattice  @xcite .    to quantify this effect , we calculated a measure of the @xmath159-symmetry of the histograms @xmath166 using formulations similar to those in refs .  @xcite . for a pure stripe phase , we have @xmath167 , while for a circular histogram @xmath168 .
the variation of @xmath169 as a function @xmath109 with different system sizes displayed in the main panel of fig .
[ fig : loop_cos ] confirms the behavior observed visually in the histograms : for the region @xmath170 , @xmath169 tends to approach unity with increasing system sizes , while in the region @xmath171 , @xmath169 appears to vanish within the system sizes available to us .
while therefore we can not exclude an intermediate different crystalline phase , it appears unlikely as there is no direct evidence either in the energy or in the structure factor of a phase transition ( see a similar situation for the square lattice qdm  @xcite ) . thus , we expect that on larger systems @xmath169 would also tend to approach unity for @xmath172
. note that the computations of @xmath169 are difficult as we need long monte carlo runs to be fully ergodic ( _ i.e. , _ in order not to get locked in one type of line configurations , in particular in the stripe phase ) , and also to ensure that statistical fluctuations are not too important when the order parameter is small .
we may speculate that these circular histograms are indications of the @xmath159 anisotropy being a dangerously irrelevant variable at the critical point @xmath164 .
although the precise characterization of the transition between the @xmath3 spin liquid and the stripe phase is an interesting issue , it is beyond the scope of the present work . at the point @xmath69 of our main interest , our qmc data ( see figs .  [
fig : loop_sf ] and  [ fig : loop_cos ] ) clearly indicate that the ground state of the effective hamiltonian is located in the stripe phase .
the resulting phase diagram obtained from these qmc results as well as the general considerations in sec .
[ sec : loop_general ] is given in fig  [ fig : loop_phase_diag ] .
in this section , we present the results of the qmc simulations for the original spin-1/2 hamiltonian ( [ eq : bfg_model ] ) at the @xmath0 plateau .
we used the stochastic series expansion ( sse ) algorithm  @xcite combined with a plaquette generalization  @xcite that is necessary to circumvent the freezing encountered in the bond formulation of the sse when @xmath173 is large  @xcite .
in addition to this difficulty , capturing ground state physics requires to reach a very low temperature - range characterized by the small energy scale @xmath174 of the effective model .
nevertheless , we have been able to simulate @xmath175 systems , with @xmath176 sites , up to @xmath177 , or slightly beyond @xmath178 for some observables . for numerical ease , we restrict ourselves to the truncated model where only the nearest - neighbor interactions ( i.e. , the black bonds in fig .
[ fig : kagome_lattice ] ) are retained in @xmath16 ( see the discussions in sec .  [
sec : mapping - to - qdm ] ) .    in the following
, we set @xmath179 as the energy scale , and use ferromagnetic values @xmath180 ( for the qmc simulations to have no sign problem ) and @xmath181 for the magnetic field ( corresponding to the middle of the @xmath0 plateau ; see fig .
[ fig : mofh ] ) . in order to locate the transition
point out of the plateau phase as a function of @xmath21 , we performed a finite - size analysis of its width .
resulting data are shown in fig .
[ fig : plateau_size ] , and the transition point is estimated to be around @xmath182 . for values @xmath183 ,
the system is thus in a planar phase ( superfluid in the bosonic language ) , and in the following we focus on the plateau region @xmath184 . in sec  [ sec : plateau_nature ] , we examine the nature of the plateau phase , and discuss what could be the scenarios for the transition to the planar phase .
sec  [ sec : plateau_renyi ] is devoted to an alternative characterization of the phase using the rnyi entanglement entropy .
plateau width .
the plateau is found to disappear around @xmath182 . ]      for magnetization plateaus , as has been mentioned in the introduction , it is useful to consider the quantity @xmath8 , where @xmath5 is the number of sites per unit cell ( @xmath185 for the kagom lattice ) .
indeed , it has been shown analytically  @xcite that when this quantity is fractional , the ground state on this plateau must be degenerate , _
i.e. _ , either in a crystalline state with magnetic superstructures or in a topological state , while a non - degenerate featureless state with a gap is forbidden .
if , on the other hand , it takes an integer value , which is the case here at the @xmath0 plateau ( see table  [ tab : effective - model - summary ] ) , all these possibilities are allowed .        guided by the results of the effective qlm in sec  [ sec : loop_model ]
, we expect the plateau phase to be in a ordered phase where the down spins form stripes .
this is indeed evident when plotting the connected @xmath186-@xmath186 spin correlation function in real space , as shown in the left panel of fig .
[ fig : qmccorr ] , where one readily sees the presence of stripes .
the corresponding phase breaks the six - fold rotation symmetry , and is three - fold degenerate .
cartoon pictures of the spin and dimer configurations in this phase are drawn on the right panel of fig .
[ fig : qmccorr ] .    in order to confirm the long - range nature of this order
, we compute the fourier transform of the above spin - spin correlations , _
i.e. _ , the diagonal spin structure factor : @xmath187 where @xmath188 is the position of spin @xmath189 , @xmath176 , @xmath190 and the average local magnetization is @xmath191 . for an ordered phase at wave - vector @xmath192 ,
@xmath193 diverges as @xmath156 , and for the three - fold degenerate solid , the three bragg peaks are located at wave - vectors @xmath194 , @xmath195 and @xmath196 .
note that this solid can be accommodated on any cluster used for the simulations .
finite - size scaling of @xmath197 as a function of @xmath115 for several values of @xmath21 in the plateau phase clearly indicates the presence of the stripe phase ( see fig .
[ fig : scalingsq ] ) . in agreement with the qlm study at @xmath69 , the order develops already for small system sizes .     inside the magnetization plateau . despite the relatively small available sizes
, we can reach the regime of saturation , confirming the striped solid anticipated from the qlm . ]
two comments are in order here .
first , @xmath193 used in fig .
[ fig : scalingsq ] is the average of the values at the three wave - vectors @xmath198 , which is valid when the simulation is perfectly ergodic .
however , as mentioned earlier , the large anisotropy @xmath199 in the parameter range of interest can create difficulties in the qmc simulations  @xcite . by monitoring the three values @xmath200 separately
, we have checked that the simulations are ergodic for the sizes @xmath201 considered here . however , for smaller values of @xmath202 or larger @xmath83 , we indeed noticed a quick loss of ergodicity , even with the plaquette algorithm , translating into very different values at the three wave - vectors .
the same behavior was observed in ref .   which used a similar qmc algorithm .
second , we sometimes found a non - monotonic behavior of the structure factor @xmath197 for larger sizes ( see fig .
[ fig : scalingsq ] ) .
again , a similar behavior was reported in ref .
where the increase of the order parameter was interpreted as the threshold where stripes start developing in the system . in practice , this complicates a proper extrapolation of the order parameter to its thermodynamic limit value .
we can however compare the order of magnitude obtained with the value for a perfectly equal superposition of the three ordered states , for which a straightforward calculation yields @xmath203 .
we conclude this section by a discussion on the possible scenarios for the quantum phase transition between the plateau ( stripe ) crystal and planar phases . since the stripe crystal breaks discrete lattice symmetries , one expects , on general grounds , that the transition will be of first order type , or that both phases coexist ( in a small region of the phase diagram ) .
this latter scenario implies the existence of an intervening supersolid phase , whose presence was found in several similar models on other frustrated lattices  @xcite . due to the limited systems
sizes available to us , we have not been able to favor one of these two scenarios .
we observed ( data not shown ) an apparent crossing in the spin stiffness , with critical exponents compatible with a continuous phase transition , and no sign of a double peak structure in the kinetic energy histograms .
however , it has been shown in similar situations that the correct nature of the phase transition may be difficult to capture . in the ( @xmath21,@xmath14 ) plane
, varying the xy interaction @xmath21 at the constant field @xmath181 corresponds to a transition point located close to the tip of plateau phase ( _ i.e. _ the tip of the insulator lobe in the bosonic picture ) , for which previous studies have revealed that such a crystal - planar phase transition may appear to be continuous while it is in fact _ weakly first - order
_  @xcite .
we believe that such a picture should also apply here , although proving this definitively would imply simulations on very large samples , which is currently out of reach for the present model .
we have shown , using a conventional approach ( _ i.e. _ by guessing a broken symmetry and then measuring the corresponding structure factor ) , that the @xmath0 plateau corresponds to a stripe crystal shown in fig .
[ fig : qmccorr ] .
it is interesting to check whether other means could detect the nature of the ground state without a priori knowledge on the broken symmetry .
an elegant approach to obtain the dominant order parameter in an unbiased way is provided by the correlation density matrix  @xcite , but it is not easily accessible within qmc simulations . instead , we will focus here on using the entanglement entropy as a means to access the nature of the ground state .    indeed , in the recent years
, a large number of studies have shown how the scaling behavior of a block entanglement entropy can give access to some ground state properties , such as the central charge in one - dimension  @xcite or the topological order underlying a given ground state wave function  @xcite .
generically , any entanglement entropy ( whether rnyi or von neumann ) will scale with the `` area '' of the boundary between the block and the rest of the system ( the so - called `` area - law '' ; see ref .   for a review ) , possibly with interesting ( universal ) sub - leading terms : @xmath204 where @xmath5 is the rnyi index , @xmath205 the length of the boundary of block a and @xmath206 a constant for instance .
topological phases can be detected using the constructions of levin - wen  @xcite or kitaev - preskill  @xcite ( see also ref .  ) , that enables one to extract a negative constant term @xmath207 , independent of @xmath5  @xcite , where @xmath208 is the topological entanglement entropy .
this was used for instance to conclude that the ground state of the bfg model at @xmath9 is indeed a @xmath3 spin liquid with @xmath209  @xcite .
the constructions of refs .
@xcite use a subtraction of terms which all scale like the area of the subsystems considered , which can be troublesome in qmc simulations which inevitably exhibit finite error bars .
computing the scaling of @xmath210 for a single block a ( without using subtraction schemes ) , is in general computationally simpler .
this is what we performed for the case of our discrete symmetry breaking ( @xmath211 ) on the @xmath0 plateau , taking for a the geometry of a cylinder block embedded in a torus sample ( see the inset of fig .  [
fig : renyi_vs_l ] ) .
we expect that for such a symmetry - broken state , the sub - leading term will be a _ positive _ constant @xmath212 , where @xmath213 is the number of degenerate ground states ( note the sign of @xmath214 ; see appendix  [ sec : appendixa ] ) .
we first remark that , for the simple type of order on the @xmath0 plateau , the precise block geometry is not very important , which allows us to consider the convenient cylinder block geometry .
second , it is quite important for this result to hold that all @xmath213 degenerate states ( quasi degenerate on a finite system but with an exponentially small splitting ) are observed in the qmc simulations .
hence , ergodicity must be ensured to extract the degeneracy , which we observe is the case if we are not too deep inside the stripe phase and system sizes are not too large .
note that such a _ positive _ correction gives a vanishing contribution to the previously mentioned geometric constructions  @xcite , as they are precisely built to solely capture long - range entanglement .
( @xmath215 ) on the @xmath0 plateau , as a function of @xmath83 , on various systems with @xmath176 sites .
simulations parameters are @xmath216 , @xmath217 and @xmath181 . in each case
, the a block corresponds to a half - system cylinder containing @xmath218 unit cells ( see inset ) .
the positive correction @xmath219 is shown by a star . ]
our numerical data are presented in fig .
[ fig : renyi_vs_l ] .
we have computed the rnyi entropies @xmath220 and @xmath221 on the @xmath0 plateau at low temperatures using the qmc extended ensemble method proposed in ref .
( with a small modification where the size of a is dynamically varied during the simulation ) .
an important first observation is that if one considers only the smallest ( four or five ) system sizes , a naive fit leads to a negative intercept and thus to an apparent topological phase .
one could perhaps interpret this unexpected finite - size effect in the following way : close enough to the transition to the superfluid phase , the superfluid correlation length can be of the order of , or larger than , the system size for small enough samples .
effectively , the system acquires a superfluid component , which we expect to contribute to @xmath222 as @xmath223 .
this form is derived in ref .   for a system that spontaneously breaks a continuous symmetry , in which case
the logarithm prefactor @xmath224 is universal ( positive and related to the number of goldstone modes ) , and @xmath225 a non - universal constant .
assuming such a form is also valid in a system which has effectively both superfluid and crystalline correlations ( _ i.e. _ , supersolid ) , it could result in fits that display a non - universal , possibly negative ( for small enough systems , since otherwise @xmath226 for large enough @xmath83 ) , constant intercept .
note that it would be interesting to check this behaviour using simpler bosonic models that exhibit supersolid behaviour .
however , when the system is large enough ( beyond @xmath227 in our simulations ) , we observe a qualitative change of behavior in @xmath222 as a function of @xmath83 .
note that this length scale coincides with the one where the stripe structure factor starts saturating , as observed in fig .
[ fig : scalingsq ] . for these larger systems , entanglement entropies
are compatible with a positive intercept @xmath228 , although the accuracy is not excellent .
for instance , fitting our @xmath220 data using ( [ eq : sq ] ) and sizes with @xmath229 , we get @xmath230 , which deviates quite significantly from @xmath231 .
indeed , as already advocated , a precise determination of @xmath206 is rather difficult due to the need of excellent ergodicity and hence large autocorrelation times for this quantity .
typically , we need more than @xmath232 measurements to get reliable data .
in search for magnetization plateaus of spin - liquid nature , we have considered the effect of a magnetic field on an anisotropic spin-1/2 model on the kagom lattice that was introduced in ref .  .
in the easy - axis limit , this model is designed to exhibit several magnetization plateaus , at magnetization per site @xmath233 ( plus saturation @xmath19 ) .
the properties of these plateaus have been investigated first by mapping the original spin model to simpler constrained models on the triangular lattice .
focusing on the @xmath0 plateau , we have investigated in detail the effective quantum loop model and mapped out its phase diagram using a reptation quantum monte - carlo algorithm . from this , we predict a 3-fold degenerate stripe phase , that breaks rotation symmetry , to appear on this magnetization plateau .
we have also observed that the @xmath0 plateau with the stripe order turns into a topological @xmath2 spin - liquid plateau when an additional repulsion among the dimer segments is added .
then , we have performed large - scale numerical simulations using the sse quantum monte - carlo algorithm on the original microscopic model .
this has confirmed the existence of a stripe phase found in the effective loop model , using both standard structure - factor measurements and its signature in entanglement entropies . in particular , we have observed that , for length scales shorter than @xmath234 , the scaling of the block entanglement entropy may be misleadingly interpreted as giving a _ negative _ intercept @xmath206 , while it becomes clearly _ positive _ for larger systems .
this should be used as a caveat in other situations .
when the xy - interaction @xmath26 increases , there is a quantum phase transition to a planar phase ( superfluid of the equivalent hardcore bosons ) .
while it appears to be continuous in our simulations , we believe that it should appear weakly first - order on larger lattices , as observed for instance in similar physical situations  @xcite .
we have not investigated in details the @xmath10 plateau , where the effective model is given by the standard quantum dimer model ( _ i.e. _ a single dimer per site ) on the triangular lattice , that has been widely studied in the past  @xcite .
in fact , the expected crystal in this phase would have a @xmath87 unit cell and an extremely small order parameter  @xcite .
given that numerical simulations of the microscopic spin model are more involved ( larger systems , many different energy scales etc . ) , it remains challenging to detect such a weak order in the direct simulations for the original bfg spin model .    before concluding ,
let us mention possible extensions of the model ( [ eq : bfg_model ] ) . in sec .
[ sec : loop_model ] , we have seen that there exists a gapped @xmath2 spin liquid phase around the rk point of the qlm ( see fig .
[ fig : loop_winding ] ) .
therefore , it would be interesting to drive the system toward the spin - liquid phase .
in fact , it is possible to consider additional interactions which would result in a non - zero interaction term @xmath51 in the effective constrained model ( [ eq : qdm_triangular ] ) .
the four - spin interaction @xmath235 [ eq .  ] or the inter - hexagon third - neighbor interaction @xmath236 [ eq .  ] , that is not included in the model , should do the job .
this could allow us to reach the @xmath3 liquid state at both @xmath10 and @xmath0 magnetization plateaus when @xmath237 , while the transition between the crystalline and topological phases could in principle be monitored by examining the rnyi entropy behavior .
this extension of the bfg model would certainly be a very interesting playground to investigate such topological phases , their detections and their transitions to conventional symmetry - breaking states .
also our construction is not restricted to spin-1/2 systems ( or to the equivalent hardcore - boson systems ) and applies to higher - spin systems as well .
in fact , we can obtain the same effective models ( qdm and qlm ) for the spin-1 version of the model with an additional single - ion anisotropy @xmath238 .    as a final remark ,
let us briefly comment on possible experimental realization of the bfg model .
because of the artificial interactions required to build the ice manifold of the bfg model , it is difficult to make direct connections between the model we considered and the existing kagome compounds @xcite
. however , it remains an interesting prospect to investigate whether such toy models could be realized experimentally in tunable artificial systems .
indeed , it was suggested very recently that the bfg interactions and hexagonal plaquettes could be reproduced using a 2d cold ion crystal , with an example given for one and two hexagons  @xcite .
_ note added : _ while finishing this manuscript , we became aware of the recent preprint by roychowdhury et al .
@xcite who recently derived the same loop model for the equivalent hardcore - boson hamiltonian , and further studied the phase diagram of the loop model with a potential term .
while we agree on the phase structure of the loop model , the position of the critical point separating the stripe from the @xmath3 liquid phase is noticeably different [ @xmath239 versus @xmath240 ] .
we acknowledge useful discussions with g. misguich , r. melko , and d. schwandt .
we also thank f. pollmann for the helpful communication on their results .
this work was performed using hpc resources from genci ( grants x2014050225 and x2015050225 ) and calmip ( grants 2014-p0677 and 2015-p0677 ) .
our qmc sse simulations are based on the code from the alps libraries  @xcite .
one of the authors ( k.t . ) was supported in part by jsps kakenhi grant no .
24540402 and no . 15k05211
let us consider @xmath213 states @xmath241 , each corresponding to one of the degenerate states in the thermodynamic limit .
this would be for instance @xmath242 and @xmath243 in an ising phase , or the three states with down spins occupying one of the three kagom sublattices ( and up spins on the other sites , see fig .  [
fig : qmccorr ] ) , in the stripe crystal phase discussed in the main text .
these @xmath213 states are always orthogonal in the thermodynamic limit .
since they are simple product states , one can readily show , for any block @xmath81 of the system which complement scales with system size , that the reduced density matrix corresponding to a cat superposition of these states , @xmath244 , is given by : @xmath245 note that this equality only holds for the reduced density matrix , and not the whole density matrix .    on any finite system , since these @xmath213 states are exponentially close in energy , and a qmc simulation would sample each @xmath246 with equal probabilities .
hence , one would compute a reduced density matrix corresponding to a mixed state : @xmath247 which is precisely the same object as eq .
( [ eq : cat ] ) .
so our qmc simulation is able to access the reduced density matrix of the cat state , provided it is ergodic .
now , each @xmath248 is excessively simple for each product state , as its spectrum contains only one non - zero eigenvalue @xmath249 since they are not entangled .
moreover , they can be diagonalized simultaneously ( they are diagonal in the @xmath186 basis in our example ) , which leads to an entanglement spectrum with non - zero eigenvalues @xmath250 with degeneracy @xmath213 . from it , the entanglement entropy ( von neumann or renyi ) is simply @xmath251 .
this results looks quite intuitive indeed and was conjectured in ref .
it was already noticed in ref .
@xcite for a cat state made of @xmath213 degenerate ising ground states .
we emphasize that the same result holds for the mixed state obtained in our qmc simulation . for more realistic wave - functions ,
the fluctuations around the ideal product state gives and area law , but the constant term will be unchanged : @xmath252    as a last remark , we would like to emphasize that in other more complex cases ( e.g. for states which are not one - site product states , such as valence - bond columnar states ) , the discussion is more involved as the constant term can now depend explicitly on the form of the cut as well as on the rnyi index @xmath5 .
74ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
 + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty , ,  and  ,  eds . ,  link:\doibase    10.1007/978 - 3 - 642 - 10589 - 0 [ _ _ ] ,  ed .
( ,  ) link:\doibase 10.1038/nature08917 [ * * , ( ) ] link:\doibase    10.1126/science.1201080 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.109.067201 [ * * ,   ( ) ] link:\doibase 10.1038/nphys2465 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.114.037201 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.61.2376 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.86.1881 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.89.137202 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.65.224412 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.66.205104 [ * * ,   ( ) ] http://link.aps.org/abstract/prb/v69/e064404 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.94.146805 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.97.207204 [ * * ,   ( ) ] link:\doibase    10.1126/science.1212207 [ * * ,   ( ) ] link:\doibase    10.1038/nphys2036 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.69.104431 [ * * ,   ( ) ] link:\doibase    10.1007/s00220 - 007 - 0342-z [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.84.1535 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.79.064412 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nature09910 [ * * ,   ( ) ] link:\doibase 10.1038/ncomms3287 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.71.224109 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.71.020401 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.74.134301 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.74.054402 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.99.127202 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.77.134421 [ * * ,   ( ) ] http://stacks.iop.org/1367-2630/13/i=8/a=085001 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.66.174405 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.66.214513 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.62.7850 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.76.224431 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.81.205102 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.73.245105 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.65.024504 [ * * ,   ( ) ] http://www.sciencedirect.com/science/article/b6wb1-47hjx86-1/2/45c94ccbbe3ebdfc944a7aa4ac9c932f [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.44.2664 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.82.4745 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.98.227202 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.80.180414 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.84.024406 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.85.134407 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.111.087203 [ * * , ( ) ] link:\doibase    10.1103/physrevb.80.212406 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.90.245143 [ * * ,   ( ) ] @noop ( ) link:\doibase    10.1103/physrevlett.99.207203 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.43.5950 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.66.046701 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.70.100410 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.74.214517 [ * * ,   ( ) ] http://stacks.iop.org/0953-8984/19/i=14/a=145203 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.95.127205 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.95.127206 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.95.127207 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.97.147202 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.97.067202 [ * * , ( ) ] link:\doibase    10.1103/physrevb.76.235107 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.79.212402 [ * * ,   ( ) ] http://stacks.iop.org/1742-5468/2004/i=06/a=p06002 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.96.110405 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.96.110404 [ * * ,   ( ) ] link:\doibase 10.1103/revmodphys.82.277 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.75.214407 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.103.261601 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.86.235116 [ * * ,   ( ) ] @noop ( ) ,   `` , ''  ( , )  chap
.  , p.   link:\doibase    10.1088/1367 - 2630/17/6/065018 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.92.075141 [ * * ,   ( ) ] http://stacks.iop.org/1742-5468/2011/i=05/a=p05001 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.86.024424 [ * * ,   ( ) ]