arxiv:cond-mat/ v1 [cond-mat.str-el] 21 Sep 2006

Transcription

1 Universal emergence of te one-tird plateau in te magnetization process of frustrated quantum spin cains arxiv:cond-mat/69555v [cond-mat.str-el] Sep 6 F. Heidric-Meisner, I. A. Sergienko, A. E. Feiguin, and E. R. Dagotto Materials Science and Tecnology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 3783, USA and Department of Pysics and Astronomy, Te University of Tennessee, Knoxville, Tennessee 37996, USA Microsoft Project Q, Te University of California at Santa Barbara, Santa Barbara, CA 936, USA (Dated: September, 6) We present a numerical study of te magnetization process of frustrated quantum spin-s cains wit S=, 3/, as well as te classical limit. Using te exact diagonalization and density-matrix renormalization tecniques, we provide evidence tat a plateau at one tird of te saturation magnetization exists in te magnetization curve of frustrated spin-s cains wit S > /. Similar to te case of S=/, tis plateau state breaks te translational symmetry of te Hamiltonian and realizes an up-up-down pattern in te spin component parallel to te external field. Our study furter sows tat tis plateau exists bot in te cases of an isotropic excange and in te easy-axis regime for spin-s=, 3/, and, but is absent in classical frustrated spin cains wit isotropic interactions. We discuss te magnetic pase diagram of frustrated spin- and spin-3/ cains as well as oter emergent features of te magnetization process suc as singularities, jumps, and even-odd effects. A quantitative comparison of te one-tird plateau in te easy-axis regime between spin- and spin-3/ cains on te one and and te classical frustrated cain on te oter and indicates tat te critical frustration and te pase boundaries of tis state rapidly approac te classical result as te spin S increases. I. INTRODUCTION Te searcfornovelquantum pasesas been a strong motivation for teoretical investigations of quantum spin models. Quantum states tat do not resemble classically ordered magnetic states may be realized as a consequence of reduced dimensionality or competing interactions. Famous examples are te Haldane pase of integer spin cains, or te spin liquid ground states of spin-/ ladder systems, bot of wic exibit a finite spin gap. 3 5 Te presence of competing interactions causes frustration and drives quantum pase transition, 6 as for instance in te frustrated spin-/ cain tat undergoes a transition from a critical, gapless pase into a spontaneously dimerized one. 7,8 Of great interest are te excitations above tese ground states, wic are reflected in experimentally accessible quantities suc as termodynamic properties. Here we are concerned wit te interplay of frustration and an additional external parameter, te magnetic field. A magnetic field can close zero-field gaps, but may as well induce gapped magnetic excitations in finite fields. 9, Te magnetic pase diagram can experimentally be mapped out by studying te magnetization process, wic teoretically may be computed by means of powerful and flexible numerical metods suc as te density matrix renormalization group tecnique (DMRG)., One of te intriguing features of te magnetization curve of low-dimensional quantum magnets is te emergence of plateaux at finite magnetic field, wic indicate te presence of massive excitations. Tis penomenon as been predicted for several low-dimensional spin models. 9,,3 9 Experimentally, it as been observed in many quasi-low dimensional magnetic materials, realizing networks of spin-/ (Refs. 3 34) as well as of spin- moments. 35 It is te purpose of tis work to study te magnetization process of frustrated quantum spin cains wit S > / and establis te existence of a plateau state at one tird of te full magnetization M. Te possible existence of magnetization plateaux in quantum spin cains was predicted by Osikawa, Yamanaka, and Affleck. 9 A notable result of tis work provides a necessary condition for te existence of plateau states at finite M: ps( M) Z. () Here, S is te spin, p denotes te period (or lengt of te elementary unit cell) of te plateau state in real space. Te teorem implies tat M always as a rational value on any plateau and it allows for translational invariance to be broken spontaneously, i.e., p can be larger tan. Note tat unfrustrated spin-s cains usually sow a smoot magnetization curve, as as been known for quite some time already, 36 implying tat competing interactions beyond a simple nearest neigbor model, anisotropies, or special geometries are responsible for te plateau formation. Indeed, teoretically, te existence of plateaux as been establised for, e.g., S=3/ cains wit onsiteanisotropy 9,6, atm=/,infrustratedand dimerized spin-/ cains at M=/4, 3 tree leg ladders at M=/3, or exotic models suc as te ortogonal dimer cain. 3 Te simplest model wit competing interactions in one dimension is a spin-/ cain wit nearest-neigbor(nn) and next-nearest-neigbor (NNN) excange interactions, also called a zigzag ladder. Only quite recently te existence of a M=/3 plateau in te magnetization curve of frustrated antiferromagnetic S=/ cains as been revealed by DMRG calculations. 4 Te plateau is accompanied by broken translational symmetry wit a period p = 3 and an up-up-down (uud) structure in te spin

