Arxiv:1008.0868V2 [Cond-Mat.Stat-Mech] 15 Feb 2011 N2 Ytm Ihcniuu Symmetry

Home , Universality class

arXiv:1008.0868v2 [cond-mat.stat-mech] 15 Feb 2011 n2 ytm ihcniuu symmetry. occur continuous describe which with systems phases to 2d ordered is in unnoticed, purpose previously our new, funda- universality, the between and fractional from connection symmetry Apart as the regarding such 9]. considerations [1, defects, mental walls understood topological domain be and new can vortices of transitions transition. terms models, phase these in KT universal- of usual Potts the some to and In addition phase Ising in with the classes, ity diagrams, to phase belonging complex transitions very have metry, oeswih hl rsrigtesame the XY preserving generalized while of which, class large models a finite analysis, and (FSS) scaling simulations size Carlo Monte extensive using study t ls.I s hrfr,itrsigt s htphase what ask un- transformation invariant the to Hamiltonians der 2d interesting for therefore, possible are is, universal- transitions It same the class. to ity belong Landau-Ginzburg- to sym- of not coarse-grained underlying Hamiltonian Violations same same Wilson the the Thus, with and [8]. systems 2d. metries for in one can possible common order transition is more it phase even first are order a universality second strong by a preempted 3d [4– in be transitions even phase and to possible sufficient of 7], class not universality are the considerations fix symmetry on purely decay functions correlation the and exponentially. 3], unbind [2, pairs temperature these critical (KT) Kosterlitz-Thouless Above the functions. alge- correlation in spin-spin resulting decaying the pairs, braically vortex-antivortex vortices, temperature, in valued joined low integer all of At are form of the model. in example XY defects, classical topological the A is system loss order. such a pseudo-long-range and the functions to of correlation leading above unbind, decaying while defects exponentially up, these temperature paired the tem- critical all undermine low the are which At order defects pseudo-long-range topological possible. the remains peratures still a functions systems, cor- decaying these relation algebraically of with spin order order as pseudo-long-range long-range walls. such the domain destroy excitations, by waves, and Goldstone controlled vortices massless as is Although such symmetries defects, continuous topological with systems h oe osdrdhr a itr fferromag- of mixture a has here considered model The niei d o dssesteagmnsbased arguments the systems 2d for 3d, in Unlike (2d) dimensional two of behavior temperature low The nttt eFıia nvriaeFdrld i rnedo Grande Rio do Federal F´ısica, Universidade de Instituto h e hs rniin eogt h dPts sn,o KT or Ising, Potts, 2d the to belong ferrom transitions the phase between different new competition an qualitatively The the simulations a from to Carlo originating lead phases Monte harmonics extensive order I higher of the the in combination transition a phase additional using Intro an to class. leads universality Hamiltonian (KT) Kosterlitz-Thouless the to ti elkonta h dX oe xiisa nsa infini unusual an exhibits model XY 2d the that known well is It e ree hssi ls fgnrlzdX models XY generalized of class a in phases ordered New θ → θ ´ boC oeoo eesnJ rno,adYnLevin Yan and Arenzon, J. Jeferson Poderoso, Fábio C. +2 π nti ae,w will we paper, this In . θ → θ Dtd coe 5 2018) 25, October (Dated: +2 π sym- hroyaiso h oe ecie yE.()for (1) Eq. by described model 0 the the concerns of question pseudo-long-range thermodynamics interesting a have An a temperature therefore, tem- low order. will, be The nematic-like 0 critical also = model. same will ∆ XY for there the pure phase the case exactly in this as at in perature transition that phase show to KT easy is It re.I ∆ If integer ferromagnetic of pseudo-long-range order. pairs a in in tempera- bound resulting low vorticity, are At vortices half-integer across half-integer aligned. [1], tures, anti-paralelly case: walls) are 1 (domain spins = strings which ∆ by new the connected which has in vortices in model presence present states The the not metastable and excitations orientation. to antiparallel 9–13] leads have [1, term spins authors second of the number of a by ied nrdcsas iycia hs bv h Tln,but line, KT the frustration above geometric phase triangular chiral the tiny latter, and a the also square introduces On the 12]. been [11, on be has lattices to behavior simulations This found with [1]. is verified class — universality Ising ferromagnetic ordered the and in pseudo-long-range nematic two the — Surprisingly, phases the [14]. be- between class phases universality transition paramagnetic KT the the transition to and the longs disor- nematic expected, completely the As a between enters phase. in- system paramagnetic further the dered On and vortex-antivortex unbind half-integer phase. pairs the nematic temperature, a of into crease melts and vanishes system vortices the half-integer between tension string aeapaedffrneo 2 to of spins (XY difference adjacent phase the coupling favors a ferromagnetic one have second usual the the while is model), term first The sas nain ne h transformation which the Hamiltonian, under nematic-like invariant purely also a is have we 0, = ∆ ,adi nain ne rotations under invariant when recovered is and model, (1) 1, XY Hamiltonian = pure the ∆ the ∆, of of symmetry value the has the of Independent ger. hr h u soe h ers egbr pn na on spins neighbors nearest 0 the over lattice, is square sum the where ei n eai-ieinteractions, nematic-like and netic = H < For u,C 55,95190ProAer S Brazil RS, Alegre Porto 91501-970 15051, CP Sul, ∆ uto fanmtcculn noteXY the into coupling nematic a of duction − q < iguieslt ls 1.I hspaper, this In [1]. class universality sing hs iga,wt diinlordered additional with diagram, phase gei n suonmtccouplings. pseudo-nematic and agnetic X h ,teHmloin q 1,hsbe stud- been has (1), eq. Hamiltonian, the 2, = ,weebt em compete. terms both where 1, ij i nt iesaig eso that show we scaling, size finite d nvraiyclasses. universality ∆cos( [∆ eodrpaetasto belonging transition phase order te < ∆ ≤ mc θ i ∆ stetmeauei asd the raised, is temperature the as , − ≤ θ j (1 + ) ,and 1, kπ/q − where , q )cos( ∆) sapstv integer. positive a is θ j k qθ θ → ≤ j i → − q θ j sa inte- an is qθ θ +2 j j +2 )] π , For . π/q (1) . 2

