# Criticality of a classical dimer model on the triangular lattice

## Abstract

We consider a classical interacting dimer model which interpolates between the square lattice case and the triangular lattice case by tuning a chemical potential in the diagonal bonds. The interaction energy simply corresponds to the number of plaquettes with parallel dimers. Using transfer matrix calculations, we find in the anisotropic triangular case a succession of different physical phases as the interaction strength is increased: a short range disordered liquid dimer phase at low interactions, then a critical phase similar to the one found for the square lattice, and finally a transition to an ordered columnar phase for large interactions. The existence of the critical phase is in contrast with the belief that criticality for dimer models is ascribed to bipartiteness. For the isotropic triangular case, we have indications that the system undergoes a first order phase transition to an ordered phase, without appearance of an intermediate critical phase.

## I Introduction

Dimer models have regained a lot of interest in the condensed matter physics community thanks to the pioneering work of Rokhsar and Kivelson (RK) (1). In this original proposal, dimers represent singlets formed by pairs of spins in Cuprate materials or frustrated antiferromagnetic systems. This idea has been followed by many other proposals giving rise, at low energy, to such kind of effective models. To cite only a few connections to the physics of dimers, related again to magnetism, the lowest-energy configurations of fully frustrated Ising magnets can generally be mapped onto dimer configurations on the dual lattice (2); (3). Other specific (quantum) dimer models have recently been derived from a spin-orbital model describing LiNiO (4), from the trimerized kagome antiferromagnet (5), or for Heisenberg antiferromagnets under applied field in the pyrochlore lattice (6). Lastly, we note that hard-core bosons or correlated fermions on frustrated lattices like the planar pyrochlore lattice can also be mapped onto dimer representations in the limit of large Coulomb repulsion (7); (8).

Although the RK model is quantum mechanical, there is a special value of the parameters of the Hamiltonian for which the ground-state is an equal weight superposition of all dimer tilings of the square lattice (RK point). Static properties at zero temperature can then be understood by the purely combinatorial problem of counting dimer tilings of the square lattice first solved in the early 60’s (9). This system turns out to be a critical model, and more precisely a conformal field theory with central charge . Subsequently to the RK work, Moessner and Sondhi (10) have shown that the same dimer model on the triangular lattice shows, instead of a single point with critical algebraic correlations, a disordered phase with short range dimer correlations - as expected from the known physics of its classical counterpart (11). From the interpretation originally given to the dimers, it is legitimate to call this phase a spin liquid.

More recently, a study performed by Alet et al. (12) on a square lattice classical dimer model showed the existence of a Kosterlitz-Thouless phase transition from a critical phase to a columnar ordered phase for the dimers. This transition is triggered by including an interaction term, consisting in the number of plaquettes doubly occupied by parallel dimers. Subsequently, a quantum mechanical dimer model can be built from this classical model proving that the single point showing criticality in the original RK Hamiltonian can be promoted to a whole critical phase (13).

An important difference between the square and the triangular lattices for quantum dimer models lies in the degeneracy of the ground state at the RK point: it is finite for the triangular lattice with periodic boundary conditions (PBC) and scales exponentially with the linear size for the square lattice with the same PBC (see Ref. (10) and Ref. (1) respectively for details). This large degeneracy of the square lattice, related to its bipartiteness, has been often designated as the responsible of the existence of critical correlations in this case. It is then natural to investigate the behavior of interacting dimer model on non-bipartite lattices interpolating between the square and the triangular lattices, and in particular try to elucidate the interplay between criticality and bipartiteness. In the non-interacting case, it is found (11) that criticality disappears immediately with the introduction of non-bipartite dimers.

In this paper we study an extension of the classical dimer model with nearest-neighbor interactions already studied elsewhere (12); (13) to the case of an anisotropic triangular lattice. Although a limiting case of this model is the square lattice, in general the lattice is non-bipartite. In Sec. II we start by reviewing the results on the square lattice, and recalling the principles of an associated field theory used to describe the critical phase. In Sec. III, we investigate the properties of the anisotropic model that interpolates between the square and triangular lattices by means of transfer matrix calculations. There are two independent paths for the interpolation: the introduction of a fugacity for diagonal bonds but also of the introduction of interactions between parallel dimers on diamond-like plaquettes (see Fig. 1). We find for a certain range of parameters the existence of a critical phase. We determine its location in the phase diagram as well as transitions from it to either liquid or ordered phases. In Sec. IV, we find that the isotropic model on the triangular lattice, does not possess a critical phase for any interaction strength but could display an ordered phase for sufficiently large interactions. Finally, Sec. V contains a discussion of the various results and conclusions.

