## Quantum spin-liquid behavior in the spin-1/2 random-bond Heisenberg antiferromagnet on the kagome lattice

Hikaru Kawamura, ∗ Ken Watanabe, and Tokuro Shimokawa

Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka 560-0043, Japan (Dated: August 26, 2014)

The effect of the quenched bond-randomness on the ordering of the S = 1 2 antiferromagnetic / Heisenberg model on the kagome lattice is investigated by means of an exact-diagonalization method. When the randomness exceeds a critical value, the ground state of the model exhibits a transition within the non-magnetic state into the randomness-relevant gapless spin-liquid state. Implications to the S = 1 2 kagome-lattice antiferromagnet herbertsmithite is discussed. /

Since the proposal by P.W. Anderson of the resonating valence bond state [1], the quantum spin-liquid (QSL) state possibly realized in certain S = 1 2 frustrated / magnets has attracted interest of researchers [2]. After a long experimental quest, several candidate materials were recently reported in geometrically frustrated magnets, including both the triangular-lattice and the kagome-lattice antiferromagnets (AFs).

Examples of the triangular-lattice AF might be S = 1 / 2 organic salts such as κ -(ET) Cu (CN) 2 2 3 [3-6] and EtMe Sb[Pd(dmit) ] 3 2 2 [7-10], which exhibit no magnetic ordering, neither regular nor random, down to a very low temperature. The QSL state of these compounds exhibits gapless behaviors characterized by, e.g. , the lowtemperature specific heat [4, 9] or the thermal conductivity [5, 8] linear in the absolute temperature T . An example of the second category, the kagome-lattice AF, might be herbertsmithite CuZn (OH) Cl 3 6 2 [11-17]. No magnetic order is observed down to 20mK in spite of the exchange coupling of 180K, while gapless behaviors are observed in various physical quantities.

Most of the recent theories pre-assumes that the system is clean enough that the quenched randomness plays a negligible role. Meanwhile, it was recently suggested in ref.[33] that the QSL of the triangular-lattice organic salts might be the randomness-induced one, the randomsinglet [34-36] or the valence-bond glass (VBG) state [37, 38], in which the randomness is self-generated for the spin degrees of freedom at low temperatures via the random freezing of the electric-polarization degrees of freedom [6, 10]. Such a randomness-induced QSL picture apparently explains many of the experimentally observed features of the triangular organic salts [33].

The origin of the QSL behavior observed in these triangular-lattice organic salts and the kagome-lattice herbertsmithite has attracted tremendous interest, and various theoretical proposals have been made. The simplest possible reference model might be the S = 1 2 AF / Heisenberg model with the nearest-neighbor (n.n.) bilinear coupling. In the triangular case, the ground state of such n.n. model is known to be the standard AF longrange order (LRO), the 120-degrees structure [18, 19]. Hence, to explain the QSL behavior observed in these organic salts, some ingredient not taken into account in the simplest Heisenberg model is required. In the kagome case, by contrast, the ground state of the n.n. model is believed to be some sort of QSL state without the magnetic LRO, although its nature is still under hot debate. Various scenarios, including the Z 2 spin liquid [20-24], the algebraic U(1) spin liquid [25-28], the chiral spin liquid [29] and the valence bond crystal (VBC) [30-32] etc. have been proposed. It is not entirely clear at the present stage, however, which of these, or any of these, applies to the experimentally observed QSL state of herbertsmithite.

In case of herbertsmithite, it has been realized from the earlier stage of research that certain amount of quenched randomness exists originated from the random substitution of magnetic Cu 2+ by non-magnetic Zn 2+ or vice versa. It was suggested initially that the substitution occurred on the kagome plane, providing the site-dilution for Cu 2+ on the kagome layer [13, 14]. More recently, it was suggested that 15% of Zn 2+ on the adjacent triangular layer were randomly substituted by Cu 2+ , keeping the kagome layer intact [16]. This might give rise to the Cu 2+ [kagome]-Cu 2+ [triangular]-Cu 2+ [kagome] exchange path, the kagome-triangular path being ferromagnetic.

