Pair potentials that are bounded at the origin provide an accurate description of the effective interaction for many systems of dissolved soft macromolecules (e.g., flexible dendrimers). Using numerical free-energy calculations, we reconstruct the eq
uilibrium phase diagram of a system of particles interacting through a potential that brings together a Gaussian repulsion with a much weaker Gaussian attraction, close to the thermodynamic stability threshold. Compared to the purely-repulsive model, only the reentrant branch of the melting line survives, since for lower densities solidification is overridden by liquid-vapor separation. As a result, the phase diagram of the system recalls that of water up to moderate (i.e., a few tens MPa) pressures. Upon superimposing a suitable hard core on the double-Gaussian potential, a further transition to a more compact solid phase is induced at high pressure, which might be regarded as the analog of the ice I to ice III transition in water.
Isotropic pair potentials that are bounded at the origin have been proposed from time to time as models of the effective interaction between macromolecules of interest in the chemical physics of soft matter. We present a thorough study of the phase b
ehavior of point particles interacting through a potential which combines a bounded short-range repulsion with a much weaker attraction at moderate distances, both of Gaussian shape. Notwithstanding the fact that the attraction acts as a small perturbation of the Gaussian-core model potential, the phase diagram of the double-Gaussian model (DGM) is far richer, showing two fluid phases and four distinct solid phases in the case that we have studied. Using free-energy calculations, the various regions of confluence of three distinct phases in the DGM system have all been characterized in detail. Moreover, two distinct lines of reentrant melting are found, and for each of them a rationale is provided in terms of the elastic properties of the solid phases.
