Coarsening Dynamics of an Antiferromagnetic Xy model on the Kagome Lattice: Breakdown of the Critical Dynamic Scaling
Abstract
We find a breakdown of the critical dynamic scaling in the coarsening dynamics of an antiferromagnetic XY model on the kagome lattice when the system is quenched from disordered states into the KosterlitzThouless (KT) phases at low temperatures. There exist multiple growing length scales: the length scales of the average separation between fractional vortices are found to be not proportional to the length scales of the quasiordered domains. They are instead related through a nontrivial powerlaw relation. The length scale of the quasiordered domains (as determined from optimal collapse of the correlation functions for the order parameter ) does not follow a simple power law growth but exhibits an anomalous growth with timedependent effective growth exponent. The breakdown of the critical dynamic scaling is accompanied by unusual relaxation dynamics in the decay of fractional () vortices, where the decay of the vortex numbers is characterized by an exponential function of logarithmic powers in time.
pacs:
64.70.qj, 64.60.Ht, 75.10.HkJuly 4, 2019
I Introduction
Thermodynamic systems quenched from a hightemperature disordered phase into a lowtemperature ordered phase exhibit characteristic growth of the length scale of the ordered domains which in typical situations can be represented as a power law , where the growth exponent depends on the dimension of the space and of the relevant order parameter as well as on the conserved or nonconserved nature of the order parameter.review ()
Phase ordering dynamics (or coarsening process) is usually accompanied by the annihilation and decay of the characteristic topological defects such as point vortices or domain walls which are generated in the initial disordered states. One of the most important notions in understanding and analyzing these coarsening processes is the socalled dynamic scaling hypothesisreview () for the equaltime spatial correlation function of the order parameter, which is closely related to the observed selfsimilarity of the coarsening systems at different time instants.
When the low temperature phase of the system is characterized by a quasilongrange order with power law decay of the spatial correlation function of the order parameter at equilibrium, then the dynamic scaling is generalized to the critical dynamic scaling. One of the well known examples is the ferromagnetic XY modelberezin72 (); kt73 () on the square lattice.
However, it should be noted that this (critical) dynamic scaling hypothesis for the phase ordering (or quasiordering) dynamics has not been proved on some general theoretical basis. Therefore, it is not clearly known what is the condition for the validity of the dynamic scaling which is usually assumed in analyses of experimental results or numerical simulations on coarsening dynamics.
In typical situations, coarsening dynamics is investigated on systems where the phase ordering is accompanied by the breaking of discrete or continuous global symmetry. However, there exist model systems that exhibit infinite ground state degeneracies of a discrete nature in addition to the usual (global) continuous symmetry. Prominent examples are geometrically frustrated spin systems such as antiferromagnetic Ising models on a triangular latticeuhkim2003 (); uhkim2007 (), antiferromagnetic XY or Heisenberg models on kagomesyozi51 () or pyrochlore lattices.moessner2006 () These systems exhibit interesting equilibrium and nonequilibrium behavior due to the geometric frustration effect, including spinglasslike relaxation dynamics without quenched disorder.
In this work we perform numerical simulations on the coarsening of an antiferromagnetic XY model on the kagome lattice (KAFXY model)huse92 (); rzchow97 (); cherepa2001 (); korshu2002 (); kolah2002 () which is one of the simplest geometrically frustrated models with infinite ground state degeneracies. Experimentally, this model can be realized in superconducting Josephsonjunction arrays or superconducting wire networkshiggins2000 (); park2001 (); xiao2002 () on a kagome lattice when a perpendicular magnetic field of half flux quantum (per plaquette) is applied on the system. We can also find physical examples in the anisotropic limit of Heisenberg antiferromagnets on the kagome lattice. It is well known that the system exhibits an infinite ground state degeneracy with finite entropybaxter70 (). The system also exhibits a finite temperature KT transitionkt73 () corresponding to the unbinding of socalled fractional vortices. That is, at low temperature below the KT transition, equilibrium of the system will be characterized by quasilongrange order of the order parameter . Analogous to the case of simple ferromagnetic XY model on a square lattice, we might expect that the coarsening dynamics of the system would exhibit a critical dynamic scaling for the equaltime spatial correlation of the order parameter.
Our simulations, however, show that the critical dynamic scaling is not obeyed very well at least for the time duration of our numerical simulations ranging several decades of time scale. That is, it was not possible to achieve a good scaling collapse for the equaltime spatial correlation functions of the order parameter. Another signature of the breakdown of the critical dynamic scaling is that the length scales corresponding to the average separation between fractional vortices are not proportional to the length scales of the growing quasiordered domains, which are instead related through a nontrivial powerlaw relation. This means that the two length scales exhibit different growth behavior in time. In terms of the decay of the fractional vortices, we also found that the fractional vortices residing in the small triangular plaquettes exhibit faster decay, while, in contrast, those vortices sitting on the larger hexagonal plaquettes exhibit much slower decay that can be fitted by an exponential of logarithmic powers in timecherepa94 ().
In addition to these features of multiple length scales in the coarsening dynamics, the time dependence of the length scale of the quasiordered domains does not exhibit a simple power law growth but rather exhibits an anomalous growth with timedependent effective growth exponent. This appears to be closely related to the unusual relaxation dynamics in the nonpowerlaw decay of the fractional () vortices. It is not clear yet whether the scaling may be restored in the asymptotic limit.
Ii The model system and simulation methods
The kagome lattice consists of cornersharing triangles (Fig. 1). In an antiferromagnetic XY model on the kagome lattice, the Hamilitonian is defined as
(1) 
where , the sum runs over all nearest neighbor pairs of sites, and denotes the angle of the planar spin at site with respect to some fixed direction in the two dimensional spin space. indicates all pairs of nearest neighbor sites in the kagome lattice.
It is easy to see that the ground states of this system have the property that for all pairs of nearest neighbors and , the angle difference satisfies which means that the sum of the three spins vectors on any unit triangle vanishes. Thus the space of the ground states of KAFXY model are equivalent (up to a global rotation) to the ground states of the three state Potts model on the kagome lattice. In Figs. 1(a)(b) are shown two examples of simple ground states with longrange order, the socalled state and state, where refers to the wave vector corresponding to the periodicity of the chirality configuration. It is easy to see, however, that in addition to these ground states with simple spatial order there also exist infinitely many ground states with no spatial order. It is well known that the system has a ground state entropy of per sitebaxter70 (). Now, if we consider the angle variable and the corresponding complex order parameter , the degenerate ground states are all completely ordered in terms of this new order parameter. And it has been shown that the system undergoes a KT transition at a finite temperature, where the spatial correlation of the order parameter exhibits an algebraic decay below the transition temperature.
As for the KT transition temperature of this system, an analytic approximation was given as cherepa2001 (); korshu2002 ()
(2) 
Numerical simulations were performed by Rzchowskirzchow97 () where he found two slightly different estimates on the transition temperature, i.e., one estimate of based on the Binder cumulants of the order parameter and another of from the helicity modulus (or correspondingly the decay exponent of the spatial correlation of the order parameter).
In the present work, the coarsening dynamics of the model system is performed via kinetic Monte Carlo methods with standard Metropolis algorithm for quenches to various temperatures near and below the KT transition. System sizes ranging up to were employed with periodic boundary conditions. In a kagome lattice, the number of the total spins is . The system is quenched from completely disordered initial states down to a given low temperature with the process of coarsening being monitored through equaltime spatial correlation functions and the decay of topological defects, etc.
In addition to quenches from disordered state to low temperatures, we also performed the socalled nonequilibrium relaxation by suddenly bringing the system from one of the ground states to some target temperatures around or below the KT transition. This method was found to be convenient for measuring the values of the critical exponent for the equilibrium spatial correlation functions.
When we let the systems evolve from random initial configurations the following quantities can be measured:

The number density of topological defects which are vortices in the KAFXY model
(3) where is the total number of vortices (both positive and negative) at time and denotes the total number of sites (i.e., spins). denotes an average over random initial configurations. We also count the separate number density of vortices residing on the hexagonal plaquettes () and those on the triangular plaquettes ().

The equaltime spatial correlation function of the order parameter .
(4) (5) 
Nonequilibrium spin autocorrelation functions
(6) (7)
Iii simulation results
We have performed dynamic Monte Carlo simulations of KAFXY model on a kagome lattice of dimensions . In order to obtain the values of the equilibrium critical exponents for the spatial decay of the equilibrium spatial correlation function of the order parameter, we have employed the socalled nonequilibrium relaxation (NER) method, where the system is suddenly brought from ground states to finite temperatures below or near the KT transition. Simulations were performed up to Monte Carlo steps which was sufficient for the equilibrium to be attained. This was confirmed by the collapse of the spatial correlations at later time stages. In Fig. 2(a) the correlation functions are displayed on a loglog scale. The values of the exponent thus determined for temperatures ranging from to are shown in Fig. 2(b) as well as in Table 1.
0.01  0.047(3)  0.065  0.254(11) 

0.02  0.082(4)  0.066  0.259(13) 
0.03  0.116(5)  0.067  0.264(12) 
0.04  0.153(8)  0.068  0.268(12) 
0.05  0.198(7)  0.069  0.280(11) 
0.06  0.235(9)  0.070  0.300(13) 
0.061  0.237(10)  0.071  0.308(12) 
0.062  0.240(11)  0.072  0.310(13) 
0.063  0.245(10)  0.073  0.330(15) 
0.064  0.250(12)  0.074  0.340(13) 
As for the generation of the initial low energy configurations, we first select one of the three phases of , for each site at random and then applied Monte Carlo procedure at very low temperature of for MC steps with the restriction that, in each of the MC procedure, a new trial phase is chosen only among the above three phases. In this way, the final states obtained were very close (in energy) to a ground state. These states were employed as the initial states of NER in order to bring the system to the equilibrium at some given target temperature.
From the graph of shown in Fig. 2(b), we find that the exponent increases almost linearly in temperature up to around , above which it starts to increase rather sharply. The linear regime can be fit approximately by . We also find that the exponent takes the value of near (Table 1). This temperature is expected to correspond to the KT transition, which apparently is a little lower than the theoreticalcherepa2001 (); korshu2002 () or the numericalrzchow97 () estimates of earlier works mentioned above.
Now, the coarsening dynamics of the model system under quench to low temperature from disordered initial states is investigated by monitoring the equaltime spatial correlation functions of the order parameter. One of the most important features of typical coarsening dynamics toward quasiordered phase is the critical dynamic scaling
(8) 
where is the scaling function and is the growing length scale. It should be noted that the above scaling ansatz is based on the existence of a single growing length scale .
Figure 3(a)(b) shows the equaltime correlation functions of KAFXY model for different time instants (from to ) at the temperature . For these data, we attempted a critical dynamic scaling. The procedure of rescaling is as follows. For a given temperature, we take the value of (determined as in Table 1) and evaluate the combination . And then, for a given time instant, the length scale is determined by the condition where the constant is chosen as or in this work. By plotting the scaled correlation functions in terms of the rescaled distance , we can check whether the critical dynamic scaling holds or not from the quality of the collapse of the scaled correlation functions.
By this procedure, we found, rather unexpectedly, that the critical dynamic scaling does not hold in the coarsening dynamics of KAFXY model, at least for the (Monte Carlo) time duration of our simulations (up to MC steps). The result of the scaling attempt is shown in Fig. 3(c) for the case of where we can see that the critical dyanmic scaling is not obeyed. We also tried different values of for the rescaled correlation with no success. It might be possible that the scaling is restored in the limit of infinitely long time.
Even though the critical dynamic scaling is not faithfully obeyed, we may still extract approximate length scale of spatial correlation in the manner described above. The length scale thus obtained exhibits a rather unusual timedependent behavior (Fig. 3(d)). That is, in the earlytime stage up to around , exhibits a slow growth, which appears approximately independent of the temperature. However, in the intermediate and late time stage of , the growth of the length scale becomes strongly dependent on the temperature. In addition, for a given temperature, no simple power law is found which is valid for the whole late time regime. Instead the local logarithmic slope exhibits a steady increase in time in the late time regime.
We can investigate the timedependence of the growth of the length scale by defining the effective local growth exponent as . Since we took the growing length scale only for discrete time instants with fixed interval in logarithmic scale, we evaluate the discrete version of the above logaritmic slope as
(9) 
where , and . Figures. 3(e)(f) show the effective growth exponent at time (for the case of (Fig. 3(e)) and (Fig. 3(f)) respectively). We can see that, in the late time regime, the growth becomes definitely faster as the temperature is higher. Also, for the temperatures , the effective local growth exponent appears to be monotonically increasing in the late time stage. At higher temperatures, we can recognize some indications from Figs. 3(e)(f) that the effective slopes tend to converge around in the latest time regime. It may be possible to interpret this as an indication of an asymptotic behavior of diffusive growth.
In order to understand the physical mechanism underlying the breakdown of the critical dynamic scaling as well as the peculiar timedependence of the length scale, we investigated the time dependence of the total number of the relevant topological defects, i.e., the fractional vortices residing on the triangular plaquettes as well as on the hexagonal plaquettes of the kagome lattice. The results are shown in Figs. 4(a)(d) where we can see that, for a given temperature, the triangular vortices are decaying much faster than the hexagonal vortices. We also find that in general (for both types of the vortices) the decay of the total vortex numbers do not exhibit a simple power law behavior valid for the whole time range. One prominent feature in terms of the temperature dependence of the decay of the fractional vortex density is that, for the vortices on the triangular plaquettes, the decay rate increases as the temperature is lowered, while on the other hand, those vortices on the hexagonal plaquettes exhibit opposite dependence on the temperature with the decay rate decreasing sharply as the temperature is lowered.
This can be interpreted as implying that there is almost no (free) energy barrier for the motion and decay of the triangular vortices but that some finite barrier exists for the hexagonal vortices. We also find an interesting feature of the hexagonal vortices (at lower temperatures) where, at initial stage, they increase a little and then start to decay. This can probably be understood as the influence of the decay of the neighboring triangular vortices with their excess energy turned over to neighboring hexagonal vortices, thus generating hexagonal vortices.
Therefore, we can conclude that, for wide range of temperatures, the relaxation of the fractional vortices residing on the hexagonal plaquettes determines the coarsening process such as the correlation length scale corresponding to the quasiordered domains. In addition, the vortex relaxations exhibit a considerable deviation from power law behavior with slowly increasing (in absolute magnitude) local logarithmic slope in the late time regime. Especially in the case of the hexagonal vortices, the decay of the vortex number density could be fit to a form with
(10)  
(11) 
with . Here, we found that typically takes values around (figure not shown hereswpark_prep ()). We found that it could also be fitted to an stretched exponential form as
(12) 
with .
It appears that the breaking of the critical dynamic scaling is related to the anomalous relaxation of the number of fractional vortices on the triangular and the hexagonal plaquettes. In order to check this possibility, we analyzed the behavior of the length scales determined by the vortex distributions as follows. We define the length scale which corresponds to the average separation between the vortices on the trianular plaquettes. Similarly, we can define the length scale corresponding to the average separation between vortices on the hexagonal plaquettes. We also define the length scales which represents the length scale corresponding to the total number density of vortices.
Figure 5 shows the three length scales at which clearly shows the faster growth of the length scale for the triangular vortices compared with that for the hexagonal vortices. In Fig. 6(a)(d) comparison is made between the length scales derived from the spatial correlation functions and the length scales derived from the total vortex densities. We see that (for both cases of and ) there exists no simple proportional relationship between the two length scales. Rather, the length scale from the spatial correlation is seen to grow faster than the vortex lenth scale . This probably implies that the vortices are not distributed evenly (statistically speaking) in the system but that the vortices are somehow distributed in a nonrandom manner such that some degree of clustering occurs. Interesting result is that the two length scales satisfy some nontrivial relationships such that with for the case of (Fig. 6(b)) and for the case of (Fig. 6(d)). This means that the vortex configuration of the system do not exhibit simple selfsimilarity at different time instants. Rather the vortices probably tend to cluster more unevenly as time passes by, leading to the correlation length scale growing faster than the length scale derived from the vortex density.
The snapshot of defect configurations for the case of is shown in Fig. 7(a)(d). We can easily recognize faster decay of those vortices on the triangular plaquettes. From the snapshots alone, however, it is not easy to detect the tendency of relative clustering of the vortices in the late time regime.
Now we turn to the autocorrelation functions of the variables and . For these two variables, we find that the simulation results on the autocorrelations are very different. For the case of , we obtain a approximately powerlaw decay behavior in the early time regime up to around . The value of the exponent ranges from (for ) to (for ). The exponent tends to increase slightly at higher temperatures. At later time, however, the autocorrelations decay faster than powerlaw relaxation of the early time stage. Due to the larger statistical fluctuations in the relaxation of the autocorrelation , we could not attempt a comparison of with the growing length scale .
In contrast, in case of the phase , the spin autocorrelaton exhibits less statistical fluctuations with a nonpower law behaviorcherepa94 () that can be reasonably fitted by
(13)  
(14) 
with the exponent and . One of the simulation results is shown in Fig. 9(ab) for the case of with , which shows a suitable fit to the above functional form with and . Note that, in the limit of , reduces to a power law of . The fitted values of implies that exhibits a considerable deviation from a powerlaw behavior. We also attempted to fit the simulation results to stretched exponential form. This resulted in rather small values for the stretching exponent of . It is interesting to note that similar behavior of relaxation in the autocorrelation function of the order parameter was reported in the coarsening dynamics of the socalled Hamiltonian XY model on a square lattice.kjkoo2006 ()
In summary, we investigated the coarsening dynamics of the antiferromagnetic XY model on a kagome lattice. We found that the critical dynamic scaling is violated. Novel and unusual coarsening dynamics may be attributable to the existence of infinite ground state degeneracies leading to nontrivial decay behavior of the density of defects residing on the triangular and hexagonal plaquettes with nontrivial growth of multiple length scales. The competition between the critical thermal fluctuations of the equilibrium and the infinite ground state degeneracy probably gives rise to the nontrivial features of the relaxation dynamics. It may be possible to collapse the equaltime spatial correlation functions through a multiscaling scheme which we haven’t tried yet. It would be also interesting to investigate the aging dynamics of this model system.wills2000 ()
References
 (1) A. J. Bray, Adv. Phys. 43, 357 (1994).
 (2) V. L. Berezinskii, Sov. Phys. JETP 34, 610 (1972).
 (3) J. M. Kosterlitz and D. Thouless, J. Phys. C 6, 1181 (1973); J. M. Kosterlitz, J. Phys. C 7, 1046 (1974).
 (4) E. Kim, B. Kim, and S. J. Lee, Phys. Rev. E 68, 066127 (2003) and references therein.
 (5) E. Kim, S. J. Lee, and B. Kim, Phys. Rev. E 75, 021106 (2007).
 (6) I. Syozi, Prog. Theor. Phys. 6, 306 (1951).
 (7) R. Moessner and A. P. Ramirez, PHYSICS TODAY, February 2006, page 24.
 (8) D. A. Huse and A. D. Rutenberg, Phys. Rev. B 45, 75367539 (1992).
 (9) M. S. Rzchowski, Phys. Rev. B 55, 1174511750 (1997).
 (10) V. B. Cherepanov, I. V. Kolokolov, and E. V. Podivilov, JETP Letters 74, 596 (2001).
 (11) S. E. Korshunov, Phys. Rev. B 65, 054416 (2002).
 (12) M. R. Kolahchi and J. P. Straley, Phys. Rev. B 66, 144502 (2002).
 (13) M. J. Higgins, Yi Xiao, S. Bhattacharya, P. M. Chaikin, S. Sethuraman, R. Bojko and D. Spencer, Phys. Rev. B 61, R894R897 (2000).
 (14) K. Park and D. A. Huse, Phys. Rev. B 64, 134522 (2001).
 (15) Y. Xiao, D. A. Huse, P. M. Chaikin, M. J. Higgins, S. Bhattacharya, D. Spencer, Phys. Rev. B 65, 214503 (2002).
 (16) R. J. Baxter, J. Math. Phys. 11, 784 (1970).
 (17) V. B. Cherepanov, arXiv:condmat/9407068 (1994).
 (18) S. Park and S. J. Lee, unpublished.
 (19) K.J. Koo, W.B. Baek, B. Kim, and S. J. Lee, Jour. Kor. Phys. Soc. 49, 1977 (2006).
 (20) A. S. Wills, V. Dupuis, E. Vincent, J. Hammann, and R. Calemczuk, Phys. Rev. B 62, R9264 (2000)