An analytical expression for the boundary free energy of the Ising square lattice with nearest- and next-nearest-neighbor interactions is derived. The Ising square lattice with anisotropic nearest-neighbor ( J x and J y ) and isotropic next-nearest-neighbor ( J d ) interactions has an order-disorder phase transition at a temperature T = T c given by the condition e −2 J x / k b T c + e − 2 J y / k b T c + e −2( J x + J y )/ k b T c (2 − e −4 J d / k b T c ) = e 4 J d / k b T c . The critical line that separates the ordered (ferromagnetic) phase from the disordered (paramagnetic) phase is in excellent agreement with series expansion, finite scaling of transfer matrix and Monte Carlo results. For a vanishing next-nearest-neighbor interaction Onsager's famous result, i.e. sinh (2 J x / k b T c )sinh (2 J y / k b T c ) = 1, is recaptured.
First, a new method proposed by Baker and Rushbrooke is used to study the simple ferromagnetic Ising model at and below the Curie temperature. Of course, the properties of the Ising model are already well known, so that the main aim here is to assess the potential and reliability of the new method, since it has wide applicability to other models which have not been otherwise studied. Between 8 and 16 coefficients of exact high-temperature expansions for fixed values of the magnetization are derived for various two- and three-dimensional lattices. A Pad\'e-approximant analysis of these expansions at the critical isotherm and magnetic phase boundary enables us to estimate the critical exponents $\ensuremath{\beta}$, ${\ensuremath{\gamma}}^{\ensuremath{'}}$, and $\ensuremath{\delta}$, and plot the spontaneous magnetization. The results are in good agreement with previous calculations. Secondly, an analysis of the exact series expansions provides no support for the conjecture that the phase boundary is a line of essential singularities. However, the same expansions strongly suggest the existence of a "spinodal" curve, whose properties are in reasonable agreement with the predictions of various heuristic arguments (based essentially upon analyticity at the phase boundary and one-phase homogeneity in the critical region). Finally, structure and a mild extension of the proven analyticity of the free energy are used to show the $\ensuremath{\Delta}\ensuremath{\le}{\ensuremath{\Delta}}^{\ensuremath{'}}$, $\ensuremath{\gamma}\ensuremath{\le}{\ensuremath{\gamma}}^{\ensuremath{'}}$.
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.
Using the parallel tempering algorithm and graphics processing unit accelerated techniques, we have performed large-scale Monte Carlo simulations of the Ising model on a square lattice with antiferromagnetic (repulsive) nearest-neighbor and next-nearest-neighbor 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, re-entrance 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, Suzuki's weak universality seems to hold.
Monte Carlo simulations are used to study the location and nature of phase boundaries for Ising square lattices with antiferromagnetic coupling
JNN between nearest neighbors and additional interactions
JNNN between next-nearest neighbors and
J3NN between third-nearest neighbors. Results in zero magnetic field are obtained for a wide range of R=
JNNN/JNN and R'=
J3NN/JNN. In addition to the c(2\ifmmode\times\else\texttimes\fi{}2) and (2\ifmmode\times\else\texttimes\fi{}1) phases, which also occur for R'=0, we find new (4\ifmmode\times\else\texttimes\fi{}4) and (4\ifmmode\times\else\texttimes\fi{}2) ordered states, for R'\ensuremath{\ne}0, which are separated from the disordered state by lines of first-order transitions. The nonuniversal critical behavior of the (2\ifmmode\times\else\texttimes\fi{}1) phase is studied using the block-distribution method and finite-size scaling. The possible existence of incommensurate phases is also explored.
The critical behaviour of the Ising model on the square lattice with ferromagnetic nearest-neighbour interactions (J) and anti-ferromagnetic next-nearest neighbour interactions (J') is discussed for small lambda =J/ mod J' mod . The singular part of the free energy is calculated to second order in lambda by perturbation theory. Near the critical temperature of the unperturbed system (J=0), this expansion is found to have a form which may be exponentiated yielding, for non-zero lambda , a non-universal critical line along which the exponents vary continuously with lambda . For small lambda , the specific heat exponent, alpha = alpha 1 lambda 2 +O( lambda 2 ) with alpha 1 approximately=1.5.