Although the solid-solid transition coexistence curves can be obtained
from simulations, we have yet to demonstrate that this transition
involves phases that are thermodynamically stable.
In particular, the melting transition transition might preempt the
solid-solid phase separation. It is therefore essential to study the
fluid-solid transition as well.
We computed the solid-fluid coexistence curve by means of
thermodynamic integration. The Helmholtz free energy of square-well
solid was calculated according to eqn.
3.2 using our
simulation data. The free energy of the fluid phase was obtained by
combining data from our simulations of a square-well fluid with the
known free energy of the hard-sphere reference fluid [37].
In the two dimensional case the fluid simulations were performed on a
system consisting of 200 disks with square well potential in a
square box. The densities ranged from =0.81, which is below
the hard-disk fluid-solid transition at
=0.87 [56],
to
=0.95, well above the melting density.
The other simulation parameters were the same as in the two
dimensional square-well solid simulation.
To simulate the fluid in three dimensions, we used a system of 108
square-well spheres in a periodic cubic box.
The density was varied from =0.9 to
=1.0 and the well
width ranged from
=0.01 to
=0.06. The simulation
parameters were equal to those chosen for the solid simulations.
For both the 2D and 3D simulations the initial random configuration
was compressed to the required density and equilibrated for 20,000
cycles before data was collected in a production run of 20,000 cycles.
To calculate the fluid-solid coexistence one needs the absolute free
energy of both the reference fluid and the reference solid phase. For
the free energy of the hard disk fluid we used a Padé
approximation proposed by Hoover and Ree [56].
The hard disk solid free energy was obtained from simulation data by
Alder et al. [55].
The free energy of the hard-sphere fluid was calculated using the
accurate Carnahan-Starling equation of state [37], whereas
the Hall equation of state [54] was used in the solid region
together with an absolute free energy value obtained from simulations
of Frenkel and Ladd [39].
Using these reference free energies and the simulated average internal
energies in eqn.
3.2 we were able to obtain the
coexistence curves for the fluid-solid transition for both the two
and three dimensional square-well models.
Figure: | Simulated ![]() ![]() ![]() ![]() ![]() |