Lennard-Jones equation of state: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
m (Added an internal link)
(→‎Melting line: Added van der Hoef equations)
Line 23: Line 23:
On the following page is the [[FORTRAN code for the Kolafa and Nezbeda equation of state]].
On the following page is the [[FORTRAN code for the Kolafa and Nezbeda equation of state]].
==Melting line==
==Melting line==
The solid and liquid densities along the melting line are given by the equations of Mastny and  de Pablo (Ref <ref>[http://dx.doi.org/10.1063/1.2753149    Ethan A. Mastny and Juan J. de Pablo "Melting line of the Lennard-Jones system, infinite size, and full potential", Journal of Chemical Physics '''127''' 104504 (2007)]</ref> Eqs. 20 and 21):
The solid and liquid densities along the melting line are given by the following equations
====van der Hoef====
van der Hoef (Ref. <ref>[http://dx.doi.org/10.1063/1.1314342 Martin A. van der Hoef "Free energy of the Lennard-Jones solid", Journal of Chemical Physics '''113''' pp. 8142-8148 (2000)]</ref> Eqs. 25 and 26):
 
:<math>\rho_{\mathrm {solid}} = \beta^{-1/4} \left[ 0.92302 - 0.09218 \beta + 0.62381 \beta^2 -0.82672 \beta^3 + 0.49124  \beta^4 -0.10847 \beta^5\right]</math>
 
and
 
:<math>\rho_{\mathrm {liquid}} = \beta^{-1/4} \left[ 0.91070 - 0.25124 \beta + 0.85861 \beta^2 -1.08918 \beta^3 + 0.63932 \beta^4 -0.14433 \beta^5\right]</math>
 
====Mastny and  de Pablo====
Mastny and  de Pablo (Ref <ref>[http://dx.doi.org/10.1063/1.2753149    Ethan A. Mastny and Juan J. de Pablo "Melting line of the Lennard-Jones system, infinite size, and full potential", Journal of Chemical Physics '''127''' 104504 (2007)]</ref> Eqs. 20 and 21):


:<math>\rho_{\mathrm {solid}} = \beta^{-1/4} \left[ 0.908629 - 0.041510 \beta + 0.514632 \beta^2 -0.708590\beta^3 + 0.428351 \beta^4 -0.095229 \beta^5\right]</math>  
:<math>\rho_{\mathrm {solid}} = \beta^{-1/4} \left[ 0.908629 - 0.041510 \beta + 0.514632 \beta^2 -0.708590\beta^3 + 0.428351 \beta^4 -0.095229 \beta^5\right]</math>  
Line 30: Line 41:


:<math>\rho_{\mathrm {liquid}} = \beta^{-1/4} \left[ 0.90735 - 0.27120 \beta + 0.91784 \beta^2 -1.16270\beta^3 + 0.68012 \beta^4 -0.15284 \beta^5\right]</math>
:<math>\rho_{\mathrm {liquid}} = \beta^{-1/4} \left[ 0.90735 - 0.27120 \beta + 0.91784 \beta^2 -1.16270\beta^3 + 0.68012 \beta^4 -0.15284 \beta^5\right]</math>


==References==
==References==

Revision as of 18:51, 23 November 2010

The equation of state of the Lennard-Jones model.

Johnson, Zollweg and Gubbins equation of state

Johnson et al [1] proposed an equation of state based on 33 parameters within a modified Benedict, Webb and Rubin equation of state, which accurately reproduces the vapor liquid equilibrium curve.

Kolafa and Nezbeda equation of state

The Kolafa and Nezbeda equation of state [2] provides us with the Helmholtz energy function: (Eq. 30):

the compressibility factor (Eq. 31)

and the internal energy (Eq. 32)

On the following page is the FORTRAN code for the Kolafa and Nezbeda equation of state.

Melting line

The solid and liquid densities along the melting line are given by the following equations

van der Hoef

van der Hoef (Ref. [3] Eqs. 25 and 26):

and

Mastny and de Pablo

Mastny and de Pablo (Ref [4] Eqs. 20 and 21):

and

References

Related reading

This page contains numerical values and/or equations. If you intend to use ANY of the numbers or equations found in SklogWiki in any way, you MUST take them from the original published article or book, and cite the relevant source accordingly.