No Arabic abstract
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.
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 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.
We investigate theoretically and experimentally the static magnetic properties of single crystals of the molecular-based Single-Chain Magnet (SCM) of formula [Dy(hfac)$_{3}$NIT(C$_{6}$H$_{4}$OPh)]$_{infty}$ comprising alternating Dy$^{3+}$ and organic radicals. A peculiar inversion between maxima and minima in the angular dependence of the magnetic molar susceptibility $chi_{M}$ occurs on increasing temperature. Using information regarding the monomeric building block as well as an {it ab initio} estimation of the magnetic anisotropy of the Dy$^{3+}$ ion, this anisotropy-inversion phenomenon can be assigned to weak one-dimensional ferromagnetism along the chain axis. This indicates that antiferromagnetic next-nearest-neighbor interactions between Dy$^{3+}$ ions dominate, despite the large Dy-Dy separation, over the nearest-neighbor interactions between the radicals and the Dy$^{3+}$ ions. Measurements of the field dependence of the magnetization, both along and perpendicularly to the chain, and of the angular dependence of $chi_{M}$ in a strong magnetic field confirm such an interpretation. Transfer matrix simulations of the experimental measurements are performed using a classical one-dimensional spin model with antiferromagnetic Heisenberg exchange interaction and non-collinear uniaxial single-ion anisotropies favoring a canted antiferromagnetic spin arrangement, with a net magnetic moment along the chain axis. The fine agreement obtained with experimental data provides estimates of the Hamiltonian parameters, essential for further study of the dynamics of rare-earths based molecular chains.
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.