Comparison of the capillary wave method and pressure tensor route for calculation of interfacial tension in molecular dynamics simulations.

We have studied the calculation of surface and interfacial tension for a variety of liquid-vapor and liquid-liquid interfaces using molecular dynamics (MD) simulations. Because of the inherently small scale of MD systems, large pressure fluctuations can cause imprecise calculations of surface tension using the pressure tensor route. The capillary wave method exhibited improved precision and stability throughout all of the simulated systems in this study. In order to implement this method, the interface was defined by fitting an error function to the density profile. However, full mapping of the interface from coordinate files produced enhanced accuracy. Upon increasing the system size, both methods exhibited higher precision, although the capillary wave method was still more reliable.