Loading…
Revisiting the Frenkel-Ladd method to compute the free energy of solids: The Einstein molecule approach
In this paper a new method to evaluate the free energy of solids is proposed. The method can be regarded as a variant of the method proposed by Frenkel and Ladd [ J. Chem. Phys. 81 , 3188 ( 1984 )] . The main equations of the method can be derived in a simple way. The method can be easily implemente...
Saved in:
| Published in: | The Journal of chemical physics 2007-10, Vol.127 (15), p.154113-154113-12 |
|---|---|
| Main Authors: | , |
| Format: | Article |
| Language: | English |
| Citations: | Items that this one cites Items that cite this one |
| Online Access: | Get full text |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Summary: | In this paper a new method to evaluate the free energy of solids is proposed. The method can be regarded as a variant of the method proposed by
Frenkel
and
Ladd
[
J. Chem. Phys.
81
,
3188
(
1984
)]
. The main equations of the method can be derived in a simple way. The method can be easily implemented within a Monte Carlo program. We have applied the method to determine the free energy of hard spheres in the solid phase for several system sizes. The obtained free energies agree within the numerical uncertainty with those obtained by
Polson
[
J. Chem. Phys.
112
,
5339
(
2000
)]
. The fluid-solid equilibria has been determined for several system sizes and compared to the values published previously by
Wilding
and
Bruce
[
Phys. Rev. Lett.
85
,
5138
(
2000
)]
using the phase switch methodology. It is shown that both the free energies and the coexistence pressures present a strong size dependence and that the results obtained from free energy calculations agree with those obtained using the phase switch method, which constitutes a cross-check of both methodologies. From the results of this work we estimate the coexistence pressure of the fluid-solid transition of hard spheres in the thermodynamic limit to be
p
*
=
11.54
(
4
)
, which is slightly lower than the classical value of Hoover and Ree
(
p
*
=
11.70
)
[
J. Chem. Phys.
49
,
3609
(
1968
)
]. Taking into account the strong size dependence of the free energy of the solid phase, we propose to introduce finite size corrections, which allow us to estimate approximately the free energy of the solid phase in the thermodynamic limit from the known value of the free energy of the solid phase with
N
molecules. We have also determined the free energy of a Lennard-Jones solid by using both the methodology of this work and the finite size correction. It is shown how a relatively good estimate of the free energy of the system in the thermodynamic limit is obtained even from the free energy of a relatively small system. |
|---|---|
| ISSN: | 0021-9606 1089-7690 |
| DOI: | 10.1063/1.2790426 |