Spin supersolid in an anisotropic spin-one Heisenberg chain

We consider an S=1 Heisenberg chain with strong exchange (Delta=J(z)/J(perpendicular)) and single-ion uniaxial anisotropy (D) in a magnetic field (B) along the symmetry axis. The low-energy spectrum is described by an effective S=1/2 XXZ model that acts on two different low-energy sectors for a finite range of fields. The vacuum of each sector exhibits Ising-like antiferromagnetic ordering coexisting with the finite spin stiffness obtained from the exact solution of the XXZ model. In this way, we demonstrate the existence of a spin supersolid phase. We also compute the full Delta-B quantum phase diagram using a quantum Monte Carlo method.     

Field-induced supersolid phase in spin-one Heisenberg models

We use numerical methods to demonstrate that the phase diagram of S=1 Heisenberg models with uniaxial anisotropy contains an extended supersolid phase. We show that this Hamiltonian is a particular case of a more general and ubiquitous model that describes the low-energy spectrum of some isotropic and frustrated spin-dimer systems. This result is crucial for finding a spin supersolid state in real magnets.     

Diverse solid and supersolid phases of bosons in a triangular lattice

We investigate the ground state of bosons with long-range interactions in the large U limit on a triangular lattice. By mapping this system to the spin-1/2 XXZ model in a magnetic field, we can apply the spin wave theory to this study. We demonstrate how to construct the phase diagrams within the spin wave theory. The phase diagrams are given in an extensive parameter region, where, besides the superfluid phase, diverse solid and supersolid phases are shown to exist in this model. Especially, we find that the phase diagram obtained in this method is consistent with the one obtained previously using numerical techniques in the Ising limit. This confirms the effectiveness of our method. We analyze the stability of all the obtained supersolids and show that they will not be ruined by the quantum fluctuations. We observe that the quantum fluctuations in the stripe supersolid phase could be enhanced by the external field. We also discuss the relevance of our result with the experiment that may be realized with ultracold bosonic polar molecules in a triangular optical lattice.     

Persistent supersolid phase of hard-core bosons on the triangular lattice

We study hard-core bosons with unfrustrated hopping (t) and nearest neighbor repulsion (U) (spin S=1/2 XXZ model) on the triangular lattice. At half filling, the system undergoes a zero temperature (T) quantum phase transition from a superfluid phase at small U to a supersolid at Uc approximately 4.45 in units of 2t. This supersolid phase breaks the lattice translation symmetry in a characteristic sqrt[3] x square root of 3 pattern, and is remarkably stable--indeed, a smooth extrapolation of our results indicates that the supersolid phase persists for arbitrarily large U/t.     

NMR evidence for the persistence of a spin superlattice beyond the 1/8 magnetization plateau in SrCu2(BO3)_{2}

We present 11B NMR studies of the 2D frustrated dimer spin system SrCu2(BO3)_{2} in the field range 27-31 T covering the upper phase boundary of the 1/8 magnetization plateau, identified at 28.4 T. Our data provide a clear evidence that above 28.4 T the spin superlattice of the 1/8 plateau is modified but does not melt even though the magnetization increases. Although this is precisely what is expected for a supersolid phase, the microscopic nature of this new phase is much more complex. We discuss the field-temperature phase diagram on the basis of our NMR data.     

Quantum and thermal transitions out of the supersolid phase of a 2D quantum antiferromagnet

We investigate the thermodynamic properties of a field-induced supersolid phase in a 2D quantum antiferromagnet model. Using quantum Monte Carlo simulations, a very rich phase diagram is mapped out in the temperature-magnetic-field plane, with an extended supersolid region where a diagonal (solid) order coexists with a finite XY spin stiffness (superfluid). The various quantum and thermal transitions out of the supersolid state are characterized. Experimental consequences in the context of field-induced magnetization plateau materials are briefly discussed.     

Supersolid hard-core bosons on the triangular lattice

We determine the phase diagram of hard-core bosons on a triangular lattice with nearest-neighbor repulsion, paying special attention to the stability of the supersolid phase. Similar to the same model on a square lattice we find that for densities rho<1/3 or rho>2/3 a supersolid phase is unstable and the transition between a commensurate solid and the superfluid is of first order. At intermediate fillings 1/3相似文献   

Density-matrix renormalization group study of the spin-Peierls instability in the antiferromagnetic Heisenberg ladder

The columnar dimerized antiferromagnetic S?=?1/2 spin ladder is numerically studied by the density-matrix renormalization-group (DMRG) method. The elastic lattice with spin-phonon coupling ?? and lattice elastic force k is introduced into the system. Thus the S?=?1?/?2 Heisenberg spin chain is unstable towards dimerization (the spin-Peierls transition). However, the dimerization should be suppressed if the rung coupling J _?? is sufficiently large, and a Columnar dimer to Rung singlet phase transition takes place. After a self-consistent calculation of the dimerization, we determine the quantum phase diagram by numerically computing the singlet-triplet gap, the dimerization amplitude, the order parameters, the rung spin correlation and quantum entropies. Our results show that the phase boundary between the Columnar dimer phase and Rung singlet phase is approximately of the form J _?? ~ \hbox{$(\frac{k}{\alpha^{2}})^{-\frac{5}{4}}$} ( k ?? 2 ) ? 5 4.     

