Outline
Comptes Rendus

Physics of integer-spin antiferromagnetic chains: Haldane gaps and edge states
Comptes Rendus. Chimie, Volume 22 (2019) no. 6-7, pp. 445-451.

Abstract

The antiferromagnetic Heisenberg spin chain with integer spin has a short-range magnetic order and an excitation energy gap above the ground state. This so-called Haldane gap is proportional to the exchange coupling constant J of the Heisenberg chain. In this study, we discuss recent results about the spin dependence of the Haldane gap and conjecture an analytical formula valid asymptotically for large spin values. We then discuss the robustness of the edge states of the spin-1 chain by studying a spin-1 ladder using the density matrix renormalization group (DMRG) algorithm. We show that the peculiar hidden topological order of the spin-1 chain disappears smoothly by increasing the ladder rung coupling without any intervening phase transition. This is evidence for the fragile character of the topological order of the spin-1 chain.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2019.05.005
Keywords: Magnetic properties, Statistical mechanics, Ab initio calculations

T. Jolicoeur 1; Olivier Golinelli 1

1 Institut de physique théorique, Université Paris-Saclay, CNRS, CEA, 91190 Gif-sur-Yvette, France
@article{CRCHIM_2019__22_6-7_445_0,
     author = {T. Jolicoeur and Olivier Golinelli},
     title = {Physics of integer-spin antiferromagnetic chains: {Haldane} gaps and edge states},
     journal = {Comptes Rendus. Chimie},
     pages = {445--451},
     publisher = {Elsevier},
     volume = {22},
     number = {6-7},
     year = {2019},
     doi = {10.1016/j.crci.2019.05.005},
     language = {en},
}
TY  - JOUR
AU  - T. Jolicoeur
AU  - Olivier Golinelli
TI  - Physics of integer-spin antiferromagnetic chains: Haldane gaps and edge states
JO  - Comptes Rendus. Chimie
PY  - 2019
SP  - 445
EP  - 451
VL  - 22
IS  - 6-7
PB  - Elsevier
DO  - 10.1016/j.crci.2019.05.005
LA  - en
ID  - CRCHIM_2019__22_6-7_445_0
ER  - 
%0 Journal Article
%A T. Jolicoeur
%A Olivier Golinelli
%T Physics of integer-spin antiferromagnetic chains: Haldane gaps and edge states
%J Comptes Rendus. Chimie
%D 2019
%P 445-451
%V 22
%N 6-7
%I Elsevier
%R 10.1016/j.crci.2019.05.005
%G en
%F CRCHIM_2019__22_6-7_445_0
T. Jolicoeur; Olivier Golinelli. Physics of integer-spin antiferromagnetic chains: Haldane gaps and edge states. Comptes Rendus. Chimie, Volume 22 (2019) no. 6-7, pp. 445-451. doi : 10.1016/j.crci.2019.05.005. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2019.05.005/

Version originale du texte intégral

1 The nonlinear σ model and the Haldane gap

It is well understood that the origin of exchange interactions between magnetic ions in an insulator is essentially of quantum mechanical nature. Indeed, it involves electron delocalization and the Pauli principle as needed for fermionic elementary constituents of matter. However, the description of magnetic properties of solids at the mesoscopic or macroscopic scale very often does not involve a fundamental way quantum mechanics. In case of ferromagnetic or antiferromagnetic chains, one can use coarse-grained magnetization vectors to describe ordered states that are perfectly classical quantities. Also excited states above ordered ground states like magnons admit a classical description, which has a wide range of applicability. This classical line of thought has been applied for a long time also in the realm of one-dimensional quantum spin systems with some success in the case of systems with rather large spin values. Therefore, it came as a surprise when Haldane [1–3] in the beginning of the 1980s made the bold claim that integer-spin antiferromagnetic chains have a gap to all excitations at odds with previous magnon-based reasonings. The picture that emerges is that integer-spin chains have only a short-range magnetic order even at zero temperature with spin–spin correlations decaying exponentially with a characteristic length, and there is a gap to all excitations. On the contrary, half-integer-spin chains have antiferromagnetic order with no characteristic scale albeit decaying as a power law and no gap to excited states. Interestingly, the first Haldane gap spin chain NENP [4] was synthesized roughly at the same time as Haldane's conjecture appeared. Since its synthesis, NENP has been an extremely good platform [5–7] to investigate Haldane gap physics. Quantitative comparisons between theory and experiments have been successful [8–10]. Detailed investigations have confirmed the existence of the Haldane gap close to 0.41J, where J is the nearest-neighbor antiferromagnetic exchange of the isotropic Heisenberg chain [11–13].