1 1

0.8 0.8

0.6 0.6 N UF T T 0.4 0.4 N UF F2 0.2 0.2

F1 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 ∆ ∆ (a) q = 3 (b) q = 8

FIG. 1: a) Schematic phase diagram for q = 3. At both extremes ∆ = 0 and 1, the order-disorder transition temperature is Tc = 0.893. There are two low temperature phases with a pseudo-long range order: generalized-nematic (three preferred spin orientations) and ferromagnetic (broken reflection symmetry). Notice the multicritical point around ∆mc ≃ 0.4. For ∆ < ∆mc, there are two transitions: paramagnetic to nematic, in the KT universality class; and nematic to ferromagnetic, in the q = 3 Potts universality class. b) Schematic phase diagram for q = 8. Besides the paramagnetic, the usual ferromagnetic, and the nematic phases, there are two new ferromagnetic phases (F1 and F2) in which spins have half-plane preferred orientations. otherwise the phase diagram retains the same topology as 0.3 on the square lattice [12]. A related model, with a similar 1.5 phase diagram, was also studied in Ref. [15]. An interesting question is whether, for q> 2, the topology of the 1 β/ν L phase diagram remains unchanged. To answer this, we 0.2 1 have explored using extensive Monte Carlo simulations m 0.5 γ/ν - L the equilibrium phase diagram of the Hamiltonian (1) 1 χ 0 with q ranging from 2 to 10. We find that for q =2, 3, 4 -4 -2 0 2 4 0.1 1/ν the topology of the phase diagram remains the same as tL L = 45 for q = 2, however, the transition between the pseudo- L = 60 L = 75 ferromagnetic and pseudo-nematic phases belongs either L = 90 0 to the 3-state Potts (q = 3) or to the 2d Ising (q =2, 4) -5 0 5 10 15 1/ν universality class, see Fig. 1a. For q ≥ 5, the topology of tL the phase diagram changes completely and new phases with a pseudo-long range order come into existence. FIG. 2: Ferromagnetic-nematic transition for q = 3. FSS analysis of the susceptibility χ1 and the order parameter m1 The simulations were performed on a square lat- (inset) for ∆ = 0.25 and different sizes L. In both cases the tice of linear size L and periodic boundary condi- data collapse is excellent using the 3-state Potts exponents. tions. Both Metropolis single-flip and the Wolff algo- The critical temperature at this point is Tc ≃ 0.365. rithm [16] were used. In accordance with the symmetry of the Hamiltonian, the possible order parame- −2 ters are mk = L | exp(ikθi)| where k = 1,...,q. Pi ized susceptibility scale with the size of the system as The corresponding generalized susceptibilities are χk = 2 2 2 m ∝ L−β/ν and χ ∝ Lγ/ν. Exactly at the transition, βL (hm i − hmki ) and the Binder cumulants are Uk = 4k 2 2 β/ν =1/8 and γ/ν = 7/4, which are the same ratios as 1 − hmki/3hmki [17, 18]. If the transition is not KT, the usual FSS can be used to get the critical exponents for the 2d Ising model. However, what distinguishes the β, γ, and ν: m = L−β/νf(tL1/ν ) and χ = Lγ/νg(tL1/ν ), KT transition from the Ising one, is the behavior of the where m is the order parameter and χ its susceptibil- order parameter and the susceptibility in the low temperature phase where they also exhibit FSS, but with ity, f and g are the scaling functions, and t = T/Tc − 1 is the reduced temperature. This FSS, however, is not non-universal critical exponents. Recall that for normal valid at the KT transition for which all of the low tem- second order phase transition, FSS exists only at the crit- perature phase is critical and the correlation length and ical point. This difference can be used to distinguish the the susceptibility are infinite [3, 19]. Nevertheless, it is KT transition from the Ising one. possible to show that at the KT transition and in the Figure 1a shows the schematic phase diagram for q = 3, low-temperature phase, the critical exponent ratios are which is topologically identical to the q = 2 case. At well defined and the order parameter and the general- high temperatures, the equilibrium state is the disordered 3

