My reference have been this paper. But they use a different method to identify the diffusion front (damage spreading), so I am not 100% sure I have done everything correctly. I used different sizes and temperature ranges to extrapolate the effective critical temperature due to finite size...
I am currently working on the 3D potts model in a temperature gradient to identify the critical temperature and critical exponent. I use the Swendsen-Wang algorithm to simulate the dynamics of the system, and I use the Hoshen-Kopelman algorithm to identify the clusters of spins.
The problem...