arxiv: v2 [cond-mat.str-el] 10 Jul 2017

Transcription

1 EP manuscript No. will be inserted by the editor Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet B. Danu 1 and B. Kumar 1 The Institute of Mathematical Sciences, C I T Campus, Chennai 6113, India. School of Physical Sciences, awaharlal Nehru University, New Delhi 1167, India. arxiv: v cond-mat.str-el] 1 ul 17 March, 17 1 Introduction Abstract. We study a frustrated spin-s staggered-dimer Heisenberg model on square lattice by using the bond-operator representation for quantum spins, and investigate the emergence of classical magnetic order from the quantum mechanical staggered-dimer singlet ground state for increasing S. Using triplon analysis, we find the critical couplings for this quantum phase transition to scale as 1/SS + 1. We extend the triplon analysis to include the effect of quintet dimer-states, which proves to be essential for establishing the classical order Néel or collinear in the present study for large S, both in the purely Heisenberg case and also in the model with single-ion anisotropy. PACS m Quantized spin models, including quantum spin frustration 75.3Kz Magnetic transitions classical and quantum 75.1.Kt Quantum spin liquids, valence bond phases and related phenomena The frustrated quantum antiferromagnets are a subject of great interest. The frustration arising due to competing interactions or lattice geometry, helped further by the quantum fluctuations, acts unfavourably towards classical magnetic order, and can give rise to exotic quantum paramagnetic states 1]. Many different materials are nown to display frustrated magnetic behaviour in a variety of ways. For instance, CuGeO 3 ], SrCu BO 3 3],CaV 4 O 9 4], and ZnCu 3 OH 6 Cl 5], are some of the very well nown materials of this ind. On the theoretical side, there is a lot of interest in exploring and understanding the exotic behaviour in variously motivated frustrated quantum spin models. The simplest way in which the frustrated antiferromagnets can exhibit quantum paramagnetism is through the formation of pairwise singlets, called dimers, in the ground state. The pioneering wor of Majumdar and Ghosh showed this exactly in a spin-1/ Heisenberg chain 6]. Similarly, there is a strong numerical support for the spin-1/ 1 - model on square lattice to have non-magnetic columnar dimer ground state in the range of.4 < / 1 <.6 7, 9,1]. The dimer ground states have been shown to occur in many other frustrated models 11 19]. A dimer singlet ground state can also be realized by taing the antiferromagnetic exchange interaction between the spins of selected pairs to be stronger than the rest. For Send offprint requests to: Quantum to classical transition in the ground state of a spin-s instance, in Fig. 1 denoting a Heisenberg model on square lattice, the staggered pattern of thic nearest-neighbour bonds is taen to be stronger compared to the other interactions and f. The spin-1/ case of this staggereddimer model without f has been studied by using the perturbative series expansion method ], exact diagonalization 1], coupled cluster method 1], and quantum Monte Carlo simulation, 3]. These studies found the staggered dimer ground state for strong enough /, which undergoes a transition to the Néel state as / becomes weaer. A case for topological order in this model has also been made on the basis of tensor networ studies 4]. While the studies concerned with dimer states focus mostly on spin-1/ problems, Darradi et al made an interesting investigation of their dependence on spin quantum number, S, in the unfrustrated f = staggereddimer model 5]. Using exact diagonalization, coupledcluster method and variational mean-field approach, they studied the S dependence of the critical coupling for transition from staggered dimer to Néel phase in this model. Based on their calculations for upto S =, they inferred that / c 1/SS + 1. The triplon analysis, based on the bond-operator representation of spins 7,8], is a nice analytical method for studying the ordering instabilities of the dimer singlet states. Using the spin-s bond-operator representation of Kumar 8], in this paper, we do mean-field triplon analysis of the frustrated staggered-dimer model on square lattice see Fig. 1. We find the ground state of this Heisenberg model to undergo continuous phase transition from the