The effect of the site-dilution randomness was analyzed by Singh by the series-expansion method [38]. Starting with the VBC picture for the non-random kagome AF [31], he argued that the random dilution induced the VBG state with gapless behaviors. In view of the work by Singh on the site-random kagome AF and of the recent experimental finding that the randomness in herbertsmithite might primarily of the bond-random-type, and also of the recent work of ref.[33] on the bond-random triangular AF, we wish to investigate the effect of the bond-randomness on the ordering of the S = 1 2 kagome / Heisenberg AF.

We consider the AF bond-random S = 1 2 Heisenberg / model on the kagome lattice whose Hamiltonian is given by

$$\mathcal { H } = \sum _ { < i j > } J _ { i j } \vec { S } _ { i } \cdot \vec { S } _ { j } - H \sum _ { i } S _ { i } ^ { z }, \quad \ \ ( 1 )$$

where /vector S i = ( S , S x i y i , S z i ) is a spin-1/2 operator at the i -th

site on the lattice, and H is the magnetic-field intensity. For the exchange coupling J ij , we assume for simplicity the random n.n. AF coupling obeying the bondindependent uniform distribution between [(1 -∆) J , (1 + ∆) J ] with the mean J . The parameter ∆ represents the extent of the randomness: ∆ = 0 corresponds to the regular system and ∆ = 1 to the maximally random system.

In herbertsmithite, the randomness for the kagomelayer Cu 2+ would come primarily from the interlayer exchange path, adding an effective ferromagnetic coupling to the original AF coupling on the kagome layer, somewhat reducing the magnitude of the AF coupling in the bond-random manner. Since an apical Cu 2+ on the triangular layer contributes to three Cu 2+ n.n. bonds on the kagome layer, 15% Cu 2+ substitution on the triangular layer affects significant amount of the kagome-layer Cu 2+ bonds. Although the substituted Cu 2+ eventually couples different kagome layers giving rise to the threedimensional (3D) coupling, it has been established experimentally that herbertsmithite behaves as a 2D system without the 3D order. Hence, we choose to consider the effect of the randomly substituted Cu 2+ somewhat indirectly, only within the 2D kagome model. Note that, in our modelling, the direct contribution from the interlayer Cu 2+ themselves is neglected.

Thus, we adopt here the simplest form of the bondrandomness, with the aim of getting insight into new qualitative features possibly induced by the bond randomness. By means of an exact diagonalization (ED) method, both the zero- and finite-temperature properties of the model are computed for finite lattices. Another reason of our choice of the randomness is a computational one, i.e. , this form of randomness turns out to yield rather stable results even for very small sizes accessible by the ED method. The ED technique meets a difficulty for certain types of randomness, e.g. , the sitedilution randomness, where small system sizes hamper a systematic study of the size dependence keeping the randomness extent (hole concentration) fixed.

The total number of spins is N = 12 -30 for T = 0, and N = 12 15 and 18 for , T &gt; 0, periodic boundary conditions being employed. Sample average is taken over 100 ( N = 12 16 18 21 24) 50 ( , , , , N = 27) and 12 ( N = 30) independent bond realizations in the T = 0 calculation, while 100 ( N = 12 15) and 20 ( , N = 18) in the T &gt; 0 calculation. In what follows, the energy, temperature and magnetic field are normalized in units of J .

We first examine the existence or nonexistence of the magnetic LRO at T = 0, by computing the squared sublattice magnetization m 2 s . We consider the three sublattices associated with the /vector q = 0 order, while for N of multiples of nine ( N = 18 27) the nine sublattices associ-, ated with the √ 3 × √ 3 order in addition. In Fig.1(a), we show the N -dependence of the ground-state m 2 s -values for several values of the randomness ∆. In the Heisen-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-b8a37dff-95fc-43c2-8ef8-28291d608bab.jpg)

-

