The two-state Ising model with nearest-neighbour bilinear exchange constant J (2) , three-spin coupling constant J (3) and four-spin coupling constant J (4) is studied using the effective-field approximation. The temperature variation in spontaneous magnetization is studied for different coordination numbers z and for various values of the parameters alpha =J (3) /J (2) and alpha '=J (4) /J (2) . It is seen that, in the absence of three-spin coupling, there exists a critical value of alpha ' (say, alpha ' c ) above which the magnetization becomes double valued at and beyond the Curie temperature T c , forming a bulge near T c which indicates the occurrence of a first-order transition. The expressions for alpha ' c for z=3, 4, 6 are derived and it is found that alpha ' c increases as z increases. In the absence of four-spin coupling it is found that first-order transitions appear for alpha >or=0. The magnetization curves are also studied for the case when all the interactions are present. The results are discussed with reference to those obtained from the molecular-field approximation.
The two-state Ising model with nearest-neighbour bilinear exchange constant J (2) , three-spin coupling constant J (3) and four-spin coupling constant J (4) is studied using the effective-field approximation. The temperature variation in spontaneous magnetization is studied for different coordination numbers z and for various values of the parameters alpha =J (3) /J (2) and alpha '=J (4) /J (2) . It is seen that, in the absence of three-spin coupling, there exists a critical value of alpha ' (say, alpha ' c ) above which the magnetization becomes double valued at and beyond the Curie temperature T c , forming a bulge near T c which indicates the occurrence of a first-order transition. The expressions for alpha ' c for z=3, 4, 6 are derived and it is found that alpha ' c increases as z increases. In the absence of four-spin coupling it is found that first-order transitions appear for alpha >or=0. The magnetization curves are also studied for the case when all the interactions are present. The results are discussed with reference to those obtained from the molecular-field approximation.
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 spin-1 Ising model with bilinear and biquadratic interactions is investigated using a generalized constant-coupling approximation with two parameters. Equations to determine these parameters are derived by minimizing the variational free energy. The averages of the dipolar and quadrupolar moments expressed in terms of these parameters are first ascertained to give exact results for a linear chain, then investigated in detail for the quadrupolar-coupled system in comparison with the molecular-field results. The pair-correlation and the heat capacity of this system are also investigated in connection with the phase transition. For the system with the bilinear and biquadratic interactions the effect of the anisotropy constant to the phase transition is discussed.
The equilibrium properties of the spin-1 Ising system Hamiltonian with arbitrary bilinear (J), biquadratic (K) and odd (L), which is also called dipolar-quadrupolar, interactions is studied for zero magnetic field in the lowest approximation of the cluster variation method. The odd interaction is combined with the bilinear (dipolar) and biquadratic (quadrupolar) exchange interactions by the geometric mean. In this system, phase transitions depend on the ratio of the coupling parameters, alpha = J/K; therefore, the dependence of the nature of the phase transition on alpha is investigated extensively and it is found that for alpha less than or equal to 1 and alpha greater than or equal to 2000 a second-order phase transition occurs, and for 1 < alpha < 2000 a first-order phase transition occurs. The critical temperatures in the case of a second-order phase transition and the upper and lower limits of stability temperature in the case of a first-order phase transition are obtained for different values of alpha calculated using the Hessian determinant. The first-order phase transition temperatures are found by using the free energy values while increasing and decreasing the temperature. Besides the stable branches of the order parameters, we establish also the metastable and unstable parts of these curves and the thermal variations of these solutions as a function of the reduced temperature are investigated. The unstable solutions for the first-order phase transitions are obtained by displaying the free energy surfaces in the form of a contour map. Results are compared with the spin-1 Ising system Hamiltonian with the bilinear and biquadratic interactions and it is found that the odd interaction greatly influences the phase transitions.
Phase transitions of the spin-one Ising model with nearest-neighbour and next-nearest-neighbour interactions are investigated for a square lattice using a molecular-field theory and a Monte Carlo method. Special attention is paid to the effects of the biquadratic next-nearest-neighbour interaction on the phase transitions of the model. Phase diagrams characterizing types of phase transition are obtained for several values of the coupling constants. The existence of the multi-critical points and the successive phase transitions are also discussed in connection with these phase diagrams.
Dynamics of a spin-3/2 Ising system Hamiltonian with bilinear and biquadratic nearest-neighbour exchange interactions is studied by a simple method in which the statistical equilibrium theory is combined with the Onsager's theory of irreversible thermodynamics. First, the equilibrium behaviour of the model in the molecular-field approximation is given briefly in order to obtain the phase transition temperatures, i.e. the first- and second-order and the tricritical points. Then, the Onsager theory is applied to the model and the kinetic or rate equations are obtained. By solving these equations three relaxation times are calculated and their behaviours are examined for temperatures near the phase transition points. Moreover, the z dynamic critical exponent is calculated and compared with the z values obtained for different systems experimentally and theoretically, and they are found to be in good agrement.
The one-dimensional spin-½ Ising model with very long-range ferromagnetic interaction and first and second neighbor antiferromagnetic interactions is solved exactly in a direct magnetic field. A study of the finite temperature behavior confirms the existence of up to four first-order phase transitions ending in classical critical points. We find a confluence of three critical lines at a ``tricritical'' point with exponents β=14 , δ=5, γ=γ′=1,α=0, and α′=12. We find also a confluence of four critical lines at a ``tetracritical'' point with exponents β=14, δ=7, γ=γ′=1, α=0, and α′=23. We show the existence of triple lines and ``quadruple'' lines, along which three or four phases are simultaneously in contact with each other, respectively. For some special values of the interaction strengths, five phases are in contact with each other at a ``quintuple'' point. We remark that our ``multiple'' lines and ``multicritical'' points are not ``true'' multiple lines and multicritical points but that they would be so in a somewhat extended model. This model satisfies a generalized form of the ordinary phase rule of Gibbs. In distinction to various models considered recently, our model does not exhibit a singular phase-boundary diameter.