2 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet staggered dimer to Néel or collinear state, and the critical couplings for both the transitions to scale precisely as 1/SS + 1. This lends support to the earlier finding for f = 5]. Interestingly, for strong enough frustration when f /, the Majumdar-Ghosh value we also find the staggered dimer phase to survive in the large S classical limit. It loos a bit unusual, and turns out to be an artefact of considering only the singlet and triplet dimer states in the minimal low-energy description. We improve the triplon analysis by including the effect of quintet states, which gives the expected classical magnetic order for large spins. We also get the classical behaviour in the large S limit by adding an infinitesimally small axial single-ion anisotropy that is ever-present 1. Through this study of the S dependence of quantum phase transitions in the staggered-dimer model, we present an interesting perspective on the approach to classical behaviour from a quantum state with increasing S, and extend the triplon analysis in a simple way to be applicable to the higher spin dimerised Heisenberg antiferromagnets. This paper is organized as follows. In Sec., we state the staggered-dimer model on square lattice, and carry out the simplest mean-field triplon analysis of its purely Heisenberg case, using bond-operator representation in the reduced space of singlet and triplet dimer states. Then, in Sec. 3, we bring quintet states into the discussion. In particular, in Sec. 3., we extend the triplon analysis to include the effect of quintets in a mean-field approach that is nearly as simple as the minimal triplon analysis, but maes significant physical improvement to the results of the calculation. We also point out that how it is further extendable to include the higher dimer states such as septets, octets and so on. In Sec. 4, we study the effect of singleion anisotropy with increasing S. We then conclude with a summary in Sec. 5. Staggered-dimer model on square lattice The Hamiltonian of the spin-s staggered-dimer model on square lattice depicted in Fig. 1 can be written as follows. Ĥ = S R,1.S R, DS R,1z + S R,z R +S R,.S R+aˆx,1 + S R,.S R+aˆx±aŷ,1 ] + f S R,1.S R+aˆx±aŷ,1 + S R,.S R+aˆx±aŷ, 1 Here, R is summed over the positions of the dimers, and the labels 1 and denote the spins of a dimer. The is the intra-dimer exchange interaction. The and f are the interactions between the spins of the nearest dimers. All of these exchange interactions are antiferromagnetic. The single-ion anisotropy, D, is taen to be a positive number. 1 In high spin magnetic materials, it is pretty common to have axial anisotropies, as has been reported in many materials such as DTN NiCl 4SCNH ] 6], NENC NiC H 8N NiCN 4] 7], CsNiCl 3 8], LiFePO 4 9] etc. Fig. 1. The staggered-dimer model on square lattice. The purely Heisenberg case D = of this model would clearly form the spin-gapped staggered-dimer singlet ground state for, f. Here, the dimerization is not spontaneous, but induced by the stronger. Besides a strong, the frustration due to competing and f also facilitates the non-magnetic dimer ground state. However, the Néel or collinear antiferromagnetic AFM order would be realized in the ground state for sufficiently strong or f. The anisotropy, D, also tends to support magnetic order. Clearly, here there is an ample scope for competition between the quantum disordered staggereddimer phase and the classical Néel or collinear ordered phases. We study this competition as a function of the spin quantum number, S, by doing the stability analysis of the staggered-dimer phase using bond operators. The Heisenberg case of Ĥ with D = is discussed in the following subsections and in Sec. 3, and the model for D is studied thereafter..1 Bond-operator representation of spin-s operators The systems where the pairs of spins dimers act as basic units, the bond-operator representation of spins has been found to be very useful in doing calculations. This approach was pioneered by Sachdev and Bhatt for studying quantum phase transitions in dimerized spin-1/ quantum antiferromagnets 7]. For the case of general spin-s, the bond operator representation was derived by Kumar 8]. The basic idea is to construct the spin eigenstates of a dimer, and associate with each one of them a distinct bosonic creation operator called bond operator. For a pair of spin-s interacting via exchange interaction, S 1 S, the suitable dimer states are the total-spin eigenstates. That is, the spin singlet: s, triplet: t m1, quintet: q m, and so on. The corresponding bond operators are defined as: s = ŝ, tm1 = ˆt m 1, qm = ˆq m. Here, denotes the vacuum of the bosonic Foc space. In order to describe the finite-dimensional spin Hilbert space f

