Pushing the Limits of EPD Zeros Method

The use of partition function zeros in the study of phase transitions is growing in the last decade mainly due to improved numerical methods as well as novel formulations and analysis. In this paper, the impact of different parameters choice for the energy probability distribution (EPD) zeros that were recently introduced by Costa et al. is explored in search for optimal values. Our results indicate that the EPD method is very robust against parameter variations and only small deviations on estimated critical temperatures are found even for large variation of parameters, allowing one to obtain accurate results with low computational cost. A proposal to circumvent potential convergence issues of the original algorithm is presented and validated for the case where it occurs.