Solution of the Ornstein–Zernike Equation in the Critical Region

A new numerical scheme for the solution of liquid state iIntegral equations using the Baxter factorization of the Ornstein–Zernike equation is proposed. For short range potentials the method yields reliable results over the whole fluid region, including the vicinity of the critical point, and opens up new possibilities for numerical study of the critical behavior of iIntegral equation approximations. To demonstrate the effectiveness of the method, numerical results are compared with the analytical solution of the mean spherical approximation for a hard-core plus Yukawa tail iInteraction potential. Accurate results for the critical exponents δ, γ, and η for this model are obtained