The modified model doesn't include enough terms to demonstrate fractionalized phase. The full model is simulated in ising_gauge_field.py