2 component parallel to te external field. Wile several iger-s cain models ave been studied in te context of magnetization plateaux, 9,5,6,,,37 te case of frustratedcainswit S > /asnotbeen addressedinsufficient detail. In particular, te existence of te M =/3 plateau in tese systems, wic we report in tis work, as not been explored. Our main result, based on exact diagonalization(ed), DMRG, classical Monte-Carlo simulations (MC), and spin-wave teory, is tat te M =/3 plateau is realized in frustrated quantum spin cains wit S > /in te caseofisotropicexcangeand in te easyaxis regime as we sow for te specific examples of S=, 3/, and. As for spin-/, 4 tis plateau always as broken translational symmetry wit period p = 3 and te uud structure. AM=/3plateauwit tis Néel-typeoforderemerges in te Ising-limit ofbot quantum and classicalfrustrated cains (see Ref. 4 anf furter references terein). Wile te classical frustrated cain does not exibit a finite magnetization plateau in te case of isotropic interactions as we sow in tis work, te M=/3 plateau is stable against quantum fluctuations. Te parameter region, owever, were suc plateau exists for te case of isotropic interactions srinks as te spin S grows. As we compute te full magnetization curves, we also obtain information on oter intriguing features suc as singularities, macroscopic jumps, and even-odd effects in te magnetization curve, and a very ric pase diagram is indeed found. We present a qualitative discussion of te magnetic pase diagram of spin- and spin-3/ cains, wic extends te picture emerging from previous studies for zero 4 46 and finite magnetic fields. 47 In contrast to te case of isotropic interactions, we find tat te M =/3 plateau region prevails in a substantial part of te magnetic pase diagram in te easy-axis regime. Its boundaries rapidly approac te classical result in te easy-axis regime. To tis end, we study te magnetic pase diagram of te classical frustrated cain as well by means of MC simulations and linear spin-wave teory, aiming at a comparison of emergent pases wit te quantum cases of S= and S=3/. In a very recent experimental study, a M=/3 plateau as been observed in te frustrated diamond spin-/ material Cu 3 (CO 3 ) (OH) (Ref. 33). Te existence of a /3 plateau as also been reported for te spin-/ trimer compound Cu 3 (P O 6 OH) (Ref. 34). Our results may be of particular relevance for materials tat ave been suggested to realize frustrated spin- cains. One example is CaV O 5 wit competing next and next-nearest neigbor antiferromagnetic interactions. 48 A promising family of materials realizing iger spin-s zigzag cains is NaR(WO 4 ) were R represents In, Sc (S=), Cr (S=/), and V (S= or S=). 49 Te excange constants in te latter case are estimated to be of te order of 8K, 49 wic may be small enoug to access a substantial part of te magnetization curve. Finally, te formation of spin- zigzag cains as recently been reported for Tl Ru O 7, for wic a Haldane gap of K is found. 5 We ope tat our results will stimulate te searc for te one-tird plateau in potential frustrated spin cain materials. Te restofte paperis organizedas follows. In Sec. II, we introduce te model and briefly discuss te numerical tecniques employed in tis work: DMRG, Lanczos, and MC simulations. In Sec. III, we present our numerical results for te magnetization process of frustrated spin cains wit S=, 3/, and and isotropic excange interactions. We demonstrate te existence of te M =/3 plateau for tese quantum spins and numerically determine te pase boundaries of te plateau region for spin and 3/ in te field vs frustration plane. We study te magnetization process of spin- cains in more detail, focusing on (i) singularities in te magnetization curve, (ii) low-field pases, and (iii) and te combined effect of frustration and onsite anisotropies. Te latter is important as most materials tat realize Haldane cains ave a substantial single-ion anisotropy. In Sec. IV, we turn to te case of anisotropic frustrated cains in te easy-axis region. Te main purpose of tis section is to demonstrate ow te M=/3 plateau observed in frustrated quantum spin-s cains is connected to te classical limit. To tis end, we first analyze a classical frustrated spin cain by means of MC simulations, wic allows us to map out te ground-state pase diagram of te classical model in a finite magnetic field. Analytical results for several pase boundaries are presented. We ten perform DMRG calculations for frustrated spin- and spin-3/ cains at (i) a fixed excange anisotropy, varying te frustration, and (ii) at a fixed frustration, varying te excange anisotropy, and determine te pase boundaries of te M=/3 state for tese two cases. A summary and conclusions are presented in Sec. IV. An Appendix is devoted to te details ofalinear spin-wave calculation for te M =/3 plateau state in te classical limit. II. MODEL AND METHODS Te Hamiltonian of a frustrated spin-s cain wit onsite anisotropy D in a magnetic field reads: H = i [J S i S i+ +J S i S i+ S z i +D(Sz i ) ], () were J > and J > are te antiferromagnetic NN and NNN excange integrals, respectively. Here, S i = (Si x,sy i,sz i ) denotes a spin-s operator acting on site i. We also include an excange anisotropy, te same for te NN and NNN interactions. Hence, in our notation, wit te usual definition of raising and lowering operators S i ± = Si x ±isy i : S i S j (S+ i S j +S i S+ j )+ Sz i Sz j. (3) We restrict our analysis to, and we set J = and D = unless stated oterwise. Te following normal-

3 3 ization for te magnetization is used [tis also applies to Eq. ()]: M = Sz total SN ; Sz total = i S z i, (4) (a) Spin, J =, J =, = N=5, PBC N=, OBC (iv) EO region: S z total = (v) (vi) were N is te total number of sites in te cain. Tus, M= at te saturation field sat. We compute te ground-state energies E (Stotal z, = ) in subspaces labeled by Stotal z on cains wit periodic boundary conditions(pbc) using te Lanczos algoritm. Te ground-state energies of substantially larger cains wit open boundary conditions (OBC) are calculated wit DMRG. Typically, we keep up to m = 4 states in our DMRG calculations. Ten, we include te Zeeman term and obtain te field-dependent ground-state energies E (S z total,) = E (S z total, = ) Sz total. (5) Te magnetization curves are constructed by solving te equations E (Stotal z, step) = E (Stotal z + s, step), wic define tose magnetic fields at wic te magnetization increases from M = Stotal z /(N S) to M = (Stotal z +s)/(n S). Steps larger tan s= may occur. In te classical limit, we treat S i as tree-dimensional unit vectors given by te polar angle Θ i and azimutal angle φ i. For eac value of, we find te ground state of H using te simulated annealing classical MC algoritm. Starting from a random configuration of angles {Θ i,φ i } and an initial inverse temperature β = β in, we proceed as follows. We perform sweeps troug te cain, and for eac site i =,... N, two steps are performed: (i) Te energy E = E({Θ i,φ i }) of te current configuration is calculated from Eq. (). (ii) New angles Θ i = Θ i+r w and φ = φ+r w are tried, were r [,] is a random number and w is te maximal cange allowed. Te new angles are accepted wit probability p = min{exp[β(e E )],}, were E = E({Θ i,φ i }) is te energy of te configuration wit Θ i and φ i replaced by Θ i and φ i. After eac sweep troug te cain, β is multiplied byte factor(β fi /β in ) /G, wereβ fi is te final inverse temperature and G is te number of sweeps. w is adjusted suc tat te acceptance ratio is kept close to 5% for efficiency. Typically, we use β in =.5, β fi =, and G = 5. Tis slow-cooling tecnique prevents te simulation from getting stuck in local minima in te energy landscape in most cases. Sometimes, close to pase boundaries, te system gets trapped in local minima, but repeating te algoritm a few times for te same set of parametersusuallysolveste problem. After β fi asbeen reaced, we fine-tune te ground-state energy by repeating steps (i) and (ii) wit β = for several sweeps until te energy cange after a sweep becomes negligible.. (i) Haldane gap (ii) (iii) M=/3 plateau (b) Spin 3/, J =, J =, = N=, PBC N=96, OBC M=/3 plateau (c) Spin, J =, J =.36, = N=9, PBC N=6, OBC M=/3 plateau FIG. : (Color online) Magnetization curve for frustrated cains wit an isotropic excange ( = ) for (a) S= (J = ), (b) S=3/ (J = ), and (c) S= (J =.36). Solid(dased) lines are DMRG(ED) results. Te orizontal dotted lines mark M=/3 and M=. III. MAGNETIZATION PROCESS OF FRUSTRATED QUANTUM SPIN CHAINS WITH ISOTROPIC EXCHANGE INTERACTIONS In tis section we present results for te magnetization curves of frustrated spin-s cains wit =. As a main result of tis paper, we provide numerical evidence tat te M =/3 plateau, previously observed in te magneti-