FIG. 1. (color online). (a) The squared sublattice magnetization m 2 s associated with the three sublattices (solid symbols) and with the nine sublattices (open symbols, for N = 18 27 , only) for several randomness ∆. The lines are linear fits of the data. (b) The spin freezing parameter ¯ versus 1 q / √ N . The lines are the linear fits using the four largestN data. The inset is a magnified view of the largeN region.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-9dddd589-9a91-46ff-87cc-056f79b76f6a.jpg)

berg model sustaining a nonzero m 2 s , the spin-wave form, m 2 s = m 2 s, ∞ + c 1 √ N , is expected for its N -dependence. As can be seen from the figure, the data extrapolated to N = ∞ on this form yields negative values for all ∆, indicating that the AF LRO is absent either in the regular or in the random case.

A nontrivial issue is whether the randomness ever induces the spin-glass (SG) order or not. To examine this possibility, we compute the SG-type spin freezing parameter ¯, defined by ¯ q q 2 = 1 N 2 ∑ i,j [ 〈 /vector S i · /vector S j 〉 2 ], 〈· · · 〉 and [ · · · ] being the ground-state expectation and the sample average, respectively, and the sum over i and j is taken over all sites of the lattice. The quantity takes a nonzero value when the spin is frozen either in a spatially periodic or random manner [39]. As can be seen from Fig.1(b), the randomness tends to enhance the ¯-value somewhat. q Yet, linear extrapolation of the data yields ¯ = 0 within q the error bar for any ∆. For ∆ = 1, although the extrapolation using all available N yields a very small nonzero value ¯ = 0 006 q . ± 0 002, the extrapolation using only four .

largest N yields ¯ = q -0 003 . ± 0 006, indicating the van-. ishing ¯. q Thus, we tend to conclude that the spin freezing does not occur for any ∆, even including ∆ = 1.

To get further information, we compute the quantities related to the scalar chirality χ c = /vector S 1 · /vector S 2 × /vector S 3 ( /vector S 1 , /vector S 2 and /vector S 3 are three corner spins on an upward triangle t on the lattice). While the chirality locally takes a nonzero value because of quantum fluctuations, the chiral freezing parameter ¯ , defined by q χ ¯ q 2 χ = ( 3 N ) 2 ∑ t,t ′ [ 〈 χ χ t t ′ 〉 2 ], vanishes in the N → ∞ limit at any ∆, indicating that the long-range chiral order does not occur, even the random one. Supplemental material gives the details. This observation further supports our conclusion above that the SG-type order does not occur for any ∆. The absence of the chiral order in the non-random case was also reported in ref.[23].

Next, we study the the finite-temperature properties of the model. In Fig.2(a), we show the temperature dependence of the specific heat per spin C (in units of Boltzmann constant) for several values of ∆. The specific heat of the random model exhibits a T -linear behavior C /similarequal γT at lower temperatures. For ∆ = 0 and 0.3, the finiteN data exhibit a lowT peak around T /similarequal 0 1 in . addition to the one at a higher temperature. Although the recent calculation for the larger lattice of N = 30 suggests that this lowT peak might go away in the N →∞ limit [40], a weaker structure like a hump or a cusp is likely to remain even in the the N →∞ limit. By contrast, our data for the random model indicates that even such a weaker lowT anomaly does not appear for ∆ ≥ 0 5. . Hence, some changeover of the lowT property seems to occur when the randomness exceeds a critical value ∆ c /similarequal 0 4. .

In the inset of Fig.2(b), we show the ground-state magnetization per spin M versus an applied field for our largest size N = 30. It has been known that the nonrandom model exhibits a plateau-like structure at the 1 / 3 of the saturation value (refer also to ref.[41]), with further additional structures at 1 / 9, 5/9 and 7/9 [42]. For ∆ &gt; ∆ , such plateau-like anomalies tend to go away, c yielding a near linear behavior similar to the one observed in the corresponding triangular model.

In Fig.2(b), we show the temperature dependence of the magnetic susceptibility per spin χ for several values of ∆. While the susceptibility goes to zero in the T → 0 limit for weaker randomness ∆ /lessorsimilar 0 3, it tends to a finite . value at ∆ = 0 5, and diverges obeying the Curie law . for stronger randomness ∆ ≥ 0 7. . The appearance of a Curie-like component suggests that a small fraction of free spins ( ∼ 4%) are generated at low temperatures.