3 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet 3 of the dimer, these bond operators in the Foc space are required to satisfy the constraint: ŝ ŝ+ˆt m 1ˆt m1 + ˆq m ˆq m + = 1, where is meant to denote the contributions from all the higher total-spin dimer states upto S. The exchange operator on the dimer can now be written as: S 1.S = SS + 1 ŝ ŝ + ˆt m 1ˆt m1 + ˆq m ˆq m + ] + ˆt m 1ˆt m1 + 3ˆq m ˆq m + ]. To write the interactions between the spins of different dimers, one needs the bond-operator representation of the individual spins. At the very least, it can constructed in the restricted space of singlet and triplet, and can also be extended further to include quintet or higher states, if necessary. The bond operator representation of S 1 and S in the restricted space of singlet and triplet has the following form 8]. SS + 1 S 1,α ± ŝ ˆt α + ˆt 3 αŝ i ɛ αβγˆt βˆt γ, 3 Here, α = x, y, z, and liewise for β and γ. The ɛ αβγ is the totally antisymmetric tensor, and the signs + and correspond to the spin-label 1 and, respectively.. Mean-field triplon analysis of the Heisenberg case The staggered-dimer Heisenberg model, Ĥ of Eq. 1 with D =, in the bond-operator representation turns into a model of interacting triplets and singlets, which is hard to mae progress with. To simplify it, we mae the following standard approximations 7, 18]. We treat the singlet bond-operators for the staggered singlet bonds of Fig. 1 in the mean-field approximation as: ŝ = ŝ = s. It now becomes a model of interacting triplet excitations with a mean singlet bacground. To eep it simple, we also neglect the interaction between triplets, and impose the local constraint on the bond-operators through an average Lagrange multiplier, λ. The resulting Hamiltonian is bilinear in triplet bond-operators, called triplons. Finally, by doing the Fourier transformation: ˆt Rα = 1 Nd ei.rˆt α, we get the following effective triplon dynamics. Ĥ t = N d SS ] λ + s λ + 1 { ] λ s SS + 1ξ ˆt αˆt α + ˆt αˆt α,α s SS + 1ξ ˆt αˆt α + ˆt αˆt α }. 4 Here, lies in the Brillouin zone of the staggerd-dimer square lattice, N d is the total number of dimers, and ξ = 3 cos xa+ 4 3 cos xa cos y a 8 f 3 cos xa cos y a. The triplon Hamiltonian, Ĥ t, can be diagonalized by the Bogoliubov transformation for bosons, ˆt α =ˆγ α cosh θ ˆγ α sinh θ, with θ = 1 tanh 1 { s SS + 1ξ /λ s SS + 1ξ ]}. The Ĥt in terms of the quasiparticle operators, ˆγ α, has the following diagonal form. Ĥ t = e N d + E ˆγ αˆγ α + 1 5,α Here, e = SS λ + s λ, and E = λλ s SS + 1ξ ] is the triplon quasiparticle dispersion. The ground state energy per dimer is given by e g λ, s ] = e + 3 E. 6 N d Minimization of the ground state energy density, e g, with respect to the mean-field parameters, λ and s, gives the following equations. s = 5 3 N d λ = + 3λSS + 1 N d λ s SS + 1ξ E 7 ξ E 8 We determine s and λ by self-consistently solving the Eqs. 7 and 8 in the space of / and f /. In the staggered-dimer phase, the quasiparticle gap is non-zero. However, for some critical values of / and f /, this gap would vanish. The contours of such critical points define the boundaries of the staggered-dimer phase. On the other side of this phase boundary, the condensation of triplons at some wavevector Q, given by E Q =, determines the antiferromagnetic order. The vanishing of the gap also fixes the value of λ at λ = s SS + 1ξ Q...1 Quantum phase diagram We set the intra-dimer interaction to 1, and vary and f as free parameters. In the present case, D =. We now describe the quantum phase diagram of the staggereddimer Heisenberg model on square lattice in the - f plane, as obtained from the mean-field triplon analysis. From our calculations, we find two different ordering wavevectors at which the triplon gap vanishes. It is Q N =, that gives Néel AFM order for sufficiently strong, and Q C =, π that gives collinear AFM order for strong enough f. We determine the phase-boundaries between the staggered-dimer and the two ordered AFM phases from the following equation. ] ξ Q 5 3 ξ Q = N d ξ Q ξ 3 SS The quantum phase diagram, thus obtained, is presented in Fig.. In Eq. 9, the interaction parameters, and f, and the spin quantum number, S, appear separately on the left and right hand side, respectively. It clearly implies that the phase boundaries for all the different values of S in the left panel of Fig. will collapse onto two lines only

4 4 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet f S=1/ S=1 S=3/ S= S= f SS Q=,π Collinear Phase Staggered Dimer Phase Q=, Néel Phase SS+1 Fig.. Left: The mean-field quantum phase diagram in the - f plane for = 1 and D =. The region above the upper phase boundaries is the collinear phase, below the lower phase boundary is the Néel phase, and between the two is the staggereddimer phase. For S=, the phase boundaries are f =.489 and f =.519 for Néel and collinear transition, respectively. Right: The quantum phase diagram in the space of rescaled parameters. one for dimer to Néel and the other for dimer to collinear transition in the space of rescaled interactions, SS + 1 and f SS + 1, as shown in the right panel of Fig.. For f =, the quantum critical point for staggereddimer to Néel transition is given by c.388/ss + 1. For =, the critical f for transition to the collinear phase is fc.377/ss + 1. For f =, we can compare the results of our calculations with those of others obtained by different methods. For S = 1/, we get c =.518, while it is around.4 from other numerical calculations,, 3]. Tae, for instance, the values of Darradi et al 5] for c.4545,.1818,.99 obtained from exact diagonalization method, and compare with c.5188,.194,.138 from our calculation for S = 1/, 1 and 3/, respectively. The agreement is quite reasonable, and seems to get better as S increases. Moreover, in the limit S, when either or f is zero unfrustrated cases, the ordered classical phases are the only possibility. However, when both and f are non-zero, the large S limit turns out to be surprising, as discussed below. In the limit S, the right hand side of Eq. 9 vanishes, which then gives two limiting phase boundaries: f =.4891 for staggered-dimer to Néel transition, and f =.519 for the transition to collinear phase. A small difference in their slopes leaves a tiny region of staggereddimer phase between the two limiting phase boundaries see Fig.. It presents a puzzling situation for a Heisenberg problem in which a quantum mechanical staggereddimer state, helped by strong frustration, survives into the classical limit. Quite possibly, this situation is an artefact of the oversimplified low-energy triplon dynamics, and can be resolved by including the corrections due to higher spin dimer-states. Below we particularly study the effect of quintets on the effective triplon dynamics. 3 Quintet corrections to triplon analysis To include quintets in the effective low-energy dynamics of the spin model, we consider the bond-operator representation in the subspace of singlet, triplet and quintet dimer-states, as given in Ref. 8]. S1, z ± ŝ ˆt N + ˆt ŝ ± s + 1 ˆt 1ˆq 1 + ˆq 1ˆt 1 + ˆt 1ˆq 1 + ˆq 1ˆt 1 + ˆq 1 ˆq 1 ˆq 1 ˆq 1 + ˆq ˆq ˆq ˆq S 1, + ± ŝ ˆt 1 N ˆt 1ŝ ± s 1 3 ˆt ˆq + ˆq ˆt ] + 1 ˆt 1ˆt 1 ˆt 1ˆt 1 ˆt 1ˆq ˆq ˆt ˆt ˆq 1 ˆq 1ˆt ˆt 1ˆq ˆq ˆt ˆt 1ˆt + ˆt ˆt 1 + ˆq ˆq 1 + ˆq 1 ˆq + ] 3 ˆq 1 ˆq + ˆq ˆq where = S + 1, = SS + 1S + 1/3, and N q = SS + 1S 1S + 1S + 3/15. Using Eqs. 1 and 11, our Heisenberg model would loo lie Ĥ = Ĥt + sĥˆtˆtˆq S + Ĥˆtˆtˆqˆq S + Ĥˆtˆtˆtˆt + Ĥˆqˆqˆqˆq, of which the first term is the triplon model wored out in the previous section. The next significant contribution in powers of S comes from Ĥˆtˆtˆq and Ĥˆtˆtˆqˆq. Below we discuss the corrections to triplon dynamics from these terms involving quintet operators. First we compute the corrections using second order perturbation theory, which is however limited in its scope. More importantly, in Sec. 3., we extend the triplon analysis to include the quintets in a