## Ii Criticality of the dimer model on the bipartite square lattice

Let us start by reviewing the main results of the square lattice dimer model. In this paper, we will only consider fully-packed dimers with a hard-core constraint: there must be one and only one dimer coming out of each site. The non-interacting model can be solved exactly on the square lattice by expressing the partition function of the system as a Pfaffian, which enables to obtain analytic expressions for the free energy and correlators, on either finite or infinite systems (9); (14). On the square lattice, dimers are found to be in a critical phase: dimer-dimer correlations decay algebraically as with relative distance , whereas the monomer-monomer (defining a monomer as a site with no dimer) correlations behave as . Classical dimers on the square lattice have also been studied within a finite-temperature model with nearest neighbors interactions, i.e. an interaction between parallel dimers on the same plaquette ( corresponding to attractive interactions) (12). Here the model is not integrable anymore and it was found numerically that the system is critical down to a finite temperature , where it gives rise to a low-temperature columnar phase with dimers aligned in columns. Monomer-monomer correlations decay to at long distances, either algebraically (in the critical phase) or exponentially (in the columnar phase). From now on, unless it is specified otherwise, we set and let the different dimer-dimer interactions ( for the purely square lattice) vary. Nevertheless, we still refer to high-T (respectively low-T) regions characterized by small (resp. large) values of the involved couplings.

The criticality encountered in the high-T phase can be understood in the framework of a height field theory. For each dimer configuration on the square lattice, a scalar height field with quantized values (multiples of ) can be defined microscopically on the dual lattice (15). The spatial variations of the height field between neighboring sites of the dual lattice are entirely determined by the presence of dimers between them, thanks to the lattice bipartiteness. The long-wavelength modes of this height correspond to a coarse-grained height field defined in continuum space, and the physics of the model is captured by the action (12); (16):

(1) |

The cosine term of this action is a locking potential that favors flat configurations ( for the square lattice), corresponding microscopically to the columnar configurations. The term, corresponding to the cost of fluctuations around these flat configurations, account for the entropy of dimer coverings. In the non-interacting case, the value of the stiffness is fixed by the exact results of Ref. (9); (14) to be (see also Ref. (17)). The renormalization of the stiffness constant accounts for the role of interaction: increases with .

As long as the cosine operator in Eq. (1) is irrelevant in the renormalization group sense, this action defines a conformal field theory with central charge and the system is critical. The dimer operator (which gives the local dimer density) is composed of a and gradient terms (18), and has a scaling dimension , corresponding to half of the exponent of dimer-dimer correlations. Similarly, the locking potential has a dimension , and becomes relevant when . Within the standard Coulomb gas description (19) of this theory, we can define operators corresponding to the insertion of a particle of electric and magnetic charges and respectively. In general, the dimension of such an electromagnetic vertex operator is . In the height model discussed here, the operator with is the dimer operator already mentioned; similarly, the monomer operator corresponds to inserting a magnetic particle with (the sign of the charge depends on the sublattice where the monomer is inserted). In the critical phase, dimer-dimer and monomer-monomer correlators decay as power laws of the distance and the decay exponents are and respectively. Going back to the height picture, the insertion of magnetic charges (for instance monomers or links breaking the bipartiteness) corresponds to dislocations and, as we show below, can be treated as perturbations in the coarse-grained field theory.

Prior to the presentation of our strategy to attack the dimer problem, we first introduce the numerical techniques used here and illustrate them on the example of the square lattice interacting dimer model. In this study we use a Transfer Matrix (TM) approach to determine the domains of existence of the different phases of the dimer model, computed on a torus of longitudinal and transverse sizes and respectively. The TM techniques allow to compute exponents of critical dimer-dimer and monomer-monomer correlations and also correlations themselves (see Appendix), on systems with large enough to be considered infinite.

For a critical phase such as the one encountered in the high-T region of the square lattice dimer model, we can compute the decay exponents of correlation functions in two ways: either by a direct inspection of the real-space decay of correlation functions (see the description in Sec..3 of Appendix), or by the leading eigenvalues of the TM (see Sec..4 of Appendix). Both methods are guided by a conformal field theory (CFT) analysis of the critical phase, and the latter is more precise to determine exponents as it allows to use efficiently topological and translation symmetries of the TM (see Sec..4 in Appendix). The calculations of correlation functions in real space offer however the advantage of providing more physical informations (for example about the correlation length) when the system is not critical. Another important quantity, namely the central charge of the CFT, can also be extracted from the finite-size scaling of the TM largest eigenvalue when the system is critical (see Sec..4 of Appendix).

