The Ising square lattice in aL×M geometry: A model for the effect of surface steps on phase transitions in adsorbed monolayers

Critical phenomena in adsorbed monolayers on surfaces are influenced by limited substrate homogeneity, such as surface steps. We consider the resulting finite-size and boundary effects in the framework of a lattice gas system with nearest neighbor attraction in aL×M geometry, with two free boundaries of lengthML, and periodic boundary conditions in the other direction (along the direction of the steps). This geometry thus models a terrace of the stepped surface, and adatoms adsorbed on neighboring terraces are assumed to be non-interacting. Also the effect of boundary fields is considered (describing the effects of missing neighbors and changed binding energy to the substrate near the boundary). Extensive Monte Carlo calculations on this model performed on a multi-transputer system are presented and analyzed in terms of phenomenological finite size scaling concepts. The fact that two scaling variables occur (/L,L/M, with being the correlation length in the bulk) is demonstrated explicitly. In the absence of boundary fields, the system forML orders nearT _c in a domain state, with domain walls running across the terrace, while at some temperature belowT _c a transition to a monodomain state occurs. This domain state slightly belowT _c is suppressed, however, by rather weak boundary fields. These results are interpreted in terms of exact theoretical predictions.     

Critical properties of XY model on two-layer Villain－ferromagnetic lattice

We investigate phase transitions of the XY model on a two-layer square lattice which consists of a Villain plane (J) and a ferromagnetic plane (I), using Monte Carlo simulations and a histogram method. Depending on the values of interaction parameters (I,J), the system presents three phases: namely, a Kosterlitz－Thouless (KT) phase in which the two planes are critical for I predominant over J, a chiral phase in which the two planes have a chiral order for J predominant over I and a new phase in which only the Villain plane has a chiral order and the ferromagnetic plane is paramagnetic with a small value of chirality. We clarify the nature of phase transitions by using a finite size scaling method. We find three different kinds of transitions according to the values of (I,J): the KT transition, the Ising transition and an XY-Ising transition with ν=0.849(3). It turns out that the Ising or XY-Ising transition is associated with the disappearance of the chiral order in the Villain plane.     

Spin-correlation function of the fully frustrated Ising model and ± J Ising spin glass on the square lattice

In this work we study the correlation function of the ground state of two-dimensional fully frustrated Ising model as well as spin glass. The Pfaffian method is used to calculate free energy and entropy as well as correlation function. We estimate the exponent of spin correlation function for fully frustrated model and spin glass. In this paper an overview of the latest results on the spin correlation function is presented.     

The affection on the nature of percolation by concentration of pentagon–heptagon defects in graphene lattice

In this paper the percolation behavior with a specific concentration of the defects was discussed on the twodimensional graphene lattice. The percolation threshold is determined by a numerical method with a high degree of accuracy. This method is also suitable for locating the percolation critical point on other crystalline structures. Through investigating the evolution of the largest cluster size and the cluster sizes distribution, we find that under various lattice sizes and concentrations of pentagon-heptagon defects there is no apparent change for the percolation properties in graphene lattice.     

Spin-correlation function of the fully frustrated Ising model and ±J Ising spin glass on a square lattice

In this work we study the correlation function of the ground state of a two-dimensional fully frustrated Ising model as well as spin glass. The Pfaffian method is used to calculate free energy and entropy as well as the correlation function. We estimate the exponent of spin correlation function for the fully frustrated model and spin glass. In this paper an overview of the latest results on the spin correlation function is presented.     

Solving the Ising model exactly on a 5×5×4 lattice using the Connection machine

We implement a recently proposed exact method for solving distrete statistical models for the 3-dimensional Ising model with open boundary conditions. Our computations were done on the Connection machine because the problem maps very naturally onto massively parallel architectures. We explicitly calculate the number of states of the system at each energy for systems of size 5×4×L _z, forL _z5. On serial or vector computers, the time for the computation scales with the volumeV likeV2 ^LLy . On the Connection Machine, the calculation can be spread across the processors. This decreases the computation requirements by a factor equal to the number of processors. We describe the method, its implementation on the Connection Machine both in PARIS and in FORTRAN, and our results. We also state the requirements for solving larger systems using this method.On leave from Physics Department, San Francisco State University, San Francisco, California, and Supercomputer Computations Research Institute, Florida State University, Tallahassee, Florida.     