5 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet 5 simple way, which nicely resolves the issue with the classical limit within triplon mean-field theory, and is also extendible in principle to include the higher spin dimerstates. 3.1 Second order perturbation theory in gapped phase The simultaneous coupling of the triplets to the singlet bacground and the quintet excitations will allow sĥˆtˆtˆq to contribute more directly to the low-energy dynamics of the dimer phase, as compared to the other terms involving quintets. Therefore, here we calculate the second order perturbative correction only due to sĥˆtˆtˆq to the ground state energy, and compute the change in the quantum critical points. We can write this term on a lattice as: Ĥ sˆtˆtˆq = s 4 m=,, i,δ ] q m] i,n ˆT i,i+δ + h.c. 1 where ˆT i,i+δ = ˆt i,α ˆt i+δ,β + ˆt i+δ,β 3]. The second order correction to ground state energy can be written as: s E = S sˆtˆtˆq 16 N q m,m h i,δ,i,δ m] ˆq i,m ˆT i,i +δ h h ˆq m] i,m ˆT i,i+δ E E h where the summation over h represents the sum over all the excited states h mj. In the -space the composite triplet operator m] δ ˆT i,i+δ can be written as, m] δ ˆT i,i+δ = 1 m] N d p, ˆT +p, eip.r ɛ. Treating the quintet excitations as local with the constraint i, m = i, m, we get E S sˆtˆtˆq = 1 N d N q 1 s m The operator following way. n p, m] n ˆT +p, ɛ.13 E E n m] ˆT +p for m =,, is defined in the ˆT ±] +p, = ˆt +px ˆt x + ˆt x ˆt +py ˆt y + ˆt y iˆt +px ˆt y + ˆt y iˆt +py ˆt x + ˆt x ˆT ±1] +p, = ˆt +pz ˆt x + ˆt x ˆt +px ˆt z + ˆt,z +iˆt +pz ˆt y + ˆt y + iˆt +py ˆt z + ˆt,z ˆT ] +p, = 3 ˆt +px ˆt x + ˆt x ˆt +py ˆt y + ˆt y + ˆt +pz ˆt z + ˆt z ] a local quintet excitations with an energy cost = + λ and two triplet excitations. After using the Bogolibov transformation of triplets in Eq. 13, and computing the non-zero matrix elements in new two triplon quasiparticle state, ˆγ 1αˆγ, such as β 1α, β p α, β. The summation term in Eq. 13 for each m can be written as, m] n ˆT +p, ɛ n p, E E n = sinh θ ɛ +pcosh θ +p sinh θ +p p, ] + sinh θ +p ɛ cosh θ sinh θ / + λ + E + E p. The ground state energy per dimer including second order correction can be written: e g λ, s ] = e + 3 E + E S sˆtˆtˆq 14 N d Minimising e g with respect to s and λ, we have calculated the critical points for transition from dimerised phase to Néel phase. We find the values of critical points as: c.153,.8,.5,.34,.5,.9] for S=1, 3/,, 5/, 3 and 5, respectively. These values approximately consistent with Darradi et al 5] calculations and lie in between their variational mean-field lie calculation and the coupled-cluster method, by which they were able to calculate the critical points only for S 3/. For comparison their values of c through CCM: c =.156 S = 1 and c =.91 S = 3/. Notably, the values of critical points from our perturabtive calculation with the inclusion of quintet states reduces approximately around - 5%, with respect to the critical values obtained from the minimal triplon mean-field theory. But one finds it difficult to generate the phase boundaries continuously in the whole parameter space using the second order perturbation theory. We, therefore, find a better and easier way to include quintets into the triplon mean-field theory. 3. Extended triplon analysis with quintet mean-fields: Improved quantum phase diagram For a neater presentation, let us first redefine the quintet operators in the real spherical-harmonics form. ˆq xy = i ˆq ˆq, ˆq x y = 1 ˆq + ˆq ˆq z = ˆq, ˆq xz = 1 ˆq 1 ˆq 1, ˆq yz = i ˆq 1 + ˆq 1 From Eq. 13 one can notice that, within the second order perturbation theory, only the non-zero matrix element will contribute to the correction energy for a state n of three quasiparticle excitations. This states will be composed of With an s describing the singlet, the Ĥˆtˆtˆq type terms being linear in ˆq s will induce some static values for the quintet operators in the mean-field approximation. Hence, we introduce five real quintet mean-field variables, q z, q x y,