The Haldane conjecture was originally based on the derivation of an effective quantum field theory using an expansion of fluctuation about local antiferromagnetic order of a spin-S chain. This derivation has been the subject of much theoretical work and is now available in streamlined form in textbooks [14]. In this section, we propose a quick presentation of the main arguments and proceed to discuss the dependence of the Haldane gap on the spin magnitude S, as there are now several numerical studies [15] extending up to S=4.

The starting point is the imaginary time representation of the partition function in terms of spin coherent states:

Z=DΩˆiexp{iSiω[Ωˆi]0βS2iΩˆi·Ωˆ+1}(1)

The spin operators Six,y,z are written in terms of the spin coherent states Six,y,zSΩˆix,y,z, where the quantities Ωˆi are classical commuting vectors of unit norm and S is the spin magnitude. In fact, the Haldane mapping is an asymptotic expansion valid for large S and its quantitative use for S=1 is not guaranteed. Numerical direct studies down to S=1 have confirmed its applicability even in this case. The first term iω[Ωˆi] in the exponential is the Berry phase contribution, which is a nontrivial function of the imaginary time evolution of the coherent state. This term is in fact responsible for the difference between integer- and half-integer-spin chains. It is irrelevant to the bulk physics of the integer-spin chain case and it is a fortunate circumstance that provided this contribution can be dropped then the field theoretical analysis of the effective model we will derive is in fact well known. From the point of view of the effective long-distance low-energy field theory, the spin–1/2 spin chain is much difficult to understand even if this is the only case where exact results are known for the eigenstates by the Bethe solution. Because we are discussing only the Haldane gap, we will drop the Berry phase from now on (in the case of integer spin the Berry phase term is responsible for the appearance of edge spins S = 1/2, which is discussed in Section 2). We now proceed to construct an effective model valid for long distance and low energies. This requires the identification of the modes that are pertinent in this limit. We assume that these are antiferromagnetic as well as ferromagnetic fluctuations as observed in all approximation schemes used to study spin-chain physics. We write the spin coherent states as

Ωˆi=()inˆi1Li2/S2+1SLi(2)
wherenˆi is a unit vector and Li is a vector perpendicular to nˆi. This way of writing the spin coherent state leads naturally to a 1/S expansion. The nearest-neighbor Heisenberg coupling is expanded as follows:
Ωˆi·Ωˆi+1=1+12(nˆinˆi+1)2+12S2(Li+Li+1)2(3)

The Berry phase contribution in the path integral leads to a coupling between the L vector and the imaginary time derivative of the order parameter field nˆi:

DLiexp{0dτdq2π(nˆ×τnˆ)q·Lq0dτdq2πχq1Lq·Lq}(4)
where the susceptibility is given by χq1=4Jcos2(qa/2), with a being the distance between spins. Progress can be made if we take the long wavelength limit q0 in the susceptibility χ14Jχ0. Then, we can integrate the ferromagnetic fluctuation vector L and obtain a simple quadratic action:
Z=Dnˆexp{120dτdx[χ0(τnˆ)2+ρs(xnˆ)2]}(5)
where we have introduced a stiffness ρs=JS2a. Now Euclidean time and space are on an equal footing, provided we rescale by the velocity c=ρs/χ0. At zero temperature, β; this field theory is in fact the well-known nonlinear σ model (NLσM). It was introduced to describe critical properties of classical Heisenberg ferromagnets at nonzero temperatures. Here, we observe the consequence of the mapping of a 1+1 dimensional quantum system onto a classical two-dimensional system. In this new language, there is a very basic fact that there is no order at any nonzero temperature for classical Heisenberg spins. In physical terms, spinwaves disorder the system beyond a finite correlation length ξ and spin correlations decay exponentially with distance. We thus infer immediately that the space correlations of the nˆ vector field decay exponentially. In addition, as we have complete equivalence between space and (Euclidean) time, the time correlations also decay exponentially:
nˆ(0,x)nˆ(τ,x)exp((cτ)ξ1)(6)

It is a classic result in Euclidean quantum mechanics that the imaginary-time decay of correlations is governed by the energy gap of the system. Thus, we conclude from Eq. 6 that the spin chain has a gap Δ=cξ1. If we rewrite the effective action in terms of the conventional NLσM, we find the action as follows:

SNLσM=12gdτdx(μnˆ)2(7)
with the coupling constant g=c/ρs=2/S. Renormalization group calculations have shown that the correlation length of this model has a known dependence upon g:
ξ=C2πgexp(2π/g)(1+O(g))(8)
where C is a pure number. Because we know the spin dependence of the velocity c=2JSa, we deduce that the gap depends upon the spin value as
ΔS/J=C0S2exp(πS)(1+O(1/S))(9)
with a prefactorC0, which is a pure number independent of the spin. This formula should be understood as being asymptotic, that is, valid when S. The case S = 2 has been studied experimentally [16] and theoretically [17]. More recent theoretical estimates are available for S = 3 [18] and even for S = 4.10 cm [15].

The pure number in factor of the σ model formula for ξ is known and is C=e/8 for a specific choice of the renormalization procedure. The renormalization procedure corresponding to the quantum lattice model cannot be easily related to the standard ways of renormalizing the σ model. Because we are not yet able to make this correspondence precise, we choose to make an extrapolation from known values of the Haldane gap up to S=4 (see Table 1). The Shanks extrapolation for three numbers a,b,c is simply calculated by the ratio (acb2)/(a+c2b) and is given in the last column of Table 1. After some trial and error, we are led to propose that the limiting value is 2exp(2)=14.778 close to 14.525. We thus conclude this section by conjecturing the asymptotic formula for the Haldane gap

ΔS/J=2S2exp(2πS)(1+O(1/S))(10)

Table 1

Gaps of the spin-S Heisenberg chain using state-of-the-art numerical results [15].

Table 1
SGap/JNLσMGap/NLσMShanks
14.104800 × 10−14.321392 × 10−29.499
28.916000 × 10−27.469771 × 10−311.936
31.002000 × 10−27.262957 × 10−413.79619.787
47.990000 × 10−45.579748 × 10−514.32014.525

It remains to be confirmed whether methods like those of ref. [19] can relate the renormalization schemes and confirm this prefactor.

2 Edge states of a Haldane ladder

Another important feature of the spin-1 chain is the existence of a hidden magnetic order, which is long range. We have already commented on the fact that ordinary spin correlations decay exponentially with a finite correlation length. However, there is nevertheless a hidden order in the system. The nearest-neighbor Heisenberg S = 1 spin chain is not a solvable system, and hence one has to use numerical methods to study its quantitative properties. However, it is possible to perturb the Hamiltonian and obtain a new model [20] called AKLT whose ground state wavefunction can be found by elementary means, although it is still not solvable. The trick is to add a biquadratic spin coupling:

HAKLT=iSi·Si+1+13(Si·Si+1)2.(11)