4 4 <S i z >/S.5.5 M=/3, J =, = Spin, J = Spin 3/, J = Spin, J =.36 Spin S J,crit,l J,crit,u / / TABLE I: Isotropic excange, =: Lower and upper critical frustrations J,crit,l/u for te formation of te M=/3 plateau for spin / (from Ref. 4), spin, 3/, and site i FIG. : (Color online) Onsite spin density S z i /S on te M=/3 plateau state of frustrated spin-s cains wit S=,3/, and and isotropic excange interactions ( = ) at J = (S= and 3/) and J =.36 (S=). We sow DMRG results for N = sites. zation curve of frustrated spin-/ cains, 4,5 exists for iger spins as well. Te boundaries of te plateau pase strongly depend on S, and te region in te J plane were te M=/3 state is realized srinks as a function of S. Our MC simulations in te classical limit sow no evidence for a magnetization plateau in te isotropic case. A. Magnetization curves of frustrated S=, 3/, and cains Te magnetization curves at an intermediate frustration strengt are presented in Fig. for te cases of spin [Fig. (a), J =], spin 3/ [Fig. (b), J =], and spin [Fig. (c), J =.36]. Te magnetizationplateau at M=/3 is clearly observed in all tree plots, and it can tus be considered a universal feature of frustrated quantum spin cains. In accordance wit te criterion Eq. () for te plateau formation, te translational invariance of te Hamiltonian is broken on te plateau, leading to a state wit period p = 3 and an uud pattern in te onsite spin density Si z parallel to te external field. Te local spin densities corresponding to te plateau state are sown in Fig.. 5 As we use OBC suited for DMRG calculations, one of te tree possible patterns(up-up-down, up-down-up, down-up-up) is selected as a unique ground state due to energetically preferred orientations of edge spins, pointing up in our case. We furter observe deviations from te perfect up-down-up pattern, selected in our case, at te boundaries: edge-spin induced spin-spin correlations decay into te bulk. Let us briefly comment on additional features of te magnetization curves presented in Fig.. In te case of spin sown in Fig. (a), we identify at least six regions, labeled (i)-(vi) in te viewgrap. First, te Haldane gap,4,4 sows up as a zero-field plateau. A specialty of open spin- cains is te presence of effective spin-/ edge spins wic ave been studied in detail bot for unfrustrated 4 and frustrated cains. 4,5 Tese edge spins couple to te so-called Kennedy triplet, wic in te termodynamic limit is degenerate wit te ground state. 53 Te existence of te effective free spin-/ edge spins emerges naturally witin te AKLT (Affleck-Kennedy-Lieb-Tasaki) description of te ground state of a spin- cain. 4,54 Witin our numerical precision, tis triplet causes a zero-field magnetization of M = /N seen in te Haldane pase. Similar edge spin effects are present in all open spin-s > / cains (see, e.g., Ref. 5 and references terein): spin-3/ cains ave spin-/ edge spins wereas spin- cains ave spin- edge spins, wic consistently result in spurious M = /N and M = /N zero-field magnetizations seen in Figs. (b) and (c). At, te frustrated spin- cain enters a second region were te magnetization initially increases wit a steep slope but ten follows = b wit b.74 for.9. Te M=/3 plateau is realized for Note, owever, tat te pase boundaries of te plateau state cannot directly be obtained from N = sites for tis frustration parameter, as finite-size effects are still significant in tis regime. Bot below and above te plateau (indicated as region iv), te magnetization increases in steps corresponding to Stotal z =. Tis may indicate binding effects between te elementary excitations in finite magnetic fields. In te ig-field region below te saturation field, we igligt te presence of a in te magnetization curve at 3.89, wic correspond to jumps in te differential susceptibility d/d. Tis splits te ig-field regionintoparts(v)and(vi). Asecond, lessobvious may exist at smaller fields.8. Te singularities are related to middle-field cusp singularities recently studied in te literature. 47,55,56 We will discuss te singularities in more detail in Sec. III C. Apart from te M =/3 plateau, te prominent features of te magnetization curve of a frustrated spin- 3/ cain wit J = are te zero-field gap and several s. At zero magnetic field, te frustrated S=3/ cain is gapless for J.9 and in a critical pase wit antiferromagnetic correlations, wile a gap opens beyond tis value. 5 Tis transition at zero field is of te