We illustrate the validity of this numerical approach in Fig. 2, where the critical phase is evidenced by the plateau of the central charge (estimated using TM largest eigenvalues for system widths up to ). The scaling dimensions of dimer and monomer operators (obtained from subleading eigenvalues) are also displayed according to the discussion above, the transition to the columnar phase is characterized by a value of the stiffness, thus a value . This criterion is used to estimate the transition temperature (estimates of can also be used in principle, but they are more sensitive to finite-size effects in that temperature range). The results of Fig. 2 lead to an estimate , in good agreement with previous results (12); (13).

## Iii Critical Phase in a Dimer Model Interpolating the Square and Triangular Lattices

In this section, we investigate whether the bipartiteness of the square lattice is a condition for the existence of the critical phase, by building a model defined on a lattice interpolating continuously between the square and the triangular lattices.

Definition of a model with lattice and interaction anisotropies — For commodity, we represent the triangular lattice as a deformed square lattice with bonds in the , and direction (see Fig. 1 and Appendix). The interpolation between both lattices is made by assigning a fugacity to the diagonal bonds. Since the elementary plaquettes of the triangular lattice are more numerous than those of the square, we need to define another parameter characterizing the interaction on diamond-like plaquettes (i.e. non-square four-site plaquettes). We finally keep the notation for the interaction between parallel dimers on a square plaquette. These elementary energy scales in the problem are illustrated in Fig. 1. For and , we have the classical dimer model on the square lattice discussed in Sec. II, while the pure triangular dimer model (that will be discussed in Sec. IV) corresponds to both and . To describe the system in function of these two anisotropy parameters, we take the approach of varying either or while keeping the other parameter constant. The different paths in the space where phase diagrams were computed numerically are represented in Fig. 3.

Description of expected phases and transitions— In the non-interacting limit () of the present model, the system is disordered with exponentially decaying correlations for any fugacity : criticality is destroyed by an infinitesimal proportion of diagonal bonds (11). This can be accounted for in the height field theory by adding a relevant term in the action Eq. (1), which can be understood by noticing that adding a single diagonal dimer on the square lattice corresponds to inserting two monomers on the same sublattice, i.e. two magnetic charges of the same sign (20). The corresponding vertex operator is written in terms of the field , dual to the height field , and leads to the new effective action:

(2) |

In the absence of interactions, this new perturbing field (the term) corresponds to the mass term of the free Majorana doublet arising in the large scale regime within the Pfaffian description of the dimer model (11). For a generic value of the stiffness, the scaling dimensions of the two perturbing terms are and respectively. The fact that the introduction of diagonal bonds is indeed relevant in the non-interacting limit is seen from the value of the stiffness in this case (since ). However, the key point of our study is that the term becomes irrelevant for whereas the locking potential is relevant only for . As we can tune the stiffness constant by modifying the interaction strength between parallel dimers ( and terms), we are hoping to reach the window, i.e. a critical phase on the triangular lattice, within the interacting dimer model defined above. The predicted phase diagram is therefore as follows: first, a high- liquid phase (for in the unperturbed model), then a transition to a critical phase () at intermediate couplings and finally a low- columnar ordered phase (for ). The dimer and monomer scaling exponents are predicted to be both equal to at the transition from the liquid to the critical phase, and and respectively at the critical-columnar transition. In both cases, transitions out of the critical phase are expected to be of Kosterlitz-Thouless type. This analysis holds only in the perturbative regime for the diagonal dimer fugacity , and we now address the question whether this scenario is realized for an arbitrary lattice anisotropy by means of numerical TM calculations.

Perturbation of the square lattice model by diagonal bonds — We first consider the model with no additional interactions and turn on the diagonal dimer fugacity (path A in Fig. 3). Our numerical estimates of the central charge and exponents and are displayed in Fig. 4 as a function of the coupling strength for the specific value . We clearly observe the emergence of a plateau witnessing a critical phase, for a wide range of couplings. The points and where and correspond roughly to the high-T limit of the plateau, as predicted above; from these criteria, the transition between the critical and liquid phase is located at , the error bar being estimated by . Similarly, the low-T end of the plateau at , corresponding to the entrance into the columnar phase, coincides with the criterion , confirming the analysis above. We use estimates of to locate the columnar transition, since as on the square lattice estimates of the monomer exponent are more affected by finite-size effects.

