LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Interfacial tension (IFT)
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Interfacial tension (IFT)


From: "Ahmed E. Ismail" <aeismail@...24...>
Date: Sat, 16 Dec 2017 01:34:45 +0000

Liquid-vapor simulations of surface tension are easy to set up: https://doi.org/10.1063/1.2209240

The basic idea is that you set up a slab of the liquid in a simulation box with the direction normal to the surfaces elongated, from which you can compute the surface tension from the difference between the pressure tensor components normal and parallel to the interface. (The formulas are provided.) Note that if long-range dispersion (PPPM/disp) is used, then there is no need to include the long-range corrections that must otherwise be calculated from the density profiles (https://doi.org/10.1080/00268979500101371).

To do liquid-liquid surface tension calculations, you can measure each of the individual liquid-gas surface tensions, and subtract those from the value obtained for the combined liquid-liquid system. (The method is explained in more detail: https::/doi.org/10.1021/ma802805y.

—AEI



On Fri, Dec 15, 2017 at 1:45 PM Giacomo Fiorin <giacomo.fiorin@...24...> wrote:
Sorry, I mistyped: it's fix print, not fix props.

On Fri, Dec 15, 2017 at 1:41 PM, Giacomo Fiorin <giacomo.fiorin@...33....24...> wrote:
Shame not to catch the chance for self-promotion:
https://doi.org/10.1063/1.4937153
(look in "parameterization")

The formula is well known, and dates back decades based on Kirkwood's work (assumes integration an orthogonal unit cell going across the interface).  The procedure for air/liquid or liquid/liquid is basically the same: use a unit cell with constant volume, and all liquids must be in equilibrium for it.  Because fluctuations of the pressure tensor are big, expect slow converge and some sensitivity to cutoffs.  Isotropic fluids (aside from the interface) will converge much better than oriented fluids like lipid membranes.

Specifically in LAMMPS, you can define the formula for gamma as an equal-style variable, and print it through a custom thermo output or fix props.

Giacomo

On Fri, Dec 15, 2017 at 1:03 PM, Robert McAllen <mm229196@...7278...> wrote:
Dear users

Does anybody know how calculating Interfacial tension (IFT) between two immiscible liquids or liquid/gas in LAMMPS?

If yes, would you please let me know an illustrative work or paper?

Sincerely



Sent with Mailtrack

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users




--
Giacomo Fiorin
Associate Professor of Research, Temple University, Philadelphia, PA
Contractor, National Institutes of Health, Bethesda, MD



--
Giacomo Fiorin
Associate Professor of Research, Temple University, Philadelphia, PA
Contractor, National Institutes of Health, Bethesda, MD
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot_______________________________________________
lammps-users mailing list
lammps-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users
--
========================
Ahmed E. Ismail
aei@...495...
aei@...496...