5 5 /S 6 4 J =, = Spin Spin 3/ Eq. (7) M=/3 plateau Ferromagnet Haldane incommensurate Double-Haldane. J FIG. 3: (Color online) Magnetic pase diagram of frustrated S= and 3/ cains wit an isotropic excange ( = ) as derived from ED and DMRG calculations of magnetization curves. Squares denote boundaries of S= pases, wile stars are for S=3/. Anarrow at J,t.75 at low fields marks te first-order transition from te (gapped) Haldane pase found for J.75 to te likewise gapped double Haldane pase for te frustrated spin- cain. 59 Kosterlitz-Touless type and is similar to te beavior of te frustrated S=/ cain. 7,8 At larger J.9, te system enters a dimerized pase, wit a finite spin gap seen in Fig. (b). For te case of >, several singularities exist in te magnetization process of spin-3/ cains bot below and above te plateau, as is indicated in Fig. (b). Te derivative of te ground-state energies E (Stotal z, = ) wit respect to M exibits discontinuous canges of te slope were te s emerge in te magnetization curve, indicating at second-order transitions. Similar to te case of spin, te magnetization increases in steps corresponding to Stotal z = 3 = S around te plateau. Tis ints at a very ric pase diagram in a magnetic field tat deserves to be explored in more detail. Te magnetization curve of frustrated spin- cains at J =.36 exibits a zero-field Haldane gap, 45 similar to te case of spin, s at ig fields, and te M=/3 plateau. B. Magnetic pase diagram of frustrated spin- and 3/ cains wit isotropic excange interactions: Te M=/3 plateau region We next discuss ow te boundaries of te plateau region depend on spin S. First, we present results for te critical frustration for te plateau formation in Table I. Second, te pase boundaries are determined as a function of frustration J for spin and spin 3/ and tese results are summarized in Fig. 3. Finally, we present a discussion of additional pases present in te case of spin S= in a separate section, Sec. IIIC. J,t Since te plateau is induced by frustration, we expect it to vanis as J, as well as for J J. In te latter case, te system separates into two decoupled unfrustrated NN cains. For te extreme quantum case of S=/, DMRG studies 4,5 report te existence of te plateau for.56 J /J.5. On te oter and, as mentioned earlier, our MC simulations for a classical frustrated spin cain indicate te absence of te M =/3 plateauinteisotropiccase =,tusweexpectteregion tat allows for te uud state to become te ground state at M=/3 to srink wit increasing spin S. We terefore determine te critical lower and upper frustration J,crit,l/u (S) for te formation of te plateau state. To tis end, we compute te onsite spin density Si z for several J, searcing for te emergence of te uud pattern. We also perform a finite-size scaling analysis of te plateau boundaries and compare cains wit open and periodic boundary conditions. Tis procedure results in estimates of te critical frustrations for spin-/ consistent wit Ref. 4. Note tat te precise determination of te critical frustrations is somtimes difficult as te gap, i.e., te widt of te plateau, is typically exponentially small close to te critical frustrations. Our results are intended to study trends and we estimate te critical frustrations listed in Table I and II to be correct wit an error of ±.. Te results for S=, 3/, and are summarized in Table I, togeter wit data from Ref. 4 for spin /. Te data collected in Table I renders S=/ a special case wit critical frustrations very different from S > /. For S > /, te critical lower J,crit,l (S) remains approximately constant, wile J,crit,u (S) J,crit,l (S) sows a tendency to decrease wit increasing S, melting te plateau region from te large J side. We next elaborate more on te pase boundaries = (J ) of te plateau state in te -J plane for S= and 3/. Our results from DMRG calculations for several J are summarized in Fig. 3, wic sows te magnetic pase diagramof frustrated S= and 3/ cains at =. Additional magnetization curves for S= and J =.,.3,.7, and are sown in Fig. 4 and in Fig. 5 for S=3/ and J =.3,.5,, and. Close to te upper critical fields J,crit,u (S), te determination of te boundaries is ampered by te existence of additional steps on te plateau. An example is sown in te inset of Fig. 4(c) for J =.7, te steps are indicated by te arrow in te plot. Suc steps ave also been seen in te case of frustrated S=/ cains, 4 and can be traced back to te open boundary conditions. First, as we ave verified for several cases, in wic suc steps are observed, we find tat te step eigt scales down to zero wit increasing system size. Consistent wit te interpretation tat te steps are boundary effects, tey are not seen if periodic boundary conditions are imposed. Finally, we compare te onsite spin density Si z on te plateau to tat of te step states, i.e., te magnetizations just below and te two magnetizations just above te plateau. Te spin densities are sown in Fig. 6 for

6 6.. (a) Spin, J =, J =, = (b) J =.3 N=5, PBC (c) J =.7 (d) J = /3 step.5 M=/3 plateau step (a) Spin 3/, J =, J =.3, = (b) J =.5 (c) J = (d) J = N=, PBC N=6, OBC M=/3 plateau FIG. 4: (Color online) Magnetization curves for frustrated S=cains wit an isotropic excange ( = ) for (a) J =, (b) J =.3, (c) J =.7, and (d) J =. Full lines are DMRG results for N = sites, dased lines are ED results for N = 5. Te inset in panel (c) is an enlarged view of te M=/3 plateau. Te squares indicate te boundaries of te plateau as seen on N = 5 sites wit PBC. J =.7andN = sites. Onteplateau(S z total = 4), we find te up-down-up pattern. For S z total = 39,4,4, i.e., just above and below te plateau, one-domain wall excitation sits in te middle of te cain. 5 In te magnetic pase diagram Fig. 3, we scale te field wit te value of S. We see tat wile te upper and lower critical frustration J,(l,u) are similar for S= and 3/, te widt and te critical field, i.e., te boundaries of te plateau region exibit a significant dependence on spin S. In particular, we empasize te non-monotonic dependence of te plateau s widt on S it is quite narrow for S=/, 4 broader for S= and 3/, and it disappears in te classical case S =. C. Magnetic pase diagram of te frustrated spin- cain: Additional pases Finally, we present a discussion of additional pases present in te case of spin S=. To our knowledge only te magnetization curve of te NN spin- cain as been computed (see, e.g., Ref. 57). Based on a field-teoretical analysis, a first suggestion for te magnetic pase diagram of frustrated S= cains as been put forward in Ref. 47 for te low- and ig-field parts. Our results basically confirm teir picture and provide additional information on te so far unexplored middle-field region were te M=/3 plateau is found. Several studies of frustrated spin cains ave focused on te zero-field pases, 4,44 47,5,58,59 wit special interest in te emergence of cirality. As is well known, unfrustrated S= cains exibit a Haldane gap at zero magnetic field, in contrast to alf-integer spin cains tat ave gapless excitations.,5, As earlier studies FIG. 5: (Color online) Magnetization curves for frustrated S=3/ cains wit an isotropic excange ( = ) for (a) J =.3, (b) J =.5, (c) J =, and (d) J =. Full lines are DMRG results for N = 6 sites, dased lines are ED results for N = sites wit PBC. ave sown, frustrated S= cains ave a gap for any J. 4,6,6 Te Haldane pase persists up to J.75, and is ten followed by anoter gapped pase, te double-haldaneregion. 4 Tistransitionisaconsequence of te frustration driving te system from te unfrustrated one-cain (J ) to te two-cain limit (J =). Te existence of a spin gap gives rise to a zero-field plateau in te magnetization curve, wic is evident from our data sown for S= and J =.,.3,,.7, [Figs. 4(a), 4(b), (a), 4(c), 4(d)]. Te zero-field magnetization of M = S z total = /N due to te Kennedy triplet is seen in te Haldane pase, but disappears at te transition to te Double-Haldane pase. 4 We terefore estimate te spin gap from te field were te total S z total increases from S z total = () to S z total = () in te Haldane(Double-Haldane) pase. Our results for te spin gap are included in te magnetic pase diagram in Fig. 3 (squares) and are in agreement wit te data of Ref. 59. Wile te magnetization curve of bot te unfrustratedcaindepicted infig. 4(a)andofcainswit large J (not sown in te figures) are mostly featureless, a more interesting beavior arises in te range of parameters close to te plateau region. Around te plateau itself we find evidence for a narrow region were te magnetization increases in steps corresponding to S z total = : states wit an odd S z total never become ground states. Suc a region exists in frustrated S=/ cains as well, bot in its anti- and ferromagnetic versions. 4,63 Te corresponding state as been called te even-odd (EO) pase. 4 Wile te EOpaseforS=/ispredominantly realized at large J, we observe te EO effects only in a narrow region around J [see Fig.(a)]. Starting from te limit of two decoupled cains, te existence of te EO regions can be motivated in a simple picture. Wen J =, spins always flip in pairs, one on

