Nanotech 2003 Vol. 3
Nanotech 2003 Vol. 3
Technical Proceedings of the 2003 Nanotechnology Conference and Trade Show, Volume 3

Nanostructures: Biological Ion Channels to Thin Oxides Chapter 10

Equilibrium Structure of Electrolyte calculated using Equilibrium Monte Carlo, Molecular Dynamics, and Boltzmann Transport Monte Carlo Simulations

Authors: T.A. van der Straaten, G. Kathawala, Z. Kuang, D. Boda, D.P. Chen, U. Ravaioli, R.S. Eisenberg and D. Henderson

Affilation: University of Illinois, United States

Pages: 447 - 451

Keywords: pair correlation function, electrolyte, comparitive simulation

We use the pair correlation function as a benchmark calculation to compare three simulation methodologies -Equilibrium Monte Carlo, Molecular Dynamics, and Boltzmann Transport Monte Carlo. The three simulation approaches handle electrostatic boundary conditions and ion-water interactions, and evaluate electrostatic forces quite differently. The system simulated is a simple homogeneous electrolyte at equilibrium at concentrations of biological interest. Ions interact with each other via the electrostatic field and by either the standard Lennard-Jones interaction potential (LJ), or by a truncated form of LJ that mimics the ionic core repulsive part of the interaction potential. Ion-water interactions are represented by the primitive model, in which the water is modeled as a uniform dielectric medium. The pair correlation function has been calculated for monovalent and divalent electrolytes using both interaction potentials (standard LJ and truncated LJ). Initial results show good comparison between the three different simulations. We also show that the Boltzmann Transport Monte Carlo used in conjunction with the Particle-Particle-Particle Mesh (P3M) technique [2] allows the electrostatic forces to be resolved on a much coarser mesh while still recovering the correct pair correlation function, resulting in a considerable reduction in computational load.

ISBN: 0-9728422-2-5
Pages: 560