6 6 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet q xy, q yz and q xz. Now the intra-dimer exchange interaction can be written as: S 1 S = SS + 1 s + { SS ] q z + q x y + q xy + q yz + q xz} + SS ] α Moreover, the constraint can be written as: s + q z + q x y + q xy + q yz + q xz + α ˆt αˆt α. 15 ˆt αˆt α = 1 16 As we are interested in the approach to classical behavior for increasing S, of all the terms in the bond-operator representation given Eqs. 1 and 11, we only eep those with the highest power of S. After approximating the singlet and quintet operators by their mean-fields, we get the following simplified representation of the spin operators. S 1x S 1y S 1z S x S y S z { s I + 1 Q { s I + 1 Q Here, I is the 3 3 identity matrix, and 3 q xx q xy q xz Q = q xy 3 q yy q yz. q xz q yz 3 q zz } ˆt x + ˆt x ˆt y + ˆt y 17 ˆt z + ˆt z } ˆt x + ˆt x ˆt y + ˆt y 18 ˆt z + ˆt z Moreover, q xx + q yy + q zz =, where q zz = q z, q xx = q z 3 q x y / and q yy = q z + 3 q x y /. The quintet-corrected effective triplon dynamics of the Heisenberg model on staggered-dimer square lattice is now given by the following Hamiltonian. H t = N d SS λ + s λ + + λ q x y + q z ] + λ 1 s,α { ɛ ˆt,α + ˆt,α,α,β I + 1 ˆt αˆt α + ˆt αˆt α ] } Q ˆt,β N + ˆt,β t αβ 19 Here, with unit lattice spacing, ɛ = cos x + cos x cos y ] 8 f cos x cos y. The triplon Hamiltonian, Ht, can be simplified by diagonalising the real symmetric matrix, s 1 Q, which is purely dependent on s and q s. Under this diagonalisation, the three component triplon vector would also rotate canonically. In terms of these new triplons, the different values of α= x, y, z do not mix with each other, and one can study H t in the same way as Ĥt of Eq. 4. Note that, due to the tracelessness of Q, its three eigenvalues would have the same general structure as its diagonal elements, q xx etc. Hence, without any serious loss of generality, we consider a simpler case with q xy = q xz = q yz =. The ground state energy per dimer, e g, of Ht can be written as: e g = e + + λ q z + 1 q x y + E x N d I+ + 1 E y + 1 E z, N d N d where the three quasiparticle dispersions are gives as follows. E,α = λ λɛ s + q αα Note that the triplon dispersions here will in general be non-degenerate. The presence of q s will pull one or more of these dispersions lower, which will eventually become gapless at some Q to cause magnetic order. We will find below that, indeed, this improved triplon dynamics exhibits stronger tendency towards ordering. Minimising the e g of Eq. with respect to λ, s, q z and q x y gives us the following mean-field equations. s = 5 q z q x y λ = + 1 N d 1 N d 1 N d λ N d + λ N d + λ N d λ ɛ s + E x q xx λ ɛ s + q yy E y λ ɛ s + q zz + 1 s + 1 s + 1 s E z q xx q yy q zz q zz = 1 λ N q 4N d q xx + s + λ 1 λ N q 4N d q yy + s + λ ɛ E y ɛ E z ɛ E x ɛ E x ɛ E y