∆ = 0

0.4 1 ) ∆ = 0.25 θ β/ν L 4 P( 0.3 m 0.5 γ/ν - ∆ = 0.7 L 4

χ 0 0.2 -10 0 10 20 30 1/ν tL L = 96 0.1 ∆ = 1 0.1 L = 112 L = 128PSfrag replacements L = 144 0 0 -20 -10 0 10 20 30 40 50 60 70 0 π/2 π 3π/2 2π 1/ν tL θ

FIG. 3: Data collapse for the susceptibility χ4 and the order parameter m4 (inset) for q = 8 at ∆ = 0.15 as the sys- FIG. 4: Angle distribution for the four phases along the ﬁxed tem crosses the border between F1 and F2 phases. The col- temperature line T = 0.16 of the q = 8 phase diagram. All lapse was obtained using the critical exponents of the 2d Ising graphs have the same vertical scale. model, γ = 7/4, ν = 1, β = 1/8 and Tc ≃ 0.058.

congruent discrete spin orientations, separated by π/4. paramagnetic (P). As the temperature is lowered, the system enters either the usual ferromagnetic (UF) or the generalized-nematic phase (N), depending on the value of ∆. Both of these order-disorder transitions belong to 1 0.1 the KT universality class. Up to the multicritical point β/ν located at ∆mc ≃ 0.4, there is a line of critical points L 0.5 1 separating the generalized-nematic from the UF phase. m γ/ν The order parameter m1 is used to distinguish between - L

1 0 the generalized-nematic and the ferromagnetic phases: χ 0.05 0 5 10 15 1/ν m1 ≃ 0 in the nematic phase and is ≈ 1 in the ferromag- tL L = 56 netic phase. Using FSS, and also the Binder cumulant, L = 64 we find that the critical points along this line belong to L = 80 L = 96 the 3-state Potts universality class. Fig. 2 shows the data 0 1 -5 0 5 10 15 20 25 30 35 collapse of m (inset) and the corresponding susceptibil- 1/ν ity χ1, for a critical point with ∆ = 0.25. The collapse tL is excellent using the critical exponents of the 2d, 3-state FIG. 5: Data collapse for the susceptibility χ1 and the order Potts model [20]: β =1/9, γ = 13/9 and ν =5/6. parameter m1 (inset) for q = 8 at ∆ = 0.35 as the system The same topology of the phase diagram persists for crosses the border between the generalized-nematic and F2 q = 4, except that the transition from the generalized- phases at Tc ≃ 0.34. The collapse was obtained using the 2d nematic to the ferromagnetic phase is once again in the Ising model critical exponents, γ = 7/4, ν = 1 and β = 1/8. universality class of the Ising model [21, 22]. For q ≥ 5, however, the topology changes dramatically. Two new Besides the order-disorder KT transitions, several new phases emerge from the T = 0, ∆ = 0 fixed point, see order-order transition lines appear in the q ≥ 5 phase di- Fig. 1b. The new phases are pseudo-ferromagnetic and agrams, Fig. 1b. The transition from the F1 to F2 phase have a broken reflection symmetry. We shall denote these is well described by the m4 order parameter, which is phases as F1 and F2. The low-temperature phase F1 has close to zero in the later phase. The data collapse for only one preferred spin orientation, while in the phase F2 several lattice sizes is shown in Fig. 3 along with the there are four preferred spin orientations with different data collapse for the susceptibility χ4. Both collapses weights, see Figs. 1b and 4. Fig. 4 exhibits histograms of were obtained using the set of critical exponents of 2d spin orientation in the low temperature phases along the Ising model: γ = 7/4, ν = 1 and β = 1/8. The tran- line T =0.16 together with a pictorial representation of sition from F2 to nematic is also continuous and is also the possible spin orientations for the q = 8 model. All in 2d Ising universality class, as can be seen in Fig. 5 in the ferromagnetic phases have a quasi-long-range order which the data collapse of the order parameter and the and a broken reflection symmetry. For F2, the distribu- corresponding susceptibility are shown. For this transition function has four significant peaks, while for F1 there tion, the relevant order parameter is m1. Finally, the is only a single narrow peak. The UF phase of the XY transition from F1 to UF and the transition from F2 to model has a broad continuous distribution of spin ori- UF are both well described by the m8 order parameter. entations. The nematic phase is characterized by eight The two transitions are found to belong to the KT uni- 4

Ising model, respectively. For q = 5, the topology of the -0.18 phase diagram changes dramatically and two new ferromagnetic phases appear. After this, up to q = 10, the maximum value explored in this work, the topology of

k -0.22 β/ν = 0.088 the phase diagram remains unchanged. It is very curi- m ous that although the systems studied here are invariant 10 β/ν = 0.122 under the continuous global symmetry θi → θi + α (for log -0.26 all spins simultaneously), for arbitrary α, the sequence of ∆ = 0.6, T = 0.32 ∆ = 0.6, T = 0.36 the phase transitions between the pseudo-ordered phases ∆ = 0.85, T = 0.16 follow the one observed for the discrete clock models: PSfrag replacements ∆ = 1, T = 0.893 -0.30 q=2 Ising; q=3 Potts; q=4 Ising; for q=5 a bifurcation 1.8 1.9 2 2.1 2.2 2.3 2.4 and a new phase transition appears [23]. If this analogy log10L persists, we expect that there should be a critical value of q above which the phase transition between F1 and F2 becomes KT [24]. This will be explored in a future FIG. 6: Scaling of the order parameter at the KT phase work. The present paper also shows a significant lack boundaries (three bottommost lines) and inside the critical of universality of 2d systems: while all the Hamiltonians phase (topmost line) for q = 8. Notice that along the KT studied in this paper have the same underlying symmetry, transition the exponent β/ν has the Ising value (1/8), while θ → θ +2π, the transitions between the different phases inside the critical phase the exponent is non universal. belong to a variety of universality classes. Furthermore, since the Hamiltonians discussed in this paper can be versality class with the critical exponent ratios, γ/ν = thought of as the leading orders in a Fourier expansion 7/4 and β/ν = 1/8. As an example, Fig. 6 shows that of a general microscopic spin-spin interaction V (θi − θj ), the work raises a troubling question: How much can we F2 is critical with respect to m8 order parameter and has a non-universal FSS characteristic of a low-temperature really deduce about the thermodynamics of 2d systems KT phase, top-most curve of Fig. 6. from the form of their coarse-grained Landau-Ginzburg- To conclude, we have studied, using extensive numer- Wilson (LGW) Hamiltonian? It is clear that the symme- ical simulations, phase diagrams of a class of general- try arguments alone are not sufficient to determine the ized XY models described by eq. (1). Previous results phase diagram of these systems, and one needs to have a show that for q = 2 besides the usual KT transition, detailed knowledge of the microscopic interactions [7]. there is also a nematic to ferromagnetic transition in the Acknowledgments. We thank D. Stariolo for discus- Ising universality class. We find that for q = 3 and 4 sions and the Brazilian agencies CNPq and Fapergs for the topology of the phase diagram remains unchanged, partial support. JJA acknowledges support from the but the ferromagnetic-generalized-nematic transition be- INCT-SC while YL acknowledges support from INCT- longs to the universality class of 3-state Potts and the FCx, and US-AFOSR/FA9550-09-1-0283.

[1] D. H. Lee and G. Grinstein, Phys. Rev. Lett. 55, 541 (2010). (1985). [14] D. G. Barci and D. A. Stariolo, Phys. Rev. Lett. 98, [2] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 200604 (2007). (1973). [15] J. Geng and J. V. Selinger, Phys. Rev. E 80, 011707 [3] J. M. Kosterlitz, J. Phys. C 7, 1046 (1974). (2009). [4] R. J. Baxter, Phys. Rev. Lett. 26, 832 (1971). [16] U. Wolﬀ, Phys. Rev. Lett. 62, 361 (1989). [5] E. Domany, M. Schick, and R. H. Swendsen, Phys. Rev. [17] D. Loison, J. Phys.: Condens. Matter 11, L401 (1999). Lett. 52, 1535 (1984). [18] M. Hasenbusch, J. Stat. Mech. p. P08003 (2008). [6] R. L. C. Vink, Phys. Rev. Lett. 98, 217801 (2007). [19] J. Tobochnik and G. V. Chester, Phys. Rev. B 20, 3761 [7] Y. Levin, Phys. Rev. Lett. 99, 228903 (2007). (1979). [8] A. C. D. van Enter and S. B. Shlosman, Phys. Rev. Lett. [20] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 89, 285702 (2002). [21] D. D. Betts, Can. J. Phys. 42, 1564 (1964). [9] S. E. Korshunov, JETP 41, 263 (1985). [22] S. K. Baek and P. Minnhagen, Phys. Rev. E 82, 031102 [10] T. J. Sluckin and T. Ziman, J. Phys. France 49, 567 (2010). (1988). [23] S. Elitzur, R. B. Pearson, and J. Shigemitsu, Phys. Rev. [11] D. B. Carpenter and J. T. Chalker, J. Phys.: Condens. D 19, 3698 (1979). Matter 1, 4907 (1989). [24] C. M. Lapilli, P. Pfeifer, and C. Wexler, Phys. Rev. Lett. [12] J.-H. Park, S. Onoda, N. Nagaosa, and J. H. Han, Phys. 96, 140603 (2006). Rev. Lett. 101, 167202 (2008). [13] L. Komendov´aand R. Hlubina, Phys. Rev. B 81, 012505