Surface tension: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
No edit summary
Line 23: Line 23:


=== Binder procedure===
=== Binder procedure===
Here, only a sketchy picture of the procedure is presented, more details can be found in Reference xxx


For given conditions of volume and temperature, the [[Helmholtz  energy function]] is computed as a function of the number of molecules:
For given conditions of volume and temperature, the [[Helmholtz  energy function]] is computed as a function of the number of molecules:
Line 34: Line 36:


Using basic thermodynamic procedures (Maxwell construction) it is possible
Using basic thermodynamic procedures (Maxwell construction) it is possible
to compute the densities of the two phases; <math> \rho_v, \rho_l </math>  
to compute the densities of the two phases; <math> \rho_v, \rho_l </math>.
 
 
Considering the thermodynamic limit for densities <math> \rho </math>  with  <math> \rho_v < \rho < \rho_l </math> the
[[Helmholtz energy function]] will be:
 
<math> A = - p_{eq} V + \mu_{eq} N + \gamma {\mathcal A}(N)  </math>.
 
where the quantities with the subindex "eq" are those corresponding to the fluid-phase equilbrium situation.


=== Explicit interfaces ===
=== Explicit interfaces ===

Revision as of 12:23, 1 August 2007

The surface tension, , is a measure of the work required to create a surface.

Thermodynamics

In the Canonical ensemble the surface tension is formally given as:

;

where

  • is the number of particles
  • is the volume
  • is the temperature
  • is the surface area
  • is the Helmholtz energy function

Computer Simulation

A review on different techniques to compute surface (interface) tension can be found in the paper by Gloor et al.

Liquid-Vapour Interfaces of one component systems

Binder procedure

Here, only a sketchy picture of the procedure is presented, more details can be found in Reference xxx

For given conditions of volume and temperature, the Helmholtz energy function is computed as a function of the number of molecules:

The calculation is usually carried out using Monte Carlo simulation

If liquid-vapour equilibrium occurs, the plot of the chemical potential, , as a function of shows a loop.

Using basic thermodynamic procedures (Maxwell construction) it is possible to compute the densities of the two phases; .


Considering the thermodynamic limit for densities with the Helmholtz energy function will be:

.

where the quantities with the subindex "eq" are those corresponding to the fluid-phase equilbrium situation.

Explicit interfaces

Mixtures

References

  1. G. J. Gloor, G. Jackson, F. J. Blas, and E. de Miguel "Test-area simulation method for the direct determination of the interfacial tension of systems with continuous or discontinuous potentials", Journal of Chemical Physics 123 134703 (2005)
  2. K. Binder "Monte Carlo calculation of the surface tension for two- and three-dimensional lattice-gas models", Physical Review A 25 pp. 1699 - 1709 (1982)