The motivation is formal, because in the real world, a superexchange theory leads to negligible higher-spin couplings. However, with the very special choice of the 1/3 coefficient for the biquadratic term, the spin operator is now the projector onto the total spin S=2 of the pair i,i+1. This can be checked by elementary means (a bit tedious). Now we proceed to exhibit a wavefunction, which is an exact zero-energy eigenstate of Hamiltonian (Eq. 11). Each local spin-1 can be viewed as the triplet state of two spin-1/2 residing at the same site i. In the case of NENP, this corresponds to microscopic physics, because spin-1/2 elementary electrons in the nickel ions are coupled by Hund's coupling with a triplet state. Now, we make singlet bonds between neighboring sites (ions) and construct the wavefunction by capturing all spins in the following way:

|ΨAKLT={|ii+1|ii+1}{|i+1i+2|i+1i+2}(12)

If we consider two neighboring spins, then the total spin of such a pair can only be 0 or 1 but not 2, because the central spins 1/2 are already engaged in a singlet state. Thus, we see that all projectors onto the S = 2 state for any pair gives simply zero and hence, this wavefunction is an exact zero-energy eigenstate of the Hamiltonian (Eq. 11). In fact, this is the exact ground state of the AKLT Hamiltonian. Although this construction may seem artificial, we note that this Hamiltonian can be viewed as a perturbation of the usual Heisenberg Hamiltonian provided “1/3 is small” meaning that the biquadratic operator does not change the physics of the system. Indeed, we are lucky and it has been shown in detail that this is the case [21]. As a consequence, reasonings based on the model AKLT state are likely to be correct for the real world dominated by the Heisenberg exchange. A very simple inference from the AKLT description is that if we consider an open chain, then there are dangling spins 1/2 at the end. Of course real materials always involve open chains and we thus can predict that local probes should detect effective spin 1/2 degrees of freedom at the end of NENP chains. The experiment has been done [22] using NMR of a zinc impurity breaking the chain and the number of NMR lines is exactly in agreement with the prediction of the spin 1/2 instead of 1 as one may guess naively. For an open chain, we expect that the two end spins 1/2 appear as almost degenerate singlet and triplet states due to a coupling through the bulk of the chain that should go exponentially to zero when the chain length increases. So the spectrum of an open chain should display quasidegenerate S = 0 and S = 1 states and, above a (Haldane) gap, an S = 2 excitation will be obtained by creating a magnon-like state [23]. If we look at the magnetization profile of the Sz=1 member of the triplet, it should display nonzero magnetization only close to the edges of the chain. This is indeed what is observed numerically for the S = 1 chain using the density matrix renormalization group (DMRG) algorithm [13,24].

We now turn to the question of the robustness of the formation of these end spin states. They are due to the special ordering pattern described by the AKLT wavefunction [20,25], and it is still present for the realistic Heisenberg exchange Hamiltonian. It is known that several physical effects can destroy the Haldane gap, such as too strong anisotropies either local on-site or exchange, and also an applied magnetic field can close the Haldane gap. If we cross a phase transition by such a mechanism, we do not expect edge states to survive because the peculiar long-range correlations of the AKLT state are destroyed. It has been noted some time ago that in fact there is no need for a sharp phase transition to kill the hidden order and the edge states. This is very clear in the case of the magnetization process of the Heisenberg S = 1 chain. Indeed, under an applied field the Haldane gap goes to zero at some critical field Hc and the magnetization remains zero up to this value. When the critical value is exceeded, the magnetization starts to increase continuously up to full saturation and the edge states disappear right at Hc. If now we consider a realistic material like NENP, then there are small but nonzero anisotropies that break the spin rotation symmetry. As a consequence, the magnetization transition [26] is now rounded and is no longer sharp. Magnetization appears immediately for infinitesimal applied fields and edge states disappear when going to saturation without any phase transition. This phenomenon has been dubbed [27] “symmetry protected topological order”. It leaves open the possibility that the hidden order and accompanying edge states may disappear even without breaking any symmetry and without any phase transition. We now show that this is the case by studying a coupled S = 1 ladder system. We thus envision a plausible molecular magnet in which neighboring Haldane chains of spins S = 1 are coupled along the rungs of a ladder by some independent exchange J. A spin Hamiltonian can thus be written as follows:

Hladder=iSi,ASi+1,A+iSi,BSi+1,B+JiSi,ASi,B(13)
where the sites along the two chains are labeled by i and A,B label the two coupled chains. We have taken the exchange along the chain as the unit of energies so the only remaining parameter is the ratio of exchange interactions along and across the chains called J. We note that recent experiments have shown that the organic molecular magnet BIP-TENO is precisely a spin-1 ladder [28,29]. However, the exchange interactions in BIP-TENO are more complex than our simple Hamiltonian (Eq. 13).

If we consider a ladder with periodic boundary conditions, then it has been established that there is no phase transition [30,31] between the J=0 decoupled chain limit and the J= with extremely strong rungs. This strong rung limit is simple to analyze: the ground state is made of singlet states involving two spins that are related by the rung coupling and the total state is just the tensor product of such singlets. The first excited state is made by breaking a rung bond from singlet to triplet, and this triplet will have a dispersion along the chain whose magnitude is given by the (relatively) small coupling along the chains. The small J limit is two weakly coupled Haldane chains and because they are gapped, they will be resilient to any small local perturbation like a small rung coupling. Even if there is no transition between these two limits, it remains unclear what happens to the edge states of an open ladder.

We first enumerate the low-lying states of this spin system starting from our knowledge of the S = 1 chain. We know that each chain will have a singlet and a triplet of low-lying states separated by the Haldane gap from the higher excitations [23]. With two chains, this means a total of 16 states that can be classified as two S = 0 singlets, three S = 1 triplets, and one S = 2 quintuplet. When adding a very small J, these will no longer be degenerated and first-order perturbation theory will lead to an order O(J) splitting of the states. Their edge nature can be revealed by computing the magnetization profile Snz along one of the coupled chains. We have performed DMRG calculations of this system to obtain the ground state wavefunction. Convergence becomes problematic when J is very small J<0.001, but results are reliable for higher values of the coupling. We use up to 1500 states per block and chain length up to 200 spins.

If we consider the low-lying Sz=0 (where Sz is the total spin) state, we observe nothing remarkable. The magnetization is uniformly vanishing and this does not change from the decoupled limit J=0 till Jreaches a very large number as expected. We now compute the ground state wavefunction in each sector Sz=1,2,3; we expect that Sz=1,2 display only edge modes whereas Sz=3 should capture a bulk magnon state. We vary the rung coupling in the Sz=1 sector in Figs. 1–3. For small J, we observe very clear edge mode oscillations of the local magnetization and when the chain is long enough, the bulk magnetization is zero (green line in the figure), see Fig. 1. If we increase J, there is coexistence of edge oscillations and bulk magnetization (see Fig. 2). Finally, at large enough Jvalues, all edge phenomena are washed out and we observe a single bump of the magnetization. We interpret this state as a single magnon state of the infinite rung limit that propagates as a particle in a box as reported in a study by White and Huse [13]. The important information is the gradual disappearance of the edge modes without any level crossing.

Fig. 1

The magnetization profile Snz for site n along the first 50 spins of an S = 1 Heisenberg ladder of 100 sites long (for a total of 200 spins). The expectation value is computed in the lowest-lying Sz=1 state, which belongs to the manifold of edge states. The rung exchange coupling is J=0.003. The green line is zero magnetization. The bulk is not magnetized, proving the edge nature of the excitations.

Fig. 2

The magnetization profile Snz in the lowest-lying Sz=1 state for site n along the S = 1 Heisenberg ladder of 50 sites long (total 100). The rung exchange coupling is now J=0.5; we observe that the magnetization is now essentially in the bulk and there are still edge oscillations of the magnetization.

Fig. 3

The magnetization profile as in Fig. 2. The rung exchange coupling is J=4; the state can be interpreted as a single magnon state clearly defined from the infinite rung limit. There is no reminder of the edge modes.

