Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.