Commensurate supersolid of three-dimensional lattice bosons

Using an unbiased quantum Monte?Carlo method, we obtain convincing evidence of the existence of a checkerboard supersolid at a commensurate filling factor 1/2 (a commensurate supersolid) in the soft-core Bose-Hubbard model with nearest-neighbor repulsions on a cubic lattice. In conventional cases, supersolids are realized at incommensurate filling factors by a doped-defect-condensation mechanism, where particles (holes) doped into a perfect crystal act as interstitials (vacancies) and delocalize in the crystal order. However, in the model, a supersolid state is stabilized even at the commensurate filling factor 1/2 without doping. By performing grand canonical simulations, we obtain a ground-state phase diagram that suggests the existence of a supersolid at a commensurate filling. To obtain direct evidence of the commensurate supersolid, we next perform simulations in canonical ensembles at a particle density ρ=1/2 and exclude the possibility of phase separation. From the obtained snapshots, we discuss its microscopic structure and observe that interstitial-vacancy pairs are unbound in the crystal order.     

Hard-Core Bosons on a Two-Dimensional Square Optical Superlattice

In this work,we theoretically study hard-core bosons on a two-dimensional square optical superlattice at T = 0.First of all,we present the mean field phase diagram of this model in terms of the chemical potential μ and the alternating potential strength △.Besides a superfluid(SF) phase at △ = 0 and a charge density wave(CDW)phase in the large △ at half filling,we demonstrate that a supersolid(SS) phase emerges in the moderate △.Then,we focus on the μ = 0,e.g.,half filling case,using large-S semi-classical spin-wave approximation to study the SS to CDW quantum phase transition.In particular,we calculate the ground-state energy and the superfluid density at the level of1/S correction.We then compare the spin-wave results with the large scale quantum Monte Carlo(QMC) simulations using the cluster stochastic series expansion(CSSE) algorithm,and find that while the spin wave method is intuitive with clear physical pictures,the quantum critical point is quite different from that of numerical results which is believed to be accurate.We suggest that as simple as it is,this model still exhibits strong quantum fluctuations near the quantum critical point beyond the power of semiclassical spin-wave approach.     

Quantum versus geometric disorder in a two-dimensional Heisenberg antiferromagnet

We present a numerical study of the spin-1/2 bilayer Heisenberg antiferromagnet with random interlayer dimer dilution. From the temperature dependence of the uniform susceptibility and a scaling analysis of the spin correlation length we deduce the ground state phase diagram as a function of nonmagnetic impurity concentration p and bilayer coupling g. At the site percolation threshold, there exists a multicritical point at small but nonzero bilayer coupling g(m)=0.15(3). The magnetic properties of the single-layer material La(2)Cu(1-p)(Zn,Mg)(p)O4 near the percolation threshold appear to be controlled by the proximity to this new quantum critical point.     

Interacting classical dimers on the square lattice

We study a model of close-packed dimers on the square lattice with a nearest neighbor interaction between parallel dimers. This model corresponds to the classical limit of quantum dimer models [D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 2376 (1988)]. By means of Monte Carlo and transfer matrix calculations, we show that this system undergoes a Kosterlitz-Thouless transition separating a low temperature ordered phase where dimers are aligned in columns from a high temperature critical phase with continuously varying exponents. This is understood by constructing the corresponding Coulomb gas, whose coupling constant is computed numerically. We also discuss doped models and implications on the finite-temperature phase diagram of quantum dimer models.     

Quantum phase diagram of the t-J(z) chain model

We present the quantum phase diagram of the one-dimensional t- J(z) model for arbitrary spin (integer or half-integer) and sign of the spin-spin interaction J(z), using an exact mapping to a spinless fermion model that can be solved exactly using the Bethe ansatz. We discuss its superconducting phase as a function of hole doping nu. Motivated by the new paradigm of high temperature superconductivity, the stripe phase, we also consider the effect the antiferromagnetic background has on the t- J(z) chain intended to mimic the stripe segments.     

Role of the spin anisotropy of the interchain interaction in weakly coupled antiferromagnetic Heisenberg chains