Repeating the same analysis for different values of the diagonal dimer fugacity, we obtain, as a function of and , the phase diagram of Fig. 5. The main features are: (i) the transition temperature to the columnar phase is essentially not affected by the diagonal dimer fugacity and remains close to the value obtained in the square lattice; (ii) as found in Ref. (11), the system with no interaction is a gapped liquid irrespective of the diagonal dimer fugacity; (iii) the dimer model (with ) comprises a critical phase even on the isotropic triangular lattice (), where it extends from up to .

These results are in full agreement with the predictions of the field-theoretical analysis developed above. We emphasize also here that the lattice bipartiteness is not a necessary condition to have criticality in a dimer model, and that this condition should rather be replaced by a condition on the existence of a stiffness constant window in an unperturbed model.

Introduction of a finite interaction anisotropy — We now introduce the interaction term on the diamond-like plaquettes. The and terms are in competition, since they tend to favor different dimer orderings. When , the square columnar configurations will eventually dominate at low enough , but the presence of the frustrating interactions will shift down the transition temperature to the ordered phase. As the algebraic correlations of the critical phase also correspond to a square-like columnar ordering at lower temperatures, we also expect the extent of the critical phase to shrink down as the strength of the interactions is increased. From the point of view of the effective action Eq. (2), we expect to decrease with the frustrating interactions.

For a sufficiently small parameter, the behavior of the system is found to be identical to the case , with a critical phase, evidenced by a plateau between the dimer liquid and the ordered phase. This is illustrated in Fig. 6 for anisotropy parameters ()=(). As expected, the transition towards the ordered phase is slightly shifted towards lower temperatures ( for while for ). The high-T boundary of the critical phase is more strongly affected by the presence of interactions: for the same lattice anisotropy, for to be compared with for .

Using the same procedure as before, the transition points and were determined in the regions of the phase diagram corresponding to paths B,C,D (with respectively and ), and E (for ) and the corresponding phase diagrams are reported in Fig. 7. As anticipated, we clearly see that the effect of adding the interaction is essentially to shift both boundaries of the critical phase towards lower temperatures: in other words, the frustration due to the terms has a net tendency to stabilize the dimer liquid and to destabilize the ordered (columnar) or quasi-ordered (critical) phase. For a given fugacity , the extent of the critical phase indeed decreases as increases. The transition between the dimer liquid and the critical phase turns out to be more affected by than the transition towards the columnar phase. We can interpret this with the following rough argument: interactions effectively increase the weight of diagonal bonds in the dominant configurations, enhancing the perturbation caused by the diagonal bonds and therefore decreasing the extent of the critical phase from the high- direction.

For large values of , a fair determination of phase boundaries becomes more difficult as finite-size effects are getting more important. For instance, the determination of with the dimer and monomer exponents criteria give values increasingly far from each other, resulting in larger error bars in the phase diagrams of Fig. 7c and d. In practice for , both our numerical results and the way of analyzing them break down due to uncontrolled finite size effects. This is actually not a surprise, as the CFT-guided analysis of the TM results (see Sec. .4 in Appendix) relies on the existence of an unperturbed critical window in the model, which is no longer present for large interactions. For , we also find that the second largest TM eigenvalue is no more found in the symmetry sector in the low-T phase. This indicates a change in nature of the lowest-energy excitations and that the scheme along which criticality was previously understood is not valid anymore. From these arguments and even if our finite-size results become less reliable, it is clear that they are no longer compatible with the existence of an extended critical phase for : this is exemplified by the absence of a plateau in the numerical estimate of the central charge. Finally, we note that in parallel to the shrinking of the critical phase observed when increasing , the columnar transition seems to be progressively shifted towards the limit as the interaction anisotropy parameter approaches (see Fig. 7d). This is again in agreement with our previous qualitative arguments on the effect of interactions.

## Iv Dimer ordering on the isotropic triangular lattice

We now turn to the more complex situation of the isotropic triangular model, where and (in this section, the isotropic interaction will only be denoted by ). The high- limit is in this case well understood: the system is in a dimer liquid phase, with exponentially decaying dimer-dimer and monomer-monomer correlation functions. This has been shown exactly for in Ref. (11), and we could confirm it for small values of the interaction strength. For intermediate values, as expected from the results of Sec. III for , we find no evidence for a critical phase in our TM calculations (in particular, no plateau).

