No Arabic abstract
Using the parallel tempering algorithm and GPU accelerated techniques, we have performed large-scale Monte Carlo simulations of the Ising model on a square lattice with antiferromagnetic (repulsive) nearest-neighbor(NN) and next-nearest-neighbor(NNN) 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, reentrance 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, Suzukis weak universality seems to hold.
We report results of a Wang-Landau study of the random bond square Ising model with nearest- ($J_{nn}$) and next-nearest-neighbor ($J_{nnn}$) antiferromagnetic interactions. We consider the case $R=J_{nn}/J_{nnn}=1$ for which the competitive nature of interactions produces a sublattice ordering known as superantiferromagnetism and the pure system undergoes a second-order transition with a positive specific heat exponent $alpha$. For a particular disorder strength we study the effects of bond randomness and we find that, while the critical exponents of the correlation length $ u$, magnetization $beta$, and magnetic susceptibility $gamma$ increase when compared to the pure model, the ratios $beta/ u$ and $gamma/ u$ remain unchanged. Thus, the disordered system obeys weak universality and hyperscaling similarly to other two-dimensional disordered systems. However, the specific heat exhibits an unusually strong saturating behavior which distinguishes the present case of competing interactions from other two-dimensional random bond systems studied previously.
We implement a new and accurate numerical entropic scheme to investigate the first-order transition features of the triangular Ising model with nearest-neighbor ($J_{nn}$) and next-nearest-neighbor ($J_{nnn}$) antiferromagnetic interactions in ratio $R=J_{nn}/J_{nnn}=1$. Important aspects of the existing theories of first-order transitions are briefly reviewed, tested on this model, and compared with previous work on the Potts model. Using lattices with linear sizes $L=30,40,...,100,120,140,160,200,240,360$ and 480 we estimate the thermal characteristics of the present weak first-order transition. Our results improve the original estimates of Rastelli et al. and verify all the generally accepted predictions of the finite-size scaling theory of first-order transitions, including transition point shifts, thermal, and magnetic anomalies. However, two of our findings are not compatible with current phenomenological expectations. The behavior of transition points, derived from the number-of-phases parameter, is not in accordance with the theoretically conjectured exponentially small shift behavior and the well-known double Gaussian approximation does not correctly describe higher correction terms of the energy cumulants. It is argued that this discrepancy has its origin in the commonly neglected contributions from domain wall corrections.
The dynamics of the one-dimensional random transverse Ising model with both nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions is studied in the high-temperature limit by the method of recurrence relations. Both the time-dependent transverse correlation function and the corresponding spectral density are calculated for two typical disordered states. We find that for the bimodal disorder the dynamics of the system undergoes a crossover from a collective-mode behavior to a central-peak one and for the Gaussian disorder the dynamics is complex. For both cases, it is found that the central-peak behavior becomes more obvious and the collective-mode behavior becomes weaker as $K_{i}$ increase, especially when $K_{i}>J_{i}/2$ ($J_{i}$ and $K_{i}$ are exchange couplings of the NN and NNN interactions, respectively). However, the effects are small when the NNN interactions are weak ($K_{i}<J_{i}/2$).
We calculate the quantum phase diagram of the {it XXZ} chain with nearest-neighbor (NN) $J_{1}$ and next-NN exchange $J_{2}$ with anisotropies $Delta_{1}$ and $Delta_{2}$ respectively. In particular we consider the case $Delta_{1}=-Delta_{2}$ to interpolate between the {it XX} chain ($% Delta_{i}=0$) and the isotropic model with ferromagnetic $J_{2}$. For $% Delta_{1}<-1$, a ferromagnetic and two antiferromagnetic phases exist. For $| Delta_{i}| <1$, the boundary between the dimer and spin fluid phases is determined by the method of crossing of excitation spectra. For large $J_{2}/J_{1}$, this method seems to indicate the existence of a second spin fluid critical phase. However, an analysis of the spin stiffness and magnetic susceptibility for $Delta_{1}=Delta_{2}=1$ suggest that a small gap is present.
We study the one-dimensional Hubbard model with nearest-neighbor and next-nearest-neighbor hopping integrals by using the density-matrix renormalization group (DMRG) method and Hartree-Fock approximation. Based on the calculated results for the spin gap, total-spin quantum number, and Tomonaga-Luttinger-liquid parameter, we determine the ground-state phase diagram of the model in the entire filling and wide parameter region. We show that, in contrast to the weak-coupling regime where a spin-gapped liquid phase is predicted in the region with four Fermi points, the spin gap vanishes in a substantial region in the strong-coupling regime. It is remarkable that a large variety of phases, such as the paramagnetic metallic phase, spin-gapped liquid phase, singlet and triplet superconducting phases, and fully polarized ferromagnetic phase, appear in such a simple model in the strong-coupling regime.