Phase diagram and critical behavior of the square-lattice Ising model with competing nearest-neighbor and next-nearest-neighbor interactions

Effect of geometrical frustration on inverse freezing

The interplay between geometrical frustration (GF) and inverse freezing (IF) is studied within a cluster approach. The model considers first-neighbor (J_{1}) and second-neighbor (J_{2}) intracluster antiferromagnetic interactions between Ising spins on a checkerboard lattice and long-range disordered couplings (J) among clusters. We obtain phase diagrams of temperature versus J_{1}/J in two cases: the absence of J_{2} interaction and the isotropic limit J_{2}=J_{1}, where GF takes place. An IF reentrant transition from the spin-glass (SG) to paramagnetic (PM) phase is found for a certain range of J_{1}/J in both cases. The J_{1} interaction leads to a SG state with high entropy at the same time that can introduce a low-entropy PM phase. In addition, it is observed that the cluster size plays an important role. The GF increases the PM phase entropy, but larger clusters can give an entropic advantage for the SG phase that favors IF. Therefore, our results suggest that disordered systems with antiferromagnetic clusters can exhibit an IF transition even in the presence of GF.

Charge-regulation phase transition on surface lattices of titratable sites adjacent to electrolyte solutions: An analog of the Ising antiferromagnet in a magnetic field

We report a charge-patterning phase transition on two-dimensional square lattices of titratable sites, here regarded as protonation sites, placed in a low-dielectric medium just below the planar interface between this medium and a salt solution. We calculate the work-of-charging matrix of the lattice with use of a linear Debye-Hückel model, as input to a grand-canonical partition function for the distribution of occupancy patterns. For a large range of parameter values, this model exhibits an approximate inverse cubic power-law decrease of the voltage produced by an individual charge, as a function of its in-lattice separation from neighboring titratable sites. Thus, the charge coupling voltage biases the local probabilities of proton binding as a function of the occupancy of sites for many neighbors beyond the nearest ones. We find that even in the presence of these longer-range interactions, the site couplings give rise to a phase transition in which the site occupancies exhibit an alternating, checkerboard pattern that is an analog of antiferromagnetic ordering. The overall strength W of this canonical charge coupling voltage, per unit charge, is a function of the Debye length, the charge depth, the Bjerrum length, and the dielectric coefficients of the medium and the solvent. The alternating occupancy transition occurs above a curve of thermodynamic critical points in the (pH-pK,W) plane, the curve representing a charge-regulation analog of variation of the Néel temperature of an Ising antiferromagnet as a function of an applied, uniform magnetic field. The analog of a uniform magnetic field in the antiferromagnet problem is a combination of pH-pK and W, and 1/W is the analog of the temperature in the antiferromagnet problem. We use Monte Carlo simulations to study the occupancy patterns of the titratable sites, including interactions out to the 37th nearest-neighbor category (a distance of √74 lattice constants), first validating simulations through comparison with exact and approximate results for the nearest-neighbor case. We then use the simulations to map the charge-patterning phase boundary in the (pH-pK,W) plane. The physical parameters that determine W provide a framework for identifying and designing real surfaces that could exhibit charge-patterning phase transitions.

Nematic phase in the J(1)-J(2) square-lattice Ising model in an external field

The J(1)-J(2) Ising model in the square lattice in the presence of an external field is studied by two approaches: the cluster variation method (CVM) and Monte Carlo simulations. The use of the CVM in the square approximation leads to the presence of a new equilibrium phase, not previously reported for this model: an Ising-nematic phase, which shows orientational order but not positional order, between the known stripes and disordered phases. Suitable order parameters are defined, and the phase diagram of the model is obtained. Monte Carlo simulations are in qualitative agreement with the CVM results, giving support to the presence of the new Ising-nematic phase. Phase diagrams in the temperature-external field plane are obtained for selected values of the parameter κ=J(2)/|J(1)| which measures the relative strength of the competing interactions. From the CVM in the square approximation we obtain a line of second order transitions between the disordered and nematic phases, while the nematic-stripes phase transitions are found to be of first order. The Monte Carlo results suggest a line of second order nematic-disordered phase transitions in agreement with the CVM results. Regarding the stripes-nematic transitions, the present Monte Carlo results are not precise enough to reach definite conclusions about the nature of the transitions.

Critical properties of the Ziff, Gulari, and Barshad (ZGB) model with inert sites

We study the model proposed by Ziff, Gulari, and Barshad to mimic the oxidation of carbon monoxide (CO) in the presence of fixed impurities distributed over the catalytic surface. Our focus is on the continuous phase transition between the active phase, where occurs the production of carbon dioxide (CO2), and the inactive phase, where all the non inert sites become filled with oxygen molecules. We employ Monte Carlo simulations to calculate the different ratios between moments of the order parameter at the critical point, as well as, we determine the critical exponents β and ν⊥ as a function of the concentration of impurities. We show that the presence of impurities over the catalytic surface changes the critical behavior of the system. The critical exponents depend on the concentration of impurities and the model does not belong to the directed percolation universality class.