7 7 <S z i > <S z i > (a) S=, J =.7, =, S z =39 (b) S z total =4 (c) S z total =4 (d) Sz total =4 5 site i 5 site i FIG. 6: (Color online) Onsite spin density S z i state for S=, J =.7, =, and N = sites (DMRG). Panels (a) (d): S z total=39 4. S z total = 4 corresponds to M=/3. eac cain. Te nontrivial question is weter a finite J preserves tis effect leading to pair-binding of excitations. Tis is indeed te case for te frustrated S=/ cain. 4,63 For spin, owever, Kolezuk and Vekua 47 conclude tat an EO pase is not possible in frustrated cains. Our DMRG results do not contradict tis prediction as we find EO effects in a region tat as not been studied in Ref. 47. To gain a qualitative understanding of te beavior in ig magnetic fields, it is instructive to consider te dispersion of a single magnon, i.e, one spin flipped wit respect to te fully polarized state. Te dispersion ǫ k, generalizing results given in Refs. 55 and 6, reads ǫ k = S[J cosk +J cosk (J +J )]. (6) Here, and trougout, k denotes te momentum. We also ave included te excange anisotropy [see Eq. (3)] for te sake of generality. Equation 6 first allows us to derive te saturation field, wic depends on S only by an overall prefactor: 6 sat = S[J J + (J +J )] for J J 4, [ ] J sat = S +J + (J +J ) for J > J 8J 4.(7) Te saturation field is plotted as a dased line in Fig.3. For = te saturation field is independent of J for J /4. For J.5, a singularity splits off te saturation field and divides te ig-field region into two pases. A second emerges at lower fields at J and ten approaces te first one. We observe te first ig-field for.5 J, were it ten merges wit te second and eventually disappears, at least on te system sizes studied ere. Just below te saturation field, were te one-magnon state is expected to give a correct description, we can understand its emergence along te lines of Ref. 55, were a similar effect as been observed for spin-/ cains. It turns out tat for J >.5 te global minimum of ǫ k at k = π splits into <S z i > <S z i >.. (a) J =, =, D=-.5 (c) D=. N=5, PBC N=6, OBC (b) D=. (d) D= FIG. 7: (Color online) Magnetization curves for te frustrated S= cain wit isotropic excange, = and J =, and a finite onsite anisotropy. (a) D =.5, (b) D =., (c) D =., (d) D =.5. DMRG results are for N = 6 sites (straigt lines); ED (dased lines) for N = 5 sites and PBC. two local minima at te position of te ig-field singularity. We can furter conclude tat te ig-field region between tis and te saturation field is incommensurate wit cosk = J /4J. Tis agrees wit Ref. 47 were tis part of te pase diagram as been suggested to exibit ciral order. We empasize tat te interpretation of te singularity based on te one-magnon dispersion Eq.(6) is justified only just below te saturation field. Indeed, even for te case of spin /, te nature of te two pases left and rigt of te singularity is currently under controversial discussion. According to Refs. 4 and 55, it signals a transition from a one-component Tomonaga-Luttinger liquid (TL) to a two-component Tomonaga-Luttinger liquid (TL) state. Tis interpretation as recently been callengedby Kolezukand Vekua. 47 Using bosonization tecniques, tey conclude tat te can be explained in terms of a renormalization of te magnetization due to an irrelevant operator in te incommensurate region neigboring te TL (or ciral pase) on te large J side. Furter studies are necessary to relate te tree featuresofte curves te plateau, singularities, and EO effects wit eac oter and sould involveadetailed study of excitations of S= and 3/ cains as well. Also, a numerical investigation of te ig-field pases including te possible existence of ciral order is a necessary step towards a complete understanding of te pases S= cains in finite magnetic fields. D. Magnetization process of frustrated S= cains wit onsite anisotropy Most materials tat realize S= cains typically also ave a significant onsite anisotropy D of te form..

