Related Experiment Video
Updated: May 18, 2026

Deciphering the Structural Effects of Activating EGFR Somatic Mutations with Molecular Dynamics Simulation
Published on: May 20, 2020
Grand-canonical Monte Carlo method for Donnan equilibria
S A Barr1, A Z Panagiotopoulos
1Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544, USA.
Abstract:
We present a method that enables the direct simulation of Donnan equilibria. The method is based on a grand-canonical Monte Carlo scheme that properly accounts for the unequal partitioning of small ions on the two sides of a semipermeable membrane, and can be used to determine the Donnan electrochemical potential, osmotic pressure, and other system properties. Positive and negative ions are considered separately in the grand-canonical moves. This violates instantaneous charge neutrality, which is usually considered a prerequisite for simulations using the Ewald sum to compute the long-range charge-charge interactions. In this work, we show that if the system is neutral only in an average sense, it is still possible to get reliable results in grand-canonical simulations of electrolytes performed with Ewald summation of electrostatic interactions. We compare our Donnan method with a theory that accounts for differential partitioning of the salt, and find excellent agreement for the electrochemical potential, the osmotic pressure, and the salt concentrations on the two sides. We also compare our method with experimental results for a system of charged colloids confined by a semipermeable membrane and to a constant-NVT simulation method, which does not account for salt partitioning. Our results for the Donnan potential are much closer to the experimental results than the constant-NVT method, highlighting the important effect of salt partitioning on the Donnan potential.
Related Concept Videos
Chemical Equilibria: Systematic Approach to Equilibrium Calculations
The first step is to identify all the chemical reactions involved, The...
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Ladder Diagrams: Complexation Equilibria
The formation constant, K1, for the formation of Cd(NH3)2+ complex from cadmium and ammonia is 3.55 × 102. Log K1 (i.e. pNH3) is 2.55, and...
Alternative Sets of Equilibrium Equations
One example of such a situation can be observed in a...
Homogeneous Equilibria for Gaseous Reactions
For gas-phase reactions, the equilibrium constant may be expressed in terms of either the molar concentrations (Kc) or partial pressures (Kp) of the reactants and products. A relation between these two K values may be simply derived from the ideal gas equation and the definition of molarity. According to the ideal gas equation:
The Response of Equilibria to the Conditions