7 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet 7.4 S=1.4 S=3/.3 Collinear Phase.3 Collinear Phase f..1 Néel Phase f..1 Néel Phase f S=5/.3..1 Collinear Phase Néel Phase f S= Collinear Phase Néel Phase Fig. 3. The quantum phase diagrams obtained from the extended triplon analysis for different values of S. The blac lines represents the phase boundaries obtained from the extended triplon analysis including quintets, and, the red lines are the phase boundaries obtained from the minimal triplon analysis shown in Fig.. The boundaries of the Néel and collinear phases merge into the Majumdar-Ghosh line, f = /, at the tricritical point filled blac circle, tc. + 1 λ N q N d q zz + s + λ q x y = 1 3λ N q 4N d q xx + s + λ 1 3λ N q 4N d q yy + s + λ ɛ E z ɛ E x ɛ. E y To generate the phase boundaries between the staggereddimer and AFM phases, we solve the above self-consistent mean-field equations for λ, s, q z and q x y, iteratively. We find the minimum ground state energy solution to have q zz and q x y =. This leads to the lowering of the dispersion, E z. The quantum phase diagrams obtained from this approach for different values of S are shown in Fig 3. The phase boundaries are generated for the ordering wavevector Q =, for the Néel phase, and Q =, π for the collinear phase, with E Qz =, which fixes the chemical potential at λ N = ɛ t. Q s + q zz Here, we find the staggered-dimer phase to shrin rapidly, as the spin quantum number S increases. For S 3/, the two critical phase boundaries of the staggered-dimer phase are found to merge at the tricritical point, tc, beyond which the Néel and collinear phases are directly separated from each other by the well-nown Majumdar-Ghosh line, f / = 1/ 31,3]. As plotted in Fig. 4, the tc tends to zero as S increases. Notably, for numerical comparison at f =, our improved critical c for S = 1 is.18141, which nearly same as the exact diagonalization estimate of c =.1818 for S=1, and for S = 3/ it is.951 in our calculation which lies in between the values of c =.99,.917 estimated from ED and CCM, respectively for S = 3/ in Ref. 5]. It is satisfying to see that our extended triplon mean-field theory is able to get the consistent physical behaviour while approaching the classical limit from a quantum state. We close this section by briefly noting that this approach for dimerized quantum antiferromagnets can be extended to include the septets, octets etc. by treating the even total-spin bondoperators in mean-field, while eeping the odd total-spin bond-operators as such. It would to lead an effective bilinear probelm in the odd spin bond-operators, which can be studied in the same manner as done here.

8 8 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet tc /SS+1 Fig. 4. The tricritical point, tc, versus 1/SS + 1. The tc is defined as the distance from the origin of the point of merger of the Néel and collinear phase boundaries in Fig Staggered-dimer model including single-ion anisotropy Let us first understand a single antiferromagnetic dimer with axial spin anisotropy. Ĥ dimer = S 1.S DS 1z + S z. 1 A non-zero D causes splitting and mixing between the total-spin states of a dimer. For instance, now the singlet mixes with the triplet as well as the quintet states. This we can see by constructing S 1z and S z in the subspace of singlet, triplet and quintet states. Thus, in the bondoperator representation, S 1,z + ŝ ˆq + ˆq ŝ + N q + ˆt ˆt ŝ ŝ ˆt 1ˆt 1 + ˆt 1ˆt 1 ± ˆt N 1ˆq 1 + ˆq 1ˆt 1 ˆt 1ˆq 1 ˆq 1ˆt 1 t S + 11S 15 ˆq 1 ˆq + 1 6S + 6S 5 ˆq 1 14 ˆq 1 + ˆq ˆq S + S + 5 ˆq 7 ˆq + ˆq ˆq. Notice the mixing between ŝ and ˆq in Eq., which in addition to the triplons, will most directly affect the dimer singlet state. Hence, we now write Ĥdimer of Eq. 1 in terms of ŝ, ˆt α and ˆq only, and ignore all the other bondoperators for simplicity. In this approximation, the Ĥdimer taes the following form. Ĥ dimer SS + 1 ŝ ŝ + ˆt αˆt α + ˆt zˆt z + ˆq ˆq +ˆt αˆt α + ˆt zˆt z + 3ˆq ˆq N D t ŝ ŝ D 3 ŝ ˆq + ˆq ŝ N D t + ˆt zˆt z D + 1 ˆt αˆt α D 1 11S + 11S 15 ˆq ˆq 3 where α = x, y. The mixing of ŝ with ˆq in Eq. 3 can be diagonalized by unitary rotation of ŝ and ˆq to new operators s and q such that ŝ = s cos θ q sin θ and ˆq = s sin θ + q cos θ. Here, tan θ = D / 3 + D D S + 11S 15 ]. The Ĥdimer in terms of the bond operators s, ˆt α and q has the following diagonal form. Ĥ dimer 1 s s + q q SS D D 11S + 11S 15 ] s s + q q D sin θ 3 D + N q ˆt 3N zˆt z D + 1 ˆt t N αˆt α 4 t 4.1 Modified bond-operator representation The eigenvalues of the diagonalized dimer Hamiltonian of Eq. 4 as a function of D are plotted in Fig. 5. As the quintet state, q, is much higher in energy, for a lowenergy effective description, we can ignore it. Moreover, in the mean-field approximation s is replaced by s. Therefore, we get Ĥ dimer 1 s + SS D D 1 11S + 11S 15 D ] sin θ 3 + SS D + 1 ]ˆt αt N ˆ α t + SS D + N q ]ˆt N zˆt z. 5 t The bond-operator representation for S 1α and S α is also required to be reconstructed for this case. From Ref. 8], we now S 1α and S α in the restricted space of { ŝ, ˆt α, ˆt z, ˆq }. After taing care of the mixing between ŝ and ˆq, and then using the mean-field approximation for s operator and ignoring the q excitations, we get the following expressions for S 1α and S α. S 1α = S α s cos θ 1 sin θ ˆt α + ˆt α 6 S 1z = S z s cos θ + sin θ ˆt z + ˆt z 7 Here, α = x, y. Using this representation, now we can do the standard triplon analysis of the full staggered-dimer model of Eq. 1 including the spin anisotropy, D.