The ground state for weaker randomness turns out to be a spin singlet corresponding to the total spin S = 0 for all evenN samples (for odd N , a doublet corresponding to S = 1 2). As the randomness gets stronger, some / fraction of samples has a spin-triplet ground state. In Fig.3, we show the ratio of evenN samples with spintriplet ground states, R , versus ∆. As the randomness

FIG. 2. (color online). The temperature dependence of (a) the specific heat per spin C , and of (b) the uniform susceptibility per spin χ , for several randomness ∆ for N = 12 (open symbols) and 18 (solid symbols). The inset of (b) represents the T = 0 magnetization curve for N = 30.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-84d13993-0350-446c-8dd9-54aba16c2353.jpg)

gets stronger beyond ∆ c /similarequal 0 4, . R becomes nonzero, and increases with ∆. For comparison, we show R of the S = 1 2 triangular-lattice model in the inset, for which / ∆ c /similarequal 0 6 separates the AF and the random-singlet states. . While the triplet tendency is enhanced in the kagome case than in the triangular case, the two data look similar. Hence, some sort of rapid change, possibly a phase transition, is likely to occur at the critical randomness ∆ c /similarequal 0 4 in the ground state of the kagome model. This . is corroborated by our data of the specific heat and the susceptibility or the magnetization of Fig.2 where the change of behavior is also observed around ∆ = ∆ c /similarequal 0 4. .

The state at ∆ &gt; ∆ has properties similar to those of c the random-singlet state identified in the corresponding triangular-lattice model, including the T -linear specific heat, the gapless susceptibility with a Curie-like component for the stronger randomness, and the near linear magnetization curve without plateau-like structures [33]. The gapless feature is also consistent with the one suggested for the site-dilution model of ref.[38]. So, we identify the rapid change observed at ∆ = ∆ c /similarequal 0 4 . as a transition from the randomness-irrelevant QSL state to

FIG. 3. (color online). The randomness ∆ dependence of the ratio of samples with triplet grounds states for various even N . The inset is the corresponding plot for the S = 1 2 / triangular-lattice model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-c674fa73-4530-42e1-9720-349735db20f6.jpg)

FIG. 4. (color online). The dynamical spin structure factor S /vector, ω ( q ) at T = 0 at the Γ point plotted versus ω for the randomness ∆ = 1, in comparison with that of the non-random model of ∆ = 0. The broadening factor η introduced as in refs.[43, 44] is set to η = 0 02. The inset is a magnified view of . the lowω region (indicated by the arrow), the ω -value transformed to that of herbertsmithite with J = 17 meV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-37dd2c26-344a-49e6-8b54-f89067a9e145.jpg)

the randomness-relevant QSL state, the random-singlet (or VBG) state. Our data are not capable of discriminating the former being either the gapped Z 2 liquid or the gapless U (1) liquid, but the observation of the vanishing chiral order seems to rule out the chiral QSL state. If the QSL state of the non-random system is the gapped Z 2 spin liquid, it would be stable against a weak randomness up to a certain point, and ∆ c would be a true critical point separating the gapped/gapless behaviors.

To further characterize each phase, we compute the dynamical spin structure factor S /vector, ω ( q ), a quantity accessible by inelastic neutron scattering, following refs.[43, 44]. The computational details are given in Supplemental material. Its ω -dependence computed for ∆ = 1 at the Γ point /vector q = (0 , 2 π √ 3 d ) ( d the lattice spacing) is shown in Fig.4, together with the one for the non-random case ∆ = 0 for comparison [44]. A similar calculation is also made for the M point /vector q = (0 , π √ 3 d ). The intensity at the

Γ-point is somewhat greater than that at the M -point, but its ω -dependence is rather similar.