Global phase diagram of a spin–orbit-coupled Kondo lattice model on the honeycomb lattice

Motivated by the growing interest in the novel quantum phases in materials with strong electron correlations and spin–orbit coupling, we study the interplay among the spin–orbit coupling, Kondo interaction, and magnetic frustration of a Kondo lattice model on a two-dimensional honeycomb lattice.We calculate the renormalized electronic structure and correlation functions at the saddle point based on a fermionic representation of the spin operators.We find a global phase diagram of the model at half-filling, which contains a variety of phases due to the competing interactions.In addition to a Kondo insulator, there is a topological insulator with valence bond solid correlations in the spin sector, and two antiferromagnetic phases.Due to the competition between the spin–orbit coupling and Kondo interaction, the direction of the magnetic moments in the antiferromagnetic phases can be either within or perpendicular to the lattice plane.The latter antiferromagnetic state is topologically nontrivial for moderate and strong spin–orbit couplings.     

Spin transport in the Néel and collinear antiferromagnetic phase of the two dimensional spatial and spin anisotropic Heisenberg model on a square lattice

We analyze and compare the effect of spatial and spin anisotropy on spin conductivity in a two dimensional S = 1/2 Heisenberg quantum magnet on a square lattice. We explore the model in both the Néel antiferromagnetic (AF) phase and the collinear antiferromagnetic (CAF) phase. We find that in contrast to the effects of spin anisotropy in the Heisenberg model, spatial anisotropy in the AF phase does not suppress the zero temperature regular part of the spin conductivity in the zero frequency limit–rather it enhances it. In the CAF phase (within the non-interacting approximation), the zero frequency spin conductivity has a finite value, which is suppressed as the spatial anisotropy parameter is increased. Furthermore, the CAF phase displays a spike in the spin conductivity not seen in the AF phase. We also explore the finite temperature effects on the Drude weight in the AF phase (within the collisionless approximation). We find that enhancing spatial anisotropy increases the Drude weight value and increasing spin anisotropy decreases the Drude weight value. Based on these studies, we conclude that antiferromagnets with spatial anisotropy are better spin conductors than those with spin anisotropy at both zero and finite temperatures.     

Exact ground states of the extended Hubbard model on the Kagomé lattice

We discuss the exact plaquette-ordered ground states of the generalized Hubbard model on the Kagomé lattice for several fillings, by constructing the Hamiltonian as a sum of products of projection operators for up and down spin sectors. The obtained exact ground states are interpreted as Néel ordered states on the bond-located electrons. We determine several parameter regions of the exact ground states, and calculate the entanglement entropy. We examine the above results by numerical calculations based on exact diagonalization and density-matrix renormalization group methods.     

Critical behavior of a stochastic anisotropic Bak–Sneppen model

In this paper we present our study on the critical behavior of a stochastic anisotropic Bak–Sneppen (saBS) model, in which a parameter α is introduced to describe the interaction strength among nearest species. We estimate the threshold fitness f _c and the critical exponent τ _r by numerically integrating a master equation for the distribution of avalanche spatial sizes. Other critical exponents are then evaluated from previously known scaling relations. The numerical results are in good agreement with the counterparts yielded by the Monte Carlo simulations. Our results indicate that all saBS models with nonzero interaction strength exhibit self-organized criticality, and fall into the same universality class, by sharing the universal critical exponents.     

Blume–Capel model on the antiferromagnetic fcc lattice

In the presence of a uniform external magnetic field, the spin-1 fcc antiferromagnetic Ising model in the four sub-lattice with the nearest-neighbor bilinear and crystal field interactions is investigated by using mean-field theory. The lattice is considered in blocks and each block contains 32 spins. While the central four spins forming a basic tetrahedron structure are allowed to fluctuate, the surrounding 28 spins take the mean-field values. The ground state phase diagram of the system is examined since it provides valuable information to make a prediction of the different phase diagram topologies. The phase diagrams are obtained on the reduced external field versus temperature, (H/|J|−kT/|J|)$(H/\left|J\right|-kT/\left|J\right|)$, plane for selected crystal field interaction values from the ground state phase diagram and it is found that the system exhibits reentrant behavior and some special points such as tricritical, bicritical and triple points.     

