arXiv:2102.11143v1 [cond-mat.str-el] 22 Feb 2021 il o antzto ( magnetization for sible oeooia cteiglnt [10–14]. length scattering nomenological the Thus, field [4]. described magnetic theory physics Tomonoga-Luttinger-liquid low-energy the a by by have separated are that phases in phases gapped states The gapless edge [3]. associated chains spin exhibit open and Oshikawa-Yamanaka- criteria, topological [2] the Affleck incompress- satisfy these dimension, phases one ible In [1]. curves field netic sue omdleetosi togycreae ma- the in correlated investigated strongly high- particularly in [16], of was context electrons interaction and [15]. [17–21] model bond-charge systems terials to the quantum used of terms, variety is these a of of One modeling the in 9.I hsapoiain h nryi opie only comprised fermions is or energy 8] the a [7, approximation, by bosons this can hard-core In critical of and gas quantum regime [9]. a high-dilute the as a of treated in vicinity be are plateau the magnons the In the at points, points 6]. critical [5, quantum extremes with chain, spin pia atcs[93] ute,gnrlcreae hop- correlated in general Further, [26–28], species [29–36]. lattices polar Hamiltonian optical typically the a effective molecules, On with bosonic an a model for is in Hubbard [22–25]. hopping boson solution space extended density-dependent exact parameter the an hand, the has other of term point this special with model bard orcint hslwi ieradotie ytaking by obtained first and The linear transition. is account the law into of this side to gapless correction the in gularity oe.Floigti prah ecnso htthe that show can we approach, effective m the this in Following potential chemical model. of role the plays which lated ( h apdpae fmgei nuaosaerespon- are insulators magnetic of phases gapped The nte ido opn emi nw as known is term hopping of kind Another h h oeo est-eedn annhpigadmagnon-ma and hopping magnon density-dependent of role The uv mgo est)peet qaero sin- square-root a presents density) (magnon curve ) ipehpigterm hopping simple (or density-dependent 1 p onaymgo est togydpnso h interactio the mo on bosonic of depends the behavior strongly from density the and magnon spin in boundary and the accordance from model good plateau, spin the ferrimagnetic the verify the from we in results terms particular, the interaction between two matrix agreement these density tative of use effects We the energy. investigate kinetic many-magn the the decreases ho increases hopping repulsion hopp (density-dependent) magnon-magnon uniform correlated The with a sion. model have Rice-Mele interactions a the is Hamiltonian particle annmgo interaction magnon-magnon , h 3 ecmaetegon-tt etrso lentn ferrim alternating of features ground-state the compare We s/S / T nue unu hs rniin nthe in transitions phase quantum induces .INTRODUCTION I. 2 c , with , 2 uecnutvt.Bsds h Hub- the Besides, superconductivity. , 5 / ermgei spin-(1/2, ferrimagnetic namgei edadtecrepnigHolstein-Primako corresponding the and field magnetic a in 2 2 s aoa´ rod ´sc eoiaeCmuainl Departa Computacional, F´ısica Te´oricaLaborat´orio e de 1 m ertrad dc¸˜od aaıa 81-0 Jo˜ao Pess Educa¸c˜ao Para´ıba, 58015-900 de da Secretaria 1 = n h nfr emnterm, Zeeman uniform the and nvriaeFdrld enmuo 06-0 Recife-PE, 50760-901 Pernambuco, de Federal Universidade ltasi the in plateaus ) ) / ,cnieigteflyplrzdmgeiaina h bos the as magnetization polarized fully the considering 2, hopping .M aSilva da M. W. n sessential is and hog phe- a through m s mag- vs. ,2 1, corre- n .R Montenegro-Filho R. R. and S hisi antcfield magnetic a in chains ) ob pn(/,52 hi bv h three-dimensional the above temperature. chain ordering 5/2) spin-(1/2, a be to iaincreo h hretase at(4-Br- salt magne- charge-transfer full V)FeCl the of the Recently, curve tization [65–70]. chains ferrimagnetic s 5,6,7–8.Frhroe iersi-aetheory spin-wave the linear Furthermore, from 73–78]. ex- 60, formal- [59, magnetic Holstein-Primakoff ism low-energy the through their mainly characterize citations, classical to the used from was spin- [72] interacting theory and linear wave state, ground ordered ranged opud a hw ob oee yspin-(1/2, by bimetallic one- modeled of the be variety to side, a shown experimental was of the compounds phase On magnetic dimensional [59–64]. cell unit strto lta) where plateau), (saturation pt order to up field. magnetic a h apdadgpespae fspin-(1/2, of phases gapless and gapped the nesadtepyiso rsrtdisltn magnetic 37]. insulating [29–31, frustrated materials to of investigated physics the are understand Hamiltonians bosonic hard-core ping iae nteiorpc[04]adaiorpc[1 3 44] systems 43, Ferrimagnetic [41, anisotropic cases. and inves- [40–42] was isotropic the and in tigated and [39]. ferromagnetic orders exhibit long-range and antiferromagnetic [38] theorem Lieb-Mattis rudsaepaedarmo ermgei spin-(1/2, ferrimagnetic of diagram phase ground-state le antcfil rsn antzto ltasat plateaus magnetization present m field magnetic plied rto 5–8 otefriantcmdl.I partic- In spin-(1/2, models. ferrimagnetic ferrimagnetic the frus- ular, to geometric ob- adding [51–58] are or tration diagrams [46–50] phase doping rich through served Furthermore, [45]. magnetization tigated 1/3 AB the the with of associated plateau states edge the of ic spin-( Since eew hwta h osenPiao Hamiltonian Holstein-Primakoff the that show we Here h rudsaeo ermgei hisstsythe satisfy chains ferrimagnetic of ground-state The = S pn emadmgo-annrepul- magnon-magnon and term pping 4 ooi oe,addslytequanti- the display and model, bosonic − sadpril statistics.. particle and ns a xeietlyivsiae 7]adshown and [71] investigated experimentally was eomlzto ru acltosto calculations group renormalization ul oaie state polarized fully neeg n h density-dependent the and energy on es utemr,w hwta the that show we Furthermore, dels. 1 h desae,ascae ihthe with associated states, edge the n n oie onais while boundaries, modified and ing / h ulbsncapoiain In approximation. bosonic full the gei his(1 chains agnetic friantcpaeu and plateau) (ferrimagnetic 2 p s s/S , et eF´ısica, de mento S ooi oesu oorder to up models bosonic ff 2 ermgei hishv long- a have chains ferrimagnetic ) aP,Brasil oa-PB, nstoi hi a eetyinves- recently was chain anisotropic ie loteatyrslst the to results exactly almost gives nvcu.Tesingle- The vacuum. on 2 Brasil / m 2 nnrplinin repulsion gnon npriua,tebehavior the particular, In ˙ S , 6]gvsgo eut for results good gives [64] S stemgeiainper magnetization the is with ) hisudra ap- an under chains ) ermgei state ferrimagnetic S = m S = hisin chains ) o S -MePy- 1 + S / 2 ) 2
S) chains in a magnetic field. We use the density matrix and on b sites: renormalization group (DMRG) [79, 80] to obtain the magnetization curves, besides local properties, from the Sz = S nb; Holstein-Primakoff Hamiltonian and the spin Hamilto- j − j b 1/2 b nian. In addition to the simple hopping and the nearest- nj nj S− = √2Sb† 1 √2Sb† 1 , (4) neighbor interacting terms, a correlated hopping term is j j − 2S ≈ j − 4S essential to obtain a good accord between the numerical results from the spin and bosonic models. In Sec. II we present the Holstein-Primakoff Hamiltonian and the an- alytical formulas for the hard-core boson approximation. In Sec. III we compare the DMRG results for the magne- tization and local properties from the Holstein-Primakoff Hamiltonian and the spin model. In Sec. IV we present a summary of the main results of the manuscript.
II. HOLSTEIN-PRIMAKOFF BOSONIC HAMILTONIAN FROM THE FULLY POLARIZED VACUUM
The alternating mixed-spin (s = 1/2,S) chain with L unit cells has the Hamiltonian
L spin s S s S z = J ( j j + j j 1) BStot, (1) H · · − − j=1 X where
z z z Stot = (sj + Sj ) (2) j X is the z-component of the total spin, and we consider the magnetic field B in the z direction, with gµ 1, B ≡ where µB is the Bohr magneton and g is the g-factor. The spin 1/2 are attached to a sites, while spin S to b sites along− the chain, and we study chains for− which FIG. 1. (a) Schematic representation of the alternating spin S =1, 3/2, 2, 5/2, as schematically shown in Fig. 1(a). model, with spin-1/2 at a sites, and spin-S at b sites, with The ground-state total spin for B = 0 is S s and has S = 1, 3/2, 2, and 5/2. The z-direction is the direction of an − z a copy in each sector in the range (S s) Stot applied magnetic field B and the superexchange coupling is J. − − ≤ ≤ −1 2 (S s), as expected from the Lieb-Mattis theorem [38], (b) Holstein-Primakoff bosonic Hamiltonian up to (S / ), − O with a ferrimagnetic (FRI) long-range ordered state. If with a hopping term t = J√sS, local potentials εa = 2SJ − a little magnetic field is applied, the ground state with and εb = 2sJ, nearest-neighbor interaction V = J, and z − Stot = (S s) is chosen. Further, the ground state has density dependent correlated hopping process X = Jps/S. − z The a sites have a hard-core constraint, na 1, while the a finite gap to spin excitations carrying a spin ∆S = ≤ constraint nb 2S is imposed on b sites, with na (nb) as the +1, which induces a magnetization plateau at mF RI = ≤ (S s). Also, a second magnetization plateau is the fully number of bosons in one a (b) site. The magnetic field B acts polarized− plateau at m = S + s. as a chemical potential µ in the bosonic model: µ = B. F P (c) The bosonic approximations investigated: t, t V ,− and fully polarized state The is an exact ground state of t X V . In the t and t V approximations, there− is a hard- the spin Hamiltonian, and we build the spin-wave the- core− constraint− on all sites;− while in the t X V model, the ory considering it as the magnon vacuum. Making the b sites can accommodate up to two magnons.− − Holstein-Primakoff mapping on a sites
z a sj = s nj , − a b a 1/2 a where nj = a†aj and nj = b†bj , we arrive in the following nj nj j j s− = √2sa† 1 √2sa† 1 , (3) spin-wave Hamiltonian j j − 2s ≈ j − 4s 3
b a SW b a nj nj = J S n s n + √sS 1 b a† 1 + H − j − j − 4S j j − 4s j ( " X nb na na nb j j √ j j+1 + bj† 1 1 aj + sS 1 aj bj†+1 1 + − 4S − 4s # " − 4s − 4S a b nj nj+1 a b + aj† 1 1 bj+1 + s nj S nj+1 + − 4s − 4S # − − ) a b 1 B S + s n n + (S− ). (5) − − j − j O j X
Dropping the classical energy of the ferromagnetic state is a paradigmatic model to the understanding of topolog- Eclass =2JLsS B S+s L, the relevant magnon Hamil- ical insulators, and can be realized by atoms in optical tonian is: − lattices [82]. The model presents the bulk-boundary cor- SW 1 respondence [83], and an interacting version was recently H = H + H + H + (S− ), (6) t X V O investigated [83] to probe the connection between topol- with HSW = SW E . ogy and particle-particle interactions [83]. The Rice-Mele H − class As sketched in Fig. 1(b), the Ht term comprises a model recovers the bulk Hamiltonian Ht for t1 = t2. hopping process and distinct local potentials on a and b However, for an open chain, the mapping of the spin sites: model requires local potentials on the boundary sites, Eq. (9). Ht = t bj†aj+1 + bj†+1aj+1 + H. c. The second term in the bosonic Hamiltonian (6) is a j X density-dependent or correlated hopping term given by a b + [(εa µ)nj + (εb µ)nj ] (7) b − − HX = X (a† + a† )n bj + H. c. , (11) j − j j+1 j X j with X with t = J√sS, hopping parameter; J s εa = 2SJ, local potential on a sites; X = . (12) − (8) 4 S εb = 2sJ, local potential on b sites; r − a µ = B. Since s = 1/2 and nj 1, a hard-core constraint must − be imposed on a sites for≤ all models considered. Hence, In an open chain, if a or b is a boundary site, the local we discard a term similar to the X term but with a and b potential is half of the above value: J variables exchanged, with X′ = 4 S/s as the correlated (boundary) hopping amplitude. As sketched in Fig. 1(b), the energy εa = SJ of the system is lowered by the hoppingp of a magnon to a (boundary) − (9) (ε = sJ. site that is already occupied. In other words, the magnon b − probability to overcome the potential barrier between a Considering the local potentials, we see that the magnon and b sites increases if there is one magnon on the b site. has a higher probability to be found on a sites, and this This term becomes relevant for higher magnon densities, probability increases with S. However, since the hopping since it is an interaction term, and X 0 as S increases. parameter t increases with √S, quantum fluctuations are The last term in the Hamiltonian→ (6), see Fig. 1(b), relevant for moderate values of S, and the magnons can is a repulsive term between magnons in nearest neighbor overcome the potential barrier between a and b sites. sites: We observe that the bulk Hamiltonian Ht is a partic- a b b ular case of the Rice-Mele model [81]: HV = V nj+1(nj + nj+1), (13) Rice-Mele j H = t2bj†aj+1 + t1bj†+1aj+1 + H. c. X j with V = J, and increases the energy for higher magnon X a b densities. For one magnon per unit cell, this term favors + (εanj + εbnj ), (10) the magnon localization on alternating a sites, since the j X local potentials (ε) favors the occupation of a sites. Thus, putting µ = 0 and considering the general case of al- ε (any density) and V (higher densities) favor magnon lo- ternating hopping: t1 (t2) for intra-cell (inter-cell) hop- calization on a sites, while quantum fluctuations (tunnel- ping. The Rice-Mele Hamiltonian was originally pro- ing between a and b sites) are favored by t (any density) posed to model the physics of electrons in polymers, but and X (higher densities). 4
In this work, we compare data from three approxima- with tions of Hamiltonian (6), as summarized in Fig. 1(c), 2 and from the spin Hamiltonian (1). The first, which we εb εa 2 ωk = − + γ identify as the t approximation, is an analytical solution s 2 k to the free hard-core model. In this approximation, all = J (S s)2 +4sS cos2(k/2), (19) sites have a hard-core constraint and can be occupied by − a b only one magnon: n 1 and n 1, for any j. This (±) (±)† (±) j j and n = α αp, where constraint implies that≤ there is not≤ an energy contribu- k k k tion from the X term, and we drop the nearest-neighbor (−) αk cos θk sin θk ak interaction V . In the second approximation, t V , we (+) = − , (20) α sin θk cos θk bk keep the hard-core constraint but add the V contribution− k to Ht. The last approximation, t X V , has the three with − − terms Ht, HX , and HV . We consider a hard-core con- 2 1 S s b cos θk = + − ; straint on a sites and a constraint nj 2 on b sites. As 2 2ωk ≤ 2 1 S s (21) we present below, the relaxation of the hard-core con- sin θ = − . ( k 2 2ωk straint on b sites and the consequent activation of the − correlated hopping term implies an excellent agreement A magnetic field (or chemical potential) adds a +B en- between the results of this approximation and the spin ergy term to the two bands. model. 0 2s = 1
vacuum A. Free hard-core magnons [t-approximation]: H -2s = -1 t inactive −µ FP and L 0 →∞ -1 -0.5 0 0.5 1 one - magnon k / π states B=2(s+S)= 0
The first approximation to the magnetization curve, a Energy / J one magnon many-magnon state, is to consider the magnons as free -2S per unit cell hard-core bosons or free fermions. This mapping is exact B=0 B=2S in the high-dilute regime of magnons, near the saturation FRI -2(S+s) -2s = -1 field. Here, we extend this approach for the full range of -1 0 1 -1 -0.5 0 0.5 1 π π the magnetization, (S s) m (S + s), and compare k / k / their results to more precise− ≤ calculations≤ considering the FIG. 2. One-magnon states from the fully polarized (FP) interaction terms. vacuum and B = 0 as a function of lattice wave-vector k. The single-magnon energies are given by the term Ht, In the t-approximation, the high energy band is inactive. For Eq. (7), in Hamiltonian (6). Using the following Fourier B = 2(s+S)= µ, the exact ground state is the FP vacuum; transforms: while there is one− magnon per unit cell in the system for B = 2S, and the ground state is ferrimagnetic (FRI). 1 ik/4 ikj aj = e− e ak; √L k (14) b = 1 e+ik/4eikj b , The magnon bands are shown in Fig. 2 for (i) B = 0; ( j √L Pk k (ii) B = 2(s + S), the saturation field, for which m = P ik/4 mFP = (S + s); and (iii) B =2S, the critical field of the where a phase e± is included to ease the calculation, ferrimagnetic plateau, for which m = mFRI = (S s), or Ht becomes: − 1 magnon per unit cell: mFP mFRI = (S +1/2) (S 1/2) = 1. As expected from− the Lieb-Mattis theorem,− − εa γk ak Ht = a† b† , (15) the magnetization in the null field is S s. In the free k k γk εb bk k − X hard-core boson or free fermion approximation, we fill the single-particle states following a Fermi distribution up to where the effective Fermi wave-vector kF . As shown in Fig. 2, if we follow this procedure the two bands should be filled γ =2t cos(k/2). (16) k for B = 0 (two magnons per unit cell) and the magnetiza- tion curves of the spin model would not be reproduced. After diagonalization, the Hamiltonian is written as We have shown in Ref. [64] that this problem can be
(−) (−) (+) (+) overcome, even for finite T , by introducing an effective Ht = ωk nk + ωk nk , (17) chemical potential µ to the upper band, in a way similar Xk to Takahashi’s solution to the ferromagnetic linear chain where the dispersion relations are [84]. In particular, µ J = 2sJ as T 0, such that the overall effect of→µ −at T −= 0 is the suppression→ ε + ε of the upper band. Hence, in the free hard-core approx- ω(±) = a b ω = J(s + S) ω , (18) k 2 ± k − ± k imation, we must consider only the lower band in the 5 calculations, as schematically indicated in Fig. 2. Thus, state for S = 3/2, 2, and 5/2. In these cases, the global for example, the energy per unit cell in the free hard-core energy minimum is reached by the algorithm if the chains approximation is written as have a b site at the left extreme and an a site at the right extreme (b-a boundaries). In the other magnetizations, k E 1 F this is irrelevant, i. e., the same state is calculated for = ω(−) (22) L L k the a-b or the b-a boundaries. We retain a maximum of k= kF X− 243 states per block and the maximum discarded weight less than 10 9. for a magnon density per unit cell n, where − For the spin model, the magnetization curves are cal- z k = πn, (23) culated from the lowest energy state for a fixed S and F B = 0: E(Sz); since for B = 0 we need only to add the Zeeman term, such that E6 (Sz)= E(Sz) BSz. In and n = mFP m. B − a gapless (non-plateau) phase, the magnetization− curves are made of finite steps in a finite-size system. Defining 1. Average local magnetizations the extreme points of these finite steps as B and B+, we have − If the chain has a magnon density per unit cell n B = [E(Sz 1) E(Sz)] (27) ± ± ± − and considering the hard-core approximation, the aver- z age magnetizations of a and b sites are given by for the step at S . In a gapless phase, B B+ as L , while in a thermodynamic-limit plateau− → state → ∞ z 1 kF a B = B+ for L . In the last case, B are quantum s = s n (a sites); − ± L k= kF k 6 → ∞ h i − − h i (24) critical fields separating the plateau insulating state from Sz = S 1 kF nb (b sites). ( L k= kF k a gapless critical Luttinger liquid phase. For the bosonic h i − P − h i P models, the magnon density per unit cell n as a function Using Eqs. (20) and discarding terms involving the upper of chemical potential µ is obtained with a similar pro- band, we obtain cedure. We calculate the lowest energy state for a fixed number of bosons N, with N = nL and µ = 0. The value a (−) 2 n = n cos θk; of the chemical potential at the extremes (µ and µ ) of h ki k (25) + nb = n(−) sin2 θ . the step at N are given by − (h ki k k µ = [E(N 1) E(N)]. (28) We, thus, have ± ± ± − A gapless phase has µ+ µ as L , while in a z 1 S s kF 1 → − → ∞ s = s ; plateau insulating phase µ+ = µ in the thermodynamic 2 −2 k= kF ω − − − k (26) 6 − Sz = S 1 + S s kF 1 , limit. The transformation from the boson to the spin ( 2 −2 k= kF ω − P − k language is performed through the following equations: after the substitution of Eqs. (21)P in Eqs. (25) and the n = m mFP, and results in Eqs. (24). − (29) B = µ. ( −
III. DMRG RESULTS FOR THE SPIN MODEL AND THE BOSONIC HAMILTONIANS B. Magnetization curves and local magnetizations
A. Methodology The magnetization curves from the spin model and the boson Hamiltonians are shown in Fig. 3(a). The models present two magnetization plateaus, one at the ferrimag- We use the density matrix renormalization group to netic magnetization mFRI = S 1/2 and the fully polar- obtain the magnetization curves and local properties of − ized state mF P = S +1/2. The saturation field B , end the spin, t V , and t X V models; and also com- F P point of the FP plateau, is obtained through the closing pare this data− with the− analytical− results from the free of the single-particle magnon gap with B, as sketched in hard-core model (t-model) in the thermodynamic limit: Fig. 2, at B = BF P = 2(S + s). The saturation field L . These approximations are summarized in Fig. BF P from any bosonic approximation is rigorously equal 1(c).→ ∞All DMRG results (boson and spin models) are ob- to its exact value since the fully polarized state, an ex- tained through the Algorithms and Libraries for Physics act eigenstate of the spin model, is the vacuum for the Simulations (ALPS) project [85] for chains with L = 128 bosonic models. unit cells, with one a site at one extreme and a b site The Oshikawa-Yamanaka-Affleck topological criteria at the other. If the system has an a site at the left [2] states that a plateau can appear in the magnetiza- extreme and a b in the right (a-b boundaries), the renor- tion curve of spin systems if malization process for the magnetization step inside the ferrimagnetic plateau becomes trapped in a metastable S m = integer, (30) u − 6
(a) 1 the magnetization curves also exhibits a plateau at n = 1. 0.9 We also observe the occupancy of the edge states of the ferrimagnetic plateau by one magnon at the indicated 0.8 magnetic fields. These edge states appears in finite-size
s 0.7 (1/2, 5/2)FRI plateau open systems associated with topological aspects [3, 45] (1/2, 2) FRI plateau of the ferrimagnetic state. m / 0.6 The non-interacting Rice-Mele model (10) does not (1/2, 3/2) FRI plateau Spin, DMRG (L=128) present edge states for uniform hopping. However, re- 0.5 →∞ t, L cently [83], it was shown that an interacting fermionic 0.4 FRI t - X - V, DMRG (L=128) (1/2, 1) plateau arrows edge states system, with a local Coulomb interaction, presents ef- 0.3 fective edge states, and that a fraction of the boundary 0 1 2 3 4 5 6 B / J charge is, in fact, related to the bulk properties. Our (b) 1 3/2 2 5/2 non-interacting Hamiltonian (7) has the modified local 1 potentials in the boundaries (9), required from the spin FRI plateau (1/2, 5/2) (n = 1) mapping, and thus localized edge states. Furthermore, our interacting model is a bosonic system having the cor-
Gapless (1/2, 2) related hopping and nearest-neighbor coulomb repulsion. Luttinger n In the Sec. IV, we study the boundary magnon den- liquid 0.5 (1/2, 3/2) sity and compare it from relevant interacting and non- interacting bosonic models. Spin, DMRG (L=128) t, L →∞ (1/2, 1) The gapless phase between the magnetization plateaus t - V, DMRG (L=128) is a Luttinger liquid phase with power-law decay of the t - X - V, DMRG (L=128) 0 transverse spin correlation functions [4] and has a dy- 0 1 2 3 4 5 6 namical exponent z = 1. B / J We notice in Fig. 3(a) that as B decreases from FIG. 3. (a) Magnetization per unit cell (m) of (1/2,S) chains, BF P , the magnetization from the free hard-core model, with S = 1, 3/2, 2, and 5/2, normalized by its saturation t-approximation, departs from that of the spin model (−) value (ms) and (b) magnon densities per unit cell n as func- at roughly half filling of the lower magnon band ω , tions of B in units of J. DMRG data for (full lines) the spin Eq. (18). At this filling, the interaction effects start to model and for ( ) the t X V approximation. We present • − − become relevant. The critical field of the ferrimagnetic in both figures the results for the t-model (dashed lines) and plateau from the free hard-core model: 2S, see Fig. 2, L . In (b), we also show the magnon density as a func- → ∞ becomes more near its value for the spin model as S in- tion of B units of J for ( ) the t V approximation, calculated creases, as also confirmed in Table I. This effect can be with DMRG for a system• with−L = 128. In the thermody- attributed to the local energy of the a sites that becomes namic limit, the ferrimagnetic plateau is observed for n = 1, while the gapless Luttinger liquid phase for 0 0.5 1 (a) Spin, DMRG(L=128) (b) TABLE I. Critical field of the ferrimagnetic plateau in units t, L→∞ (1/2, 5/2) of J for the spin-(1/2, S) chains in the free hard-core boson t - X - V, DMRG (L=128) 0.25 (1/2, 2) approximation t-approximation and the spin model. − 0.9 2S β > a z (1/2, 3/2) S t-approx.: 2S spin model: β 0 > / S b β z (1/2, 1) 1 2 1.76 0.13 (a) (1/2, 1) chain (b) (1/2, 5/2) chain (c) edge states 0.08 spin: cell n = 8/128 = 0.0625 0.08 spin: cell n = 8/128 = 0.0625 1 (1/2, 1) 0.06 0.06 0.8 0.6 cell: spin 0.04 0.04 spin: a site a site: spin spin: a site t - X - V 0.4 0.02 spin: b site FRI b site: spin 0.02 → spin: b site t - X - V 0.2 0 0 0 0.15 0.15 FRI+1 1 0.8 0.1 0.1 δ 0.6 (1/2, 3/2) n = 16/128 = 0.12500 n = 16/128 = 0.12500 0.05 0.05 0.4 0.2 0 0 0 1 0.5 0.5 0.8 0.6 (1/2, 2) n = 72/128 = 0.56250 n = 72/128 = 0.56250 0.4 0.2 0 0 0 1 0.8 0.8 Average magnon densities along the chain 0.8 0.6 0.6 0.4 n = 116/128 = 0.90625 0.4 n = 116/128 = 0.90625 Density change along the chain 0.6 (1/2, 5/2) 0.2 0.2 0.4 0 0 0.2 0 64 128 0 64 128 0 l l 0 l 10 FIG. 5. Average magnon densities along (a) the (1/2, S = 1) chain and (b) the (1/2, S = 5/2) with L = 128 unit cells. Magnon densities at ( ) a sites, na,l , at ( ) b sites, nb,l , and ( ) total magnon density, ncell,l , for unit cell l from the spin model and (N) the t• X V happroximation,i • bothh calculatedi with• DMRG for systems withh Li = 128. The following total average magnon densities− are− shown: n = 8/128, 16/128, 72/128, and 116/128, from top to bottom, corresponding, respectively, to the magnetizations (a) m = 184/128, 176/128, 120/128, and 76/128; and (b) m = 376/128, 368/128, 312/128, and 268/128. (c) Average magnon distribution change δ n F RI+1→FRI,l = n N=L,l n N=L−1,l along the chain as a function of unit cell l shows the presence of the edge states forh thei four chains studh iied. (d)−h Thei average magnon probability density n along a chain with L = 128 at (circles) a and (diamonds) b sites for (white symbols) N = L 1 and (red symbols) N = L bosonsh i from the t X V approximation of the (1/2,1) chain. (e) Sketch of the ground states for− N = L 1 and N = L bosons. Average magnon− densities− are shown as grey circles, fluctuations between a and b sites are indicated by− a blue stripe, and we use a red filled curve to represent the localized orbital in the boundary unit cell. For the two higher magnon densities, there is an excel- alize the spatial extent of the edge state, we consider the lent agreement between the spin model and the t X V magnon distribution change, δ n F RI+1 F RI , between a approximation for the two chains. − − total number of magnons N =hLi 1 [Sz→= L(S s) + 1] and N = L [Sz = L(S s)]: − − In Fig. 5(c), we show the magnon distribution in − the edge localized state occupied by one magnon for the δ n F RI+1 FRI,l = n N=L,l n N=L 1,l. (33) four chains studied. To calculate it, we notice that the h i → h i − h i − edge state appears between the two magnetization steps: We notice in the data shown in Fig. 5(c) that a tiny dis- Sz = L(S s) and Sz = L(S s) + 1, which will join crepancy is observed between the results for the t X V - in the thermodynamic− limit, and− make the S s ferri- approximation and the spin model in the case of the− (1/2,− magnetic plateau, see Fig. 3(a) and (b). Thus,− to visu- 1) chain, while for the other chains the agreement is excel- 9 ≠ lent. The hole added to the many-magnon state at N = L (a) Bound. pot. Bulk pot. (b) Bound. pot. = Bulk pot. 138 138 becomes well localized in the boundary cell, mainly on t - X - V t - X - V 136 t - X - V 136 t - X - V the boundary a site, thus implying the presence of this 134 with 134 with z X=0, V=0 X=0, V=0 empty orbital in the N = L 1 [S = L(S s)] state. 132 t, finite L 132 t, finite L − − The localization in the boundary increases with S as ex- N 130 130 pected from the increasing of the bulk gap. In fact, the 128 128 126 B = 0.06 J 126 fluctuations of the magnon density in the boundary be- 124 124 come neglegible as S , and the boundary hosts one 122 122 magnon in this limit.→ We ∞ can attibute this behavior to 0 0.5 1 1.5 2 2.5 0 0.5 1 1.5 2 2.5 B / J B / J the increasing potential barrier between bulk and edge, (c) Bound. ≠ Bulk (d) Bound. = Bulk Eqs. (8) and (9), respectively, as S . 1 1 → ∞ As an example of the data used to make the results 0.8 0.1 0.8 l , shown in Fig. 5(d), we exhibit in Fig. 5(d) the average f N → 0.6 0.6 magnon density along the chain, n , for a and b sites, i h i N in the t X V approximation of the (1/2,1) chain. 0.4 0.4 − − is shown in Figs. 6(c) and (d), respectively, where Nf With the help of Fig. 6, we discuss the relevance of the (Ni) is the highest (lowest) value of N between the two difference between the local potential in the boundaries magnon density steps separated by the occupancy of the and the bulk sites on the boundary magnon densities, edge state. as well as the importance of interactions in it. We use For the bulk different from the boundaries, Fig. 6(a), the parameters of the spin-(1/2, 1) chain, but focus on we see that the edge state is very robust, appearing even the bosonic Hamiltonian, for three models in finite-size in the t approximation. This does not occur for the t chains: (1) the interacting model considering the map- approximation if the boundary is equal to the bulk, Fig. ping to the spin system (X =0, V = 0) - the t X V 6(b), which is an expected result for the Rice-Mele model model; (2) the corresponding6 non-interacting6 model,− i.− e., [83]. However, we notice in Fig. 6(a) that, in comparison the t X V model with X = 0 and V = 0, which has with the spin model, the plateau width is very underesti- a hard-core− − constraint on a sites and the b sites can host mated in the t approximation, since we are not discarding up to 2 bosons; and (3), the t approximation, which has the single-particle states in the upper band shown in Fig. a hard-core constraint on a and b sites. One of the sets, 2. We observe that this upper band was discarded for the Figs. 6(a) and (c), shows data for the boundary poten- L results shown in Figs. 3 and 4, due to the in- tials distinct from the bulk potentials, with the values in troduction→ ∞ of the effective chemical potential connected Eqs. (9) and (8), respectively. The other set, Figs. 6(b) with the Takhashi’s constraint [64]. Considering also the 10 spin model, the error in the plateau width for t X V densities favor magnon localization on the a sites. On interacting case is very tiny, with a lower critical− field− the other hand, quantum fluctuations, magnon tunneling 0.06J. Also in Fig. 6(a), bulk and boundaries dis- between a and b sites, are favored by t, at any density, ≈tinct, we note that the edge state in the non-interacting and X for higher densities. We use the density matrix t X V model appears only for two bosons over the renormalization group to investigate the spin model and N−= 127− state, Fig. 6(a), with one of the bosons oc- the bosonic Hamiltonians t V and t X V in finite- cupying an extended state and the other an edge state, size open systems, while for− the t-approximation− − we have as shown in Fig. 6(c). This shows that the distinction considered its analytical solution. We compare the mag- between boundaries and bulk is not sufficient to a good netization and magnon densities per unit cell as a func- representation of the edge state in the spin model, the tion of a magnetic field, average bulk density and local interactions are essential for it. In the case of the t ap- densities along the chains from the spin model and the proximation, the particle statistics provides a change in bosonic approximations. From the ferrimagnetic plateau the single particle quantum states that mimic the effect (one magnon per unit cell) to saturation (empty chain), of the X and V interactions in the t X V model. In the t X V -approximation is in excellent agreement Fig. 6 (b), where the boundary is equal− to− the bulk, we with the− spin− model, while the t and t V results de- show that edge states appear in the t X V model in part from that of the spin model as the magnon− density the interacting and non-interacting cases.− As− mentioned increases. This, thus, shows the relevance of both interac- above, the edge state of the non-interacting case does tion terms, magnon-magnon repulsion and the correlated not appear in the t approximation. Further, the bound- hopping term X, and its associated particle fluctuations, ary density, Fig. 6 (d), in the interacting t X V case to describe the many-body system near the ferrimagnetic appears in the right extreme and has a value− − approxi- plateau. The edge state associated with the insulating mately equal to that in Fig. 6 (c). From these results, ferrimagnetic plateau is well reproduced by the bosonic we can assert that the interactions X and V , and the par- model. In particular, we have shown that the magnon ticle statistics, are essential to understand the boundary boundary densities are strongly dependent on the inter- magnon density. actions and the particle statistics. The use of the fully polarized state as the vacuum en- ables a better understanding of the underlying quantum V. SUMMARY AND DISCUSSION processes in the spin Hamiltonian, compared to a ferri- magnetic vacuum, which is not an exact eigenstate of the We have investigated the Holstein-Primakoff bosonic spin Hamiltonian. Our results also suggest that, beyond Hamiltonian, up to order s/S, of the alternating spin- hopping and magnon-magnon interaction, the density- (s = 1/2, S) chains in a magnetic field, and considering dependent hopping term increases the range of magne- the fully polarized as the vacuum.p Three bosonic Hamil- tizations for which effective bosonic models can make a tonians were considered: the first has only a hopping good description of the physical data from general spin term (t-approximation) and distinct local potentials (ε) systems having ions with spins higher than 1/2. on the spin-1/2 (a site) and the spin-S (b site) sites , this approximation is a Rice-Mele Hamiltonian with bound- aries different from the bulk and uniform hopping term. ACKNOWLEDGMENTS The second one has the terms of the t-approximation plus a magnon-magnon repulsion V (t-V -approximation); and We acknowledge support from Coordena¸c˜ao de Aper- the last approximation shows the t and V terms, and a fei¸coamento de Pessoal de N´ıvel Superior (CAPES), Con- density-dependent correlated hopping term X (t X V - selho Nacional de Desenvolvimento Cient´ıfico e Tec- approximation). In the t X V approximation,− b −sites nol´ogico (CNPq), and Funda¸c˜ao de Amparo `aCiˆencia e can accommodate up to two− magnons,− while in the oth- Tecnologia do Estado de Pernambuco (FACEPE), Brazil- ers a hard-core constraint is imposed on a and b sites. ian agencies, including the PRONEX Program which is The local potentials, at any density, and V for higher funded by CNPq and FACEPE, APQ-0602-1.05/14. [1] T. Giamarchi, C. R¨uegg, and O. Tchernyshyov, B 90, 035150 (2014); H.-P. Hu, C. Cheng, H.-G. Luo, Bose–Einstein condensation in magnetic insulators, Na- and S. Chen, Topological incommensurate magnetization ture Physics 4, 198 (2008). plateaus in quasi-periodic quantum spin chains, Scientific [2] M. Oshikawa, M. Yamanaka, and I. Affleck, Magneti- Reports 5, 8433 (2015). zation plateaus in spin chains: “haldane gap” for half- [4] T. Giamarchi, Quantum Physics in One Dimension (Ox- integer spins, Physical Review Letters 78, 1984 (1997). ford University Press, 2004); R. Chitra and T. Gia- [3] H. Hu, C. Cheng, Z. Xu, H.-G. Luo, and S. Chen, marchi, Critical properties of gapped spin-chains and lad- Topological nature of magnetization plateaus in period- ders in a magnetic field, Physical Review B 55, 5816 ically modulated quantum spin chains, Physical Review (1997). 11 [5] S. Sachdev, Quantum Phase Transitions (Cambridge [25] C. Vitoriano and M. D. Coutinho-Filho, Fractional University Press, 2001). Statistics and Quantum Scaling Properties of the Hub- [6] M. Vojta, Quantum phase transitions, Reports on bard Chain with Bond-Charge Interaction, Physical Re- Progress in Physics 66, 2069 (2003). view Letters 102, 146404 (2009). [7] I. Affleck, Bose condensation in quasi-one-dimensional [26] M. Maik, P. Hauke, O. Dutta, M. Lewenstein, and J. Za- antiferromagnets in strong fields, Physical Review B 43, krzewski, Density-dependent tunneling in the extended 3215 (1991). Bose–Hubbard model, New Journal of Physics 15, 113041 [8] E. S. Sørensen and I. Affleck, Large-scale numerical ev- (2013). idence for Bose condensation in the S =1 antiferromag- [27] S. Baier, M. J. Mark, D. Petter, K. Aikawa, L. Chomaz, netic chain in a strong field, Physical Review Letters 71, Z. Cai, M. Baranov, P. Zoller, and F. Ferlaino, Extended 1633 (1993). Bose-Hubbard models with ultracold magnetic atoms, [9] A. M. Tsvelik, Field-theory treatment of the Heisenberg Science 352, 201 (2016). spin-1 chain, Physical Review B 42, 10499 (1990). [28] S. Fazzini, A. Montorsi, M. Roncaglia, and L. Barbi- [10] K. Okunishi, Y. Hieida, and Y. Akutsu, δ-function Bose- ero, Hidden magnetism in periodically modulated one di- gas picture of S = 1 antiferromagnetic quantum spin mensional dipolar fermions, New Journal of Physics 19, chains near critical fields, Physical Review B 59, 6806 123008 (2017). (1999). [29] R. Bendjama, B. Kumar, and F. Mila, Absence of single- [11] J. Lou, S. Qin, T.-K. Ng, Z. Su, and I. Affleck, Finite-size particle bose-einstein condensation at low densities for spectrum, magnon interactions, and magnetization of S bosons with correlated hopping, Physical Review Letters = 1 Heisenberg spin chains, Physical Review B 62, 3786 95, 1 (2005). (2000). [30] K. P. Schmidt, J. Dorier, A. L¨auchli, and F. Mila, Single- [12] I. Affleck, W. Hofstetter, D. R. Nelson, and particle versus pair condensation of hard-core bosons U. Schollw¨ock, Non-Hermitian Luttinger liquids and with correlated hopping, Physical Review B 74, 174508 flux line pinning in planar superconductors, Journal of (2006). Statistical Mechanics: Theory and Experiment 2004, [31] K. P. Schmidt, J. Dorier, A. M. L¨auchli, and F. Mila, Su- P10003 (2004). persolid Phase Induced by Correlated Hopping in Spin- [13] I. Affleck, Luttinger liquid parameter for the spin-1 1/2 Frustrated Quantum Magnets, Physical Review Let- Heisenberg chain in a magnetic field, Physical Review ters 100, 090401 (2008). B 72, 132414 (2005). [32] M. Eckholt and J. J. Garc´ıa-Ripoll, Correlated hopping [14] L. Vanderstraeten, F. Verstraete, and J. Haegeman, Scat- of bosonic atoms induced by optical lattices, New Journal tering particles in quantum spin chains, Physical Review of Physics 11, 093028 (2009). B 92, 125136 (2015). [33] H. C. Jiang, L. Fu, and C. Xu, Pair superfluid and su- [15] A. Dobry and A. Aligia, Quantum phase diagram of the persolid of correlated hard-core bosons on a triangular half filled Hubbard model with bond-charge interaction, lattice, Physical Review B 86, 045129 (2012). Nuclear Physics B 843, 767 (2011). [34] O. J¨urgensen, F. Meinert, M. J. Mark, H.-C. N¨agerl, and [16] S. Kivelson, W.-P. Su, J. R. Schrieffer, and A. J. Heeger, D.-S. L¨uhmann, Observation of Density-Induced Tunnel- Missing bond-charge repulsion in the extended Hubbard ing, Physical Review Letters 113, 193003 (2014). model: Effects in polyacetylene, Physical Review Letters [35] K. Biedro´n, M.L¸acki, and J. Zakrzewski, Extended Bose- 58, 1899 (1987). Hubbard model with dipolar and contact interactions, [17] J. Hirsch, Bond-charge repulsion and hole superconduc- Physical Review B 97, 245102 (2018). tivity, Physica C: Superconductivity and its Applications [36] D. Johnstone, N. Westerberg, C. W. Duncan, and 158, 326 (1989). P. Ohberg,¨ Staggered ground states in an optical lattice, [18] J. E. Hirsch and F. Marsiglio, Superconducting state Physical Review A 100, 043614 (2019). in an oxygen hole metal, Physical Review B 39, 11515 [37] J. Stasi´nska, R. W. Chhajlany, O. Dutta, and M. Lewen- (1989). stein, Effects of extended correlated hopping in a bose- [19] G. I. Japaridze and A. P. Kampf, Weak-coupling phase bose mixture, arXiv:1912.13359 (2019). diagram of the extended Hubbard model with correlated- [38] E. Lieb and D. Mattis, Ordering energy levels of inter- hopping interaction, Physical Review B 59, 12822 (1999). acting spin systems, Journal of Mathematical Physics 3 [20] J. Vidal and B. Dou¸cot, Strongly correlated hopping and , 749 (1962). many-body bound states, Physical Review B 65, 045102 [39] G.-S. Tian, Coexistence of the ferromagnetic and antifer- (2001). romagnetic long-range orders in the generalized antiferro- [21] A. Anfossi, P. Giorda, and A. Montorsi, Entanglement magnetic heisenberg model on a bipartite lattice, Journal in extended Hubbard models and quantum phase transi- of Physics A: Mathematical and General 27, 2305 (1994). tions, Physical Review B 75, 165106 (2007). [40] E. P. Raposo and M. D. Coutinho-Filho, Quantum crit- [22] L. Arrachea and A. A. Aligia, Exact Solution of a Hub- ical properties of ferrimagnetic hubbard chains, Physical bard Chain with Bond-Charge Interaction, Physical Re- Review Letters 78, 4853 (1997); Field theory of ferri- view Letters 73, 2240 (1994). magnetic hubbard chains, Physical Review B 59, 14384 [23] A. Schadschneider, Superconductivity in an exactly solv- (1999). able Hubbard model with bond-charge interaction, Phys- [41] F. C. Alcaraz and A. L. Malvezzi, Critical behaviour of ical Review B 51, 10386 (1995). mixed heisenberg chains, Journal of Physics A: Mathe- [24] F. Dolcini and A. Montorsi, Finite-temperature proper- matical and General 30, 767 (1997). ties of the Hubbard chain with bond-charge interaction, [42] B. Gu, G. Su, and S. Gao, Thermodynamics of Physical Review B 66, 075112 (2002). spin- 1 2 antiferromagnet-antiferromagnet-ferromagnet 12 and ferromagnet-ferromagnet-antiferromagnet trimer- [58] K. Sekiguchi and K. Hida, Partial Ferrimagnetism in S ized quantum Heisenberg chains, Physical Review B 73, = 1/2 Heisenberg Ladders with a Ferromagnetic Leg, an 134427 (2006); S.-S. Gong, S. Gao, and G. Su, Ther- Antiferromagnetic Leg, and Antiferromagnetic Rungs, J. modynamics of spin-1/2 tetrameric Heisenberg antiferro- Phys. Soc. Jpn. 86, 084706 (2017). magnetic chain, 80, 14413 (2009); S.-S. Gong, W. Li, [59] S. K. Pati, S. Ramasesha, and D. Sen, A density ma- Y. Zhao, and G. Su, Magnetism and thermodynamics trix renormalization group study of low-energy exci- of spin-(1/2, 1) decorated Heisenberg chain with spin-1 tations and low-temperature properties of alternating pendants, 81, 214431 (2010). spin systems, Journal of Physics: Condensed Matter [43] T. Sakai and S. Yamamoto, Critical behavior of 9, 8707 (1997); Low-lying excited states and low- anisotropic Heisenberg mixed-spin chains in a field, Phys- temperature properties of an alternating spin-1–spin-1/2 ical Review B 60, 4053 (1999). chain: A density-matrix renormalization-group study, [44] L. M. Ver´ıssimo, M. S. S. Pereira, J. Streˇcka, and M. L. Physical Review B 55, 8894 (1997). Lyra, Kosterlitz-Thouless and Gaussian criticalities in a [60] K. Maisinger, U. Schollw¨ock, S. Brehmer, H. J. Mikeska, mixed spin-(1/2, 5/2, 1/2) branched chain with exchange and S. Yamamoto, Thermodynamics of the (1, 1/2) fer- anisotropy, Physical Review B 99, 134408 (2019). rimagnet in finite magnetic fields, Physical Review B 58, [45] R. R. Montenegro-Filho, F. S. Matias, and M. D. R5908 (1998). Coutinho-Filho, Topology of many-body edge and ex- [61] A. S. F. Ten´orio, R. R. Montenegro-Filho, and M. D. tended quantum states in an open spin chain: 1/3 Coutinho-Filho, Quantum phase transitions in alter- plateau, Kosterlitz-Thouless transition, and Luttinger nating spin-1/2–spin-5/2 heisenberg chains, Journal of liquid, Physical Review B 102, 035137 (2020). Physics: Condensed Matter 23, 506003 (2011). [46] A. M. S. Macˆedo, M. C. dos Santos, M. D. Coutinho- [62] J. Streˇcka and T. Verkholyak, Magnetic Signatures of Filho, and C. A. Macˆedo, Magnetism and phase separa- Quantum Critical Points of the Ferrimagnetic Mixed tion in polymeric hubbard chains, Physical Review Let- Spin-(1/2, S) Heisenberg Chains at Finite Temperatures, ters 74, 1851 (1995). Journal of Low Temperature Physics 187, 712 (2017). [47] R. R. Montenegro-Filho and M. D. Coutinho-Filho, [63] J. Streˇcka, Breakdown of a Magnetization Plateau in Fer- Doped ab2 hubbard chain: Spiral, nagaoka and rimagnetic Mixed Spin-(1/2,S) Heisenberg Chains due to resonating-valence-bond states, phase separation, and a Quantum Phase Transition towards the Luttinger Spin luttinger-liquid behavior, Physical Review B 74, 125117 Liquid, Acta Physica Polonica A 131, 624 (2017). (2006). [64] W. M. da Silva and R. R. Montenegro-Filho, Magnetic- [48] G. Sierra, M. A. Mart´ın-Delgado, S. R. White, D. J. field–temperature phase diagram of alternating ferrimag- Scalapino, and J. Dukelsky, Diagonal ladders: A class netic chains: Spin-wave theory from a fully polarized vac- of models for strongly coupled electron systems, Phys. uum, Physical Review B 96, 214419 (2017). Rev. B 59, 7973 (1999). [65] M. Hagiwara, K. Minami, Y. Narumi, K. Tatani, and [49] R. R. Montenegro-Filho and M. D. Coutinho-Filho, Mag- K. Kindo, Magnetic Properties of a Quantum Ferrimag- netic and nonmagnetic phases in doped ab2 t-j hubbard net: NiCu(pba)( D 2 O ) 3 · 2 D 2 O, Journal of the chains, Phys. Rev. B 90, 115123 (2014). Physical Society of Japan 67, 2209 (1998). [50] K. Kobayashi, M. Okumura, S. Yamada, M. Machida, [66] M. Hagiwara, Y. Narumi, K. Minami, K. Tatani, and and H. Aoki, Superconductivity in repulsively interacting K. Kindo, Magnetization Process of the S = 1/2 and 1 fermions on a diamond chain: Flat-band-induced pairing, Ferrimagnetic Chain and Dimer, Journal of the Physical Phys. Rev. B 94, 214501 (2016). Society of Japan 68, 2214 (1999). [51] K. Hida, Magnetic Properties of the Spin-1/2 [67] O. Kahn, Y. Pei, M. Verdaguer, J. P. Renard, and Ferromagnetic-Ferromagnetic-Antiferromagnetic Trimer- J. Sletten, Magnetic ordering of manganese(II) cop- ized Heisenberg Chain, J. Phys. Soc. Jpn. 63, 2359 per(II) bimetallic chains; design of a molecular based (1994). ferromagnet, Journal of the American Chemical Society [52] K. Takano, K. Kubo, and H. Sakamoto, Ground states 110, 782 (1988). with cluster structures in a frustrated Heisenberg chain, [68] A. Gleizes and M. Verdaguer, Ordered magnetic bimetal- Journal of Physics: Condensed Matter 8, 6405 (1996). lic chains: a novel class of one-dimensional compounds, [53] R. R. Montenegro-Filho and M. D. Coutinho-Filho, Journal of the American Chemical Society 103, 7373 Frustration-induced quantum phase transitions in a (1981). quasi-one-dimensional ferrimagnet: Hard-core boson [69] M. Verdaguer, A. Gleizes, J. P. Renard, and map and the tonks-girardeau limit, Physical Review B J. Seiden, Susceptibility and magnetization of 78, 014418 (2008). CuMn(S2C2O2)2 7.5H2O, Physical Review B 29, [54] N. B. Ivanov, Spin models of quasi-1d quantum ferrimag- 5144 (1984). · nets with competing interactions, Condens. Matter Phys. [70] P. J. Van Koningsbruggen, O. Kahn, K. Nakatani, Y. Pei, 12, 435 (2009). J. P. Renard, M. Drillon, and P. Legoll, Magnetism of A- [55] T. Shimokawa and H. Nakano, Frustration-Induced Fer- copper(II) bimetallic chain compounds (A = iron, cobalt, rimagnetism in S = 1/2 Heisenberg Spin Chain, J. Phys. nickel): one- and three-dimensional behaviors, Inorganic Soc. Jpn. 80, 043703 (2011). Chemistry 29, 3325 (1990). [56] S. C. Furuya and T. Giamarchi, Spontaneously magne- [71] H. Yamaguchi, T. Okita, Y. Iwasaki, Y. Kono, N. Ue- tized Tomonaga-Luttinger liquid in frustrated quantum moto, Y. Hosokoshi, T. Kida, T. Kawakami, A. Mat- antiferromagnets, Physical Review B 89, 1 (2014). suo, and M. Hagiwara, Experimental realization of Lieb- [57] F. Amiri, G. Sun, H. J. Mikeska, and T. Vekua, Ground- Mattis plateau in a quantum spin chain, Scientific Re- state phases of a rung-alternated spin- 12 Heisenberg lad- ports 10, 9193 (2020); H. Yamaguchi, Y. Iwasaki, der, Physical Review B 92, 1 (2015). Y. Kono, T. Okita, A. Matsuo, M. Akaki, M. Hagiwara, 13 and Y. Hosokoshi, Low-energy magnetic excitations in [80] S. R. White, Density matrix formulation for quantum the mixed spin-(1/2, 5/2) chain, Physical Review B 102, renormalization groups, Physical Review Letters 69, 10.1103/PhysRevB.102.060408 (2020). 2863 (1992). [72] Y. Noriki and S. Yamamoto, Modified Spin-Wave The- [81] M. J. Rice and E. J. Mele, Elementary Excitations of a ory on Low-Dimensional Heisenberg Ferrimagnets: A Linearly Conjugated Diatomic Polymer, Physical Review New Robust Formulation, J. Phys. Soc. Jpn. 86, 034714 Letters 49, 1455 (1982). (2017). [82] N. R. Cooper, J. Dalibard, and I. B. Spielman, Topologi- [73] S. Brehmer, H.-J. Mikeska, and S. Yamamoto, Low- cal bands for ultracold atoms, Reviews of Modern Physics temperature properties of quantum antiferromagnetic 91, 015005 (2019). chains with alternating spins and, Journal of Physics: [83] Y.-T. Lin, D. M. Kennes, M. Pletyukhov, C. S. Weber, Condensed Matter 9, 3921 (1997). H. Schoeller, and V. Meden, Interacting Rice-Mele model: [74] S. Yamamoto and T. Fukui, Thermodynamic properties Bulk and boundaries, Physical Review B 102, 085122 of heisenberg ferrimagnetic spin chains: Ferromagnetic- (2020). antiferromagnetic crossover, Physical Review B 57, [84] M. Takahashi, Few-dimensional heisenberg ferromagnets 14008 (1998). at low temperature, Physical Review Letters 58, 168 [75] S. Yamamoto, S. Brehmer, and H.-J. Mikeska, Elemen- (1987). tary excitations of heisenberg ferrimagnetic spin chains, [85] B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, Physical Review B 57, 13610 (1998). J. Freire, S. Fuchs, L. Gamper, J. Gukelberger, E. Gull, [76] S. Yamamoto, T. Fukui, K. Maisinger, and S. Guertler, A. Hehn, R. Igarashi, S. V. Isakov, U. Schollw¨ock, Combination of ferromagnetic and D. Koop, P. N. Ma, P. Mates, H. Matsuo, O. Parcol- antiferromagnetic features in heisenberg ferrimagnets, let, G. Paw lowski, J. D. Picon, L. Pollet, E. Santos, Journal of Physics: Condensed Matter 10, 11033 (1998). V. W. Scarola, U. Schollw¨ock, C. Silva, B. Surer, S. Todo, [77] N. B. Ivanov, Magnon dispersions in quantum Heisen- S. Trebst, M. Troyer, M. L. Wall, P. Werner, and S. Wes- berg ferrimagnetic chains at zero temperature, Physical sel, The ALPS project release 2.0: open source software Review B 62, 3271 (2000). for strongly correlated systems, J. Stat. Mech.: Theory [78] S. Yamamoto, Bosonic representation of one-dimensional Exp. 2011 (05), P05001. heisenberg ferrimagnets, Physical Review B 69, 064426 [86] T. Sakai and S. Yamamoto, Coexistent quantum (2004). and classical aspects of magnetization plateaux in [79] S. R. White, Density-matrix algorithms for quantum alternating-spin chains, Journal of Physics Condensed renormalization groups, Physical Review B 48, 10345 Matter 12, 9787 (2000). (1993). > (1/2, 1) b a S t-approx. spin model (1/2, 2)