The size-dependence of the random system turns out to be rather mild. The gapless behavior is observed both at the Γ and the M points, without any sharp structure. At ∆ = 1, the observed intensity exhibits a mild peak at ω = 0, with a very broad, almost flat plateau extending to higher ω . The ω = 0 peak is not pronounced, the peak height being about twice the intensity at the plateau part. In the non-random case of ∆ = 0, by contrast, the lowenergy peak is much more enhanced, the peak intensity being an order of magnitude greater [44]. The observed features of the data, a broad ω /similarequal 0 peak and almost flat intensity extending to higher ω , especially no gap appearing at any /vector q , are qualitatively consistent with the recent inelastic neutron-scattering data on single crystals [17]. Since the experimental data now available are limited to the low-energy region of ω /lessorsimilar 10meV ( J /similarequal 17meV), we show in the inset a magnified view of the experimentally accessible lowω range. The data of the random system appears to resemble the experimental one much more than the one of the non-random system. This observation lends support to the view that the randomness is certainly playing a significant role in herbertsmithite.

The remaining question is whether herbertsmithite lies at ∆ &lt; ∆ or at ∆ c &gt; ∆ . c If the former is the case, the QSL state might be a mild modification of that of the non-random system, though the translational symmetry is lost at the Hamiltonian level as soon as ∆ becomes nonzero. If the latter is the case, on the other hand, the QSL state of herbertsmithite would be the randomsinglet state. Although the direct contribution from the interlayer Cu 2+ and the effect of the DzaloshinskiiMoriya interaction needs to be considered for further comparison, the gapless features observed experimentally in various quantities, the single-crystal S /vector, ω ( q ) data in particular, seems to favor the latter scenario, i.e. , the QSL state of herbertsmithite be the random-singlet state. Further experiments at higher fields as well as inelastic neutron measurements at higherω are desirable.

In summary, we investigated both the zero- and finitetemperature properties of the spin-1/2 bond-random Heisenberg AF on the kagome lattice, to find that the model exhibited beyond a critical randomness a gapless QSL behavior, which was attributed to the randomsinglet (or VBG) state. The result might consistently explain the recent experimental data on herbertsmithite.

The authors are thankful to ISSP, the University of Tokyo for providing us with CPU time. This study is supported by Grants-in-Aid for Scientific Research No. 25247064. The calculation was performed by use of TITPACK Ver.2.