9 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet 9 E -4 s E t t -8 1,-1 q 1,-1 q,- q D Fig. 5. Energy levels splittings of dimer bond with the single ion anisotropy for S=1. 4. Triplon analysis and magnetic order for D The triplon Hamiltonian for the Ĥ of Eq. 1, including single-ion anisotropy has the following form. { Ht D = e D N d + 1 λ α s φ α ξ ˆt αˆt α + h.c,α 1 s φ α ξ ˆt αˆt α + h.c} { + 1 λ z s φ z ξ ˆt zˆt z + h.c 1 s φ z ξ ˆt zˆt z + h.c}. 8 Here, ξ = cos x a + 4 cos x a cos y a 8 f cos x a cos y a, and the parameters φ α and φ z as functions of S, and D are defined as, φ α=x,y = cos θ sin θ and φ z = cos θ sin θ. Similarly, the spin anisotropy induced split effective chemical potentials are: λ α=x,y = λ D λ D +, and e D = SS λ + D +D + N q + 1, and, λ z = s + λ D D 1 11S + 11S 15 D sin θ 3 ]. After Bogoliubov diagonalization, we get Ĥt D = e D N d + E α ˆγ αˆγ α + 1,α + E z ˆγ zˆγ z Here, E α = λ α λ α s φ α ξ for α = x, y and E z = λz λ z s φ z ξ are the quasiparticle dispersions. The ground state energy per dimer is e D g λ, s ] = e D + 1 E α + 1 E z. 3 N d N d,α Minimizing, e D g, with respect to λ and s gives the following mean-field equations. { s = 5 1 λ α s φ α ξ } λ z s φ z ξ 31 N d E α E z,α λ = + D + D 1 11S + 11S 15 + D sin θ λ α s φ α ξ + 1 λ z s φ z ξ 3 N d E α N d E z,α To generate the boundaries of the gapped staggereddimer phase, we trac the quasiparticle gap by solving Eqs self-consistently for λ and s. Below we discuss the quantum phase diagram thus generated. 4.3 Quantum phase diagram for D The closing of triplon gap at Q =, in the Néel phase, and at Q =, π in the collinear phase occurs through E Qz =, which fixes λ z at λ z = s φ z ξ Q. This is similar to the dispersions of the Heisenberg case in Sec. 3.. For positive D, the E α always remains gapped for α = x, y. Recall that, for D =, in the absence of quintet corrections we always found the dime phase, even for S classical limit. But now, for D >, we find that, above some critical value, D c at f =.5, the Néel and collinear phase boundaries merge into a single line for any value of S, and the staggered-dimer phase completely vanishes. In Fig. 6, we show this for S = 1, S = 3/, S = and S = 5. Thus, a sufficiently strong D will always mae the dimer phase disappear for any S 1. We estimate the critical values of D for first few values of S, which are D c.6633 for S = 1, D c.1563 for S = 3/, D c.7 for S =, and D c.94 for S = 5. We note that for large S, the critical D c becomes negligibly small, and the ground state is always a classical. In fact, D c 1/SS + 1 for large S, as estimated from the data for upto S = 15 in Fig. 7. This calculation clearly shows that, in the limit S, an infinitesimal amount of single-ion anisotropy would suffice to completely suppress the quantum mechanical dimer phase, and to fully establish classical order. 5 Conclusion: Here, we have studied the quantum phase transitions in a frustrated staggered-dimer model for arbitrary spin S on square lattice by doing triplon mean-field theory, without and with single-ion anisotropy. First, we performed

