Few years ago, application of the mean field Bethe scheme on a given system was shown to produce a systematic change of the system intrinsic symmetry. For instance, once applied on a ferromagnet, individual spins are no more equivalent. Accordingly a new loopwise mean field theory was designed to both go beyond the one site Weiss approach and yet preserve the initial Hamitonian symmetry. This loopwise scheme is applied here to solve the triangular antiferromagnetic Ising model. It is found to yield Wannier's exact result of no ordering at non-zero temperature. No adjustable parameter is used. Simultaneously a non-zero critical temperature is obtained for the triangular Ising ferromagnet. This simple mean field scheme opens a new way to tackle random systems. PACS . 75.25.+z Spin arrangements in magnetically ordered materials (including neutron and spin-polarized electron studies, synchrotron-source X-ray scattering, etc.) -- 05.50.+q Lattice theory and statistics (Ising, Potts, etc.) -- 75.50.-y Studies of specific magnetic materials 1
Motivated by a recent proposal of a Bethe approximation for the triangular Ising antiferromagnet [Phys. Rev. B 56, 8241 (1997)], which seems to predict a disordered phase at any temperature in zero field, we analyze in some detail several mean-field-like approximations for this model, namely, the Bethe approximation itself, the cluster variation method, and the hard-spin mean-field theory. We show (i) that the disordered phase predicted by the Bethe approximation is unphysical at low enough temperature because of a negative entropy; (ii) how the results of the cluster variation method (namely, zero-temperature entropy and critical temperature of the spurious transition) converge to the exact ones for increasing cluster size; (iii) that it is possible to construct a cluster variation approximation which yields a disordered phase which is stable down to zero temperature; (iv) a few, so far unknown, zero-temperature results (entropy and internal energy) of the hard-spin mean-field theory.
Motivated by a recent proposal of a Bethe approximation for the triangular Ising antiferromagnet [Phys. Rev. B {\bf 56}, 8241 (1997)], which seems to predict a disordered phase at any temperature in zero field, we analyze in some detail several mean-field like approximations for this model, namely the Bethe approximation itself, the cluster variation method and the hard-spin mean-field theory. We show: (i) that the disordered phase predicted by the Bethe approximation is unphysical at low enough temperature because of a negative entropy; (ii) how the results of the cluster variation method (namely, zero temperature entropy and critical temperature of the spurious transition) converge to the exact ones for increasing cluster size; (iii) that it is possible to construct a cluster variation approximation which yields a disordered phase which is stable down to zero temperature; (iv) a few, so far unknown, zero temperature results (entropy and internal energy) of the hard-spin mean-field theory.
We use a recently developed quantum simulation approach to study the properties of a three-dimensional Ising model consisting of S = 1/2 quantum spins localized at the sites of a simple cubic lattice. We assume nearest-neighbor interaction between spins with an exchange interaction that can be either ferromagnetic or antiferromagnetic. It is found that the computational method quickly converges towards the expected equilibrium spin configurations. The resulting spontaneous magnetization curves corresponding to the two types of magnetic interactions under consideration were found to be almost identical to the ones obtained via quantum mean field theory at all temperatures. The derived total energies, total free energies, magnetic entropies and specific heats per mole of spins show no sizeable differences from known theoretical values. Furthermore, the results of the simulations for two different 3D Ising systems containing 4×4×4 and 20×20×20 spins localized at the sites of a simple cubic lattice were found to be almost identical to each other. This finding suggests that the self-consistent algorithm approach of the current simulation method allows one to obtain the physical bulk properties of a large magnetic system by relying on simulations of a much smaller spin system sample. Therefore, the method presently considered appears to be not only very accurate as gauged by comparison to mean field theory, but also able to greatly increase the speed of simulations.
We study a triangular Ising model subject to nearest-neighbor pair interactions
K, an external magnetic field
H, and triangular three-spin interactions
P.
K>0(<0) corresponds to ferromagnetic (antiferromagnetic) coupling. The phase diagram and thermodynamics are invariant under $(K,H,P)\ensuremath{\leftrightarrow}(K,\ensuremath{-}H,\ensuremath{-}P)$. A simple position-space renormalization group is formulated, based on three interpenetrating three-spin cells, which preserves the three-sublattice symmetry of the antiferromagnetic phase. The complete (
K,H,P) phase diagram consists of five interconnecting transition surfaces (two second-order and three first-order), which bound the four single-phase regions (\ensuremath{\uparrow}\ensuremath{\uparrow}\ensuremath{\uparrow}, \ensuremath{\uparrow}\ensuremath{\uparrow}\ensuremath{\downarrow}, \ensuremath{\uparrow}\ensuremath{\downarrow}\ensuremath{\downarrow}, and \ensuremath{\downarrow}\ensuremath{\downarrow}\ensuremath{\downarrow}). The edges of these surfaces form loci of second-order, critical-endpoint, and Baxter-Wu transitions. When
P=0, this system is known to provide a useful model of the order-disorder transition in monolayer films adsorbed on a hexagonal substrate. Calculation of the magnetization
M(K,H) and heat capacity
cM(K,H) of the Ising model allows quantitative comparison with experimental data for helium on graphite. The observed dependence
Tc(n) of the transition temperature on coverage is well fit with a single effective antiferromagnetic pair interaction. The corresponding calculated heat capacities (no additional adjustable parameters) are in remarkable agreement with experimental data near and below
Tc.
We use an effective field model which explicitly incorporates correlations to determine the critical temperatures for Ising ferromagnets on the common lattices. We specifically retain only the first nontrivial odd correlation function, i.e., the triplet correlation function, while neglecting higher order correlations. Within this approximation, assuming the equivalence of triplet correlations involving both central and noncentral sites in the nearest neighbor cluster yields values for the critical temperature below the Bethe approximations for z⩾4. Likewise, assuming the inequivalence of triplet correlations involving the central and noncentral sites yields values of the critical temperature only slightly higher than the Bethe approximation for z⩾4. For lattices with four- and six-fold coordination this equivalence and inequivalence of triplet correlation functions is sufficient to cause the critical temperatures to approach their two- and three-dimensional values. For instance, for quadratic lattices we obtain values for the critical temperature of 2.5284 and 2.9796 as compared to the exact and series results of 2.2692 and 2.7044 for the simple quadratic and diamond lattices, respectively, and the Bethe approximation of 2.8854.
A two-dimensional Ising model on a particular topological structure, a ``triangles-in-triangles'' Kagom\'e structure, which is called a triangular-Kagom\'e lattice, is investigated. In this model, we consider two interaction parameters
Jaa and
Jab, corresponding to the interactions among spins on the inner sites and between spins on the inner and vertex sites of the triangular in the Kagom\'e lattice, respectively. By summing over all spins at inner sites in the partition function, we arrive at the partition function of the Kagom\'e lattice. The effective interaction of the corresponding Kagom\'e lattice is always ferromagnetic, even for antiferromagnetic
Jaa. The critical properties of the system depend only on ${J}_{aa}∕\ensuremath{\mid}{J}_{ab}\ensuremath{\mid}$. When ${J}_{aa}∕\ensuremath{\mid}{J}_{ab}\ensuremath{\mid}>\ensuremath{-}1$, the system has long range order at low temperature. However, when ${J}_{aa}∕\ensuremath{\mid}{J}_{ab}\ensuremath{\mid}<\ensuremath{-}1$, the partition function of the triangular-Kagom\'e lattice can be related to that of the Kagom\'e lattice with effective interaction at a temperature higher than its critical temperature. Therefore, the system will be in paramagnetic phase at all temperatures for ${J}_{aa}∕\ensuremath{\mid}{J}_{ab}\ensuremath{\mid}<\ensuremath{-}1$. The phase diagram for ${J}_{aa}∕\ensuremath{\mid}{J}_{ab}\ensuremath{\mid}$ is given exactly.