Baxter solution to the O-Z equation near the critical point (Q919785)

The paper deals with the solution of the Ornstein-Zernike equation \((1)\quad h(r)=c(r)+\rho \int dr'h(| \vec r-\vec r|)c(r'),\) which is used in the theory of homogeneous classic liquids. To solve (1) some kind of equation is needed which relates h(r) and c(r) to the known interparticle pair potential, u(r), such as, e.g. (2) \(c(r)=(h(r)+1)[1- \exp (\beta u(r))].\) The authors comment the usual method of solving (1), (2) using an iterative procedure after transforming (1) with the Fourier transform, and find it inadequate in the neighborhood of the critical point. They report their results in solving numerically a system obtained from (1), (2) by Baxter's method and show that their method gives the most suitable results in the critical region.