8 8. J =., J =, =, MC results, PBC N=4 N=3 N= FIG. 8: (Color online) Magnetization curve of a frustrated classical cain ( = ) for J = obtained by means of MC simulations (circles). Te orizontal, dotted line marks M=/3. Te inset sows an enlarged view of te region around.75, illustrating te size dependence. D i (Sz i ) [see also Eq. ()]. We terefore study te effect of bot frustration J and a finite D on te magnetization process of frustrated S= cains. Te problem of a spin- cain wit frustration and onsite anisotropy D as previously attracted attention as well, owever, for a sligtly different model Hamiltonians tan our Eq. (), including biquadratic terms of te form (S i S j ). 5,37 Our DMRG and ED results for te magnetization process of frustrated S= cains wit a finite D are presented in Fig. 7 for fixed J = and =. A negative D does not affect te existence of te M=/3 plateau state. We ave verified te existence of te /3 plateau up to D =. Wile te widt of te plateau is at D = and J =, it increases steadily wit D <. At D =.5, we obtain.95. Tus, a negative D stabilizes te M=/3 plateau state. As an example, te magnetization curve for D =.5 is sown in Fig. 7(a). A positive onsite anisotropy D, owever, destroys te plateau state wic disappears for D. At larger D, an anomaly around M=/3 emerges, as is sown in Fig. 7(d) for te case of D =.5. Suc a beavior at finite D is expected since negative D favors te states wit a maximal projection S z i on eac site (suc as te uud state on te M=/3 plateau). On te contrary, a positive D plays te role of an easy-plane XY anisotropy and tus suppresses te plateau formation. Let us summarize te main results of tis Sec. III on te magnetization process of frustrated cains wit an isotropic excange. First, we ave presented numerical results for S=, 3/, and sowing tat a M=/3 plateau is realized in te magnetization curve of frustrated spin cains wit isotropic excange, wile a classical frustrated model does not support tis state. Second, we ave presented a magnetic pase diagram for frustrated S= and 3/ cains wit a special focus on te boundaries of te M =/3 plateau region. A tendency is found tat as S > / grows, te widt of te plateau becomes more narrow. Moreover, singularities exist in te magnetization curves of spin- cains for J /J.5, wic is at least close to te saturation field due to te emergence of two minima in te dispersion relation of te one-magnon excitations above te fully polarized state, similar to te case of S=/. 55 Tis effect furter indicates te presence of an incommensurate region for fields > above te ig-field singularity. Te inclusion of an onsite anisotropy term in te Hamiltonian of frustrated spin- cains stabilizes te M =/3-plateau wen te onsite anisotropy is negative (D < ), but te plateau is quickly destroyed by a positive D. Te magnetization curves of frustrated spin-3/ cains exibit a ric beavior wit several singularities and te M=/3 plateau. Around te plateau and at small J, te magnetization increases in steps of S = S. It is beyond te scope of tis work to fully map out te pase diagram of spin-3/ cains but interesting results are expected to emerge from future studies. IV. THE MAGNETIZATION PROCESS OF ANISOTROPIC FRUSTRATED SPIN CHAINS IN THE EASY-AXIS REGIME > In tis section we compare our results for frustrated spin-s cains wit S > / to a classical frustrated spin cain. Beyond a critical excange anisotropy c (J ), te classical frustrated cain supports a M =/3 plateau. Using linear spin-wave teory summarized in te Appendix, we derive several expressions for te pase boundaries of te classical model and in particular compare te pase boundaries of te plateau region to te cases of S=/, and3/cains. Te quantumcasesarestudied by means of ED and DMRG. We find tat, as S increases, te classical pase boundaries are rapidly approaced. Deviations between te quantum cases and te classical case are most significant close to te critical frustration for plateau formation. Te case of S=/ is clearly singled out: ere, te critical frustration and te pase boundaries are very different from spin S / cains and te classical model. Finally, we briefly discuss additional pases of frustrated spin- cains in te easy-axis regime >. A. Classical frustrated cain We perform MC simulations on cains of up to N = 36 sites wit periodic boundary conditions wic allow us to determine te ground-state at finite magnetic fields as outlined in Sec. II. Knowing te states gives us a andle on analytical expressions for teir energies. Te closely related problem of a classical Ising-like

9 9 /S J =, = Spin / Lines: Classical spin Spin Spin 3/ Eq. (7) CN N: C: canted DN: Ferromagnet P: M=/3 plateau J FIG. 9: (Color online) Magnetic pase diagram (field vs frustration J ) of a frustrated classical cain at =. Te pase boundaries of te classical model are given by Eqs. (a) (4). Tey are denoted by straigt lines. Te pases are labeled by capital letters (see Sec. IVB): N is a Néel pase, C is a canted pase wit periodicity tree, DN is te double-néel pase, P te M=/3 plateau state, CN is a canted Néel-like pase, and IC is an incommensurate region. Numerical estimates of te pase boundaries of te M=/3 plateau region for spin /,, and 3/ are denoted by circles, squares, stars, respectively. Dased lines are guides to te eye. Te dot-dased line is te saturation field sat from Eq. (7). Te dotted line scematically indicates te pase boundary between pases CN and IC, wic we ave not determined analytically. Heisenberg model on a triangular lattice was studied by Miyasita 4 and later on compared to S=/ Isinglike Heisenberg model on te same lattice. 64 Similar to te Ising limit and spin /, 38 te M=/3 plateau was found and results for te magnetic pase diagram were reported. Here we focus on te one-dimensional case only, but we allow J J in our study, wile Ref. 4 concentrated on equal excange constants on all sides of te triangles. As an example for our MC calculations, we sow te magnetization curve for J = and = in Fig. 8. At least five pases can be distinguised. First, a zerofield plateau exists, corresponding to a gap of.. Te ground state as an up-up-down-down (uudd) structure in te S z component, wile at small J, a Néel state (updown) is realized. LargerJ empasizeste double-cain caracter of our model. An uudd pattern results in Néel states on eac of te single cains in te limit of J. Next, separated by a first-order transition te system enters a region were te magnetization increases linearly wit. Tis pase as periodicity p = 3 and is a precursor of te plateau state: all spins lie in te same plane, one points exactly in te z direction wile te oter two are canted away from te +z direction by angles θ and θ suc tat IC cosθ = ( ) +. (8) + J +J /S J =, J = Ferromagnet P: M=/3 plateau Spin Spin 3/ Eq. (7) C H DN FIG. : (Color online) Magnetic pase diagram (field vs excange anisotropy ) of a frustrated classical cain at J = obtained by MC simulations. Te pase boundaries of te classical model are given by Eqs. (a) (4). Tey are denoted by straigt lines. Results for te pase boundaries of te M=/3 plateau region for spin and spin 3/ from DMRG calculations are included (squares: spin, stars: spin 3/). Te dased line is te saturation field sat from Eq. (7). Capital letters (C, DN,andP)refer topases sowninfig. 9 and are discussed in te text. H denotes a elical pase. 65 Te dot-dased line indicates te transition between pases H and C. Dotted lines are pase boundaries tat are only indicated scematically. As is increased, θ vanises continuously, so tat te transition into te plateau state is second order. Te M=/3 plateau exists for J + J = Te system ten enters a ig-field region wit a complicated dependence of M on. Witin te error of our simulations te pase transition to tis state from te plateau state is always second-order. Te ig-field state is probably incommensurate wit te lattice period. Tis is illustrated by te cain size dependence of te magnetization curve and its small anomalies sown in te inset of Fig. 8. Cains of finite lengt wit PBC cannot accommodate a truly incommensurate state. Tis picture is consistent wit our analytical spin-wave calculations. From Eq. (6) we obtain tat, for te present value of J = and = sat +, te fully saturated state is instable against fluctuations wit cosk = J /4J =.35. On te oter and, te plateau state is instable against te fluctuations wit cos(3k) = as sown in te Appendix. Tis dependence of k on empasizes te incommensurability of te ig-field pase. In tis work, we ave not attempted to study te k() dependence or te nature of tis state in furter detail. B. Te magnetic pase diagram of te classical model We derive analytical expressions for te boundaries of te plateau of te classical model as well as of several surroundingpasesasafunctionof andj. Teresults

