Molecular dynamics simulation of surface tension
The main goal of this work is to give a systematic estimation of surface tension for planar and spherical cases and related characteristics using the molecular dynamics (MD). For planar case, a thin film is considered. For spherical case, a droplet and a bubble are considered. Distributions of density and pressure tensor components are obtained through an averaging process during MD trajectory. The surface tension for planar case is obtained directly from pressure tensor. Base on the mechanical arguments for the force and moment balances and the mechanical stability called as “mechanical route” bubble and droplet surface tensions are calculated together with the corresponding radii. Also dependence of surface tension on cut-off radius, temperature and size of the system is investigated.