- ∗ kawamura@ess.sci.osaka-u.ac.jp
- [1] P.W. Anderson: Mater. Res. Bull. 8 (1973) 153.
- [2] L. Balents: Nature 464 (2010) 199.
- [3] Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato and G. Saito: Phys. Rev. Lett. 91 (2003), 107001.
- [4] S. Yamashita, Y. Nakazawa, M. Oguni, Y. Oshima, H. Nojiri, Y. Shimizu, K. Miyagawa and K. Kanoda: Nature Physics 4 (2008) 459.
- [5] M. Yamashita, N. Nakata, Y. Kasahara, T. Sasaki, N. Yoneyama, N. Kobayashi, S. Fujimoto, T. Shibauchi and Y. Matsuda: Nature Physics 5 (2009) 44.
- [6] M.A. Jawad, I. Terasaki, T. Sasaki, N. Yoneyama, N. Kobayashi, Y. Uesu and C. Hotta: Phys. Rev. B 82 (2010) 125119.
- [7] T. Itou, A. Oyamada, S. Maegawa, M. Tamura and R. Kato: Phys. Rev. B 77 (2008) 104413.
- [8] M. Yamashita, N. Nakata, Y. Senshu, M. Nagata, H.M. Yamamoto, R. Kato, T. Shibauchi and Y. Matsuda: Science 328 (2010) 1246.
- [9] S. Yamashita, T. Yamamoto, Y. Nakazawa, M. Tamura and R. Kato: Nature Commun. 2 (2011) 275.
- [10] M.A. Jawad, N. Tajima, R. Kato, and I. Terasaki: Phys. Rev. B 88 (2013) 075139.
- [11] M.P. Shores, E.A. Nytko, B.M. Bartlett, and D.G. Nocera: J. Am. Chem. Soc., 127 (2005) 13462.
- [12] J.S. Helton, K. Matan, M.P. Shores, E.A. Nytko, B.M. Barlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. Chung, D.G. Nocera, and Y.S. Lee: Phys. Rev. Lett. 98 (2007), 107204.
- [13] M.A. de Vries, K.V. Kamenev, W.A. Kockelmann, J. Sanchez-Benitez, and A. Harrison: Phys. Rev. Lett. 100 (2008) 157205.
- [14] P. Mendelse and F. Bert: J. Phys. Soc. Jpn. 79 (2010) 011001.
- [15] J.S. Helton, K. Matan, E.A. Nytko, B.M. Barlett, Y. Qiu, D.G. Nocera, and Y.S. Lee: Phys. Rev. Lett. 104 (2010) 147201.
- [16] D.E. Freedman, T.H. Han, A. Prodi, P. M¨ uller, Q.-Z. Huan, Y.-S. Chen, S.M. Webb, Y.S. Lee, T.M. McQueen, and D.G. Nocera: J. Am. Chem. Soc., 132 (2010) 16185.
- [17] T.-H. Han, J.S. Helton, S. Chu, D.G. Nocera, J.A. Rodriguez-Rivera, C. Broholm and Y.S. Lee: Nature, 492 (2012) 406.
- [18] B. Bernu, P. Lecheminant, C. Lh¨illier u and L. Pierre: Phys. Rev. B 50 (1994) 10048.
- [19] L. Capriotti, A.E. Trumper and S. Sorella: Phys. Rev. Lett. 82 (1999) 3899.
- [20] S. Sachdev: Phys. Rev. B 45 (1992) 12377.
- [21] Y.-M. Lu, Y. Ran and P.A. Lee: Phys.. Rev. B 83 (2011) 224413.
- [22] S. Yan, D.A. Huse and S.R. Whilte: Science, 332 (2011) 1173.
- [23] S. Depenbrock, I.P. McCulloch, U. Sch¨ ollwock: Phys. Rev. Lett. 109 (2012), 067201.
- [24] H.-C. Jiang, Z. Wang, and L. Balents: Nature Physics 8 (2012), 902.
- [25] M.B. Hastings: Phys. Rev. B 63 (2000) 014413.
- [26] Y. Ran, M. Hermele, P.A. Lee, and X.-G. Wen: Phys. Rev. Lett. 98 (2007) 117205.
- [27] H. Nakano and T. Sakai: J. Phys. Soc. Jpn. 80 (2011) 053704.
- [28] Y. Iqbal, F. Becca, S. Sorrela, and D. Poilblanc: Phys. Rev. B 87 (2013) 060405.
- [29] L. Messio, B. Bernu, and C. Lh¨illier: u Phys. Rev. Lett. 108 (2012) 207204.
- [30] J.B. Marston and C. Zeng: J. Appl. Phys. 69 (1991) 5962.
- [31] R.R.P. Singh and D.A. Huse: Phys. Rev. B 76 (2007) 180407; R.R.P. Singh and D.A. Huse: Phys. Rev.B 77 (2008) 144415.
- [32] G. Evenbly and G. Vidal: Phys. Rev. Lett. 104 (2010), 187203.
- [33] K. Watanabe, H. Kawamura, H. Nakano and T. Sakai: J. Phys. Soc. Jpn. 83 (2014) 034714.
- [34] C. Dasgupta and S.-K. Ma: Phys. Rev. B 22 (1980) 1305.
- [35] R.N. Bhatt and P.A. Lee: Phys. Rev. Lett. 48 (1982) 344.
- [36] D.S. Fisher: Phys. Rev. B 50 (1994) 3799.
- [37] M. Tarzia and G. Biroli: Europhys. Lett. 82 (2008) 67008.
- [38] R.R.P. Singh: Phys. Rev. Lett. 104 (2010) 177203.
- [39] D.X. Viet and H. Kawamura: Phys. Rev. Lett. 102 (2009) 027202.
- [40] S. Suigiura and A. Shimizu: Phys. Rev. Lett. 111 (2013) 010401.
- [41] H. Nakano and T. Sakai, to appear in J. Phys. Soc. Jpn. (arXiv:1408.4538).
- [42] S. Nishimoto, N. Shibata and C. Hotta: Nature Communications 4 (2013) 459.
- [43] E.R. Gagliano and C.A. Balseiro: Phys. Rev. Lett. 59 (1987) 2999.
- [44] A.M. L¨ auchli and C. Lh¨illier: u arXiv:0901.1065 (2009).

