Albert Reiner, Gerhard Kahl:
Abstract: Continuing our investigation into the numerical properties of the Hierarchical Reference Theory, we study the square well fluid of range λ from slightly above unity up to 3.6. After briefly touching upon the core condition and the related decoupling assumption necessary for numerical calculations, we shed some light on the way an inappropriate choice of the boundary condition imposed at high density may adversely affect the numerical results; we also discuss the problem of the partial differential equation becoming stiff for close-to-critical and sub-critical temperatures. While agreement of the theory's predictions with simulational and purely theoretical studies of the square well system is generally satisfactory for λ greater than about 2, the combination of stiffness and the closure chosen is found to render the critical point numerically inaccessible in the current formulation of the theory for most of the systems with narrower wells. The mechanism responsible for some deficiencies is illuminated at least partially and allows us to conclude that the specific difficulties encountered for square wells are not likely to resurface for continuous potentials.
The published text is available from the journal's web site, or as a local copy in the following formats:
Alternatively, it is available as e-print cond-mat/0208471 from arXiv.org.
Note that the bzip2-compressed version may be sent with the incorrect type text/plain due to a misconfiguration of the server.
The article is Copyright 2002 American Institute of Physics. It may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics.
The article appeared in J Chem Phys and may be found at the address given above.
Sun Aug 25 16:48:52 MEST 2002