In quasi-one-dimensional(q1D) quantum antiferromagnets, the complicated interplay of intrachain and interchain exchange couplings may give rise to rich phenomena. Motivated by recent progress on field-induced phase transitions in the q1D antiferromagnetic(AFM) compound YbAlO_3, we study the phase diagram of spin-1/2 Heisenberg chains with Ising anisotropic interchain couplings under a longitudinal magnetic field via large-scale quantum Monte Carlo simulations,and investigate the role of the spin anisotropy of the interchain coupling on the ground state of the system. We find that the Ising anisotropy of the interchain coupling can significantly enhance the longitudinal spin correlations and drive the system to an incommensurate AFM phase at intermediate magnetic fields, which is understood as a longitudinal spin density wave(LSDW). With increasing field, the ground state changes to a canted AFM order with transverse spin correlations. We further provide a global phase diagram showing how the competition between the LSDW and the canted AFM states is tuned by the Ising anisotropy of the interchain coupling.     

Phase separation in supersolids

We study quantum phase transitions in the ground state of the two dimensional hard-core boson Hubbard Hamiltonian. Recent work on this and related models has suggested "supersolid" phases with simultaneous diagonal and off-diagonal long range order. We show numerically that, contrary to the generally held belief, the most commonly discussed "checkerboard" supersolid is thermodynamically unstable. Furthermore, this supersolid cannot be stabilized by next-near-neighbor interaction. We obtain the correct phase diagram using the Maxwell construction. We demonstrate that the "striped" supersolid is thermodynamically stable and is separated from the superfluid phase by a continuous phase transition.     

面心立方格子上自旋S=1的Ising反铁磁体的基态研究

本文发展了文献[1]中提出的利用费密变换求Ising模型基态的方法,根据由费密变换后,激发态能量大于零求得相边界,提出了直接比较各真空态能量来求解各基态的相边界。采用此方法,研究了具有近邻对相互作用的最一般的自旋S=1的Ising反铁磁体模型,严格地求得它在面心立方格子上的基态能量和零温相图。同时,还利用多相点方程,将参数k,l平面分成60个小区域,它们分别对应于各个拓扑不等价的零温相图。 关键词：     

Extended Bose-Hubbard model with pair hopping on triangular lattice

We study systematically an extended Bose-Hubbard model on the triangular lattice by means of a meanfield method based on the Gutzwiller ansatz. Pair hopping terms are explicitly included and a three-body constraint is applied. The zero-temperature phase diagram and a variety of quantum phase transitions are investigated in great detail. In particular, we show the existence and the stability of the pair supersolid phase.     

Finite-temperature phase diagram of hard-core bosons in two dimensions

We determine the finite-temperature phase diagram of the square lattice hard-core boson Hubbard model with nearest neighbor repulsion using quantum Monte Carlo simulations. This model is equivalent to an anisotropic spin-1/2 XXZ model in a magnetic field. We present the rich phase diagram with a first order transition between a solid and superfluid phase, instead of a previously conjectured supersolid and a tricritical end point to phase separation. Unusual reentrant behavior with ordering upon increasing the temperature is found, similar to the Pomeranchuk effect in 3He.     

A renormalization-group study of an Ising spin-glass with annealed vacancies

Systems driven and characterized by fluctuations in density and magnetization can be realistically modeled using the Blume–Emery–Griffiths model; a spin-1 Ising model with bilinear, biquadratic, and crystal-field interactions. In this study, renormalization-group techniques are used on an exactly solvable system in which frustration is present due to competing ferromagnetic and antiferromagnetic interactions. Thus, this calculation models a spin-glass system with annealed vacancies. To determine the effects of these competing bilinear interactions, an exactly solvable frustrated hierarchical model has been constructed, similar to those introduced to study spin glasses [S.R. McKay, A.N. Berker, S. Kirkpatrick, Phys. Rev. Lett. 48 (1982) 767]. Phase diagrams have been calculated for a series of planes of constant biquadratic coupling while varying the temperature and concentration of annealed vacancies in the system. In addition, a phase diagram was produced for constant concentration of annealed vacancies as the biquadratic coupling (i.e. clustering bias) was varied. Each phase diagram reveals three qualitatively unique basins of attraction, each corresponding to a phase distinguished by a unique renormalization-group trajectory. The sink of each trajectory is interpreted to determine the nature of each phase: dense paramagnetic, dilute paramagnetic and spin-glass.     

Spin nematics and magnetization plateau transition in anisotropic kagome magnets

We study S=1 kagome antiferromagnets with an isotropic Heisenberg exchange J and strong easy-axis single-ion anisotropy D. For D>J, the low-energy physics can be described by an effective S=1/2 XXZ model with antiferromagnetic Jz approximately J and ferromagnetic J perpendicular approximately J2/D. Exploiting this connection, we argue that nontrivial ordering into a "spin-nematic" occurs whenever D dominates over J, and discuss its experimental signatures. We also study a magnetic field induced transition to a magnetization plateau state at magnetization 1/3 which breaks lattice translation symmetry due to ordering of the Sz and occupies a lobe in the B/Jz-Jz/J perpendicular phase diagram.