## Supplemental Material

## Chirality-related quantities

We compute both the local scalar-chirality amplitude ¯ and the chiral freezing parameter χ ¯ . q χ The former is defined by

$$\bar { \chi } = \sqrt { \frac { 3 } { N } \left [ \sum _ { t } \langle \chi _ { t } ^ { 2 } \rangle \right ] }, \\ \intertext { a over } t \text{ is taken over all } N / 3$$

where the sum over t is taken over all N/ 3 upward triangles on the kagome lattice. The chiral freezing parameter is defined by

$$\bar { q } _ { \chi } = \frac { 3 } { N } \sqrt { \left [ \sum _ { t, t ^ { \prime } } \langle \chi _ { t } \chi _ { t ^ { \prime } } \rangle ^ { 2 } \right ] }. \\ \text{sum over } t \text{ and } t ^ { \prime } \text{ is taken over all } \Lambda$$

where the sum over t and t ′ is taken over all N/ 3 upward triangles.

These quantities are computed by the ED method at T = 0. The largest system size is limited to N = 24

-

-

FIG. 5. (color online). The T = 0 chiral freezing parameter normalized by its local amplitudes, ¯ q χ /χ ¯ , 2 is plotted versus 1 / √ N for several randomness ∆. The lines are linear fits of the data. In the inset, the local scalar-chirality amplitude ¯ χ at T = 0 is plotted versus 1 /N .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0606v2-figure-6a8dbcb4-d9ee-466b-8590-2df10dc0d7df.jpg)

for the chirality-related quantities due to the computational load. The computed ¯ and ¯ χ q χ normalized by its local amplitudes, ¯ q χ /χ ¯ , are shown in the figure. 2 As can be seen from its inset, the local scalar-chirality amplitude becomes nonzero for any ∆, even including ∆ = 0, indicating that the spin tends to take noncoplanar configuration in the ground state at the local level because of quantum fluctuations. While the scalar chirality takes a nonzero value locally, our data indicates that the longrange scalar chirality, both the uniform and the glasstype random ones, tends to zero for N →∞ irrespective of the ∆-value.

## Dynamical spin structure factor

The dynamical spin structure factor S ( q , ω ) at T = 0 is defined by

$$\underset { \dots } { i t s o f } \underset { \dots } { \hat { \chi } } \quad S ( q, \omega ) = - \frac { 1 } { \pi } \text{Im} \left [ \langle S ^ { z } ( - q ) \frac { 1 } { \omega - ( \mathcal { H } - E _ { 0 } ) + i \eta } S ^ { z } ( q ) \rangle \right ], \\ \quad. \quad \underset { \dots } { \dots } \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where S z ( q ) is the Fourier transform of the z -component of the spin, and E 0 is the ground-state energy. We then compute S ( q , ω ) by the ED method combined with a continued fraction technique [see refs.[43, 44]]. The broadening factor η is set to η = 0 02. .
Hikaru Kawamura Ken Watanabe Tokuro Shimokawa
Quantum spin-liquid behavior in the spin-1/2 random-bond Heisenberg antiferromagnet on the kagome lattice
The effect of the quenched bond-randomness on the ordering of the $S=1/2$
antiferromagnetic Heisenberg model on the kagome lattice is investigated by
means of an exact-diagonalization method. When the randomness exceeds a
critical value, the ground state of the model exhibits a transition within the
non-magnetic state into the randomness-relevant gapless spin-liquid state.
Implications to the S=1/2 kagome-lattice antiferromagnet herbertsmithite is
discussed.