Dynamic approach to finite-temperature magnetic phase transitions in the extended J1-J2 model with vacancy order

The recently discovered iron-based superconductors A(y)Fe(2-x)Se(2) (A=K, Rb, Cs, Tl) show a long-range antiferromagnetic order with an unexpectedly high transition temperature T(N)~550 K and a unique √5×√5 vacancy order. Taking the extended J(1)-J(2) model as a minimal model, we investigate the finite-temperature magnetic phase transitions in a square lattice with a √5×√5 vacancy superstructure by using large-scale Monte Carlo simulations. By the parallel tempering technique, the block spin checkerboard and stripe antiferromagnetic states are detected to be the ground states for three representative sets of model parameters. The short-time dynamic approach is applied to accurately determine the critical temperature as well as the static and dynamic exponents. Our results indicate that the dramatic enhancement of the critical temperature as observed in experiments should be mainly due to a combination of the vacancy order and the block lattice contraction.

General approach for studying first-order phase transitions at low temperatures

By combining different ideas, a general and efficient protocol to deal with discontinuous phase transitions at low temperatures is proposed. For small T's, it is possible to derive a generic analytic expression for appropriate order parameters, whose coefficients are obtained from simple simulations. Once in such regimes simulations by standard algorithms are not reliable; an enhanced tempering method, the parallel tempering-accurate for small and intermediate system sizes with rather low computational cost-is used. Finally, from finite size analysis, one can obtain the thermodynamic limit. The procedure is illustrated for four distinct models, demonstrating its power, e.g., to locate coexistence lines and the phase density at the coexistence.

Scaling behavior of a square-lattice Ising model with competing interactions in a uniform field

Transfer-matrix methods, with the help of finite-size scaling and conformal invariance concepts, are used to investigate the critical behavior of two-dimensional square-lattice Ising spin-1/2 systems with first- and second-neighbor interactions, both antiferromagnetic, in a uniform external field. On the critical curve separating collinearly ordered and paramagnetic phases, our estimates of the conformal anomaly c are very close to unity, indicating the presence of continuously varying exponents. This is confirmed by direct calculations, which also lend support to a weak-universality picture; however, small but consistent deviations from the Ising-like values η=1/4, γ/ν=7/4, β/ν=1/8 are found. For higher fields, on the line separating row-shifted (2×2) and disordered phases, we find values of the exponent η very close to zero.

The effects of the next-nearest-neighbour density-density interaction in the atomic limit of the extended Hubbard model

We have studied the extended Hubbard model in the atomic limit. The Hamiltonian analysed consists of the effective on-site interaction U and the intersite density-density interactions W(ij) (both nearest-neighbour and next-nearest-neighbour). The model can be considered as a simple effective model of charge ordered insulators. The phase diagrams and thermodynamic properties of this system have been determined within the variational approach, which treats the on-site interaction term exactly and the intersite interactions within the mean-field approximation. Our investigation of the general case taking into account for the first time the effects of longer-ranged density-density interaction (repulsive and attractive) as well as possible phase separations shows that, depending on the values of the interaction parameters and the electron concentration, the system can exhibit not only several homogeneous charge ordered (CO) phases, but also various phase separated states (CO-CO and CO-nonordered). One finds that the model considered exhibits very interesting multicritical behaviours and features, including bicritical, tricritical, and critical end points and isolated critical points.

Ordering in spatial evolutionary games for pairwise collective strategy updates

Evolutionary 2×2 games are studied with players located on a square lattice. During the evolution the randomly chosen neighboring players try to maximize their collective income by adopting a random strategy pair with a probability dependent on the difference of their summed payoffs between the final and initial states assuming quenched strategies in their neighborhood. In the case of the anticoordination game this system behaves like an antiferromagnetic kinetic Ising model. Within a wide region of social dilemmas this dynamical rule supports the formation of similar spatial arrangement of the cooperators and defectors ensuring the optimum total payoff if the temptation to choose defection exceeds a threshold value dependent on the sucker's payoff. The comparison of the results with those achieved for pairwise imitation and myopic strategy updates has indicated the relevant advantage of pairwise collective strategy update in the maintenance of cooperation.

Square lattice gases with two- and three-body interactions revisited: a row-shifted (2x2) phase

Monte Carlo simulations have been used to study the phase diagrams for square Ising-lattice-gas models with two- and three-body interactions for values of interaction parameters in a range that has not been previously considered. We find unexpected qualitative differences as compared with predictions made on general grounds.