Whereas we always found a low-T columnar ordering up to now, the low- behavior of the isotropic triangular model reveals more complex. Indeed, from the isotropy of the interactions and since all dimer fugacities are equal, it appears that configurations that minimize the energy are more numerous. At first glance, columnar ordered ground-state configurations emerge ( for the number of ways to put a dimer on a link connected to a given site and for the possible directions of a dimer column, given an orientation of dimers). However, as first remarked in Ref. (10), the number of configurations that minimize the energy is much larger as one can easily create zero-energy defects by just translating a line of dimers or flipping all dimers along a column. These two types of moves from a given columnar ground-state are represented in Fig. 8. The almost extensive ground-state degeneracy generated by these costless defects clearly changes the picture for the low- behavior of the isotropic triangular dimer model. Indeed, at , strictly speaking, the system is no longer ordered with respect to a local order parameter such as the columnar one defined for the square lattice. However, it is still possible that at finite , the thermal fluctuations select a specific ordering pattern: this would be an illustration of the “order by disorder” effect (21); (22).

In fact, for the quantum dimer model on the triangular lattice, Moessner and Sondhi found in perturbation theory (with respect to the quantum kinetic terms) that indeed (quantum) fluctuations select the ordered columnar states (10). Noticing that in our classical model, thermal fluctuations play the same role as quantum kinetic fluctuations in the quantum dimer model (12) since both essentially count the number of flippable plaquettes in a given state, it is likely that thermal fluctuations (instead of quantum fluctuations) trigger a similar order-by-disorder scenario here. From general arguments (22), we expect that if it is indeed the case, the transition temperature to the columnar phase should be quite low. It could also happen that, for our specific model, this transition never occurs and the systems stays in a liquid phase down to .

We now try to settle this issue with the help of numerical TM calculations. As the critical phase disappeared, we can no longer use the CFT-based analysis of the previous sections and we have to resort to more standard thermodynamical means of detecting the hypothetic transition to a low- ordered phase. We first computed dimer-dimer correlations by TM iterations (see Appendix). In Fig. 9, connected correlation functions are displayed as a function of dimer-dimer distance for the example of a torus for interaction strength up to 5. As exemplified by the log-linear scale, they are clearly short-ranged with an exponential decay (plus an oscillating part depending on the parity of ). We obtain an estimate of the correlation length by fitting the correlations to an exponential decay . The resulting curve , displayed in the inset of Fig. 9, shows a well-pronounced peak at a finite value of , which could correspond to an ordering temperature. We also find that the peak positioned at for the system size , shifts towards higher interaction strength for .

Another thermodynamical insight is given by the specific heat per site , which is also accessible to the TM calculations (see Sec. .2 of Appendix for details). Fig. 10 displays the specific heat as a function of for different system widths . For each value of , the specific heat peaks at a temperature close to the ones of the curves. As increases, the peak in sharpens and shifts to lower temperatures. We discuss below the finite-size scaling of the position of the peak but already note at this stage that the clear narrowing of the specific heat peak with increasing system size, conjugated with the absence of a power-law envelop, is very suggestive of a first-order phase transition.

Having possible signs of a phase transition, we now try to find an order parameter for the low- phase. As already discussed, this a not a simple task because of the zero-energy modes that are responsible for the large ground-state degeneracy. The specific geometry of the TM is of help here to characterize the appearance of long-range order. Consider a very long cylinder with PBC in the small direction. With this geometry, the ground-states having dimers perpendicular to the long direction (horizontal dimers) are much more numerous than ground-states with only vertical or diagonal dimers. Indeed, there are horizontal ground-states (corresponding to the line-shifting modes, see Fig. 8a). The same line-shifting modes provide only a much smaller () degeneracy for vertical or diagonal ground-states. The other family of low-energy modes (column-flipping modes) give a small degeneracy for the horizontal ground-states (see Fig. 8b), and a degeneracy for both vertical and diagonal ground-states. Consequently, in the limit of large , we expect horizontal line-shifting modes to predominate.

If there is long-range order in the system, the previous analysis indicates that the long cylinder geometry induces a preferential ordering in the horizontal dimers (for entropic reasons). Consequently, we consider the average occupation of horizontal bonds on a given line of the lattice as a good indicator of a possible phase transition. At infinite temperature, occupations of all bonds are all likely and equal to . In contrast, from the arguments above, at very low , the probability for a horizontal bond to be occupied on an infinitely long cylinder is . Defining

we expect the “order parameter” to vanish in the low-coupling limit and to saturate to 1 at large enough coupling. Even if, strictly speaking, this above argument is not rigorous (because and are both finite in our computations), we expect to reflect a true physical behavior. is shown in Fig. 11 for lattice widths and we clearly observe the predicted behavior (saturation to at small and to at large couplings). As for the specific heat data, we observe that the shape of is strongly influenced by the system size, with similar trends: as increases, the temperature region of ordering (where is maximal) is getting narrower and is shifted towards lower temperatures.