10 are summarized in te magnetic pase diagrams in Figs. 9 and. Atlowfields,wefindtreepases: anéelstateatsmall J (region N in Fig. 9), te canted commensurate pase C, and te double-néel state DN for J 4 wit an uudd pattern. According to Eq. (8), te magnetization of te C pase at = is M( = ) = (cosθ ) = 3 3( +), (9) wic vanises only in te isotropic limit. Terefore, as a non-trivial consequence of te interplay between frustration and excange anisotropy, we observe a finite zero-field magnetization in a model wit purely antiferromagnetic interactions. Te boundaries between te commensurate states mentioned above are given by: N C = { (J +J )( )+ 3 } (J +J )(+ )[ (3J 5J ) (J +J )], (a) DN C = { (J +J )( )+ 3 } (J +J )(+ )[ (3J J ) (J +J )]. (b) For te lower boundaries of te M =/3 plateau we find: N P = (J J ); C P = J +J ; DN P = (J J ). (a) (b) (c) Note tat te critical field C P does not depend on. Above te plateau, we ave identified two pases. At small J, a canted Néel state is realized (labeled CN in Fig. 9), separated by a first-order transition line N CN = (J J )( )[ (J +J )+(J J )] () from te gapped Néel state N. Te P CN transition is also first-order. It occurs at P CN = 3 [J J + (J +J )+ J J (+ )+J ( )+J ( + )]. (3) For larger J and magnetizations M > /3, te situation becomes more complicated since incommensurability arises, as discussed above. For te range of parameters tat we study, te P-IC pase transition is always second order. Hence, te pase boundary can be derived from an analysis of spin-wave instabilities around te plateau state (see Appendix), P IC = J +J [ + (4 +4 7)]. (4) Equations (a) and (3) result in an expression for te miminal frustration J,crit ( ) required for te formation of te M=/3 plateau state in te vicinity of = : J,crit J = + ( ) +5. (5) Tis implies J,crit.79 for te classical frustrated cain at =. To conclude tis section, we present a cut troug te pase diagram at fixed J = in Fig., now plotting te magnetic field /S vs te excange anisotropy. Te double-néel pase DN vanises at.485, and below tis value of, a elical pase H emerges at zero magnetic field (see, e.g., Ref. 65 for a discussion of elical order in te isotropiccase). Te M=/3plateau only exists above a critical value of te excange anisotropy:.6. We obtain tis result by comparing te energy of te elical pase H, wic accordingto our MC simulations extends up to te plateau at small, to te energy of te plateau pase P. Note tat a similar plot as been presented in Ref. 4 for te triangular lattice wit equal couplings J on eac triangle. A pase C as not been reported in tat work: at J = J, we consistently also only find te double-néel pase below te plateau, separated by a first-order transition from te plateau [see Fig. 9]. We empasize tat we ave not yet fully explored all regions of te magnetic pase diagrams in Figs. 9 and of te classical frustrated cain. For instance, an analytical expression for te transition between pases CN and IC still needs to be found. Moreover, additional pases may arise in te limit J J. In tis parameter region,

11 .. (a) Spin, J =, J =, = (b) J =. N N=, OBC N=5, PBC (c) J = (d) J = C step jump S z total =S N DN step S z total =S S z total =S (a) Spin 3/, J =, J =, = (b) J =. N (c) J = (d) J = N=, PBC N=6, OBC step N DN P S z total =S= FIG. : (Color online) Magnetization curves of frustrated spin- cains wit an anisotropic excange ( = ) for (a) J =, (b) J =., (c) J =, and (d) J =. DMRG results (straigt lines) are for N = 6 sites, te dased lines are ED results (PBC). Te capital letters stand for: Néel pase N, canted pase C, double-néel pase DN. FIG. : (Color online) Magnetization curves for frustrated spin-3/ cains wit an anisotropic excange ( = ) for (a) J =, (b) J =., (c) J =, and (d) J =. DMRG results (straigt lines) are for N = 6 sites, te dased lines are ED results (PBC). Te capital letters stand for: Néel pase N, canted pase C, double-néel pase DN. we expect te plateau to vanis as te cains decouple. C. M=/3 plateau in te magnetization process of spin- and spin-3/ cains: Comparison to te classical case We now turn to te magnetization process of frustrated quantum spin cains wit S= and 3/ in te easy-axis regime >. Numerical results for are presented in Figs. and for J =,.,, and =. We briefly discuss te magnetization process and emergent pases starting wit te low-field part. Te focus ere is on to wat extent te S= and 3/ cains resemble te pases found for te classical cain. Results on te zerofield pase diagramoffrustratedspin- cainsand > can be found in Refs. 44 and 46. At a finite > and for spin-, essentially five pases ave been identified: a Néel state at small J and >.8, te Haldane pase, te Double Haldane pase and te Double-Néel pase for.95 on te large J side (see Ref. 44 and furter references terein). Our results reported below are in agreement wit tis picture. We now concentrate on te example of =. First, at small J, bot te S= and 3/ cains are in te Néel pase N, see for instance Figs. (a) and (a). Tis pase is gapped and terefore, a zero-field plateau exists in te magnetization curves. Te Néel pase terminates at J,N C.9 for S= and J,N C.39 for S=3/. For larger J, an intermediate region C follows, wic is gapped in te case of S= [see Fig. (c)]. Tis region must contain bot te Haldane and te Double- Haldane pase, wit a transition at J. 44 Te commensurabilty and te existence of a gap in te case of S=3/in tis regionstillneeds tobe clarified. Te widt in J of tis region at zero field becomes narrower as S grows, approacing te boundaries of te pase C of te classical model (see Fig. 9). At J,C DN. for S= and J,C DN.94 for S=3/, te systems undergoes anoter first-order transition and enters te double-néel pase DN, wic is gapped for bot values of S (see Figs. (d) and (d), J = ). At small J and >, te N pase terminates at a first-order transition, indicated by macroscopic jumps in te magnetization curves of bot S= and 3/ cains[see te case of J = in Figs. (a) and (a)]. Te critical frustrations for te formation of te M =/3 plateau state are J,crit for S= and J,crit. for S=3/. Interestingly, an EO region ( Stotal z = ) exists around te plateau at small J in te case of S= [see Fig. (b) and (c) for J =. and J = ] and extends far below te plateau. For S=3/, we find a corresponding region wit Stotal z = 3. An example is sown in Fig. (b) for J =.. We are tus led to conjecture tat frustrated spin-s cainsgenerallyrealizeregionswit Stotal z = S. We aveverifiedtat tisis te caseforspin S=aswell. Te M =/3 plateau becomes very broad wit increasing J, as is evident from te case of J = depicted in Figs. (c) and (c). Similar to te case of isotropic cains, boundary induced steps appear on te plateau. We ave not determined te upper critical frustrations for te plateau formation, but we did verify tat te plateau disappears for large J in bot te cases of S= and 3/. For magnetizations M > /3, we mention tat singularities are again found for J >.5, see te examples in Figs. (c) and (c). We proceed wit a comparison of te critical frustrations J of te M=/3 plateau state as a function of S, listed in Table II. For S > /, a clear trend towards te classical result is seen. Numerical results for te pase

