Editing Wandering interface method

Jump to navigation Jump to search
Warning: You are not logged in. Your IP address will be publicly visible if you make any edits. If you log in or create an account, your edits will be attributed to your username, along with other benefits.

The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.

Latest revision Your text
Line 1: Line 1:
The '''wandering interface method''' (WIM) is a general, versatile and robust procedure for the calculation of surface [[Helmholtz energy function |free energies]]
The '''wandering interface method''' (WIM) is a general, versatile and robust procedure for the calculation of surface [[Helmholtz energy function |free energies]]. It can be employed to calculate liquid-vapor [[surface tension]] and solid-fluid interfacial tension alike. It does not require the explicit evaluation of the [[Virial equation of state | virial]] and may be employed for both continuous and discontinuous [[Idealised models |model potentials]].
<ref>[http://link.aps.org/doi/10.1103/PhysRevE.75.061609  L.G. MacDowell and P. Bryk "Direct Calculation of Interfacial Tensions from Computer Simulation: Results for Freely Jointed Tangent Hard Sphere Chains", Physical Review E '''75''' 061609 (2007)]</ref>.
It can be employed to calculate liquid-vapor [[surface tension]] and solid-fluid interfacial tension alike. It does not require the explicit evaluation of the [[Virial equation of state | virial]] and may be employed for both continuous and discontinuous [[Idealised models |model potentials]].


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({\mathcal 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({\mathcal A})</math>.
Line 9: Line 7:
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 Monte Carlo]] scheme. This provides a [[Markov chain |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:
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 Monte Carlo]] scheme. This provides a [[Markov chain |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({\mathcal A}) \propto -\gamma {\mathcal A} </math>
:<math> k_BT\ln P({\mathcal A}) \propto -\gamma {\mathcal A} </math>
Hence, the surface tension, <math>\gamma</math> may be readily evaluated as the slope of <math>P({\mathcal A})</math> on a logarithmic plot.
Hence, the surface tension, <math>\gamma</math> may be readily evaluated as the slope of <math>P({\mathcal A})</math> in a logarithmic plot.


Since the system will tend to stretch continuously so as to achieve a vanishing surface area, 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:
Since the system will tend to stretch continuously so as to achieve a vanishing surface area, 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:
Line 26: Line 24:


Several other choices are possible and have been tested since then.
Several other choices are possible and have been tested since then.
==References==
==References==
<references/>
#[http://link.aps.org/doi/10.1103/PhysRevE.75.061609  L.G. MacDowell and P. Bryk "Direct Calculation of Interfacial Tensions from Computer Simulation: Results for Freely Jointed Tangent Hard Sphere Chains", Physical Review E '''75''' 061609 (2007)]
[[category: Computer simulation techniques]]
[[category: Computer simulation techniques]]
Please note that all contributions to SklogWiki are considered to be released under the Creative Commons Attribution Non-Commercial Share Alike (see SklogWiki:Copyrights for details). If you do not want your writing to be edited mercilessly and redistributed at will, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource. Do not submit copyrighted work without permission!

To edit this page, please answer the question that appears below (more info):

Cancel Editing help (opens in new window)