In order to have a better understanding of this ordering phenomenon, we have computed the order parameter away from the purely isotropic case, for a diagonal dimer fugacity but for values of lower than . The evolution of with is shown on Fig.12 for different interaction anisotropies . As soon as is larger than , the mean occupation of horizontal links is at zero temperature, which is consistent with the fact that, as in the purely square columnar state, only configurations have the minimal energy. Nevertheless, the ordering occurring when increases is characterized by larger values of in an intermediate temperature range. As increases, the maximal value of gets closer to , which is the value found in the ordered phase discussed previously in the isotropic case. Notice that the extent of the temperature range where approaches increases with the length of the system (see curves in the Fig. 12), which supports the scenario of order by disorder: the asymptotic value can be understood by the number of lowest-energy defects (which are lowest-energy configurations in the isotropic case). With interaction anisotropy, these defects are allowed only at finite temperatures, and proliferate down to lower temperatures as decreases. This clearly shows that the ordered phase in the isotropic triangular lattice model differs both qualitatively and quantitatively (by the number of lowest-energy configurations) from that of its anisotropic version, and is a good example of the “order by disorder” scenario.

At this stage, we have three indications (from the correlation length, the specific heat and ) that the isotropic triangular interacting dimer model could order at low temperatures. In all cases, the transition temperature was seen to decrease quite consequently with system size. To check whether the system could order at finite in the thermodynamic limit, we perform a finite-size scaling of these effective transition temperatures. In the insert of Fig. 11, the finite-size temperature transitions and corresponding to respectively the maxima of the specific heat and of , are plotted as a function of inverse transverse system size . A precise finite-size scaling form is difficult to determine (due noticeably to large error bars), but all reasonable finite size dependence (e.g. linear or quadratic in ) lead to a finite value of in the thermodynamic limit. A rough estimate can be made with the help of linear interpolation (dashed line on the figure) and gives (in units of ).

Our numerical results therefore seem to be consistent with a finite temperature ordering of the isotropic lattice model, probably triggered by an order by disorder mechanism (21); (22). The sharp behavior of both specific heat and “order parameter” at the transition (see Fig. 11), as well as the absence of any criticality behavior in the central charge, suggests that this transition is first order. We finish by noting that the samples used in the computations are of relatively moderate size, and that it is still possible that the extrapolated transition temperature actually vanishes with larger samples available. To have a better understanding of this transition with the presently available system sizes, one could also study this model on samples with other geometries and adapted order parameters, and check if extrapolations of finite-size temperatures on different geometries give the same in the thermodynamic limit. Another insight could be given by Monte Carlo simulations for this model. We expect however these simulations to be difficult as the presence of the low-energy modes (and the corresponding large degeneracy of the ground-state) will certainly induce ergodicity and freezing problems in the Monte Carlo process (23). Such investigations are beyond the scope of the present paper.

## V Conclusions

To summarize, we have constructed a simple classical interacting dimer model on a lattice that interpolates between square and triangular lattices. This is of particular interest since the two limiting models singularly exhibit very different behaviors at infinite temperature (critical and short-ranged phases for the square and triangular lattices respectively). Since the topology of the triangular lattice can be simply obtained from the square lattice by adding one extra diagonal bond on each square plaquette, we have introduced a fugacity parameter for these extra bonds. Similarly, the extra bonds leading to new local interactions in diamond-like plaquettes, we have also considered dimer interactions on these plaquettes, differing from those on the square plaquettes.

This anisotropic dimer model has been investigated in great details using simple considerations as well as numerical transfer matrix techniques (on strips with up to 14 sites wide) supplemented by predictions from conformal field theory. A very rich phase diagram has been obtained with, in particular, a novel intermediate behavior where the critical phase is now restricted to a finite intermediate temperature range and does not extend, as for the pure square lattice model, up to infinite temperature. At high-temperature, the critical phase is hence replaced by a liquid dimer phase as in the isotropic triangular lattice. It is of interest that such a behavior appears immediately for arbitrary small fugacity of diagonal bonds so that one can view the typical behavior of the square lattice as a “singular” limit. However, we note that the critical phase itself survives in an extended vicinity of the square lattice model (although its extension in temperature jumps abruptly) and this stability is therefore not related to the bipartiteness of the lattice. While the critical phase disappears progressively with increasing interaction anisotropy, the ordered phase survives for any range of anisotropy, as well as the high-T liquid dimer. Criticality in this model hence requires to have sufficiently anisotropic interactions. These findings are reproduced in the schematic phase diagram of Fig. 13 for the case of isotropic dimer fugacity . Lastly, we have devoted special attention to the limiting case of the isotropic dimer model on the triangular lattice where we found evidence of a direct (likely first-order) transition between the high-T liquid phase and a low-T ordered phase triggered by an order-by-disorder mechanism.