New Baxter phase in the Ashkin–Teller model on a cubic lattice

The mean field theory results are obtained from the Bogoliubov inequality for the spin-1/2 Ashkin–Teller model on a cubic lattice for different cluster sizes. The phase diagram, magnetization and free energy are obtained. From those expressions we observed a new phase in the model. Denoted in the course of this work by Baxter⁽²⁾ this new phase presents $〈S〉\ne 〈\sigma 〉\ne 0$. The phase transitions between the Baxter⁽²⁾ and the others well known phases for the model are studied and classified.     

Nonlocal field correlation functions on a lattice in the HP1 σ model

Connected two-point field-strength correlation functions are measured on a lattice in the quaternionic projective σ model within pure SU(2) theory. The correlation lengths extracted from exponential fits for these correlation functions, λ₁⁻¹ = 1.40(3) GeV and λ⁻¹ = 1.51(3) GeV, are found to be in good agreement with the results of other known calculations. The dependence of bilocal functions on the connector shape is also studied.     

Phase transitions of “extended” Ising models on a square lattice

The ferromagnetic square lattice Ising spin system is dynamically coupled to another set of Potts variables . We show that the usual Ising phase transition is universally preserved, but the transition temperatureT _c is shifted upwards. To investigate the transition and to calculateTM _c we use both a method by Müller-Hartmann-Zittartz as well as a systematic expansion about the soluble Ising limit. The comparison shows that the MHZ-formula forT _c is presumably a very accurate fit to the correct transition temperature. The results are relevant for special cases of more generalq-state models, for instance the Ashkin-Teller model.Work performed within the research program of the Sonderforschungsbereich 125 Aachen-Jülich-Köln     

Projected-wave-function study of the spin-1/2 Heisenberg model on the Kagomé lattice

We perform a Gutzwiller projected-wave-function study for the spin-1/2 Heisenberg model on the Kagomé lattice to compare energies of several spin-liquid states. The result indicates that a U(1)-Dirac spin-liquid state has the lowest energy. Furthermore, even without variational parameters, the energy turns out to be very close to that found by exact diagonalization. We show that such a U(1)-Dirac state represents a quantum phase whose low-energy physics is governed by four flavors of two-component Dirac fermions coupled to a U(1) gauge field. These results are discussed in the context of recent experiments on ZnCu(3)(OH)(6)Cl(2).     

Computer simulation study of the two-dimensional nematic lattice model based on the modified Gruhn–Hess pair potential model

A modified Gruhn–Hess pair potential model, where the potential parameters related to the elastic constants are different from the original model, was investigated with the aid of Monte Carlo (MC) simulation and Mean Field (MF) theory based upon a two-dimensional nematic lattice model. The model produces a ground state in perfect nematic order, where particles are aligned in the lattice plane. Both the MF predictions and the simulation results for the second-rank ordering tensor show that the system is biaxial in the low-temperature region, with a positive primary order parameter and the main director aligned along the lattice axis. A transition to uniaxial order takes place at higher temperatures with a negative primary order parameter and the director is orthogonal to the lattice plane. This orientational order survives up to temperatures higher than the transition temperature of the three-dimensional lattice model, possibly at all finite temperatures. MF predictions agree qualitatively with simulation but, in quantitative terms, the transition temperature is overestimated by 52%.     

Effect of the geometry of workpiece on polishing velocity in free annular polishing

Base on Coulomb friction model, the workpieces with different geometry rotating in free annular polishing are simulated. From simulation, the following conclusions are drawn. The angular velocity of workpiece is higher than that of polishing pad if the ring rotates uncontrolled in free annular polishing. The circular workpiece can synchronize with polishing pad through controlling the rotation of ring, which depends on the radii of ring and workpiece, the friction coefficients of polishing pad-workpiece and ring-workpiece, and the angular velocity of polishing pad. The workpiece with sharp corner cannot contact with the ring contiguously, which causes the contact state changing and the angular velocity of workpiece fluctuating ceaselessly, and this type of workpiece should be controlled with clamp to rotate synchronistically with the polishing pad.