Wandering interface method: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 3: Line 3:
The method is based on a clear and appealing physical basis: since an interface has usually an extra energy cost, a system will try to stretch in order to minimize the interfacial area. The original feature of WIM is to allow the system's interfacial area to fluctuate freely and extract the surface free energy from the related surface area probability distribution <math> P(A)</math>.  
The method is based on a clear and appealing physical basis: since an interface has usually an extra energy cost, a system will try to stretch in order to minimize the interfacial area. The original feature of WIM is to allow the system's interfacial area to fluctuate freely and extract the surface free energy from the related surface area probability distribution <math> P(A)</math>.  


During the course of the simulation, trial deformations of the simulation box at constant volume are attempted and either accepted or rejected according to the usual metropolis scheme. This provides a markovian random walk of the interfacial area. Collecting the surface area probability distribution, one can readily calculate the surface free energy, since, according to statistical mechanics, one can showv that:
This goal is achieved in practice by considering a canonical ensemble, with the usual constant N, V, T variables, where the box shape as described by the surface area,
:<math> k_BT\ln P(A) \propto \gamma A </math>
<math>A</math> is considered as an extra ensemble variable.
During the course of the simulation, trial deformations of the simulation box at constant volume are attempted and either accepted or rejected according to the usual metropolis scheme. This provides a markovian random walk of the interfacial area. Collecting the surface area probability distribution, one can readily calculate the surface free energy, since, according to statistical mechanics, one can show that:
:<math> k_BT\ln P(A) \propto -\gamma A </math>
Hence, the surface tension, <math>\gamma</math> may be readily evaluated as the slope of <math>P(A)</math> in a logarithmic plot.


:<math> \ln P(A) = \gamma /k_BT </math>
Since the system will tend to stretch continuously so as to achieve a vanishing surface area, in practice one must constraint the extent of deformations in order to bracket the surface area fluctuations within a convenient interval. This is easily achieved in practice by considering a modified canonical ensemble where the microstate probability distribution is given by:
:<math> P(U,A) \propto \exp(-\beta U + W(A) )</math>.
where <math>W</math> is a  biasing function chosen so as to constraint the fluctuations within some convenient interval.

Revision as of 22:07, 10 May 2009

The Wandering Interface Method (WIM) is a general, versatile and robust procedure for the calculation of surface free energies. It can be employed to calculate liquid-vapor surface tensions and solid-fluid interfacial tensions alike.

The method is based on a clear and appealing physical basis: since an interface has usually an extra energy cost, a system will try to stretch in order to minimize the interfacial area. The original feature of WIM is to allow the system's interfacial area to fluctuate freely and extract the surface free energy from the related surface area probability distribution .

This goal is achieved in practice by considering a canonical ensemble, with the usual constant N, V, T variables, where the box shape as described by the surface area, is considered as an extra ensemble variable. During the course of the simulation, trial deformations of the simulation box at constant volume are attempted and either accepted or rejected according to the usual metropolis scheme. This provides a markovian random walk of the interfacial area. Collecting the surface area probability distribution, one can readily calculate the surface free energy, since, according to statistical mechanics, one can show that:

Hence, the surface tension, may be readily evaluated as the slope of in a logarithmic plot.

Since the system will tend to stretch continuously so as to achieve a vanishing surface area, in practice one must constraint the extent of deformations in order to bracket the surface area fluctuations within a convenient interval. This is easily achieved in practice by considering a modified canonical ensemble where the microstate probability distribution is given by:

.

where is a biasing function chosen so as to constraint the fluctuations within some convenient interval.