Spin Glass Phase Exists in the Random Weak Disorder for the Villain Model

In this work we have studied non random Villain model by introducing simple defects to calculate degeneracies of the first excited states using Pfaffian approach through a perturbation theory. The distributions of excitations of the ground states are displayed graphically. The results are indicated that spin glass occurs in the weak disorder for the Villain model. At the concentration of defect bonds 0.03 p = , the distribution behaves in the same manner as for 0.5 p = for different sizes of lattice. The latest result of the spin glass is presented in this paper.


Introduction
In this paper we have presented computational results for excitations of ground states for the Villain model [1] by introducing simple defects.In the Villain model, we have replaced ferromagnetic bonds by antiferromagnetic bonds randomly and calculated first excitations in the concentrations of defect bonds p = 0.01 to 0.5 for the whole random system.The initial model is fully frustrated [2] with highly degeneration of ground states and the replacement of ferromagnetic bonds with antiferromagnetic bonds reduces frustration [3].The spin glasses are disorder system [3] which is defined as the random combinations of ferromagnetic and antiferromagnetic bonds.
We have expected the Villain model to have a spin glass at zero temperature in the aforementioned limit of weak disorder.In this study, the square lattice L L × is wounded cylindrically and the periodic boundary condition is in one dimension.The energy gap is 4J in the mentioned system if L is even.Otherwise it is 2J [4] [5].In this work, we propose a simple picture of spin glass phase in the weak disorder that has evaluated from power law distribution of first excitations of ground states.
An outline of this paper is as follows.In Section 2 we introduce formalism and discuss with model and method how one can calculate first excitations.In Section 3 we present results from numerical simulation of the Villain model with antiferromagnetic exchange randomness and in Section 4 we draw concluding remarks.

Formalism
The frustration of two dimensional Ising model with interaction in one dimension is studied in this paper extensively.The two dimensional fully frustrated Ising model, or Villain model [1], consists of Ising spins on a square lattice with nearest neighbor bonds 1 ij J = with periodic boundary conditions.The Hamiltonian can be written as where 1 i σ = ± are Ising spins on a square lattice, L L × as shown in Figure 1.The method is used in the calculation basing on the Pfaffian [5] through a degenerate state perturbation theory [6]- [8].The degeneracies of the excitated states are calculated from an expression of the eigenvalues with disorder subspace.We write the partition function for the two dimensional Ising model as where ij J is the nearest-neighbour bond strength for sites ij and the product is over all bonds on the N-site lattice.The Matrix D [9]- [15] is Hermitian, pure imaginary, and of order 4N.The eigenvalue of Matrix D can be written ( ) where r is an integer and X is a real number.At T = 0 there is a degeneracy at 0 ε = equal to the number of frustrated plaquettes.The expression the ground-state degeneracy are written a where i M is the degeneracy of the ith excitated states.The entropy of the bimodal spin glass Ising model can be expressed in terms of the degeneracies as ( ) The specific heat of the system can be written as

Results and Discussion
We have carried out a numerical calculation to lift degeneracies of excited states at zero temperature.The disorder is confined to a frustrated square patch with periodic boundary conditions in one dimension.Square lattices of size L L × were considered over a range of concentrations of negative bond [ ] 0.01, 0.5 .p ∈ As p for the Villain model is zero, the concentration of negative bonds increases from p = 0.01 to 0.5.We have shown that there is a mechanism from the Villain phase to spin glass phase, which is characterized by the excitations of ground states.These excitations can exist for the disordered spin glass on the finite system.We have calculated In Figure 3 the distribution peak increases with the increase of lattice size and it's value scale is 3  L .But in Figures 4-7, the distribution peak goes higher for some L and after that it decreases with increasing value of L. This behavior of distributions tells us that spin glass phase occurs in random weak disorder.The aforementioned system is symmetrical within the range 0.      where it was 20,000.

Conclusion
We have carried out a numerical study of the Villain model, showing that there is an observation that indicates a zero temperature spin glass phase with weak disorder.We have used a computer code that is able to collect numerical data for simulation.The Pfaffian approach is used to collect data for first excitations of ground states, based on perturbation theory.It has been observed clearly that the spin glass phase occurs at concentration of defect bonds 0.03 0.5.p ≤ ≤

Figure 1 .
Figure 1.The Villain model on L L × lattice; along the horizontal direction, bonds are ferromagnetic while along vertical direction, bonds are ferromagnetic and antiferromagnetic alternately.

2 .
and several concentrations for Villain model.The distribution of excitations for p = 0.01 are the most likely value in the scale 2 ln L L as shown in Figure It has been shown that the distributions of degeneracy of the first excited state are fat-tailed for p = 0.03 as well as for p = 0.5 [4] not developing sharp peaks when L is increased.Spin glass behavior occurred in low concentrations (p = 0.03) as shown in Figure 4.

Figure 2 .
Figure 2. At concentration of defect bonds 0.01 p = , the distribution of ( ) 2 1 0 ln M M L L

Figure 3 .
Figure 3.At concentration of defect bonds 0.02 p = , the distribution of ( ) 3 1 0 M M L along

Figure 4 .
Figure 4.At concentration of defect bonds 0.03 p = , the distribution of ( ) 3 1 0 M M L along

Figure 5 .
Figure 5.At concentration of defect bonds 0.1 p = , the distribution of ( ) 3 1 0 M M L along

Figure 6 .
Figure 6.At concentration of defect bonds 0.2 p = , the distribution of ( ) 3 1 0 M M L along y-axis

Figure 7 .
Figure 7.At concentration of defect bonds 0.5 p = , the distribution of ( ) 3 1 0M M L along y-axis