MOLECULAR-DYNAMICS SIMULATION OF IONS IN CHARGED CAPILLARIES

Electroneutrality and chemical potential balance in conventional model s for ions in charged capillaries are discussed. The Poisson-Boltzmann (PB) equation and molecular dynamics (MD) simulation of ions inside a n infinitely long charged pore show that inconsistencies exist. The PB solution shows that chemical potentials of individual ions are the sa me and electroneutrality is obeyed provided that the correct boundary condition is used at the wall. The MD simulations show that for a neut ral system, the chemical potentials of the ions are different and henc e they cannot be in equilibrium with a bulk solution. Modifications of existing models and set-ups of simulations are necessary to model the behaviour of ions in capillaries in equilibrium with a bulk solution.