Non-equilibrium hybrid insertion/extraction through the 4th dimension in grand-canonical simulation

Luc Belloni 1, *
* Corresponding author
1 LIONS - Laboratoire Interdisciplinaire sur l'Organisation Nanométrique et Supramoléculaire
NIMBE UMR 3685 - Nanosciences et Innovation pour les Matériaux, la Biomédecine et l'Energie (ex SIS2M)
Abstract : The process of inserting/deleting a particle during grand-canonical Monte-Carlo (MC) simulations is investigated using a novel, original technique: the trial event is made of a short nonequilibrium molecular dynamics (MD) trajectory during which a coordinate w along a 4$^{th}$ dimension is added to the particle in the course of insertion/deletion and is forced to decrease from large values down to zero (for insertion) or increased from 0 up to large values (for extraction) at $imposed$ $\nu_w$ velocity. The probability of acceptation of the whole MC move is controlled by the chemical potential and the external work applied during the trajectory. Contrary to the standard procedures which create/delete $suddenly$ a particle, the proposed technique gives time to the fluid environment to $relax$ during the gradual insertion/extraction before the acceptation decision. The reward for this expensive trial move is a gain of many orders of magnitude in the success rate. The power and wide domain of interest of this hybrid “H4D” algorithm which marries stochastic MC and nonequilibrium deterministic MD flavors are briefly illustrated with hard sphere, water, and electrolyte systems. The same approach can be easily adapted in order to $measure$ the chemical potential of a solute particle immersed in a fluid during canonical or isobaric simulations. It then becomes an efficient application of the Jarzynski theorem for the determination of solvation free energy.
Document type :
Journal articles
Complete list of metadatas

Cited literature [26 references]  Display  Hide  Download

https://hal-cea.archives-ouvertes.fr/cea-02183138
Contributor : Serge Palacin <>
Submitted on : Monday, July 15, 2019 - 10:33:05 AM
Last modification on : Wednesday, July 17, 2019 - 1:25:15 AM

File

Hybrid_4D_JCP_V0.pdf
Files produced by the author(s)

Identifiers

Citation

Luc Belloni. Non-equilibrium hybrid insertion/extraction through the 4th dimension in grand-canonical simulation. Journal of Chemical Physics, American Institute of Physics, 2019, 151 (2), pp.021101. ⟨10.1063/1.5110478⟩. ⟨cea-02183138⟩

Share

Metrics

Record views

57

Files downloads

59