10 1 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet f.4.3. D=. D=.3 D=.5 y/x=.5 S=1 f.4.3. D=. D=.1 D=.13 y/x=.5 S=3/.1.1 f D=. D=.1 D=.1 y/x=.5 S= f D=. D=.1 y/x=.5 S= Fig. 6. The quantum phase diagram with the single ion anisotropy for S = 1 top-left, for S = 3/ top-right, for S = bottom-left, and for S = 5 bottom-right. The boundaries of the Néel and collinear phases merge into a single line at f =.5 for strong enough single-ion anisotropy D D c tc D c /SS+1 Fig. 7. Critical values of D c as a function of spin quantum number S from S = 1 upto S = 15. The tc and D c as function of 1/SS + 1 for merging two phase boundaries shown in subfigure. Both tc and D c tends to zero as we approaches classical limit. the minimal bond-operator mean-field theory in the reduced Hilbert space of singlet and triplets, without singleion anisotropy. Through this, in the purely Heisenberg case, we found two ordered phases, viz. Néel and collinear states, sandwiching the spin-gapped staggered-dimer phase. The critical couplings for these transitions are found to scale precisely as 1/SS + 1. Puzzlingly, here we also found that this quantum staggered-dimer phase, helped by strong frustration, persists even in the limit S which is nown to be a fully classical ordered state. This contradiction leads us to revise the singlet-triplet meanfield calculation by including the quintet states. To include quintet states first we performed a second order perturbation theory and present the improve critical points for transition from dimer to Néel phase. Next, we discuss the quintet participation in the Hamiltonian through a minimal singlet-triplet-quintet mean field theory. We find that this extended mean-field analysis is sufficient to consistently establish the classical AFM order for larger S. In fact through this, we find the Néel and collinear phases merge into a single line, f = /, at the tricritical point tc 1/SS + 1. In our case, f = / is the levelcrossing line between two classical AFM phases. Hence, we call the merger point, tc, as a tricritical point because it is a point where the phase boundaries of t two continuous phase transitions meet with a line of first order transition. Here, it is dimer to Néel, dimer to collinear, and, Néel to collinear phase transition. Thus, in general, infinitesimally small and f would be sufficient to establish classical order in the limit of very large S. In the presence of single-ion anisotropy, we first rewored the

11 B. Danu, B. Kumar: Quantum to classical transition in the ground state of a spin-s quantum antiferromagnet 11 spin-s bond-operator representation to include the effects of singlet-quintet mixing on a dimer. Then, we used it to carry out the standard triplon analysis to wor out the quantum phase diagram. Interestingly, for any S 1, we always find a critical value of the anisotropy, D c, for which the dimer phase between the Néel and collinear phases vanishes. We also find that D c 1/SS + 1 for large S. Thus, an infinitesimally small D would be sufficient to ill the quantum mechanical dimer phase and establish classical order in the limit of very large S. Through this model calculation, we have presented an interesting perspective on the change of behaviour from quantum to classical in frustrated spin systems with increasing S. It would be nice to chec this in other models, and if possible, using cold atoms in optical lattices 33 36]. The methods that we have developed here would, in general, be useful to dimerized quantum antiferromagnets of larger spins. 3. R. Ganesh et al., Phys. Rev. B. 84, P. Chandra et al., Phys. Rev. B. 38, B. Danu et al., Phys. Rev. B. 94, A. V. Gorshov et al., Phys. Rev. Lett. 17, S. Gopalarishnan et al., Phys. Rev. Lett. 17, Simon et al., Nature 47, oseph W. Britton et al., Nature 484, Acnowledgments B. D. thans R. Ganesh for useful discussions. B. K. acnowledges financial support under the UPE-II and DST- PURSE programs of NU, New Delhi. References 1. C. Lacroix, P. Mendels and F. Mila, Springer Series in Solid-State Sciences, Introduction to Frustrated Magnetism, 164, 11.. M. Hase et al., Phys. Rev. Lett. 7, Kageyama et al., Phys. Rev. Lett. 8, Taniguchi et al.,. Phys. Soc. pn. 64, Helton et al., Phys. Rev. Lett. 98, C. K. Majumdar and D. K. Ghosh. Math. Phys. 1, S. Sachdev and R. N. Bhatt, Phys. Rev. B 41, B. Kumar, Phys. Rev. B 8, D. Poilblanc et al., Phys. Rev. B 43, Oitmaa and Z. Weihong, Phys. Rev. B 54, B. S. Shastry and B. Sutherland, Physica B. 18, D.. Klein,. Phys A: Math. Gen. 15, I. Bose and P. Mitra, Phys. Rev. B 44, B. Kumar, Phys. Rev. B 66, B. Danu, B. Kumar and R. V. Pai, EPL 1, B. Kumar, Phys. Rev. B 71, A. Gelle et al., Phys. Rev. B 77, R. Kumar and B. Kumar, Phys. Rev. B 77, R. Kumar et al., Phys. Rev. B 8, R. P. Rajiv Singh et al., Phys. Rev. Lett. 61, S.E. Krüger et al., Phys. Rev. B 61, S. Wenzel et al., Phys. Rev. Lett. 11, F.-. iang, Phys. Rev. B 85, Ching-Yu Huang et al., Phys. Rev. B 84, R Darradi et al.,. Phys.: Condens. Matter 17, S. A. Zvyagin et al., Phys. Rev. Lett. 98, M. Orendáč et al., Phys. Rev. B 5, W.. L. Buyers et al., Phys. Rev. Lett. 56, iying Li et al., Phys. Rev. B. 73, 4416.