We expect that the same scheme applies for the lowest-lying Sz=2 state. Its behavior is displayed in Figs. 4–7. We again observe clear edge modes exponentially localized at the end of the ladder with no bulk magnetization at small J. There is a smooth crossover to a regime that we interpret as a two-magnon state (see Fig. 7).

Fig. 4

The magnetization profile Snz for site n along the S = 1 Heisenberg ladder of 100 sites long. The expectation value is computed in the lowest-lying Sz=2 state, which belongs to the manifold of edge states. The rung exchange coupling is J=0.003; there is zero bulk magnetization and only edge modes.

Fig. 5

The magnetization profile Snz as in Fig. 4. The rung exchange coupling is J=0.015; there is now coexistence of a small bulk magnetization and edge modes.

Fig. 6

The magnetization profile Snz as in Fig. 5. The rung exchange coupling is J=0.5; the bulk excitations display a two-rounded-peak structure.

Fig. 7

The magnetization profile Snz as in Fig. 6. The rung exchange coupling is J=4; edge modes have disappeared completely and we observe two magnons with characteristics of a particle-in-a-box wavefunction.

The situation should be different for the Sz=3 sector because in the decoupled limit J=0, one is forced to excite at least one magnon in one of the chains. This magnon will be delocalized and its total magnetization should spread all over the chain. This is exactly what we observe in Figs. 8 and 9. The bulk is always magnetized even when J0. In addition to this magnon contribution, there are also edge oscillations. When increasing the coupling, the magnon contribution becomes clearer (see Fig. 9). Finally, for very large rung coupling, Fig. 10, we are left with only rung magnon excitations and no edge modes. This is a gradual phenomenon with no physical discontinuity.

Fig. 8

The magnetization profile Snz in the lowest-lying Sz=3 state. The rung exchange coupling is J=0.003; although there are edge oscillations, we note that the bulk magnetization is nonzero. This is not a pure edge state.

Fig. 9

The magnetization profile Snz in the lowest-lying Sz=3 state. The rung exchange coupling is J=0.5; edge oscillations are still present but the three-peak structure appears.

Fig. 10

The magnetization profile Snz in the lowest-lying Sz=3 state. The rung exchange coupling is J=4; we observe a pure three-magnon state.

In principle, NMR on well-chosen chain-breaking impurities should be able to measure these physical effects.

3 Conclusion

Antiferromagnetic spin chains are a fascinating playground for states of magnetic matter that do not fit into the usual symmetry-breaking views of magnetism. The crafting of molecular magnets has allowed chemistry and physics to join forces and investigate experimentally the Haldane conjecture that predicts a fundamental difference between integer- and half-integer-spin chains. In the case of integer-spin chains, the so-called Haldane weakens with the spin value, and we have proposed here a conjecture for its asymptotic behavior for large spin S. Another very special feature of the Haldane gap state is the existence of long-range correlations that lead to edge modes, which have been observed in experiments [22]. These correlations are fragile and can be destroyed even without crossing any phase transition. We have given an explicit mechanism for the destruction of the hidden order by showing that an S = 1 spin ladder interpolates between a haldane-style phase and a simple tensor product phase with no kind of long-range order. The crossover we describe should be accessible to NMR measurements, provided one finds the appropriate molecular magnet. In this respect, the recent studies of BIP-TENO [28,29] or other magnets [32] show that it is indeed feasible to study spin-1 ladders and observe the interesting crossover between the Haldane phase and the rung singlet phase.

Acknowledgments

