Improved grand canonical sampling of vapour-liquid transitions

Nigel B. Wilding

Research output: Contribution to journalArticle

5 Citations (Scopus)
56 Downloads (Pure)

Abstract

Simulation within the grand canonical ensemble is the method of choice for
accurate studies of first order vapour-liquid phase transitions in model
fluids. Such simulations typically employ sampling that is biased with respect
to the overall number density in order to overcome the free energy barrier
associated with mixed phase states. However, at low temperature and for large
system size, this approach suffers a drastic slowing down in sampling
efficiency. The culprits are geometrically induced transitions (stemming from
the periodic boundary conditions) which involve changes in droplet shape from
sphere to cylinder and cylinder to slab. Since the overall number density
doesn't discriminate sufficiently between these shapes, it fails as an order
parameter for biasing through the transitions. Here we report two approaches to
ameliorating these difficulties. The first introduces a droplet shape based
order parameter that generates a transition path from vapour to slab states for
which spherical and cylindrical droplet are suppressed. The second simply
biases with respect to the number density in a tetragonal subvolume of the
system. Compared to the standard approach, both methods offer improved
sampling, allowing estimates of coexistence parameters and vapor-liquid surface
tension for larger system sizes and lower temperatures.
Original languageEnglish
Article number414016
Number of pages7
JournalJournal of Physics: Condensed Matter
Volume28
Issue number41
Publication statusPublished - 22 Aug 2016

Fingerprint

sampling
Vapors
vapors
Sampling
Liquids
liquids
slabs
Energy barriers
Free energy
Surface tension
liquid surfaces
Phase transitions
Boundary conditions
Temperature
interfacial tension
liquid phases
Fluids
simulation
free energy
boundary conditions

Cite this

Improved grand canonical sampling of vapour-liquid transitions. / Wilding, Nigel B.

In: Journal of Physics: Condensed Matter , Vol. 28, No. 41, 414016, 22.08.2016.

Research output: Contribution to journalArticle

@article{36d3e57dc91d46bba46cf4dd672a9eb4,
title = "Improved grand canonical sampling of vapour-liquid transitions",
abstract = "Simulation within the grand canonical ensemble is the method of choice foraccurate studies of first order vapour-liquid phase transitions in modelfluids. Such simulations typically employ sampling that is biased with respectto the overall number density in order to overcome the free energy barrierassociated with mixed phase states. However, at low temperature and for largesystem size, this approach suffers a drastic slowing down in samplingefficiency. The culprits are geometrically induced transitions (stemming fromthe periodic boundary conditions) which involve changes in droplet shape fromsphere to cylinder and cylinder to slab. Since the overall number densitydoesn't discriminate sufficiently between these shapes, it fails as an orderparameter for biasing through the transitions. Here we report two approaches toameliorating these difficulties. The first introduces a droplet shape basedorder parameter that generates a transition path from vapour to slab states forwhich spherical and cylindrical droplet are suppressed. The second simplybiases with respect to the number density in a tetragonal subvolume of thesystem. Compared to the standard approach, both methods offer improvedsampling, allowing estimates of coexistence parameters and vapor-liquid surfacetension for larger system sizes and lower temperatures.",
author = "Wilding, {Nigel B.}",
year = "2016",
month = "8",
day = "22",
language = "English",
volume = "28",
journal = "Journal of Physics: Condensed Matter",
issn = "0953-8984",
publisher = "IOP Publishing",
number = "41",

}

TY - JOUR

T1 - Improved grand canonical sampling of vapour-liquid transitions

AU - Wilding, Nigel B.

PY - 2016/8/22

Y1 - 2016/8/22

N2 - Simulation within the grand canonical ensemble is the method of choice foraccurate studies of first order vapour-liquid phase transitions in modelfluids. Such simulations typically employ sampling that is biased with respectto the overall number density in order to overcome the free energy barrierassociated with mixed phase states. However, at low temperature and for largesystem size, this approach suffers a drastic slowing down in samplingefficiency. The culprits are geometrically induced transitions (stemming fromthe periodic boundary conditions) which involve changes in droplet shape fromsphere to cylinder and cylinder to slab. Since the overall number densitydoesn't discriminate sufficiently between these shapes, it fails as an orderparameter for biasing through the transitions. Here we report two approaches toameliorating these difficulties. The first introduces a droplet shape basedorder parameter that generates a transition path from vapour to slab states forwhich spherical and cylindrical droplet are suppressed. The second simplybiases with respect to the number density in a tetragonal subvolume of thesystem. Compared to the standard approach, both methods offer improvedsampling, allowing estimates of coexistence parameters and vapor-liquid surfacetension for larger system sizes and lower temperatures.

AB - Simulation within the grand canonical ensemble is the method of choice foraccurate studies of first order vapour-liquid phase transitions in modelfluids. Such simulations typically employ sampling that is biased with respectto the overall number density in order to overcome the free energy barrierassociated with mixed phase states. However, at low temperature and for largesystem size, this approach suffers a drastic slowing down in samplingefficiency. The culprits are geometrically induced transitions (stemming fromthe periodic boundary conditions) which involve changes in droplet shape fromsphere to cylinder and cylinder to slab. Since the overall number densitydoesn't discriminate sufficiently between these shapes, it fails as an orderparameter for biasing through the transitions. Here we report two approaches toameliorating these difficulties. The first introduces a droplet shape basedorder parameter that generates a transition path from vapour to slab states forwhich spherical and cylindrical droplet are suppressed. The second simplybiases with respect to the number density in a tetragonal subvolume of thesystem. Compared to the standard approach, both methods offer improvedsampling, allowing estimates of coexistence parameters and vapor-liquid surfacetension for larger system sizes and lower temperatures.

UR - http://iopscience.iop.org/article/10.1088/0953-8984/28/41/414016

M3 - Article

VL - 28

JO - Journal of Physics: Condensed Matter

JF - Journal of Physics: Condensed Matter

SN - 0953-8984

IS - 41

M1 - 414016

ER -