Monte Carlo study of the Yukawa coupled two-spin Ising model

We consider a particular four state spin system composed of two Ising spins ( s_x, σ_x) with independent hopping parameters κ ₁ κ ₂, coupled by a bilinear Yukawa term, ys_xσ_x. The Yukawa term is solely responsible for breaking the global Z₂ × Z₂ symmetry down to Z2. This model is intended as an illustration of general coupled Higgs system where scalars can arise both as composite and elementary excitations. For the Ising example in 2d, we give convincing numerical evidence of the universality of the two-spin system with the one-spin Ising model, by Monte Carlo simulations and finite size scaling analysis. We also show that as we approach the phase transition, universality arises by having a single correlation length that diverges.