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.
Using extensive quantum Monte Carlo simulations, we study a minimum model of interacting hard-core bosons on a square lattice with hoppings of different lengths, featuring nearest-neighbor hopping (
t1), anisotropic next-nearest-neighbor hopping (${t}_{2}^{\ensuremath{'}}$), and nearest-neighbor repulsion (
V1). The paradigmatic checkerboard supersolid (CSS) phase emerges as ${t}_{2}^{\ensuremath{'}}$ turns on, with the solid order being characterized by ordering vector $\mathbf{q}=(\ensuremath{\pi},\ensuremath{\pi})$. This serves as a rare example of the CSS phase which is obtained by doping a checkerboard solid and harbors spontaneously broken gauge and translational symmetries. A first-order supersolid-to-superfluid transition is observed. Moreover, we find a solid order-by-thermal disorder behavior together with a superfluid-to-solid transition upon increasing temperature. The underlying picture of the order-by-disorder phenomenon is figured out within the framework of the entropy effect.
The phase diagrams of Ising antiferromagnets in a magnetic field
H are investigated for various values of the ratio
R between nearest- and next-nearest-neighbor interaction. While meanfield approximations and the existing real-space renormalization-group treatments yield phase diagrams which are sometimes even qualitatively incorrect, accurate results are obtained from Monte Carlo calculations. For
R<0 only an antiferromagnetically ordered phase exists. Its transition to the disordered phase is first order for temperatures below the tricritical point (
Tt,
Ht). For $R\ensuremath{\rightarrow}0$ also ${T}_{t}\ensuremath{\rightarrow}0$. For
R=0 we find very good agreement with the results of M\"uller-Hartmann and Zittartz. For
R>0 and
H1<H<H2 a new phase with anomalous high ground-state degeneracy is found (two sublattices have only one-dimensional order). These sublattices undergo order-disorder transitions at
T=0, such that for
T>0 one is left with a "superantiferromagnetic" phase. At low temperatures in this phase a pronounced tendency is observed to form a simpler (2 \ifmmode\times\else\texttimes\fi{} 2) superstructure but with many antiphase domain boundaries. For $R\ensuremath{\rightarrow}\frac{1}{2}$ and
H<H1 the regime of the antiferromagnetic phases goes to zero temperature, while for
R>21 the superantiferromagnetic phase exists also for
H<H1. The order-disorder transition associated with this phase seems to have non-Ising critical exponents which vary as a function of
R and
H. Estimates for the exponents lead us to suggest that Suzuki's "weak universality" is valid. The behavior of the model at
T=0 is related to known results on hard-core lattice gases. It is shown that it is useful to interpret the transitions at
T=0 as generalized percolation transitions. Since the model may have applications to adsorbate phases in registered structures at (100) surfaces of cubic crystals, the transcription of our results to temperature-coverage phase diagrams and adsorption isotherms is discussed in detail, and possible experimental applications are mentioned.
Using the parallel tempering algorithm and graphics processing unit accelerated techniques, we have performed large-scale Monte Carlo simulations of the Ising model on a square lattice with antiferromagnetic (repulsive) nearest-neighbor and next-nearest-neighbor interactions of the same strength and subject to a uniform magnetic field. Both transitions from the (2x1) and row-shifted (2x2) ordered phases to the paramagnetic phase are continuous. From our data analysis, re-entrance behavior of the (2x1) critical line and a bicritical point which separates the two ordered phases at T=0 are confirmed. Based on the critical exponents we obtained along the phase boundary, Suzuki's weak universality seems to hold.
The Ising square lattice model with nearest-neighbor (nn) interactions (
J1) is one of the few exactly solvable models [1]. Adding next-neareast- neighbor (nnn) interactions (
J2) or a magnetic field (or both) leads to the non solvability of the model and only some approximate solutions are possible. In this brief report we will review some results obtained within effective field theory. We will show that besides second-order transitions there are also lines of first-order transitions and the coordinates of tricritical points are calculated.
We present a transfer-matrix scaling study of a hard-square lattice-gas model with an exclusion shell extending up to second-nearest neighbors, and with a uniaxial finite repulsion active at third-nearest-neighbor separations. The phase diagram, as a function of temperature and chemical potential, shows three distinct regions bounded by two merging critical lines. At low temperature the system undergoes a transition to an ordered phase whose c(4\ifmmode\times\else\texttimes\fi{}2) structure is characterized---within this calculation---by the asymptotic degeneracy of the four largest eigenvalues of the transfer matrix. At high temperature this phase is reached from the disordered state via an intermediate, degenerate phase which, on average, exhibits long-range order in one dimension only. These two transition regimes are found to be separated by a multicritical point which is likely to arise from the competition between interaction anisotropy and thermal-disordering effects. Evidence of anisotropic critical behavior is found in the correlation lengths exponent ratios, ${\ensuremath{\nu}}_{?}$/${\ensuremath{\nu}}_{\ensuremath{\perp}}$, which appear to vary continuously along both transition lines. However, the accuracy of the predictions we make for the thermodynamics of the model is still limited by sensitive finite-size effects which are presumably emphasized by the relatively extended range of the interactions.