A mixed alchemical and equilibrium dynamics to simulate heterogeneous dense fluids: Illustrations for Lennard-Jones
1Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712, USA.
Abstract:
An algorithm to efficiently simulate multi-component fluids is proposed and illustrated. The focus is on biological membranes that are heterogeneous and challenging to investigate quantitatively. To achieve rapid equilibration of spatially inhomogeneous fluids, we mix conventional molecular dynamics simulations with alchemical trajectories. The alchemical trajectory switches the positions of randomly selected pairs of molecules and plays the role of an efficient Monte Carlo move. It assists in accomplishing rapid spatial de-correlations. Examples of phase separation and mixing are given in two-dimensional binary Lennard-Jones fluid and a DOPC-POPC membrane. The performance of the algorithm is analyzed, and tools to maximize its efficiency are provided. It is concluded that the algorithm is vastly superior to conventional molecular dynamics for the equilibrium study of biological membranes.
More Related Videos
06:37Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
Related Concept Videos
Dynamic Equilibrium
Mixtures of Acids
A Mixture of a Strong Acid and a Weak Acid
In a mixture of a strong acid and a weak acid, the strong acid dissociates completely and becomes a source of almost all the hydronium ions...
Mixtures of Acids
In a strong and weak acid mixture, the strong acid dissociates completely and becomes a source of almost all the hydronium ions present in the solution. In contrast, the weak acid shows...
The Fluid Mosaic Model
Solution Equilibrium and Saturation
Free Energy and Equilibrium
Recall that Q is the numerical value of the mass action...
