Surface tension: Difference between revisions
Carl McBride (talk | contribs) 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 11: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
- 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)
- 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)