Critical behavior of the two-dimensional thermalized bond Ising model

A suitably modified Wolff single-cluster Monte Carlo simulation has been performed to investigate the critical behavior of a two-dimensional Ising model with temperature-dependent annealed bond dilution, also known as the thermalized bond Ising model, which is intended to simulate the thermal excitations of electronic bond degrees of freedom as in covalently bonded network liquids. A finite-size scaling analysis of the susceptibility and the fourth-order cumulant, results in a reliable estimation of the critical exponents in the thermodynamic limit. The exponents are found to be consistent with those predicted by the Fisher renormalization relations, despite the well known violations of the renormalization relations when approximate methods such as real space renormalization group are employed to investigate two-dimensional Ising model with annealed bond dilution, and the temperature variation of the bond concentration in thermalized bond model system.