We acknowledge stimulating exchange interactions with M. Verdaguer and J. P. Renard over the years. We thank Gregoire Misguish for help with DMRG ITensor library (ITensor Library, Version 2.1.1, http://itensor.org). We also thank TOPMAT workshop of PSI2 project funded by the IDEX Paris-Saclay ANR-11-IDEX-0003-02 where part of this work has been performed. We also thank Ph. Lecheminant for careful reading of our manuscript and K. Totsuka for useful discussion.


References

[1] F.D.M. Haldane Ground State Properties of Antiferromagnetic Chains with Unrestricted Spin: Integer Spin Chains as Realisations of the O(3) Non-linear Sigma Model (arXiv:1612.00076 [cond-mat]ArXiv: 1612.00076) | arXiv

[2] F.D.M. Haldane Phys. Rev. Lett., 50 (1983) no. 15, pp. 1153-1156 https://link.aps.org/doi/10.1103/PhysRevLett.50.1153 | DOI

[3] F.D.M. Haldane Phys. Lett. A, 93 (1983) no. 9, pp. 464-468 http://www.sciencedirect.com/science/article/pii/037596018390631X | DOI

[4] A. Meyer; A. Gleizes; J.J. Girerd; M. Verdaguer; O. Kahn Inorg. Chem., 21 (1982) no. 5, pp. 1729-1739 | DOI

[5] J.P. Renard; M. Verdaguer; L.P. Regnault; W.A.C. Erkelens; J. Rossat-Mignod; W.G. Stirling Europhys. Lett., 3 (1987) no. 8, p. 945 http://stacks.iop.org/0295-5075/3/i=8/a=013 | DOI

[6] L.P. Regnault; J. Rossat-Mignod; J.P. Renard J. Magn. Magn. Mater., 104–107 (1992), pp. 869-870 http://www.sciencedirect.com/science/article/pii/030488539290400I | DOI

[7] J.P. Renard; M. Verdaguer; L.P. Regnault; W.A.C. Erkelens; J. Rossat-Mignod; J. Ribas; W.G. Stirling; C. Vettier J. Appl. Phys., 63 (1998) no. 8, p. 3538 https://aip.scitation.org/doi/10.1063/1.340736 | DOI

[8] O. Golinelli; T. Jolicoeur; R. Lacaze Phys. Rev. B, 46 (1992) no. 17, pp. 10854-10857 https://link.aps.org/doi/10.1103/PhysRevB.46.10854 | DOI

[9] O. Golinelli; T. Jolicoeur; R. Lacaze Phys. Rev. B, 45 (1992) no. 17, pp. 9798-9805 https://link.aps.org/doi/10.1103/PhysRevB.45.9798 | DOI

[10] O. Golinelli; T. Jolicoeur; R. Lacaze J. Phys. Condens. Matter, 5 (1993) no. 9, p. 1399 http://stacks.iop.org/0953-8984/5/i=9/a=024 | DOI

[11] I. Affleck J. Phys. Condens. Matter, 1 (1989) no. 19, p. 3047 http://stacks.iop.org/0953-8984/1/i=19/a=001 | DOI

[12] O. Golinelli; T. Jolicoeur; R. Lacaze Phys. Rev. B, 50 (1994) no. 5, pp. 3037-3044 https://link.aps.org/doi/10.1103/PhysRevB.50.3037 | DOI

[13] S.R. White; D.A. Huse Phys. Rev. B, 48 (1993) no. 6, pp. 3844-3852 https://link.aps.org/doi/10.1103/PhysRevB.48.3844 | DOI

[14] A. Auerbach Interacting Electrons and Quantum Magnetism, Springer New York, New York, NY, 1994 | DOI

[15] S. Todo; H. Matsuo; H. Shitara Parallel Loop Cluster Quantum Monte Carlo Simulation of Quantum Magnets Based on Global Union-Find Graph Algorithm (arXiv:1810.07485 [cond-mat]ArXiv: 1810.07485) | arXiv

[16] G.E. Granroth; M.W. Meisel; M. Chaparala; T. Jolicoeur; B.H. Ward; D.R. Talham Phys. Rev. Lett., 77 (1996) no. 8, pp. 1616-1619 https://link.aps.org/doi/10.1103/PhysRevLett.77.1616 | DOI

[17] U. Schollwoeck; O. Golinelli; T. Jolicoeur Phys. Rev. B, 54 (1996) no. 6, pp. 4038-4051 https://link.aps.org/doi/10.1103/PhysRevB.54.4038 | DOI

[18] S. Todo; K. Kato Phys. Rev. Lett., 87 (2001) no. 4 https://link.aps.org/doi/10.1103/PhysRevLett.87.047203 | DOI

[19] T. Jolicoeur; J. Le Guillou Mod. Phys. Lett. B, 05 (1991) no. 08, pp. 593-596 https://www.worldscientific.com/doi/abs/10.1142/S021798499100071X | DOI

[20] I. Affleck; T. Kennedy; E.H. Lieb; H. Tasaki Phys. Rev. Lett., 59 (1987) no. 7, pp. 799-802 https://link.aps.org/doi/10.1103/PhysRevLett.59.799 | DOI

[21] U. Schollwoeck; T. Jolicoeur; O. Golinelli Phys. Rev. B, 53 (1996) no. 6, pp. 3304-3311 https://link.aps.org/doi/10.1103/PhysRevB.53.3304 | DOI

[22] M. Hagiwara; K. Katsumata; I. Affleck; B.I. Halperin; J.P. Renard Phys. Rev. Lett., 65 (1990) no. 25, pp. 3181-3184 https://link.aps.org/doi/10.1103/PhysRevLett.65.3181 | DOI

[23] T. Kennedy J. Phys. Condens. Matter, 2 (1990) no. 26, p. 5737 http://stacks.iop.org/0953-8984/2/i=26/a=010 | DOI

[24] U. Schollwoeck Ann. Phys., 326 (2011) no. 1, pp. 96-192 http://www.sciencedirect.com/science/article/pii/S0003491610001752 | DOI

[25] M. den Nijs; K. Rommelse Preroughening transitions in crystal surfaces and valence-bond phases in quantum spin chains, Phys. Rev. B, Volume 40 (1989) no. 7, pp. 4709-4734 https://link.aps.org/doi/10.1103/PhysRevB.40.4709 | DOI

[26] K. Katsumata; H. Hori; T. Takeuchi; M. Date; A. Yamagishi; J.P. Renard Phys. Rev. Lett., 63 (1989) no. 1, pp. 86-88 https://link.aps.org/doi/10.1103/PhysRevLett.63.86 | DOI

[27] Y. Fuji Phys. Rev. B, 93 (2016) no. 10, p. 104425 https://link.aps.org/doi/10.1103/PhysRevB.93.104425 | DOI

[28] K. Katoh; Y. Hosokoshi; K. Inoue; T. Goto J. Phys. Soc. Jpn., 69 (2000) no. 4, pp. 1008-1011 https://journals.jps.jp/doi/10.1143/JPSJ.69.1008 | DOI

[29] T. Sakai; N. Okazaki; K. Okamoto; K. Kindo; Y. Narumi; Y. Hosokoshi; K. Kato; K. Inoue; T. Goto Phys. B Condens. Matter, 329–333 (2003), pp. 1203-1204 http://www.sciencedirect.com/science/article/pii/S0921452602021233 | DOI

[30] D. Allen; D. Senechal Phys. Rev. B, 61 (2000) no. 18, pp. 12134-12142 https://link.aps.org/doi/10.1103/PhysRevB.61.12134 | DOI

[31] S. Todo; M. Matsumoto; C. Yasuda; H. Takayama Phys. Rev. B, 64 (2001) no. 22, p. 224412 https://link.aps.org/doi/10.1103/PhysRevB.64.224412 | DOI

[32] C. Mennerich; H.-H. Klauss; M. Broekelmann; F.J. Litterst; C. Golze; R. Klingeler; V. Kataev; B. Bchner; S.-N. Grossjohann; W. Brenig; M. Goiran; H. Rakoto; J.-M. Broto; O. Kataeva; D.J. Price Phys. Rev. B, 73 (2006) no. 17, p. 174415 https://link.aps.org/doi/10.1103/PhysRevB.73.174415 | DOI


Comments - Policy