Two-dimensional Ising model with Einstein site phonons

We consider a simple Ising magnetic model in two dimensions with Einstein site phonons and study it using Monte Carlo simulations that take into account both degrees of freedom simultaneously. In nonfrustrated systems, like the square lattice with ferromagnetic and antiferromagnetic interactions, we find that the coupling of the magnetic to the elastic degrees of freedom gradually lowers the magnetic ordering transition until it is completely suppressed at a critical value of the coupling constant. Above this the system suffers a simultaneous magnetic and structural transition into a dimerized state with lower crystalline symmetry and ferromagnetic clusters antiferromagnetically aligned. In the case of the Kagomé lattice with antiferromagnetic interactions, which is frustrated, we show that a similar ordered state takes place when the coupling constant is above a critical value.