Critical behavior of the restricted primitive model revisited

Reassessment of the critical temperature and density of the restricted primitive model of an ionic fluid by Monte Carlo simulations performed for system sizes with linear dimension up to $\text{L/σ=34}$ and sampling of ${\text{∼10}}^9$ trial moves leads to $T_c^{*}\text{=0.049\hspace{0.05em}17±0.000\hspace{0.05em}02}$ and ${\rho }_c^{*}\text{=0.080±0.005.}$ Finite size scaling analysis based in the Bruce–Wilding procedure gives critical exponents in agreement with those of the three-dimensional Ising universality class. An analysis similar to that proposed by Orkoulas et al. [Phys. Rev. E 63, 051507 (2001)], not relying on an a priori knowledge of the universality class, leads to an inaccurate estimate of $T_c^{*}$ and to unexpected behavior of the specific heat and value of the critical exponent ratio γ/ν.