Abstract
We study the ground state phase diagram of a two dimensional mixed-spin system of coupled alternating spin-1 and 1/2 chains with a stripe supersolid phase. Utilizing different analytical and numerical approaches such as mean field approximation, cluster mean field theory and linear spin wave theory, we demonstrate that our system displays a rich ground state phase diagram including novel stripe supersolid, solids with different fillings and super-counterfluid phases, in addition to a stripe solid with half filling, superfluid and Mott insulating phases. In order to find a minimal mixed-spin model for stripe supersolidity, in the second part of the paper we consider two kinds of mixed-spin system of coupled alternating spin-1 and 1/2 chains with (i) anisotropic nearest neighbor interactions, (ii) anisotropic hoppings and study their ground state phase diagrams. We demonstrate that, for the systems with uniform hoppings, the repulsive intra-chains interactions are necessary for stripe supersolidity. In this case the minimal two dimensional mixed-spin model is a system of spin-1 and spin-1/2 XXZ chains, interacting via Ising Hamiltonian. In the case of anisotropic hoppings, a system of coupled Ising chains is the minimal model.
Similar content being viewed by others
Introduction
Supersolids are characterized by the coexistence of diagonal solid and off-diagonal superfluid long-range orders1,2,3,4. Combination of these two apparently antithetical properties has attracted the attentions of both experimentalists and theorists, and searching for this exotic phenomenon has become one of the main subjects of condensed matter and cold atoms physics5,6,7,8. Since the pioneering work of Jaksch et al., in describing the dynamics of an ultracold dilute gas of bosonic atoms in optical lattices with a Bose-Hubbard model9, lots of efforts have been devoted to search for supersolid phases in Bose-Hubbard models, experimentally and theoretically on one dimensional (1D) chains10,11,12,13, two dimensional (2D)13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36 lattice structures, bilayer systems of dipolar lattice bosons37,38 and three dimensional (3D) cubic lattices17,39,40,41.
Possible representations of bosonic operators with spin operators42, and bosonic statistics of magnetic excitations43,44 make spin systems another appropriate ground for searching various supersolid phases. It has been shown that 1D spin-1 chains45,46,47, 2D frustrated spin-1/248,49,50,51,52,53,54,55,56,57,58,59 and spin-160,61,62,63 models in an external magnetic field and 3D spin models64,65 possess different kinds of stripe, checkerboard and star supersolid phases. In these systems finite transverse magnetization implies the off-diagonal long-range superfluid order, the diagonal long-range solid orders are specified by longitudinal staggered magnetization, and the states with magnetization plateaux indicate the Mott insulating phases43.
In spite of many studies on uniform spin systems, the supersolidity of mixed-spin systems has not been addressed so far. Mixed-spin systems, or quantum ferrimagnets which are composed of different spins, mostly of two kinds, are a special class of spin models where their universality class is completely different from uniform spin models66,67,68,69. Ferrimagnets, which occur rather frequently in nature, are somehow between the antiferromagnets and the ferromagnets. Their lowest energy band is gapless which shows a ferromagnetic behavior while there is a finite gap to the next band above it which has the antiferromagnetic properties. It is the acoustical and optical nature of excitations which is the result of two different types of spin in each unit cell.
Recently, we have studied a 2D frustrated ferrimagnetic spin model, originating from an inhomogeneous 2D bosonic system, composed of two kinds of hard-core (a) and semi-hard-core (b) bosons, respectively with the nilpotency conditions: \({({a}_{i}^{\dagger })}^{2}=0\) and \({({b}_{i}^{\dagger })}^{3}=0\), and shown that the model on a square lattice with nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions displays the checkerboard supersolid phase70 which is not observed in the 2D uniform spin-1/2 system on square lattices with short-range interactions22,71,72. Actually the interactions between spins with different sizes decrease the quantum fluctuations and cause the stabilization of the checkerboard supersolid order. Aside from the mentioned theoretical study, the spin supersolid phase has also been recently realized in the MnCr2S4 mixed-spin compound with spinel structure73. In this paper, we introduce a different system of coupled alternating spin τ = 1 and \(\sigma =\frac{1}{2}\) chains (CAS) (see Fig. 1 and Eq. (1)) and show that our CAS system possesses a stripe supersolid (STS) phase, characterizing by the coexistence of stripe solid (ST) and superfluid (SF) orders. We investigate the ground state phase diagram of the CAS model using different analytical and numerical approaches such as mean field (MF) approximation, cluster mean field theory (CMFT) and linear spin wave theory (LSWT). Competition between NN and NNN interactions causes the system to undergo various first- and second-order phase transitions, and different solids, Mott insulators (MI), SF and super-counterfluid (SCF) to appear in the ground state phase diagram of the model, in addition to the STS. By studying the behavior of spin wave excitations, we investigate the stability of MF orders and demonstrate that, except at the superfluid-supersolid transition lines, overall quantum fluctuations are small in our CAS system and the MF predictions concerning the stability of phases are reliable.
In the second part of this paper, we look for a minimal mixed-spin CAS system, possessing an stable supersolid phase. In this respect, we consider two kinds of anisotropic CAS model: (i) a CAS system with anisotropic NN interactions where the intra-chains and inter-chains NN interactions are not the same, and (ii) a CAS system with anisotropic hopping energies in which the intra-chains and inter-chains hoppings are different. By obtaining the CMFT ground state phase diagrams of these systems, we demonstrate that the appearance of the STS order strongly depends on the amounts of intra-chains NN interaction. By studying the behavior of spin wave excitations, and also the behavior of diagonal and off-diagonal order parameters by CMFT with larger cluster sizes, we investigate amount of quantum fluctuations and consequently the stability of the STS phase in all anisotropic CAS systems. Our results indicate that the repulsive intra-chains NN interactions are necessary for the emergence of the STS phase. Finally, based on our achievements, we present a minimal mixed-spin CAS model with stable supersolid phase in the ground state phase diagram.
Our CAS model could be related to the mixed-valance iron polymer74 in the spin system and ladder-like optical lattices75 in the bosonic systems. Also our model could be realized in coupled one dimensional optical lattices76 by alternatively changing the optical depth.
This paper is organized as follows. In Sec. 2 we introduce our CAS system and present the bosonic counterpart of this model. In Sec. 3 we obtain the ground state phase diagram of the CAS model using MF approximation, CMFT and LSWT. The linear spin wave dispersions and number of excitation modes, are also presented in this section. The anisotropic CAS models with anisotropic NN interactions, and anisotropic hoppings are investigated in Sec. 4. In Sec. 5 we bring our methods and explain MF, CMFT and LSWT in detail. Finally, we will summarize our results and give the concluding remarks in Sec. 6.
Our Model
Coupled alternating spin 1 and \(\frac{{\bf{1}}}{{\bf{2}}}\) chains
Let us consider a 2D system of coupled alternating spin τ = 1 and \(\sigma =\frac{1}{2}\) chains (CAS), describing by the following Hamiltonian:
with
where α = x, y and z. The Hamiltonians H σ , H τ , and H στ include the intra-chains and inter-chains interactions, respectively. The parameters \({{\mathscr{J}}}^{x,y}\)(=−2J) and \({{\mathscr{J}}}^{z}\)(=V1) are the NN interactions, V2 denotes the inter-chains NNN interactions, and h is a magnetic field along z direction. The magnetic field h is proportional to applied magnetic fields as h = g σ μBB σ = g τ μBB τ , where μB is the Bohr magneton, g σ and g τ are the g-factors, respectively, for spins-1/2 and spins-1, and B σ and B τ are the external magnetic fields applied to the subsystems with spins 1/2 and 1, respectively. Throughout this paper we consider g σ B σ = g τ B τ , and study the effects of a uniform h on the ground state phase diagram of the system. Our CAS system is schematically shown in Fig. 1.
The Hamiltonian in Eq. (1) possesses the translational symmetry of the 2D lattice with the translational vector \(a\hat{x}+2a\hat{y}\) as well as the rotational U(1) symmetry. Spontaneously breaking of these symmetries by varying the model parameters, causes the system to experience various first- and second-order phase transitions. Different diagonal and off-diagonal long-range orders appear in the ground state phase diagram of the above model which will be discussed in next sections.
Bosonic counterpart of the CAS chains
In bosonic language, the systems of coupled uniform spin-1/2 chains are equivalent to the systems of coupled 1D optical lattices, containing hard-core bosons. These systems possess different superfluid and solid phases76, in the presence of intra-chains hoppings, repulsive intra-chains and attractive inter-chains interactions, but no supersolid phase is formed due to the hard-core nature of the bosons.
The three-body-constrained bosons (or semi-hard-core bosons) with the nilpotency condition \({({b}_{i}^{\dagger })}^{3}=0\), may remove this problem. This nilpotency condition signifies that one can put up two b particles on each lattice site. Recent studies on the coupled 1D optical lattices containing two kinds of boson, a hard-core and a semi-hard-core boson, show that the system displays different MIs and SF orders77, in the presence of intra-chains hoppings and inter-chains interactions, but still solid and supersolid phases are absent in the ground state phase diagram. Indeed in the absence of intra-chains repulsive interactions the translational symmetry of the system preserves and consequently no solidity occurs in the system. So, considering repulsive intra-chains interactions on the semi-hard-core bosons can result in a supersolid order in the system of coupled 1D optical lattices.
Here, by using the relations between semi-hard-core bosons and spin-1 operators70, we map our CAS model in Eq. (1) to a bosonic system of coupled alternating hard-core and semi-hard-core bosonic lattices. We demonstrate that the presence of repulsive intra-chains interactions is sufficient for the appearance of solid and supersolid phases.
Let us consider a system of coupled alternating 1D optical lattices with hard-core and semi-hard-core bosons, a and b, which interact via the following Hamiltonian:
where H a (H b ) contains interactions between a (b) bosons and ta (tb), Va (Vb) and μa (μb) are the hopping energy, the interaction and the chemical potential in a (b) sublattice. H ab gives the interactions between a and b bosons with tab the hopping energy and Vab, and V2 the interaction energies. \({a}_{i}^{\dagger }\) (a i ) and \({b}_{j}^{\dagger }\) (b j ) are respectively the creation (annihilation) operators of a and b particles at sites i and j, on a 2D square lattice. The particles a are canonical hard-core bosons and satisfy the canonical commutation relations. The number of a bosons at site i is \({n}_{i}^{a}={a}_{i}^{\dagger }{a}_{i}\), and the nilpotency condition for these bosons is \({({a}_{i}^{\dagger })}^{2}=0\). The b particles are semi-hard-core bosons and satisfy the following statistics algebra42,70:
where \({n}_{i}^{b}\) is the number of b bosons which possesses the relation \({({n}_{i}^{b})}^{\dagger }={n}_{i}^{b}\).
In order to obtain the bosonic Hamiltonian in Eq. (3) we have used the following linear spin-boson transformations between a bosons and spin-\(\frac{1}{2}\)78, and between b bosons and the spin-1 operators42:
Since these spin-boson transformations are isomorphic, all symmetries and physical properties of the CAS system (1) and the bosonic system (3) are identically the same. The bosonic Hamiltonian in Eq. (3) is transformed to the mixed-spin Hamiltonian in Eq. (1) by defining the following relations:
Ground State Phase Diagram
In order to obtain the ground state phase diagram of the CAS model, first we use a MF approximation to investigate the system classically and then using LSWT we take into account the effects of quantum fluctuations around MF ground state. Moreover, utilizing a generalized CMFT we obtain the modified ground state phase diagram of the CAS model.
In MF approximation, we first divide the system into four sublattices, the subsystem with spins σ into A and B, and the subsystem with spins τ into C and D, and then approximate the local spins’ averages with MF order parameters. The four-sublattice structure is expected to be emerged due to the NN and NNN interactions. By defining the MF order parameters: \(\langle {\sigma }_{{i}_{A}}^{\alpha }\rangle ={m}_{A}^{\alpha }\), \(\langle {\sigma }_{{i}_{B}}^{\alpha }\rangle ={m}_{B}^{\alpha }\), \(\langle {\tau }_{{i}_{C}}^{\alpha }\rangle ={M}_{C}^{\alpha }\), and \(\langle {\tau }_{{i}_{D}}^{\alpha }\rangle ={M}_{D}^{\alpha }\), where 〈…〉 denotes the expectation value on the MF ground state, the Hamiltonian in Eq. (1) is readily simplified to a single particle MF Hamiltonian (see Eq. 11). The expectation values of spin operators on the ground state of the MF Hamiltonian (sublattices’ magnetizations) are given in terms of other MF order parameters (see Sec. 5.1 for detail). Various kinds of long-range diagonal and off-diagonal orders are defined by these magnetizations. In the Table 1, we have defined all possible phases, appearing in the MF ground state phase diagrams.
We have plotted in Fig. 2-top the MF ground state phase diagrams of the CAS system, for the two different strengths of frustration: V2/V1 = 0.4 and 0.6. We have also illustrated the schematic pictures of various orders, at the bottom of Fig. 2. For small values of J/V1, independent of the strengths of frustrations, the ground state phase diagram is symmetric with respect to the J = 0 line. Far from J = 0 line, the system, however, behaves differently for \(J > 0\) and J < 0 regions. For large values of \(J > 0\) the U(1) symmetry of the system breaks spontaneously, and the SF long range order emerges in the system where each boson is spread out over the entire lattice, with long-range phase coherence. For J < 0, instead of SF phase, the SCF phase appears in the phase diagram, where the transverse components of the spins lie in opposite directions (see Fig. 2, the schematic picture of SCF). This phase is characterized by a transverse staggered magnetization and a longitudinal magnetization79. Although SCF is not a superfluid, but as we will show by means of LSWT, its excitation spectrum is identically the same as the SF phase. In this phase due to the transverse staggered magnetization the easy plane U(1) symmetry reduces to the Z2 one80 and the translational symmetry is also broken. But, it is not a kind of solid since this phase possesses a gapless excitation and the longitudinal staggered magnetization is zero. This phase is not seen in the ground state phase diagram of the 2D mixed-spin system with staggered arrangement of spin-1/2 and spin-170. In the staggered 2D mixed-spin system the ground state phase diagram is completely symmetric with respect to J = 0.
By decreasing |J|, at small magnetic field, aside from the U(1) symmetry (which is completely broken in \(J > 0\) region and is decreased to the Z2 symmetry in J < 0 one), the translational symmetry of the system also breaks and a phase transition occurs from SF and SCF to the STS phase in which both diagonal and off-diagonal long range orders coexist in the system. The STS-SF and STS-SCF phase transitions are of first- or second-order which are attributed to the behavior of the low energy spin wave excitations. Indeed, the abrupt and smooth changes of the low energy excitations close to a transition point results in the discontinuous and continuous variations of the diagonal and off-diagonal order parameters.
In addition to the SF, SCF and STS phases, various kinds of stripe solids: ST(3/6), bST(4/6) and aST(5/6), respectively with fillings 3/6, 4/6, and 5/6 also appear in the phase diagram of the CAS system. The schematic pictures of these phases are plotted in the bottom of Fig. 2. Moreover, in order to show the fillings of these phases we have plotted in Fig. 3 the bosonic version of them.
In these solid phases, depending on the values of magnetic field, the translational symmetries of both subsystems or one of them break spontaneously, (in bST(4/6) the translational symmetry of the subsystem with spin τ and in aST(5/6) the translational symmetry of the subsystem with spin σ). The stripe solid orders with fillings 4/6 and 5/6 are the characteristics of our mixed-spin CAS system and are not seen in the phase diagram of uniform XXZ spin-1/2 models. For V2/V1 = 0.6, around J = 0, instead of bST(4/6), the system prefers to be in the MI(4/6) phase at moderate magnetic field where both the translational and U(1) symmetries are preserved in the system. Actually, for larger values of V2/V1, the V2 interactions try to make the spins τ and σ antiparallel, such that the translational symmetry of both subsystems preserves at moderate magnetic field. This behavior that the translational symmetry does not break even at large interactions, is the characteristic of the two-component systems with inter-components interaction, which has also been seen in the staggered mixed-spin system at V2/V1 < 0.570.
In order to obtain the more accurate phase diagrams, we employ CMFT with different cluster sizes. CMFT is an extension of the standard MF approximation in which clusters of multiple sites are used as an approximate system instead of single sites22,33,70,81. Recently we have generalized CMFT for the staggered 2D mixed-spin system70. Treating exactly the interactions within the cluster and including the interaction of spins outside the cluster as an effective field, one can partially take into account fluctuations around classical ground state as well as the effects of correlations of particles. Similar to our MF analysis, we assume a background with four-sublattice structure and embed a cluster of N C sites into this background. Instead of treating the many-body problem in the whole system, we consider the following effective cluster Hamiltonian:
where the interactions within clusters are given by H C while the interactions of spins inside the clusters with the rest of the system are included in \({H}_{\bar{c}}\). The Hamiltonian H C is given by Eq. (1) where the summations run over sites i, j ∈ C, and the Hamiltonian \({H}_{\bar{c}}\) is given in terms of the effective fields \({\overrightarrow{h}}_{i}^{eff}\) and \({\overrightarrow{g}}_{i}^{eff}\) as:
These effective fields are given in terms of CMFT ground state magnetizations, which are computed self-consistently (see Sec. 5.2 for detail).
In Fig. 2-middle we have plotted CMFT phase diagrams for the two different strengths of frustration: V2/V1 = 0.4 and 0.6, using cluster of four sites [CMFT-(2 × 2)]. Comparison between MF and CMFT phase diagrams shows that, due to the effects of quantum fluctuations the bST(4/6) and a ST(5/6) solid phases respectively convert to the ST(4/6) and ST(5/6) ones in which the translational symmetry of both subsystems break. It is also seen that except for the STS-SF transition lines there is no considerable changes in the MF phase diagrams in the presence of quantum fluctuations. The slight deviations of the STS-SF transition lines, for \({V}_{2}/{V}_{1} > 0.5\), at large magnetic field, are attributed to the large amount of quantum fluctuations at these boarders. In order to see the behavior of quantum fluctuations in each phases, we utilize LSWT and study the variations of spin waves’ number in all sublattices. Using Holstein-Primakoff (HP) transformations (see Eqs (17) and (18)) in Sec. 5.3), the Hamiltonian in Eq. (1) transforms to the following spin wave Hamiltonian:
where E0 is the classical MF energy, Hk is a square matrix in Fourier space, consisting the coefficients of bilinear terms, and ψk is a vector in terms of HP bosonic creation and annihilation operators. Dimensions of ψk and Hk depend on the number of sublattices in the MF ordered phases (see Sec. 5.3 for detail). Paraunitary diagonalization82 of Hk yields the excitation spectra in each phase, as well as the HP bosons’ number.
Amount of quantum fluctuations in different phases is given by HP bosons’ numbers n A , n B , n C and n D , respectively in the sublattices A, B, C and D. As it is seen from Fig. 4, in our mixed-spin model overall quantum fluctuations are not strong enough to destroy the MF orders and the MF predictions are reliable. However, at the second-order SF-STS and first-order SF-aST(5/6) transition lines, they are not negligible and we should take them into account for reaching to the accurate ground state phase diagrams. Also, since n A (n C ) is not equal to n B (n D ), the bST(4/6) and aST(5/6) solids convert to the ST(4/6) and ST(5/6) ones, respectively (see Fig. 2, CMFT-(2 × 2) phase diagram). It should be noticed that although quantum fluctuations break the translational symmetry of both subsystems but the fillings do not change. Moreover, in the presence of quantum fluctuations, part of the second-order STS-SF transition line below the ST(5/6) phase transforms to a first order one. This means that the MF prediction concerning the kind of transition order at this region is not correct, and more precise approaches should be employed to obtain the STS-SF critical and tricritical points.
We also investigate the behavior of both diagonal and off-diagonal order parameters considering clusters with larger sizes in CMFT. Employing clusters of eight spins [CMFT-(2 × 4)], we have computed the sublattices longitudinal and transverse magnetizations for different values of h and J (see Fig. 5). Since quantum fluctuations are strong around first order transition lines, some modifications around these lines are expected. These are clearly seen by comparison of the phases boarders of CMFT-(2 × 2) and CMFT-(2 × 4) in Fig. 5. The behavior of the order parameters indicates that the STS phase around ST(5/6) solid becomes narrower in the presence of quantum fluctuations. It seems that these regions tend to be disappeared when we use CMFT with larger cluster sizes. The STS phase appeared at small magnetic fields is however stable. The instability of STS phase at larger magnetic field can be explained as follows. The appearance of the STS phase is in fact the result of the competition between the staggered magnetization along z direction (diagonal order) and the total transverse magnetization (off-diagonal order). At large magnetic field, in the presence of quantum fluctuations, both order parameters decrease but the staggered magnetization is more sensitive and diminishes around the phase boarder. So the STS region at larger magnetic field is decreased in CMFT-(2 × 4).
The CMFT becomes exact when the cluster size goes to infinity. In practice, we are faced with computational limitations due to the increasing of clusters’ sizes and can not consider clusters of large sizes. We thus should employ other techniques such as quantum Monte Carlo simulations to obtain the exact phase diagram of the CAS model, at large magnetic field.
We have also plotted in Fig. 6, the spin wave excitation spectra in all phases of the CAS system. Number of excitation modes and their behavior depend on the number of sublattices as well as their longitudinal and transverse magnetizations. According to the translational symmetry of the CAS system, the primitive vectors in the SF and MI phases are \({\overrightarrow{a}}_{1}=a\hat{x}\) and \({\overrightarrow{a}}_{2}=a\hat{y}\), and there are two excitation modes in the system. However, when the translational symmetry breaks in x direction, as in the different solids, STS and SCF phases, the primitive vectors are \({\overrightarrow{a}}_{1}=2a\hat{x}\) and \({\overrightarrow{a}}_{2}=2a\hat{y}\) and the first Brillouin zone is folded in x direction (see Fig. 6-top). In these phases there exist four excitation modes in the system.
In STS and SF phases, as a result of the continues U(1) symmetry breaking, a gapless Goldstone mode with a roton-like minimum Δ appears in the excitation spectra (see Fig. 6). The slope of the line connecting this roton minimum to the origin in the critical velocity of the superfluid. The critical velocity in the STS phase depends inversely on the values of magnetic field. The critical velocity of the STS at lower magnetic field is larger than the one in the higher field. For J < 0, in SCF phase, the spins in each subsystem are antiparallel in xy plane. In this phase the translational symmetry breaks which results in a non-zero transverse staggered magnetization. Due to the translational symmetry breaking, four excitation modes appear in the energy spectrum. Moreover, since the U(1) symmetry also decreases to Z2 one in this phase the low energy excitation is gapless with linear dispersion around k = (0, 0), and the roton minimum is folded back to the origin.
Anisotropic CAS Models
In the previous section we obtained the ground state phase diagram of the isotropic CAS model in which the inter-chains and intra-chains NN interactions (V1) and hopping energies (J) were the same. In this section we consider two anisotropic CAS models: (i) a CAS system with different inter- and intra-chains NN interactions, and (ii) a CAS system with different inter- and intra-chains hopping energies, and investigate the effects of these anisotropies on the stability of the STS order appeared in the ground state phase diagram of the system. Following, we study these two anisotropic systems, separately.
Anisotropic CAS with different NN interactions
Let us consider the intra-chains and inter-chains NN interactions to be respectively V1 and \({V}_{1}^{\text{'}}\).
In the absence of intra-chains interaction, at V1 = 0, each 1D lattices are described by a XY Hamiltonian and no supersolid order appears in the ground state phase diagram of the system. In this case, when the inter-chains interactions are attractive (\({V}_{1}=\mathrm{0,}\,{V}_{1}^{\text{'}} < 0\)), at h = 0 the system is in the ST(3/6) phase and a phase transition occurs to the SF phase at \(J/|{V}_{1}^{\text{'}}|\,\approx \,0.23\). In the presence of magnetic field, the system enters the Full phase (not shown). When the inter-chains interactions are repulsive (\({V}_{1}=\mathrm{0,}\,{V}_{1}^{\text{'}} > 0\)) the system displays MI(4/6), in addition to the SF order, but no supersolidity occurs in the phase diagram of the model (see Fig. 7). This is due the fact that in the absence of intra-chains interactions, the translational symmetry of the chains preserves and the system is always MI(4/6) or SF, below the saturation field. It should be noticed that the MI(4/6) phase, appeared in the ground state phase diagram of the isotropic CAS model (see Fig. 2), is a result of competition between NNN interaction and magnetic field. However, in the mentioned anisotropic CAS model this phase emerges at h = V2 = 0 where the intra-chains interactions are absent (see left panel of Fig. 7).
In the absence of inter-chains NN interactions, at \({V}_{1}^{\text{'}}=0\), the anisotropic CAS system is composed of coupled 1D spin-1 and spin-1/2 XXZ models in a longitudinal magnetic field. In this system, when the intra-chains interactions are attractive (\({V}_{1}^{\text{'}}\) = 0, V1 < 0), at h = 0 the system is in the MI(4/6) phase and there is a phase transition to the SF order at J/|V1| ≈ 0.26, and in the presence of magnetic field, at h ≠ 0, trivially the system is in the fully polarized phase (not shown). However, when the intra-chains interactions are repulsive (\({V}_{1}^{\text{'}}=\mathrm{0,}\,{V}_{1} > 0\)), due to the breaking of the translational and U(1) symmetries, independent of V2, the system exhibits the STS phase. In addition to the STS phase, the ST(3/6) and ST(4/6) solids, and also SF orders also appear in the phase diagram of the system (see upper panels of Fig. 8). Comparison of the phase diagrams of V2 = 0 and V2 ≠ 0 in Fig. 8 shows that the presence of the NNN interactions V2 decreases the STS and ST(4/6) regions. It is surprising that in the absence of V2 the STS phase emerges in the phase diagram even at zero magnetic field. The amount of quantum fluctuations in this anisotropic CAS system is plotted in Fig. 8 (lower panels). As it is seen the fluctuations in STS phase are not strong and we expect the STS phases appeared in the CMFT phase diagrams to be stable in the presence of quantum fluctuations. This achievement is also confirmed by our CMFT-(2 × 4) results. According to our CMFT-(2 × 4) data (not shown), for \({V}_{1}^{\text{'}}\) = V2 = 0 although the SF-STS and STS-ST(3/6) critical points shift to the lower values of magnetic fields, but the STS region does not become narrower (see green stars in Fig. 8).
Above results indicate that, the presence of finite repulsive intra-chains interactions are necessary for the STS and ST phases to be emerged in the anisotropic CAS phase diagram. To check this idea we have also plotted in Fig. 9 the V1−\({V}_{1}^{\text{'}}\) phase diagram of the anisotropic CAS system, for the two different values of NNN interactions: V2 = 0, and 0.4 at h = 1 and J = 0.2. These figures show that in the presence of attractive intra-chains interactions, the two Mott insulating phases, MI(4/6) and Full, also appear in the phase diagram in addition to the SF order. However, there is no signature of STS and ST phases in this region in both values of NNN interaction. The STS and ST phases emerge in the V1 > 0 region, independent of the strengths of \({V}_{1}^{\text{'}}\) and V2, where the translational symmetry breaks in the presence of repulsive intra-chains interactions. The small amount of quantum fluctuations in these regions is the reason of the stability of all orders.
From Fig. 9 it is also seen that in the isotropic CAS model (\({V}_{1}={V}_{1}^{\text{'}}\)) it is impossible to find the STS phase in the absence of NNN interactions. Previous studies on the 2D Bose-Hubbard model with three-body-constrained bosons, show solid and superfluid phases in the presence of isotropic NN interactions and hopping terms83, but no supersolidity occurs in this system. Actually, NNN interactions can stabilize the supersolid phases in this model.
Anisotropic CAS with different hoppings
In this subsection we consider an anisotropic CAS model with different inter- and intra-chains hopping energies and investigate the effects of this anisotropy on the ground state phase diagram of the system. Let us consider the intra- and inter-chains hopping energies to be respectively J and J′. Suppose one of the hopping energies J or J′ to be zero. Our CMFT results show that the presence of inter-chains hoppings together with the repulsive intra-chains NN interactions are sufficient for the appearance of STS phase (see Fig. 10). The behavior of HP bosons’ number shows that in the STS phase the quantum fluctuations are exactly zero in the sublattices B and C, but considerable in A and D. This means that we should expect some modifications on the sublattices’ magnetizations in the presence of fluctuations. However, surprisingly we see that the order parameters computed by CMFT-(2 × 2) and -(2 × 4) are exactly the same (see Fig. 10) which is an indication of the stability of the STS phase in this system.
Minimal model
Based on the above results we conclude that, in the absence of inter-chains hoppings, the minimal mixed-spin CAS model for the supersolidity is given by the following Hamiltonian:
where, \({H}_{XXZ}^{\sigma }\) and \({H}_{XXZ}^{\tau }\) are respectively the spin-1/2 and spin-1 XXZ Hamiltonians. \({H}_{Z}^{\sigma \tau }\) is Ising Hamiltonian and couples the spin-1/2 and spin-1 chains. In bosonic language, the presence of off-site intra-components repulsive interactions and hopping energies together with the inter-components repulsive or attractive interactions, are sufficient to find the STS phase in the two-component hard-core Bose-Hubbard model. On the other hand, in the presence of NN inter-chains hoppings, the minimal mixed-spin CAS model for the supersolidity is, instead, spin-1 and spin-1/2 Ising chains interacting via a XY Hamiltonian. In two-component bosonic system inter-component hoppings could be possible by including different hyperfine states of the same atoms. So, in the case of possible inter-component hoppings, repulsive off-site intra-component interactions are necessary to find the STS phase.
Methods
In this section we explain the details of MF, CMFT and LSWT.
Mean field approximation
In the MF approach, we approximate local spins’ averages with MF order parameters. By defining the following MF order parameters: \(\langle {\sigma }_{{i}_{A}}^{x,y,z}\rangle ={m}_{A}^{x,y,z}\), \(\langle {\sigma }_{{i}_{B}}^{x,y,z}\rangle ={m}_{B}^{x,y,z}\), \(\langle {\tau }_{{i}_{C}}^{x,y,z}\rangle ={M}_{C}^{x,y,z}\), and \(\langle {\tau }_{{i}_{D}}^{x,y,z}\rangle ={M}_{D}^{x,y,z}\), where 〈…〉 denotes the expectation value on the MF ground state, the Hamiltonian in Eq. (1) is readily simplified to the MF Hamiltonian:
with
where J and V1 are NN hopping and interaction, respectively. Diagonalization of HMF leads to four coupled equations for the MF order parameters. Solving these equations self-consistently gives the sublattices’ magnetizations.
Cluster mean field theory
In CMFT we divide the lattice into clusters of N C sites. In CMFT-(2 × 2) each cluster includes N C = 4 spins, two spins 1 and two spins 1/2 (see Fig. 11). The effective magnetic fields seen by each spin inside a 2 × 2 cluster are given by:
where the order parameters \({m}_{A,B}^{x,y,z}\) and \({M}_{C,D}^{x,y,z}\) are the averages of spin operators on the CMFT-(2 × 2) ground state, which are computed self-consistently.
In CMFT-(2 × 4) each cluster includes N C = 8 spins, four spins 1 and four spins 1/2 (see Fig. 12). The effective magnetic fields seen by each spin inside a 2 × 4 cluster are given by:
where the order parameters \({m}_{A,B}^{x,y,z}\) and \({M}_{C,D}^{x,y,z}\) are given by
Here, the expectation values are on the CMFT-(2×4) ground state, which are computed self-consistently.
Linear spin wave theory
Let us perform the following local rotations:
where \( {\mathcal R} \) is the rotation matrix and ϕ i , θ i , φ i , and ϑ i are the Euler angles given by \(\cos \,{\theta }_{i}=\langle {\sigma }_{i}^{z}\rangle /\sigma \), \(\tan \,{{\varphi }}_{i}=\langle {\sigma }_{i}^{y}\rangle /\langle {\sigma }_{i}^{x}\rangle \), \(\cos \,{\vartheta }_{j}=\langle {\tau }_{j}^{z}\rangle /\tau \) and \(\tan \,{{\phi }}_{j}=\langle {\tau }_{j}^{y}\rangle /\langle {\tau }_{j}^{x}\rangle \). Here, \(\langle {\sigma }_{i}^{\alpha }\rangle \) and \(\langle {\tau }_{i}^{\alpha }\rangle \) are the sublattices’ magnetizations computed by means of MF theory. Using the following Holstein-Primakoff (HP) transformations:
and
the Hamiltonian in Eq. (1) is expressed in terms of HP bosons as:
where E0 is the classical MF energy, Hk is a square matrix consists of the coefficient of bilinear terms, and ψk is a vector in terms of HP bosonic creation and annihilation operators. The matrix Hk has the following general form:
where the dimensions of A, B, ψk and Hk depend on the number of sublattices in the MF ordered phases. Paraunitary diagonalization82 of Hk yields the excitation spectra in each phase.
Now we obtain the spin wave Hamiltonian Hk and the vector ψk for all phases. In all solids, supersolid and SCF phases, according to the translational symmetry breaking, the classical background has a four-sublattice structure and four HP bosons should be employed to attain the excitation spectra of these phases, (see Fig. 3). In this respect, we consider a general background and divide the subsystem with spin σ to two sublattices with HP bosons \(\hat{a}\) and \(\hat{b}\), and the subsystem with spin τ to two sublattices with bosons \(\hat{c}\) and \(\hat{d}\). The vector ψk and the matrices A and B are given by
where
with
Here, α and β are 1, 2, 3, 4, and m and n are the sublattices’ labels: a, b, c and d. The coefficients \({f}_{mn}^{\beta }\), \({g}_{mn}^{\alpha }\) and e m are given in terms of θ m , θ n and ϕ m , ϕ n as follows:
where S m and S n are the spins of sublattices m and n, respectively.
In the MI(4/6), Full and SF phases where the translational symmetry of the original Hamiltonian is preserved, the classical background has a two-sublattice structure and the excitations of these phases are achieved by defining the two HP bosons: \(\hat{a}\) and \(\hat{d}\) for the sublattices with the spin σ and τ respectively. The Fourier transformation of the bilinear term of the LSW Hamiltonian is written by:
and A and B matrices are given by
where
Summary and Conclusion
To summarize, in the present paper, employing three analytical and numerical approaches, MF approximation, CMFT with different cluster sizes and LSWT, we have studied the ground state phases of a 2D mixed-spin system of coupled alternating spin chains described by the spin Hamiltonian in Eq. (1). Our study, indicates that the CAS system displays a rich ground state phase diagram including STS and SCF phases in addition to the different solids, SF and MI phases. We have also considered two kinds of anisotropic CAS model, (i) CAS model with different intra-chains and inter-chains NN interactions and (ii) CAS model with different intra-chains and inter-chains hoppings, and investigated the effects of these anisotropies on the ground state phases. We have demonstrated that the emergence of the STS phase strongly depends on the strength of intra-chains NN interactions and hopping energies. We have shown that, for the systems with the inter-chains interactions, the repulsive intra-chains hoppings are necessary and sufficient for stripe supersolidity. In this case the minimal two dimensional mixed-spin model is a system of spin-1 and spin-1/2 XXZ chains, interacting via an Ising Hamiltonian. However, in the presence of inter-chains hoppings, the STS phase emerges even in the absence of intra-chains interactions, and a system of coupled Ising chains is the minimal model.
Our mixed-spin model is equivalent to a bosonic system of hard-core and semi-hard-core bosons and could be realized in coupled one dimensional optical lattices by alternatively changing the optical depth. Study of temperature phase diagram as well as thermodynamic properties of the CAS system and also study of the ground state phase diagram with other approaches are left for future work.
References
Matsuda, H. & Tsuneto, T. Off-Diagonal Long-Range Order in Solids. Progress of Theoretical Physics Supplement 46, 411–436 (1970).
Andreev, A. F. & Lifshitz, I. M. Quantum theory of defects in crystals. Physics-Uspekhi 13, 670–670 (1971).
Thouless, D. J. The flow of a dense superfluid. Annals of Physics 52, 403–427 (1969).
Liu, K.-S. & Fisher, M. E. Quantum lattice gas and the existence of a supersolid. Journal of Low Temperature Physics 10, 655–683 (1973).
Kim, E. & Chan, M. H. W. Probable observation of a supersolid helium phase. Nature 427, 225–227 (2004).
Kim, E. & Chan, M. H. W. Observation of superflow in solid helium. Science 305, 1941–1944 (2004).
Li, J.-R. et al. Stripe phase with supersolid properties in spin–orbit-coupled Bose–Einstein condensates. Nature 543, 91–94 (2017).
Leonard, J., Morales, A., Zupancic, P., Esslinger, T. & Donner, T. Supersolid formation in a quantum gas breaking a continuous translational symmetry. Nature 543, 87–90 (2017).
Jaksch, D., Bruder, C., Cirac, J. I., Gardiner, C. W. & Zoller, P. Cold bosonic atoms in optical lattices. Physical Review Letters 81, 3108 (1998).
Batrouni, G. G., Hebert, F. & Scalettar, R. T. Supersolid Phases in the One-Dimensional Extended Soft-Core Bosonic Hubbard Model. Physical Review letters 97, 087209 (2006).
Mathey, L., Danshita, I. & Clark, C. W. Creating a supersolid in one-dimensional Bose mixtures. Physical Review A 79, 011602 (2009).
Burnell, F. J., Parish, M. M., Cooper, N. R. & Sondhi, S. L. Devil’s staircases and supersolids in a one-dimensional dipolar Bose gas. Physical Review B 80, 174519 (2009).
Goral, K., Santos, L. & Lewenstein, M. Quantum Phases of Dipolar Bosons in Optical Lattices. Physical Review letters 88, 170406 (2002).
Kovrizhin, D. L., Venketeswara Pai, G. & Sinha, S. Density wave and supersolid phases of correlated bosons in an optical lattice. EPL (Europhysics Letters) 72, 162 (2005).
Sengupta, P., Pryadko, L. P., Alet, F., Troyer, M. & Schmid, G. Supersolids versus Phase Separation in Two-Dimensional Lattice Bosons. Physical Review letters 94, 207202 (2005).
Scarola, V. W. & Das Sarma, S. Quantum Phases of the Extended Bose-Hubbard Hamiltonian: Possibility of a Supersolid State of Cold Atoms in Optical Lattices. Physical Review letters 95, 033003 (2005).
Yi, S., Li, T. & Sun, C. P. Novel Quantum Phases of Dipolar Bose Gases in Optical Lattices. Physical Review letters 98, 260405 (2007).
Danshita, I. & Sa de Melo, C. A. R. Stability of Superfluid and Supersolid Phases of Dipolar Bosons in Optical Lattices. Physical Review letters 103, 225301 (2009).
Capogrosso-Sansone, B., Trefzger, C., Lewenstein, M., Zoller, P. & Pupillo, G. Quantum Phases of Cold Polar Molecules in 2D Optical Lattices. Physical Review letters 104, 125301 (2010).
Danshita, I. & Yamamoto, D. Critical velocity of flowing supersolids of dipolar Bose gases in optical lattices. Physical Review A 82, 013645 (2010).
Mila, F., Dorier, J. & Schmidt, K. P. Supersolid Phases of Hardcore Bosons on the Square Lattice: Correlated Hopping, Next-Nearest Neighbor Hopping and Frustration. Progress of Theoretical Physics Supplement 176, 355–374 (2008).
Yamamoto, D., Masaki, A. & Danshita, I. Quantum phases of hardcore bosons with long-range interactions on a square lattice. Physical Review B 86, 054516 (2012).
Wessel, S. Phase diagram of interacting bosons on the honeycomb lattice. Physical Review B 75, 174301 (2007).
Gan, J. Y. et al. Extended Bose-Hubbard model on a honeycomb lattice. Physical Review B 75, 214509 (2007).
Wessel, S. & Troyer, M. Supersolid Hard-Core Bosons on the Triangular Lattice. Physical Review letters 95, 127205 (2005).
Heidarian, D. & Damle, K. Persistent Supersolid Phase of Hard-Core Bosons on the Triangular Lattice. Physical Review letters 95, 127206 (2005).
Melko, R. G. et al. Supersolid Order from Disorder: Hard-Core Bosons on the Triangular Lattice. Physical Review letters 95, 127207 (2005).
Boninsegni, M. & Prokof’ev, N. Supersolid Phase of Hard-Core Bosons on a Triangular Lattice. Physical Review letters 95, 237204 (2005).
Hassan, S. R., de Medici, L. & Tremblay, A.-M. S. Supersolidity, entropy, and frustration: t-t’-V model of hard-core bosons on the triangular lattice. Physical Review B 76, 144420 (2007).
Sen, A., Dutt, P., Damle, K. & Moessner, R. Variational Wave-Function Study of the Triangular Lattice Supersolid. Physical Review letters 100, 147204 (2008).
Pollet, L., Picon, J. D., Buchler, H. P. & Troyer, M. Supersolid Phase with Cold Polar Molecules on a Triangular Lattice. Physical Review letters 104, 125302 (2010).
Zhang, X.-F., Dillenschneider, R., Yu, Y. & Eggert, S. Supersolid phase transitions for hard-core bosons on a triangular lattice. Physical Review B 84, 174515 (2011).
Yamamoto, D., Danshita, I. & Sa de Melo, C. A. R. Dipolar bosons in triangular optical lattices: Quantum phase transitions and anomalous hysteresis. Physical Review A 85, 021601 (2012).
Isakov, S. V., Wessel, S., Melko, R. G., Sengupta, K. & Kim, Y. B. Hard-Core Bosons on the Kagome Lattice: Valence-Bond Solids and Their Quantum Melting. Physical Review letters 97, 147202 (2006).
Bonnes, L. & Wessel, S. Generic first-order versus continuous quantum nucleation of supersolidity. Physical Review B 84, 054510 (2011).
Kuno, Y., Suzuki, K. & Ichinose, I. Superfluidity and solid order in a two-component Bose gas with dipolar interactions in an optical lattice. Physical Review A 90, 063620 (2014).
Trefzger, C., Menotti, C. & Lewenstein, M. Pair-Supersolid Phase in a Bilayer System of Dipolar Lattice Bosons. Physical Review letters 103, 035304 (2009).
Trefzger, C., Alloing, M., Menotti, C., Dubin, F. & Lewenstein, M. Quantum magnetism and counterflow supersolidity of up–down bosonic dipoles. New Journal of Physics 12, 093008 (2010).
Yamamoto, K., Todo, S. & Miyashita, S. Successive phase transitions at finite temperatures toward the supersolid state in a three-dimensional extended Bose-Hubbard model. Physical Review B 79, 094503 (2009).
Xi, B., Ye, F., Chen, W., Zhang, F. & Su, G. Global phase diagram of three-dimensional extended Boson Hubbard model: A continuous-time quantum Monte Carlo study. Physical Review B 84, 054512 (2011).
Ohgoe, T., Suzuki, T. & Kawashima, N. Commensurate Supersolid of Three-Dimensional Lattice Bosons. Physical Review letters 108, 185302 (2012).
Batista, C. D. & Ortiz, G. Algebraic approach to interacting quantum systems. Advances in Physics 53, 1–82 (2004).
Giamarchi, T., Ruegg, C. & Tchernyshyov, O. Bose–Einstein condensation in magnetic insulators. Nature Physics 4, 198–204 (2008).
Kreisel, A., Hasselmann, N. & Kopietz, P. Probing anomalous longitudinal fluctuations of the interacting bose gas via Bose-Einstein condensation of magnons. Physical Review letters 98, 067203 (2007).
Peters, D., McCulloch, I. P. & Selke, W. Spin-one Heisenberg antiferromagnetic chain with exchange and single-ion anisotropies. Physical Review B 79, 132406 (2009).
Rossini, D., Giovannetti, V. & Fazio, R. Spin-supersolid phase in Heisenberg chains: A characterization via matrix product states with periodic boundary conditions. Physical Review B 83, 140411 (2011).
Sengupta, P. & Batista, C. D. Spin supersolid in an anisotropic spin-one heisenberg chain. Physical Review letters 99, 217205 (2007).
Ng, K.-K. & Lee, T. K. Supersolid Phase in Spin Dimer XXZ Systems under a Magnetic Field. Physical Review letters 97, 127204 (2006).
Chen, P., Lai, C.-Y. & Yang, M.-F. Field-induced spin supersolidity in frustrated S = 1 2 spin-dimer models. Physical Review B 81, 020409 (2010).
Laflorencie, N. & Mila, F. Quantum and thermal transitions out of the supersolid phase of a 2D quantum antiferromagnet. Physical Review letters 99, 027202 (2007).
Schmidt, K. P., Dorier, J., Lauchli, A. M. & Mila, F. Supersolid Phase Induced by Correlated Hopping in Spin-1/2 Frustrated Quantum Magnets. Physical Review letters 100, 090401 (2008).
Guo, L.-F., Chen, Q.-H. & Li, P. Stripe phases in a frustrated spin-1/2 dimer Heisenberg model. International Journal of Modern Physics B 28, 1450143 (2014).
Thomson, S. J. & Kruger, F. Bose and Mott glass phases in dimerized quantum antiferromagnets. Physical Review B 92, 180201 (2015).
Murakami, Y., Oka, T. & Aoki, H. Supersolid states in a spin system: Phase diagram and collective excitations. Physical Review B 88, 224404 (2013).
Albuquerque, A. F., Laflorencie, N., Picon, J.-D. & Mila, F. Phase separation versus supersolid behavior in frustrated antiferromagnets. Physical Review B 83, 174421 (2011).
Wierschem, K. & Sengupta, P. Columnar antiferromagnetic order and spin supersolid phase on the extended Shastry-Sutherland lattice. Physical Review letters 110, 207207 (2013).
Picon, J.-D. et al. Mechanisms for spin supersolidity in S = 1/2 spin-dimer antiferromagnets. Physical Review B 78, 184418 (2008).
Ng, K.-K. & Yang, M.-F. Field-induced quantum phases in a frustrated spin-dimer model: A sign-problem-free quantum Monte Carlo study. Physical Review B 95, 064431 (2017).
Momoi, T. & Totsuka, K. Magnetization plateaus of the Shastry-Sutherland model for SrCu 2 (BO 3) 2: Spin-density wave, supersolid, and bound states. Phy. Rev. B 62, 15067 (2000).
Sengupta, P. & Batista, C. D. Ground state and thermal transitions in field-induced spin-supersolid phase. Journal of Applied Physics 103, 07C709 (2008).
Sengupta, P. & Batista, C. D. Field-induced supersolid phase in spin-one Heisenberg models. Physical Review letters 98, 227201 (2007).
Toth, T. A., Lauchli, A. M., Mila, F. & Penc, K. Competition between two-and three-sublattice ordering for S = 1 spins on the square lattice. Physical Review B 85, 140403 (2012).
Su, L., Wierschem, K. & Sengupta, P. Magnetic phases in the S = 1 Shastry-Sutherland model with uniaxial anisotropy. Physical Review B 89, 245432 (2014).
Ueda, H. T. & Momoi, T. Nematic phase and phase separation near saturation field in frustrated ferromagnets. Physical Review B 87, 144417 (2013).
Selke, W. Multicritical points in the three-dimensional XXZ antiferromagnet with single-ion anisotropy. Physical Review E 87, 014101 (2013).
Trumper, A. E. & Gazza, C. Antiferromagnetically coupled alternating spin chains. Physical Review B 64, 134408 (2001).
Abouie, J. & Langari, A. Cumulant expansion for ferrimagnetic spin (S1, s2) systems. Physical Review B 70, 184416 (2004).
Abouie, J., Ghasemi, S. A. & Langari, A. Thermodynamic properties of ferrimagnetic spin chains in the presence of a magnetic field. Physical Review B 73, 014411 (2006).
Langari, A., Abouie, J., Asadzadeh, M. Z. & Rezai, M. Phase diagram of the XXZ ferrimagnetic spin-(1/2, 1) chain in the presence of transverse magnetic field. Journal of Statistical Mechanics: Theory and Experiment 2011, P08001 (2011).
Heydarinasab, F. & Abouie, J. Inhomogeneous hard-core bosonic mixture with checkerboard supersolid phase: Quantum and thermal phase diagram. Physical Review B 96, 104406 (2017).
Hebert, F. et al. Quantum phase transitions in the two-dimensional hardcore boson model. Physical Review B 65, 014513 (2001).
Batrouni, G. G. & Scalettar, R. T. Phase Separation in Supersolids. Physical Review letters 84, 1599–1602 (2000).
Tsurkan, V. et al. Ultra-robust high-field magnetization plateau and supersolidity in bond-frustrated MnCr 2 S 4. Science Advances 3, e1601982 (2017).
Zheng, Y.-Z. et al. Spin-frustrated complex, [FeIIFeIII (trans–1,4–cyclohexanedicarboxylate)1.5]∞: interplay between single-chain magnetic behavior and magnetic ordering. Inorganic chemistry 48, 2028–2042 (2009).
Li, X., Zhao, E. & Liu, W. V. Topological states in a ladder-like optical lattice containing ultracold atoms in higher orbital bands. Nature communications 4, 1523 (2013).
Safavi-Naini, A., Capogrosso-Sansone, B. & Kuklov, A. Quantum phases of hard-core dipolar bosons in coupled one-dimensional optical lattices. Physical Review A 90, 043604 (2014).
Singh, M., Mondal, S., Sahoo, B. K. & Mishra, T. Quantum phases of constrained dipolar bosons in coupled one-dimensional optical lattices. Physical Review A 96, 053604 (2017).
Matsubara, T. & Matsuda, H. A lattice model of liquid helium, I. Progress of Theoretical Physics 16, 569–582 (1956).
Kuklov, A. B. & Svistunov, B. V. Counterflow superfluidity of two-species ultracold atoms in a commensurate optical lattice. Physical Review letters 90, 100401 (2003).
Chen, P. & Yang, M.-F. Quantum phase transitions in a two-species hard-core boson Hubbard model in two dimensions. Physical Review B 82, 180510 (2010).
Yamamoto, D. Correlated cluster mean-field theory for spin systems. Physical Review B 79, 144427 (2009).
Colpa, J. H. P. Diagonalization of the quadratic boson Hamiltonian. Physica A: Statistical Mechanics and its Applications 93, 327–353 (1978).
Chen, Y.-C., Ng, K.-K. & Yang, M.-F. Quantum phase transitions in the attractive extended Bose-Hubbard model with a three-body constraint. Physical Review B 84, 092503 (2011).
Author information
Authors and Affiliations
Contributions
F.H. performed the calculations. J.A. wrote the manuscript. All authors approved the final manuscript. J.A. supervised the research.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Heydarinasab, F., Abouie, J. Spin supersolid phase in coupled alternating spin chains. Sci Rep 8, 7955 (2018). https://doi.org/10.1038/s41598-018-26169-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-26169-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.