since they have been shown to possess resonating valence bond ( rvb ) phases on the triangular,@xcite kagome@xcite and other ( non - bipartite ) lattices,@xcite quantum dimer models ( qdm ) have been one of the main paradigms in the field of quantum spin liquids .
these models , where the hilbert space is spanned by hard - core dimer coverings of the lattice , are expected to capture the phenomenology of quantum antiferromagnets where the wave function is dominated by short - range valence bond configurations . on the triangular lattice
, the simplest qdm is defined by the hamiltonian : @xmath2 where the sum runs over all plaquettes ( rhombi ) including the three possible orientations .
the kinetic term , controlled by @xmath3 , flips the two dimers on every flippable plaquette , i.e. , on every plaquette with two parallel dimers , while the potential term controlled by the interaction @xmath4 describes a repulsion ( @xmath5 ) or an attraction ( @xmath6 ) between nearest - neighbor dimers .
the rvb phase is now relatively well understood .
the first result goes back to rokhsar and kivelson@xcite who showed that , for @xmath7 ( the rokhsar - kivelson or rk point ) , the ground state is the sum of all configurations with equal amplitudes : @xmath8 since then , dimer - dimer correlations have been shown to be short - ranged in a range of parameters below @xmath7 , @xcite and the excitation spectrum to be gapped.@xcite besides , it has topologically degenerate ground states on non - simply connected clusters.@xcite this degeneracy is not related to a standard symmetry breaking . indeed , there is no local order parameter,@xcite but only @xmath0 topological order.@xcite    in qdm , the nature of the phase transition from a liquid to a solid is a long standing problem .
it goes back to jalabert and sachdev@xcite ( see also @xcite ) who studied a three - dimensional frustrated ising model related to the square lattice qdm .    in a more general context , senthil and
fisher@xcite showed that models of mott insulators can be cast in the form of a @xmath0 gauge theory . in this language
, the transition from a fractionalized insulator to a conventionally ordered insulator appears to be a condensation of @xmath0-vortices ( dubbed _ visons _ ) . using a duality
relation,@xcite they showed that such transitions correspond to an ordering transition in a _ frustrated _ ising model in transverse field .
as we will see , this applies to the present qdm .    building on a mapping between qdm s at @xmath9 and ising models in a transverse field , moessner _
_ @xcite have developed a landau - ginzburg approach and suggested that the transition could be continuous and in the three - dimensional @xmath10 universality class .
numerical evidence in favor of this scenario has been obtained with green s function quantum monte carlo by ralko _
_ , who have shown that , at the transition point , the static form factor of the crystal decreases to zero on the crystal side of the transition , while the dimer gap also decreases to zero on the liquid side.@xcite more recently , the vison spectrum has also been numerically determined using green s function quantum monte carlo , with the conclusion that a soft mode indeed develops at the transition.@xcite    in spite of these results , a simple picture for the wave functions of the visons and the evolution of their spectrum is still missing . one difficulty is that , like any vortex , these excitations can not be created by operators which are local ( in the dimer variables ) . to attack this problem
, we follow two strategies .
first of all , we look at the problem in the context of a @xmath0 gauge theory on the triangular lattice . as usual
, this theory can be mapped onto a dual ising model.@xcite the duality transforms the nonlocal excitations of the gauge theory into local excitations , which we study using a semiclassical approach .
this simple @xmath11 expansion already captures most aspects of the confinement - deconfinement transition ( soft mode and condensation ) of the @xmath0 gauge theory . secondly ,
building on the explicit form of the excitations of the @xmath0 gauge theory , we construct single vison wave functions for the qdm and study the properties of the system at the rk point in the variational hilbert space spanned by these states .
these `` variational visons ''  living on the triangular _ plaquettes _ and experiencing a flux @xmath12 emanating from each _ site _ of the lattice  turn out to be linearly independent .
the associated dispersion is in good qualitative agreement with that obtained in monte carlo simulations and the value of the gap at the rk point ( @xmath13 ) , obtained with a single variational parameter , has the correct order of magnitude ( monte carlo simulations@xcite give 0.089 ) .
improving quantitatively further these variational results would require more adjustable parameters to account in more details for the local ( dimer - dimer , etc ) correlations in the vicinity of the core of the vortex , a task which has not been carried out here .
the connection between qdms and @xmath0 gauge theories has already been discussed from different perspectives ( see in particular refs .  ) , and is rooted in the existence , in both families of models , of ising - like degrees of freedom subjected to local constraints ( hard - core constraints for the dimers , gauss law for the gauge theory ) .    in this section ,
we review two known mappings : i ) from a @xmath0 lattice gauge theory on the hexagonal lattice to the triangular lattice qdm ( valid in a particular limit , sec .
[ ssec : z2 ] ) and ii ) from the @xmath0 theory to its dual frustrated ising model ( sec.[ssec : dim ] ) .
the latter ising model is then studied using a semiclassical approximation in secs .
[ sssec : cpd ] and [ sssec : sca ] .
since the @xmath0 gauge theory - qdm mapping is formally justified only in the confined phase of the gauge theory , the relevance of the results obtained in the rvb ( deconfined ) phase of the qdm is not guaranteed _
a priori_. some reasons why this is expected to be the case are discussed in the next section , when we build on the results obtained for the excitations of the @xmath0 gauge theory in its deconfined phase to construct elementary excitations in the rvb phase of the qdm .
is the set of the three bonds forming the edges of the plaquette @xmath14 .
2 ) @xmath15 $ ] is the set of the six ( fat ) bonds emanating from the triangular site @xmath16 .
3 ) @xmath17 denotes the bonds forming a `` zigzag '' string extending ( to the left ) from the triangular plaquette @xmath14 to an edge of the lattice.,height=132 ]    to write down a @xmath0 gauge theory analogous to a qdm , the starting point is to define pauli matrices @xmath18 on the bonds @xmath19 of the triangular lattice such that @xmath20 if bond @xmath19 is occupied by a dimer and @xmath21 if it is empty , while @xmath22 changes the state of bond @xmath19 .
since each kinetic term of the qdm flips the dimers around a rhombus , it corresponds to a product of four @xmath22 around this rhombus .
this would , however , not keep the structure of the simple @xmath0 gauge theories , in which the kinetic term acts on an elementary plaquette of the lattice , a crucial ingredient to get a simple ising model by duality.@xcite an alternative is to consider the more standard @xmath0 gauge theory which , in its hamiltonian formulation , is defined by : @xmath23 where @xmath14 runs over the sites of the dual honeycomb lattice , and @xmath24 are the three bonds forming the triangular plaquette around site @xmath14 ( see fig .  [ fig : def ] ) .
the hallmark of this model is to have local conserved quantities .
indeed , @xmath25 } \tau^x_{l[a]}]=0,\ ] ] where @xmath16 is a site of the triangular lattice , and the product over @xmath15 $ ] runs over the six links emanating from @xmath16 ( see fig .
[ fig : def ] ) .
this allows one to define different sectors according to whether @xmath26 } \tau^x_{[a]}$ ] is equal to @xmath21 or @xmath27 .
since in the qdm the number of dimers emanating from a given site is exactly equal to 1 , it is clearly better to consider the sector where @xmath28 } \tau^x_{l[a]}=-1\ ] ] for all @xmath16 since this forces the number of dimers emanating from a site to be odd ( defining an _ odd ising gauge theory _ in the terminology of ref .  ) .
then , the true constraint is recovered in the limit @xmath29 if @xmath30 since , in the ground state , the number of dimers is then minimal .
a _ bona fide _
qdm is then recovered if @xmath31 is treated with degenerate perturbation theory .
since @xmath31 changes the number of dimers , its effect vanishes to first order . to second order
however , one recovers exactly the qdm of eq .
( 1 ) with @xmath9 and @xmath32 .
so , if @xmath33 , the @xmath0 gauge theory maps onto the qdm at @xmath34 . at finite @xmath35 the model can be viewed as a  soft - dimer " model , where 1 , 3 or 5 dimers may touch a given site .
as we shall see below , the limit @xmath33 lies deep inside the confined phase of the gauge theory , and from previous work on the qdm , it is known that at @xmath34 , the model is in a valence bond crystal phase .
it would , of course , be very interesting to connect the two models away from this limit , when the gauge theory is in its deconfined phase and the qdm is in the dimer liquid phase .
an interesting step in this direction is provided by higher order perturbation theory .
indeed , to fourth order in @xmath36 , a repulsion between dimers @xmath37 is generated .
however , other terms of the same order , involving dimer shifts along loops of length 6 , are also generated ( see appendix [ sec:4th ] ) .
so , a _ rigorous
_ mapping does not extend beyond the @xmath34 point .
as usual , this @xmath0 gauge theory is best analyzed by mapping it onto a dual ising model in a transverse field .
this can be achieved by introducing spin-@xmath38 operators @xmath39 on the dual honeycomb lattice : @xmath40 where @xmath17 represents all bonds cutting a straight path ( say , horizontal , see fig .
[ fig : string ] or [ fig : def ] ) starting at @xmath14 ( we implicitly assume a finite lattice with open boundary conditions ) .
combined with the constraint , this definition implies that @xmath41 for two neighboring sites @xmath42 separated by the bond @xmath19 .
@xmath43 is such that each hexagon has exactly one @xmath44 bond ( see fig .
[ fig : hex ] ) . in terms of these spin operators ,
the hamiltonian is the fully frustrated ising model ( ffim)@xcite discussed by moessner and sondhi:@xcite@xmath45 gauge theories appearing in the description of mott insulators . ] @xmath46 note that choosing other paths to define @xmath47 leads to other signs for @xmath48 , which , however , are always such that an odd number of minus signs appear around each hexagon , leading to the same ising model up to a gauge transformation .
the present choice leads to the smallest unit cell ( 4 honeycomb sites ) .     connecting the boundary of the system to triangles @xmath14 and @xmath49 .
dotted segment : path @xmath50 . in this example
@xmath50 crosses one dimer of the reference configuration @xmath51 ( ellipses ) and , therefore , @xmath52 ( eq .  [ eq : epsilon]),height=143 ]    in the following , we study the phase diagram and the excitations of this model within a semiclassical ( large @xmath53 ) approximation .
first , we determine the classical ground state of the model as a function of @xmath36 . to this end
, we replace the spin-@xmath38 operators by classical 3-component vectors of unit length .
since the @xmath54 component does not appear in the hamiltonian , it is clear that in the ground state the magnetization must lie in the @xmath55-@xmath56 plane . to find the lowest energy solution , we use the following parametrization : @xmath57 with @xmath58 .
the corresponding energy is @xmath59    initially defined for nearest neighbors ( eq .  [ eq : sst ] ) , the coefficients @xmath48 have been upgraded to a matrix @xmath60 by setting all other elements to 0 .
this matrix describes the motion of a particle on a honeycomb lattice with 4 sites per unit cell and a flux @xmath12 per hexagonal plaquette .
it reduces to a 4@xmath614 matrix after fourier transformation .
the eigenvalues associated with the momentum @xmath62 are@xcite @xmath63    } } .&\nonumber\end{aligned}\ ] ] likewise , we denote by @xmath64 the column vector of components @xmath65 .    for @xmath66 , we expect the @xmath67 s to be small ( or zero ) .
we can , therefore , expand eq .
[ eq : ecla ] to quadratic order : @xmath68 the largest eigenvalue of @xmath60 being @xmath69 , we find that , as long as @xmath70 , the energy is minimized by @xmath71 and all spins point in the @xmath55 direction .    at @xmath72 , all the real eigenvectors of @xmath60 for the eigenvalue @xmath69 ( satisfying @xmath73 ) minimize eq .
[ eq : ecla2 ] .
let us choose four complex vectors which form a basis of the subspace associated with the eigenvalue @xmath69:@xcite @xmath74 \exp\left(\frac{i\pi}{6}x +             i\frac{\pi}{2}y\right),\nonumber \\
v_2(x , y)&=&v_4^*=\left[\begin{array}{c }             e^{i\pi/12 } \\
e^{-5i\pi/6}\\             1/f\\             e^{-5i\pi/12}/f            \end{array}\right ] \exp\left(i\frac{5\pi}{6}x +      i\frac{\pi}{2}y\right),\nonumber
\\            f&=&2\sin(5\pi/12),\end{aligned}\ ] ] where the 4 entries of the vector refer to the 4 sites in the unit cell ( numbered as in fig .
[ fig : hex ] ) and @xmath75 are the ( bravais ) coordinates of the unit cell .
these eigenvectors correspond to the four points labelled b in the rectangular brillouin zone of fig .
[ fig : bril ] .    , the other bonds have @xmath76 .
the four sites of the unit cell ( small blue triangles ) are labelled from 1 to 4 and the bravais vector are @xmath55 and @xmath54.,height=143 ]    the most general real eigenvector @xmath64 can be parametrized by three angles @xmath77 , @xmath78 , and @xmath79 and a normalization factor @xmath80 : @xmath81 .
\label{eq : rho_angles}\ ] ] to find the ground state , the energy has to be minimized as a function of these 3 parameters .
the analysis is made easier when one realizes that the second and fourth moments of the spin deviations are independent of the three angles .
indeed , @xmath82 replacing @xmath64 by eq .
[ eq : rho_angles ] in the expression for the energy ( eq .
[ eq : ecla ] ) , we get up to order @xmath83 : @xmath84 minimizing with respect to @xmath85 gives : @xmath86 for @xmath87 the energy is minimized for @xmath88 . for @xmath89 , we have to expand eq .
[ eq : ecla ] to the 6@xmath90 order to find the angles @xmath77 , @xmath78 and @xmath79 which minimize the classical energy . indeed , eq .
[ eq : elambda ] shows that , _ up to order @xmath83 , the energy is independent of the three angles_. at this order , the energy is constant and minimum on a 3-dimensional sphere , a consequence of the @xmath10 symmetry discovered by moessner and sondhi.@xcite one , therefore , has to go to the next order to find the actual minima . at 6@xmath90 order in @xmath80 ,
the energy is minimized when @xmath91 is minimum .
a numerical investigation shows that the solutions ( for @xmath36 close to but below @xmath69 ) are 48-fold degenerate and can be deduced from each other by symmetry operations ( 12 translations and 4 point - group operations ) .
this confirms the result obtained previously on symmetry grounds.@xcite    motivated by the ffim / dimer model correspondence , we are interested in the average `` dimer density '' : @xmath92 for all pairs @xmath93 of nearest - neighbor honeycomb sites . in the classical limit , this may be approximated by @xmath94 where @xmath64 ( eq .  [ eq : rho_angles ] ) is the classical ground state .
close to the transition at @xmath95 , @xmath64 scales as @xmath96 ( eq .  [ eq : lambda2 ] ) and the `` dimer density '' , thus , shows small deviations about @xmath38 .
to visualize the `` dimerization '' pattern in the vicinity of @xmath97 , the appropriate quantity is , therefore , a relative rescaled `` dimer '' density defined by @xmath98 , and plotted in fig .
[ fig : vbc ] for one of the 48 ground states .
the obtained pattern is highly reminiscent of the @xmath1 vbc observed in the triangular lattice qdm .
in particular , the 48 spin configurations give only 12 different `` dimer '' patterns , with a smaller unit cell containing 12 sites of the triangular lattice .
it is also interesting to notice that @xmath67 vanishes in the triangles located inside the large diamonds ( marked with a dot in fig .
[ fig : vbc ] ) . from eq .
[ eq : sigma_rho ] , the corresponding sites have a magnetization pointing exactly in the @xmath55 direction , which means a total absence of `` vison '' and local correlations identical to that of the @xmath99 paramagnetic ground state .
the `` dimer '' density @xmath100 is _ uniform _ inside these diamonds , in qualitative agreement with the intriguing observation made in ref .  .     in one of the 48 classical ground states of the ffim in transverse field , for a transverse field @xmath101 just below @xmath97
thin red bonds represent @xmath102 ( corresponding to the _ highest _ `` dimer density '' @xmath103 ) .
blue bonds are for @xmath104 , that is , a _ lower _ dimer density , with a width proportional to @xmath105 .
@xmath106 takes only 4 different values : 0 , -0.259 , -0.518 , and -0.776 .
, width=257 ]      if @xmath99 , the ground state is the fully polarized state @xmath108 .
the first ( degenerate ) excited state is obtained by flipping one spin at some arbitrary triangle . to first order in @xmath35 ,
the spectrum is obtained by diagonalizing the perturbation @xmath109 in this subspace .
this amounts to diagonalizing @xmath60 , and leads to the eigenvalues @xmath110 , a result already obtained previously.@xcite    to go further , beyond the limit @xmath111 , we perform a @xmath11 semiclassical expansion .
we generalize the spin-1/2 hamiltonian to an arbitrary value @xmath53 of the spin : @xmath112 ( which reduces to eq .
[ eq : imtf ] when @xmath113 ) .
spin deviations away from the @xmath55 directions are represented using holstein - primakoff bosons . to leading order in @xmath11 : @xmath114
one truncates the hamiltonian to quadratic order in @xmath115 , and diagonalizes it through a bogoliubov transformation : @xmath116 for the @xmath117 operators to be bosonic creation operators , the matrices @xmath118 and @xmath4 must satisfy @xmath119 . it is convenient to introduce a unitary matrix @xmath120 which transforms @xmath60 into a diagonal matrix @xmath121 : @xmath122 in this new basis , we can look for diagonal solutions for @xmath118 and @xmath4 : @xmath123 since @xmath60 is real , we can choose @xmath124 .
after some algebra , one finds that the @xmath125 and @xmath126 which diagonalize @xmath127 are given by @xmath128 the energies of the bogoliubov excitations are given by @xmath129 as expected , the spectrum becomes gapless at @xmath72 and we recover the localized spin flip energy @xmath130 when @xmath131 .
the bogoliubov spectrum is shown fig .
[ fig : sw ] for a few values of @xmath101 at or above @xmath97 .
one clearly sees that the dispersion is linear around the gapless points at the transition .
( eq .  [ eq : spinwave ] ) for the path a@xmath132b@xmath132c@xmath132a in the brillouin zone ( see fig .  [
fig : bril ] ) and different values of @xmath36:3 ( top ) , 2.75 , 2.5 , and @xmath133 ( bottom ) . at the critical point @xmath72 ,
the spectrum is linear around @xmath134 and @xmath135.,width=264 ]     and @xmath136 ) of the hexagonal lattice shown in fig .
[ fig : hex ] .
a , b , and c are the high - symmetry points used in fig .
[ fig : sw ] .
the large hexagon is the brillouin zone of the underlying triangular lattice .
, width=200 ]    this behavior is remarkably similar to that found numerically by ralko _
et al._@xcite for the qdm by green s function quantum monte carlo . in the next section ,
we build on this resemblance to develop a variational approach to the vison spectrum of the qdm .
note that , as stated above , the @xmath0 gauge theory and the qdm can only be rigorously mapped onto each other deep into the confined ( resp .
vbc ) phase , and the resemblance between their spectra in the deconfined ( resp .
rvb ) phases might seem fortuitous .
a somewhat deeper connection will be described in the next section .
in the ising model , elementary excitations for large enough @xmath36 are spin flips , induced by @xmath47 , which delocalize and get dressed under the effect of @xmath109 . in the equivalent dual gauge theory ,
this excitation is produced by the nonlocal string operator @xmath137 . by analogy , it is natural to define a `` point '' vison creation operator for the qdm by@xcite @xmath138 where the dimer operator @xmath139 is defined by @xmath140 if bond @xmath19 is occupied , and @xmath141 otherwise .
again , we consider here a finite lattice with open boundary conditions .
let us first see to which extent this operator is the analog of the vortex creation operator in ising gauge theories .
in gauge theories , the wilson loop operator defined by @xmath142 plays a central role since it allows one to distinguish ( in the absence of matter field , as here ) the deconfined and confined phases depending on whether its ground state expectation value ( or flux ) tends to zero exponentially with the perimeter of the domain @xmath120 ( deconfined ) or with the area of the domain ( confined ) . for a @xmath0 gauge theory @xmath143 and the flux going through @xmath120 measured by @xmath144 can only take two values @xmath145 . in that respect ,
an interesting property of the operator @xmath137 is that it changes the flux between -1 and + 1 if the site @xmath14 is inside the domain .
then , @xmath146 commutes with the wilson loop and does not change the flux unless @xmath14 and @xmath49 sit on _ opposite sides _ of the boundary .
since the deconfined phase can be accessed from the @xmath147 limit , the ground state can be thought of as an ising paramagnet ( @xmath148 and @xmath149 everywhere ) `` dressed '' perturbatively by successive applications of @xmath150 , where @xmath14 and @xmath49 are nearest neighbors . since @xmath150 only changes the flux for pairs @xmath151 _ across the boundary _ , the expectation value of the wilson loop operator behaves according to a _ perimeter
a similar perturbative argument can be used to derive the area law in the confined phase .
most of these standard @xmath0 gauge theory results apply to the qdm , with one difficulty however .
the wilson loop operator can not be defined in the same way since flipping the dimer occupation along a loop will often lead to an unphysical state that violates the condition of having exactly one dimer emanating from each site .
different ways to overcome this problem can be envisaged .
one possibility is to accept that the wilson loop operator gives zero when applied to states that are non - flippable along the loop .
the property @xmath143 is lost , but @xmath144 has eigenvalues @xmath141 , @xmath21 , or @xmath27 , and eigenstates with eigenvalues @xmath21 or @xmath27 are still interchanged when a vison operator ( eq .  [ eq : bare_vison ] ) is applied inside the domain .
so confinement or deconfinement of visons is still expected to lead to area or perimeter laws .    alternatively ,
starting from the observation that , when it does not give zero , the wilson loop operator just shifts the dimers along the contour , one can try to define a ( complicated ) flux operator in the dimer space which shifts the dimer along a `` fattened '' contour which `` adapts '' locally to the dimer configuration it acts upon .
although tricky to define in practice , this operator would have the advantage that its square is still equal to 1 , preserving the manifest @xmath0 structure of the theory .    another way to underline the deep connection between the two models
is to consider the ising model as a `` soft - dimer '' model and introduce the projection operator onto the `` hard - core '' dimer hilbert space , defined by : @xmath152 where the operator @xmath153 counts the number of frustrated bonds ( dimers ) around the plaquette @xmath154 : @xmath155 by construction , @xmath156 on all the ising configurations which correspond to a valid  hard - core " dimer covering ( of the triangular lattice ) and @xmath157 otherwise .
now , @xmath158 commutes with @xmath159 since all terms in @xmath158 contain an even number of @xmath160 operators .
this means that @xmath158 conserves the total flux , so that a spin state with -1 flux becomes a dimer wave function with an odd number of visons after projection .
equation  [ eq : bare_vison ] defines the simplest operator which changes the flux , that is , which creates a vison .
so we may consider @xmath161 as a first variational approximation to the true lowest eigenstate of @xmath162 in the @xmath27 flux sector ( @xmath163 is the total number of dimer coverings ) .
the ground state @xmath164 is a zero - energy eigenstate of @xmath162 , implying that the expectation value @xmath165 of the kinetic energy term exactly compensates the expectation value @xmath166 of the potential energy term .
indeed , @xmath167 where @xmath168 is the total number of rhombi ( @xmath169 the number of sites ) and @xmath170 the probability to have 2 parallel dimers on a given rhombus in the classical dimer problem ( with uniform measure over all dimer configurations ) . using the pfaffians method,@xcite one finds in the thermodynamic limit @xmath171    in @xmath172 the expectation value of the potential energy term is the same as in @xmath164 .
in fact , all the observables which are diagonal in the dimer basis commute with @xmath160 and thus have the same expectation value in the ground state @xmath164 and in the trial wave function @xmath172 .
the increase of the energy is only kinetic . because the kinetic terms corresponding to the three diamonds around @xmath14 anticommute with @xmath47 , their expectation values has changed sign .
thus we find @xmath173      one can lower the energy by constructing plane waves . to compute the associated dispersion relation we have to evaluate the following matrix elements : @xmath174 we begin with the overlap matrix @xmath175 ( @xmath176 and @xmath177 are defined in eq .
[ eq : bare_vison ] ) which simplifies to @xmath178 where the _ local _ operator @xmath179 counts the number of dimers across some path @xmath50 connecting the triangles @xmath14 and @xmath49 ( see fig .
[ fig : string ] ) , and @xmath180 is equal to that number in the reference configuration @xmath51 , chosen with all the dimers horizontal ( fig .
[ fig : string ] ) .
this follows from two simple properties : @xmath181 is independent of the configuration @xmath182 , and @xmath183 .
we finally write @xmath184 where @xmath185 represents the average with equal weight over all dimer coverings .
the average of any such diagonal observable can be computed using the pfaffian of a ( modified ) kasteleyn matrix.@xcite in the present case , the `` string '' observable @xmath186 is coded in the kasteleyn matrix by changing the signs of the matrix elements corresponding to bonds crossed by @xmath50 , _
i.e. _ , setting some bond fugacities to @xmath27 . to evaluate numerically such an expectation value , we construct the modified kasteleyn matrices corresponding to a large enough triangular lattice , in which the path @xmath50 is embedded .
finite - size effects decay exponentially with the system size , so that lattices with @xmath187 sites ( with periodic boundary conditions ) can safely be used to evaluate @xmath188 with high accuracy up to distances @xmath189 between triangles @xmath14 and @xmath49 .
the matrix elements @xmath188 , plotted in fig .
[ fig : sij ] as a function of the distance between @xmath14 and @xmath49 , decay exponentially ( a result anticipated by read and chakraborty@xcite ) .
these quantities have already been evaluated by ioselevich _
et al._@xcite using a classical monte carlo sampling .     between two point - vison states
( defined in eq .
[ eq : point_vison ] ) as a function the distance @xmath100 between @xmath14 and @xmath49 .
the ( blue ) square correspond to eq .
[ eq : sdef ] and the three ( red ) crosses correspond to eq .
[ eq : sija ] with @xmath190 .
the calculations are done using an exact evaluation of the pfaffians on a @xmath191-site lattice with periodic boundary conditions .
dashed line : guide to the eye corresponding to an exponential decay with the dimer - dimer correlation length @xmath192.@xcite , width=264 ]     ( eqs.[eq : sdef ] and [ eq : sija ] ) , for different values of @xmath193 ( 0 , -0.2 , -0.5 and -0.8 ) , along the path a@xmath132b@xmath132c@xmath132a in the brillouin zone ( see fig .
[ fig : bril ] ) .
the bottom panel is a zoom on the lowest eigenvalues of @xmath194 . in these calculations ,
the matrix elements @xmath195 are neglected for triangles @xmath14 and @xmath49 at distance @xmath196 ( 102@xmath90 neighbor on the hexagonal lattice ) . up to this distance , the finite - size lattice used in the calculation ( @xmath187 sites ) gives practically the infinite volume limit for @xmath195 .
the value of @xmath197 used here is large enough to ensure a good convergence of the spectrum since the curves obtained with a smaller truncation distance ( @xmath198  shown here with the same colors ) are almost superposed with that for @xmath199 , except for the bottom of the spectrum at @xmath200 and @xmath190 .
the overlap spectrum turns out to be gapped for @xmath201 and @xmath202 ( and probably at @xmath200 and @xmath190 too ) , indicating the linear independence of the vison states.,width=415 ]    by construction , a product like @xmath203 ( where @xmath204 form a closed loop of triangles ) is equal to the parity of the number of sites enclosed in the loop .
thus , the signs of the matrix elements @xmath188 are similar to those of the hopping amplitude of a particle moving on the hexagonal lattice and subjected to a magnetic field corresponding to half a flux quantum per hexagon . in such a case
, the magnetic unit cell has to be doubled compared to the original lattice cell .
since the original unit cell contains one triangular site and two triangular plaquettes , the magnetic one contains four triangular plaquettes and thus four sites of the hexagonal lattice . in
fourier space , @xmath205 is a @xmath206 matrix , as the matrix @xmath60 discussed previously .
the same is also true for the hamiltonian matrix elements described below .    the spectrum of the overlap matrix @xmath205 is plotted in fig .
[ fig : sk ] ( for @xmath207 describing a representative path in the brillouin zone ) . as an important result , the eigenvalues are _ strictly positive _ for all @xmath207 . to our knowledge , it is the first time that the _ linear independence _ of point vison states is proved .
let us now turn to the matrix elements of the hamiltonian .
we wish to transform eq .
[ eq : hdef ] into an expression which can be evaluated using the pfaffians , that is the expectation value of a diagonal observable in the dimer basis .
first , we write @xmath162 as a sum of projectors @xmath208 and expand eq .
[ eq : hdef ] into @xmath209 with the notation @xmath210 .
because of the double sum @xmath211 , this does not yet have the form of a diagonal observable amenable to an evaluation with pfaffians . to go further we note that @xmath212 vanishes if @xmath213 and/or @xmath214 is not flippable around the rhombus @xmath215 .
in addition , @xmath216 if @xmath213 and @xmath214 differ anywhere outside @xmath215 .
so we can restrict the sum to pairs of configurations @xmath217 which differ by a dimer flip in @xmath215 and which are identical elsewhere on the lattice : @xmath218 \nonumber\\          & & \hat\pi_{r_0 }              \left [ ( -1)^{n(c , j)}| c \rangle   + ( -1)^{n(\bar c , j ) } | \bar c \rangle \right]\end{aligned}\ ] ] where @xmath219 if @xmath220 is flippable on rhombus @xmath215 , and @xmath221 otherwise .
@xmath220 and @xmath222 only differ inside @xmath215 , so the signs @xmath223 and @xmath224 are the same if @xmath14 is not inside @xmath225 , and are opposite if @xmath226 .
let us note @xmath227 if @xmath226 and @xmath228 otherwise .
this leads to @xmath229          \hat\pi_{r_0 }           \left [ | c \rangle   + \delta_{j,{r_0 } } | \bar c \rangle \right ] .
\label{eq : hij0}\end{aligned}\ ] ] using the explicit form of @xmath230 , we get @xmath231          \hat\pi_{r_0 }           \left [ | c \rangle   + \delta_{j,{r_0 } } | \bar c \rangle \right ]          = \frac{1}{2}(1-\delta_{i,{r_0 } } ) ( 1-\delta_{j,{r_0 } } )          \nonumber\\\end{aligned}\ ] ] and , finally , @xmath232 the last expression is the average of a diagonal observable and can , thus , be evaluated using pfaffians : @xmath233 where we used an operator notation for the `` flippability '' @xmath234 . to get a nonzero contribution
, there must be at least one rhombus @xmath215 containing @xmath14 and @xmath49 .
so @xmath235 if @xmath14 and @xmath49 are not nearest neighbors .
using modified kasteleyn matrices in a way similar to that described for the overlap matrix @xmath236 , the nonzero matrix elements can be calculated .
we find @xmath237 for @xmath238 , @xmath239 when @xmath14 and @xmath49 are first neighbors .
although the matrix elements of @xmath162 are simple , the band structure is , however , not that of a simple tight - binding hamiltonian , because of the non - orthogonality of the present variational vison states . in fourier space , @xmath205 and @xmath240 are @xmath206 matrices .
the spectrum is obtained by solving the generalized eigenvalue problem @xmath241 .
the results are shown in fig .
[ fig : var_disp ] ( @xmath201 curve ) .
the minimum of @xmath242 is found at @xmath243 , in agreement with the monte carlo results.@xcite the corresponding energy is @xmath244 , which is significantly lower than the energy ( 0.56 ) found above for a completely localized point vison @xmath172 .
the bottom of this variational point vison band is , however , still high compared to ivanovs@xcite estimate ( 0.089 ) of the gap in the vison sector .
the next section provides an improved family of variational states .    ) dispersion relation for different values of @xmath193 , along the path a@xmath132b@xmath132c@xmath132a in the brillouin zone ( see fig .
[ fig : bril ] ) .
the curve for @xmath201 corresponds to the variational energies of point - vison states ( eq .  [ eq : point_vison ] ) .
the absolute minimum is @xmath245 , found at @xmath246 for a value @xmath247 of the variational parameter @xmath193 .
the exact value of the gap ( 0.089)@xcite is marked as a dotted horizontal line . in this calculation , the matrix elements @xmath195 are neglected for triangles @xmath14 and @xmath49 at distance @xmath196 . this value is large enough to ensure a perfect convergence of the spectrum below @xmath248 , as can be checked from the fact that the dispersion curves obtained with a smaller truncation distance ( @xmath198  shown here with the same colors ) are practically identical
. some higher energy states , however , are not fully converged.,width=264 ]      the vison wave function of eq .
[ eq : point_vison ] , differs from the rk wave functions only through minus signs .
it has the `` accidental '' property that dimer - dimer correlations are the same as in the ground state . as a consequence ,
the excitation energy of such a state is carried only by the kinetic part of the hamiltonian .
clearly , some _ local _ reweighting of the dimer configurations in the vicinity of the vortex core would allow an optimized balance between the potential and kinetic costs of the excitation , and would lead to an improved variational wave function .
this amounts to `` dressing '' locally the initial point - vison state by some _ even _  but fluctuating  number of additional point visons .    as a simple improvement of the vison wave function
, we introduce a variational parameter @xmath193 to reweight the configurations depending on their `` flipability '' at the core of the vison .
this gives the following vison state : @xmath249 where @xmath250 if the dimer configuration @xmath220 is flippable around one of the three rhombi @xmath251 , @xmath252 , and @xmath253 containing the triangle @xmath14 , and @xmath254 otherwise .
as for the point vison of eq .
[ eq : point_vison ] , the vison states of eq .
[ eq : dressed_vison ] are not orthogonal and we have to evaluate their overlaps : @xmath255 repeating the transformation leading to eq .
[ eq : saverage ] , we get @xmath256 which is an expectation value for a diagonal operator that we evaluate using pfaffians . in addition to the sign changes due to @xmath186
, some entries of the kasteleyn matrix have to be modified to incorporate the flippability operators @xmath257 .
more precisely , counting only the coverings which satisfy @xmath258 is done by `` isolating '' the rhombus @xmath259 , that is , by switching to zero in the kasteleyn matrix the 14 bonds which connect the sites of rhombus @xmath259 to their neighbors outside @xmath259 .    beyond some distance between the vison cores @xmath14 and
@xmath49 ( fourth neighbor on the hexagonal lattice ) , no rhombus can touch simultaneously both triangles . in that case
, it can be shown that @xmath260 is independent of @xmath193 .
the eigenvalues of the overlap matrix @xmath261 are displayed in fig .
[ fig : sk ] for a few selected values of @xmath193 .
although a full convergence as a function of the truncation distance @xmath197 ( see caption of fig .  [
fig : sk ] ) has not been obtained , we believe that there is a finite gap for all the values of @xmath193 shown here , and that the dressed vison states are linearly independent .    the evaluation of the hamiltonian matrix elements @xmath262 for dressed vison can still be done using pfaffians , but the algebraic manipulations are slightly more lengthy than the previous ones and we refer the reader to appendix  [ app : hija ] .
we find nonzero matrix elements @xmath263 up to ( and including ) the 9@xmath90 neighbor ( compared to first neighbor for point vison ) .
the results for the variational dispersion relation are shown fig .
[ fig : var_disp ] .
the qualitative shape of the lowest band is almost unchanged compared to the point vison states ( @xmath201 ) , except for an almost uniform shift which lowers the gap .
the minimum of @xmath264 is found at @xmath246 for a variational parameter @xmath265 , and the corresponding energy is @xmath245 , about @xmath266 higher than the exact value .
the fact that simple wave functions like those of eq .
[ eq : point_vison ] or eq .
[ eq : dressed_vison ] reproduce qualitatively the shape of the exact dispersion relation is presumably due to the fact that the dimer - dimer correlation length is rather small at the rk point ( of order of one lattice spacing@xcite ) . as a consequence , the exact vison states only differ from eq .
[ eq : point_vison ] at short distances from the vortex core , and the long - distance part ( string ) , responsible for the flux @xmath12 per hexagon , is essentially exact .
this is , of course , no longer true away from the rk point and in the direction of the crystal , where the correlation length rapidly grows .    in an attempt to extend this variational approximation _
away _ from the rk point , we computed the gap of the dressed vison states in perturbation theory to first order in @xmath267 .
however , at this order , the gap turns out to close very slowly and does not lead to a meaningful estimate for the critical @xmath4 at the liquid / crystal transition .
this failure is closely related to the remark above : the size of the core of the vison presumably grows rapidly away from the rk point , a feature which can not be accounted for with the present variational states .
in this paper , we have developed two simple approaches to describe the vison excitations of the qdm on the triangular lattice .
the first one is based on a soft - dimer version of the model , which exactly takes the form of a @xmath0 gauge theory .
we have shown that a semiclassical spin - wave approximation to the spectrum of the dual ising theory captures the important fact that the disappearance of the rvb liquid is due to a vison condensation .
measures the presence ( -1 ) or absence ( + 1 ) of a vison , and @xmath160 creates or annihilates a vison ( see eq .
[ eq : psigmaz ] for instance ) .
the classical ground state turns out to have @xmath268 in the dimer liquid phase and @xmath269 in the crystal phase ( sec .
[ sssec : cpd ] ) .
so , the vison creation / annihilation operator acquires a finite expectation value at the transition , the usual signature of a particle condensation .
this is almost identical to the standard hard - core - boson @xmath270 spin-@xmath38 correspondence .
if the particle density is represented by @xmath271 ( as here for visons ) , magnetic long - range order in the @xmath56 ( or @xmath54 ) direction ( for the spin variables ) is equivalent to bose condensation ( off - diagonal long - range order in the boson operators ) .
an important difference is , however , that the number of visons is _ not _ conserved . only their _ parity _
is conserved .
accordingly , the vison condensed phase spontaneously breaks a discrete ( @xmath0 ) gauge symmetry ( @xmath160 , which is not gauge invariant , acquires a finite expectation value in the crystal phase ) , and not a continuous [ @xmath272 one .
consequently , the condensed phase is gapped and does not have a goldstone mode . in the dimer model ,
the spectrum is , therefore , gapless only _ at _ the transition . ]
it also reproduces the qualitative shape of the vison spectrum in the disordered phase and , more importantly , at the transition ( linear spectrum at the correct points of the brillouin zone ) , as well as the spatial pattern of the ordered crystalline state .
the second approach is a variational approximation to the vison wave functions at the rk point of the ( hard - core ) qdm . it reproduces semiquantitatively the vison dispersion relation , and provides a simple picture for the vison wave functions which goes beyond the naive point - vison approximation .    beyond the problem of vison dispersion and condensation
, we expect these approaches to be useful for other problems related to vison excitations in qdm , in particular , their mutual interaction and their interaction with vacancies and mobile holes .
this is left for future investigations .
we are very grateful to f.  becca , m.  ferrero , d.  ivanov , v.  pasquier , and a.  ralko for numerous insightful discussions .
fm acknowledges the financial support of the swiss national fund , of manep , and of the rgion midi - pyrnes through its chaire dexcellence pierre de fermat program , and thanks the ipht ( saclay ) and the laboratoire de physique thorique ( toulouse ) for their hospitality .
the goal of this appendix is to derive an effective hamiltonian for the model of section ii.a defined by the hamiltonian @xmath273 in the limit @xmath274 to fourth order in @xmath36 . for simplicity ,
let us define the unperturbed hamiltonian @xmath275 and the perturbation @xmath276 .
since the hilbert space of the model is restricted by the constraint that the number of dimers starting from a given site is odd , the ground state manifold of the unperturbed hamiltonian for positive @xmath277 consists of all states having exactly one dimer emanating from each site .
let us denote by @xmath278 the projector onto the hilbert space generated by these states , and by @xmath53 the resolvent defined by @xmath279 it is easy to check that @xmath4 changes the parity of the number of dimers .
indeed , the term of @xmath4 acting on the triangle @xmath14 transforms the states with 0 and 3 ( resp . 1 and 2 ) dimers around the triangle @xmath14 into each other .
this implies that @xmath280 , and more generally that the effective hamiltonian only contains even powers of the perturbation .
thanks to the property @xmath280 , the fourth order contribution reduces to 3 terms , and the effective hamiltonian up to fourth order reads : @xmath281    up to a constant , this effective hamiltonian can be written as a qdm acting on 4 and 6-site plaquettes .
the 4-site hamiltonian has the form of the regular rk model with @xmath282 and @xmath37 .
the 6-site hamiltonian only consists of kinetic terms that flip the dimers around the three possible types of 6-site plaquettes shown in fig .
[ fig : boucles6 ] , with amplitude @xmath283 for types ( 1 ) and ( 2 ) , and with amplitude @xmath284 for type ( 3 ) .    .
, width=188 ]
for the dressed vison states ( eq .  [ eq : dressed_vison ] ) , the equivalent of eq .
[ eq : hij0 ] is @xmath285\nonumber \\          & & \hat\pi_{r_0 }           \left [ ( 1+\alpha\hat f_j)| c \rangle+ \delta_{j,{r_0 } } ( 1+\alpha\hat f_j)| \bar c
\rangle \right ] \label{eq : hija}\end{aligned}\ ] ] using the explicit form of the projector @xmath286 , we get @xmath287 unlike eq .
[ eq : hij0b ] , both the coverings @xmath220 and @xmath222 enter the expression .
this does not have the form of a diagonal observable and the terms @xmath288 ( with @xmath289 ) need to be eliminated to allow for an evaluation with pfaffians . by inspecting the possible relative positions of two rhombi @xmath225 and @xmath259 ,
one arrives at the following two relations :    * if two rhombi @xmath225 and @xmath259 have 0,1,3 or 4 sites in common , @xmath290 for any pair of configurations @xmath217 which differ by a flip around the rhombus @xmath225 . *
it @xmath225 and @xmath259 have 2 sites in common , let us call @xmath291 the bond of @xmath259 which does not touch @xmath225 .
then we have @xmath292 , where @xmath293 if @xmath291 is occupied by a dimer of @xmath220 , and 0 otherwise .
we may combine the two cases above into some compact notation @xmath294 , valid for any pair of configurations @xmath217 which differ by a flip around the rhombus @xmath225 .
accordingly , we define an operator @xmath295 for each triangle @xmath14 and rhombus @xmath225
. the matrix element of eq .
[ eq : hija ] is now expressed as the expectation value of a diagonal operator : @xmath296 this expression has to be expanded into a polynomial in @xmath297 and @xmath298 operators before each term can be evaluated thanks to the pfaffian of an appropriate kasteleyn matrix . as before the @xmath186 introduces some sign changes , and each @xmath297 ( or @xmath298 ) operator requires to `` isolating '' the corresponding rhombus ( or bond ) by switching to zero the corresponding matrix elements .
several tens of terms typically appear for each pair of triangles @xmath299 , and an automated treatment by computer had to be coded to obtain the hopping amplitudes .
the results of fig .
[ fig : var_disp ] represent several hundreds of cpu hours using the software maple to generate all the correlators and evaluate the corresponding pfaffians on a finite lattice with @xmath187 sites .