Pair correlations on two one-dimensional Ising models with next-nearest-neighbor interactions are calculated exactly. For an appropriate range of values of the interaction energies, each of these models exhibits a disorder point, at which the nature of the short-range order changes abruptly, and at which there is a cusp in the graph of range of order vs. temperature. The zero-field susceptibilities are calculated exactly, and are found to have quite smooth temperature dependence.
Pair correlations on two one-dimensional Ising models with next-nearest-neighbor interactions are calculated exactly. For an appropriate range of values of the interaction energies, each of these models exhibits a disorder point, at which the nature of the short-range order changes abruptly, and at which there is a cusp in the graph of range of order vs. temperature. The zero-field susceptibilities are calculated exactly, and are found to have quite smooth temperature dependence.
The existence of a disorder point TD is established for some one- and two-dimensional Ising lattices with antiferromagnetic next-nearest-neighbor interactions. Within the disordered phase, the decay of axial next-nearest-neighbor pair correlations with increasing spin separation is positive monotonic exponential below TD and oscillatory with exponential envelope above TD. A general definition of a disorder point is formulated, and a method of estimating TD described.
For various Ising models two approaches are discussed, one is that of simulating lattices, also called gauging on exact equations, the other is that of calculating analytical expressions for the boundary free energy of Ising lattices. The first approach allows to conjecture a solution for some Ising models, that have sofar not been solved, once some exact partial result for the problem is known. The second approach aims at furnishing such a partial result in the form of a condition for the critical temperature. An example of such a result was recently given for the 2D Ising square lattice with nearest and next-nearest-neighbor interactions. The critical line that separates the ordered (ferromagnetic) phase from the disordered (paramagnetic) phase showed good agreement in the moderate and strong nearest neighbor coupling limit with several results obtained by Monte Carlo, transfer matrix and series expansion results. We extend the discussion of the critical line, finding an excellent fit, now also in other points, like the Padé point, as well as cusp behavior at the Onsager point where the lattice decouples into two 2D square lattices with only nearest-neighbor interaction. Combination of this result with a geometrical argument in the simulation approach leads to a critical exponent 2−2≈0.5858, comparable to the exponent 4/7≈0.5714 found from renormalization arguments.
Exact, closed expressions are obtained for the perpendicular and parallel susceptibilities of the one-dimensional spin-1/2 Ising model (i) with nearest neighbor interactions in a uniform magnetic field and (ii) with next-nearest neighbor interactions but without the magnetic field, in the thermodynamic limit. Graphs of the susceptibilities as a function of temperature for various magnetic field strengths in (i), and as a function of different ratios of the interaction strengths in (ii), are discussed.
The method of developing exact power-series expansions for the partition function ZN and related thermodynamic functions for the Ising model valid below the critical point is generalized to include exchange interactions between first-, second-, and third-neighbor pairs. Expansions of the spontaneous magnetization M0(T) and zero field susceptibility χ0(T) are derived through to sixth order of perturbation for the s.q. lattice, and through to fifth order of perturbation for the Δ′r, b.c.c., s.c., and f.c.c. lattices, when interactions J1SizSjz and J2SkzSlz are present between first- and second-neighbor spins, respectively (second-neighbor model). These expansions have also been obtained for the case where interactions of equal magnitude (J1 = J2 = J3) are present between first-, second-, and third-neighbor pairs (third-equivalent-neighbor model); here expansions through to fifth order of perturbation are obtained for the s.q., Δ′r, b.c.c., and s.c. lattices and through to fourth order for the f.c.c. lattice. The Padé approximant·procedure is employed to discuss the effects of an extended but finite range of interaction on the behavior of M0(T) and χ0(T) for T→Tc− as characterized by the critical exponents β and γ′, respectively. For the second-equivalent-neighbor model lattices, it is found that 0.122 ≤ β ≤ 0.134 in two dimensions, and that 0.308 ≤ β ≤ 0.328 in three dimensions; from which it is concluded that β remains unchanged from its value in the nearest-neighbor model. The corresponding limits for γ′ in three dimensions are 1.18 ≤ γ′ ≤ 1.28; from this and the results for the b.c.c. lattice in particular, it is concluded that γ′ is probably 54 and hence the transition in χ0 is symmetrical about Tc (γ′ = γ). A repetition of this analysis for the third-equivalent-neighbor model three-dimensional lattices shows a marked shift in the estimated range of β and γ′; the results are 0.345 ≤ β ≤ 0.365, and 1.01 ≤ γ′ ≤ 1.14. In each of the above cases, the corresponding high-temperature (T > Tc) expansions of χ0(T) obtained previously have been analyzed to yield estimates of the critical exponent γ. The over-all results and in particular the estimates of γ for the s.q. and b.c.c. lattices suggest that this index is unaffected by extending the range of interaction, and that if γ is a rational fraction then it is the same fraction for the n.n. model and second- and third-equivalent neighbor models. Finally the high-temperature expansions of ZN in zero field, and of χ0(T) for the second-neighbor model are used to examine the dependence of the critical temperature Tc, the critical energy (E∞ − Ec)/kTc, and the critical entropy (S∞ − Sc)/k on the relative strengths of J1 and J2 for values of J2/J1 in the range 0 to 1. It is found that the variation of the critical point is well represented by Tc(α)=Tc(0)[1+m1α],where α = J2/J1 and lies in the range 0 ≤ α ≤ 1; and Tc(0) is the critical temperature of the nearest-neighbor model. The values of m1 are 0.61, 2.47, 0.84, 1.45, and 1.35 for the f.c.c., s.c., b.c.c., s.q., and Δ′r lattices, respectively. All these calculations are compared with the corresponding results for the Heisenberg model.
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.
The direct correlation function of the one-dimensional lattice-gas (Ising) model with nearest-neighbor and next-nearest-neighbor interactions is calculated exactly. It is shown that, depending on the strengths and signs of the coupling constants, the direct correlation function can either have a finite range equal to that of the interactions, or decay exponentially in a monotonic or oscillatory fashion. At the critical point, which occurs at zero temperature, the direct correlation function is found in all cases to have exactly the range of the interactions, while its values become unbounded; this is in contrast to earlier nonrigorous results. The second moment of the direct correlation function is found to diverge at the critical point, in the neighborhood of which it becomes proportional to the correlation length of the density (spin) fluctuations, confirming a prediction of scaling theory in one dimension.