###### Acknowledgements.

We thank F. Becca, S. Capponi and C. Castelnovo for fruitful discussions, the Agence Nationale de la Recherche (France) for support and IDRIS (Orsay, France) for computer time.## Transfer matrix analysis of the anisotropic dimer model

Other works (12); (13) already discussed the construction of the TM for an interacting dimer model. We give a few technical details related to the triangular lattice and to the observables studied in the body of the text.

### .1 1D and 2D dimer configurations

The TM allows to treat a 2D system with periodicity in the transverse direction () as a succession of 1D rows of lengths . In a row at abscissa (the axis being the direction of propagation of the TM), a dimer configuration is specified by the dimer occupation of horizontal (oriented along ) links at abscissa and of vertical and diagonal links connecting sites at abscissa to sites at , at relative positions and respectively. The number of such configurations, which is the size of the TM, is (This number can be evaluated as where is a one-dimensionnal transfer matrix defined between the dimer configurations allowed on a triangle, taking into account the hard-core constraint).

The topology of the triangular lattice is such that on a -wide cylinder, the parity of the number of vertical or diagonal dimers linking two successive rows is conserved along the axis (assuming that is even). This defines two topological sectors (even and odd), and the TM is block-diagonal with corresponding blocks and respectively. We also use the invariance by translation along the axis to reduce the size of T. Translation invariance allows to consider, instead of 1D configurations, only their representatives in a symmetry sector (defined by the transverse wave vector ), from which other configurations are obtained by translations along . Note that the TM elements take into account not only representatives but also all configurations connected to a given representative. These two invariances allow to reduce considerably the TM size in various symmetry and topological sectors. We give in Table 1 the size of the different sectors useful in the following for up to . Values of non-zero TM elements, i.e. between compatible 1D-configurations, depend on the number of doubly occupied plaquettes of both types and of diagonal links occupied in these configurations.

### .2 Partition function and correlation functions

With the transfer matrix T constructed, we have access to statistical properties of the model on a torus by means of the partition function, which reads

Internal energy and specific heat — From the TM elements, one also gets the internal energy of the system:

We computed the specific heat per site by a numerical differentiation of with respect to the inverse of the coupling strength. For the calculations of Sec. IV, was computed exactly for -wide systems (with so that does not depend on ). For wider systems , was estimated by searching the leading TM eigenvector and evaluating from it the energy per row (the consistency of this method, valid in the limit of infinitely long systems, was checked for first).

Correlation functions — Correlation functions are computed by means of TM iterations. Dimer-dimer correlations, for two dimers with a specified orientation and relative position r, are obtained as the ratio over the partition function of a modified partition function taking into account only the 2D-configurations where these two dimers are present. Mathematically, a projector P onto the corresponding 1D configurations is inserted in the matrix product in at positions corresponding to both dimers:

where is the projection onto the axis of the relative position . Even if we do not discuss them in the body of the text, monomer-monomer correlations are also accessible via TM calculations. In the TM language, the insertion of a monomer in a given row requires to replace the TM by a modified matrix T’ which forbids the occupation of links emerging from this site (this involves a change of topological sector). For both dimer-dimer and monomer-monomer correlations, we checked that in the non-interacting case, the exact results of Ref. (11) are recovered.

### .3 Analysis of real-space correlations

When the system is in a critical phase, conformal transformation techniques (24) can be used to analyze the real-space decay of correlation functions on a -wide infinite cylinder. Consider the correlation of two dimers (or monomers) located at a relative distance along the cylinder axis. If the exponent of the infinite-plane power-law correlations is , correlations on the cylinder scale as:

which can be approximated, for distances much larger than , by an exponential decay with a correlation length . We checked our calculations on the non-interacting case on the square lattice for a cylinder. We found a correlation length , close to the expected (thermodynamical limit) value . The small discrepancy can be attributed to the correction present in the correlation function (14). Repeating the same calculations for monomer correlations, we extracted from the exponential decay a critical exponent , in good agreement (for such a small value) with the exact value .

### .4 CFT analysis of largest TM eigenvalues

The CFT analysis is also useful to compute the central charge , as well as the scaling dimension of the dimer and monomer operators, from the size dependence of the leading TM eigenvalues (12). The central charge is estimated by the -dependence of the largest eigenvalue of T, which is directly connected to the free energy per site , in the limit of an infinitely long system:

In practice, we add a term to the previous fitting expression (13), which is justified by symmetry reasons and improves the agreement between raw data and the fitted expression of . Estimating thus requires to perform a fit of with at least 3 values of . The size of the TM is a limitation to the number of sizes accessible; nevertheless, one can use the already mentioned topological and translation invariance to adress larger sizes. Indeed, one can assert by symmetry reasons and check numerically that the leading eigenvalue of the TM is found in the even sector. The scaling dimensions of dimer and monomer exponents are determined from the largest eigenvalues in other symmetry sectors (12):

This allows to determine and (with ) from the power method (25), requiring less memory than for a full diagonalization. To minimize finite-size effects on and , we use different estimates for these quantities. Each estimate is the result of a fit with either (for ) or (for ) consecutive values of . In this way one obtains estimates and (for ). The finite-size effects on these estimates become smaller when increasing ; in order to obtain more reliable values we make extrapolations of these estimates, assuming that they scale as (13):

Such estimates and their extrapolations are displayed in Fig. 2, 4 and 6. In general, the extrapolations were done with estimates and , while for extrapolations we used the estimates , and . The difference between the extrapolation and the closest estimate (the one with the largest ) gives the order of magnitude of the error.

### References

- D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 2376 (1988).
- M.E. Fisher, J. Math. Phys. 7, 1776 (1966).
- R. Moessner, S.L. Sondhi, Phys. Rev. B 68, 054405 (2003).
- F. Vernay, A. Ralko, F. Becca and F. Mila, Phys. Rev. B74, 054402 (2006).
- M. E. Zhitomirsky, Phys. Rev. B71, 214413 (2005).
- D.L. Bergman, R. Shindou, G.A. Fiete, L. Balents, Phys. Rev. Lett. 96, 097207 (2006); D.L. Bergman, G.A. Fiete, L. Balents, Phys. Rev. B 73, 134402 (2006).
- A. Sen, K. Damle and T. Senthil, preprint cond-mat/0701476.
- D. Poilblanc, K. Penc and N. Shannon, preprint cond-mat/0702367.
- P.W. Kasteleyn, Physica 27, 1209 (1961); H.N.V. Temperley and M.E. Fisher, Phil. Mag. 6, 1061 (1961); M.E. Fisher, Phys. Rev. 124, 1664 (1961).
- R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001).
- P. Fendley, R. Moessner, S.L. Sondhi, Phys. Rev. B 66, 214513 (2002).
- F. Alet et al., Phys. Rev. Lett. 94, 235702 (2005); Phys. Rev. E 74, 041124 (2006).
- C. Castelnovo, C. Chamon, C. Mudry, and P. Pujol, Ann. Phys. 322/4, 903 (2007); see also D. Poilblanc et al., Phys. Rev. B74, 014437 (2006).
- M. Fisher, J. Stephenson, Phys. Rev. 132,1411 (1963).
- H.W.J. Blöte and H.J. Hilhorst, J. Phys. A 15, L631 (1982).
- J. Kondev and C.L. Henley, Phys. Rev. B 52, 6628 (1995); Nucl. Phys. B 464, 540 (1996); C. Zeng and C.L. Henley, Phys. Rev. B 55, 14935 (1997); R. Raghavan, C.L. Henley and S. L. Arouh, J. Stat. Phys. 86, 517 (1997).
- C.L. Henley, J. Stat. Phys. 89, 483 (1997).
- E. Fradkin et al., Phys. Rev. B 69, 224415 (2004).
- B. Nienhuis, in Phase transitions and Critical Phenomena, edited by C. Domb and J. L. Lebowitz (Academic, London, 1987), Vol. 11.
- C. Castelnovo, C. Chamon, C. Mudry, and P. Pujol, Phys. Rev. B 73, 144411 (2006).
- J. Villain, R. Bidaux, J.-P. Carton and R. Conte, J. Phys. 41 1263, (1980)
- E.F. Shender, P.C.W. Holdsworth in Fluctuations and order: a new synthesis, edited by M.M. Millonas (Springer-Verlag, 1996)
- C. Castelnovo, P. Pujol and C. Chamon, Phys. Rev. B 69, 104529 (2004).
- J. Cardy, in Fields, strings and critical phenomena, edited by E. Brezin and J. Zinn-Justin (Les Houches, 1988).
- J.H. Wilkinson, in The algebraic eigenvalue problem (Oxford University Press, 1988).