12 Spin S J,crit/J / 5 3/...79 TABLE II: Easy-axis regime, = : Lower critical J,crit for te formation of te M=/3 plateau for spin /,, 3/, and (DMRG results) and te classical limit S. boundaries of te plateau state for bot S= (squares) and S=3/ (stars) are included in Figs. 9 ( vs J at = ) and ( vs at J = ). It turns out tat te pase boundaries in Fig. 9 converge rapidly wit increasing S towards te classical result. For instance, on te scale of Fig. 9, te classical result for te transition C P from te pase C to te plateau cannot be distinguised from te quantum cases. Te data points for S= on te boundary C P quantitatively deviate from te classical result, but te data points for S=3/ are already very close to te classical result for J. Deviations are strongest in bot te limit of small and large J. Wit respect to Fig. ( vs at J = ), we note tat sizable deviations between te plateau boundaries of te classical cain and te quantum cases only appear at small excange anisotropies close to =. Wile te spin- cain realizes te plateau in te isotropic limit =, te plateauonly opens in ofspin-3/cains beyond a critical < crit,3/. < crit,class. We can finally conclude tat te magnetic pase diagram of frustrated S= and 3/ cains in te easy-axis regime strongly resembles te pases found in te classical limit. Examples are te Néel pase N, te double- Néel pase DN, and most importantly, te M=/3 plateau state. In te latter case, te similarity is not only of qualitative nature, but te pase boundaries of te plateau rapidly approac te classical result wit increasing S. We empasize our findings on increases of te magnetization in steps corresponding to S total = S wic may indicate interesting binding effects between excitations. V. SUMMARY In tis work we ave studied te magnetization process of frustrated spin cains. Our numerical results using DMRG and ED tecniques reveal a complex and ric beavior of te magnetization process of S=, 3/, and cains. We find jumps, zero- and finite-field plateau states, as well as singularities. Our primary result is te M=/3 plateau state wit broken translational invariance and an up-up-down pattern in te spin component parallel to te field. Te existence of tis plateau indicates te presence of a gap in te spectrum of excitations at finite magnetic fields. We ave numerically estimated te pase boundaries of te plateau for S= and 3/ cains wit isotropic excange interactions. An additional onsite anisotropy term D(S z i ), necessary for te description of many Haldane materials, stabilizes te plateau, wen D is negative but rapidly destroys te plateau state wen D is positive. Te M=/3 plateau state can be followed down from te Ising-limit by decreasing te ratio /J for fixed S. For tat reason and for its Néel-type of order, it as been classified as a classical state, 4,6 but quantum fluctuations seem to stabilize it in te case of isotropic excange: a classical frustrated cain does not sow a plateau ere. We ave studied te classical frustrated spin cain in detail, wic realizes te M =/3 plateau in te easy-axis regime. Our work provides several analytical results for te magnetic pase diagram of te classical frustrated cain including te pase boundaries oftem=/3plateaufor >. Inteeasy-axisregime, te pase diagrams of frustrated S= and 3/ cains do not only qualitatively resemble te classical result, but also te pase boundaries of te M =/3 plateau rapidly approac te classical result wit increasing S. Te extreme quantum case of S=/ is singled out, te critical frustrations are different from S > / and te classical cases bot in te isotropic ( = ) and te easy-axis regime. Our conjecture is tat te plateau exists for all frustrated spin-s cains in te easy-axis regime, as our results for S=, 3/, and te classical case suggest. Our results on te emergence of te M=/3 plateau may be of relevance to te materials mentioned in te introduction tat are suggested to realize zigzag spin- cains: CaV O 5 (Ref. 48), NaR(WO 4 ) (Ref. 49), and Tl Ru O 7 (Ref. 5). Moreover, te M=/3 plateau can be expected to emerge in te magnetization process of all materials tat realize a spin-s Heisenberg model on a triangular lattice, as is well known for spin / (see, e.g., Ref. 66). Indeed, plateaux at bot M=/3 and M =/3 ave recently been reported for te S= bilinearbiquadratic Heisenberg model on te triangular lattice. 8 As by-products to te plateau study, and in particular close to te plateau region, we ave identified additional pases. We ave obtained te pase diagram of te frustrated spin- cain wit isotropic interactions, extending te results of oter studies. 47 We igligt anoter interesting finding: all frustrated spin-s cains seem to realize regions were te magnetization increases in steps corresponding to S total = S. Extensions of our work will comprise a full caracterization of different pases in terms of correlation functions and a more detailed study of excitations. A timely subject is te emergence of cirally ordered states in frustrated quantum magnets, wic is investigated bot experimentally 48,67 and teoretically ,47 We ope tat our work will stimulate future researc activities in tese directions. Acknowledgments - It is a pleasure to tank D.C. Cabra, A. Honecker, U. Scollwöck, and T. Vekua for