Early molecular dynamics methods for osmotic pressure, introduced by Murad and Powles [1], utilized semipermeable membranes within the simulation to separate solvent and solution regions. Luo and Roux [2] refined this approach by employing flat-bottom potentials as virtual semi-permeable walls to simplify the process and avoid complications from direct membrane modeling. These techniques require multiple simulations to explore the concentration dependence of the osmotic equation-of-state. To streamline this, Milner et al. [3] developed the osmotic force balance method, applying an external harmonic potential to solutes, thereby enabling easier assessment of osmotic pressure gradients. They successfully applied this to study substances like NaCl and mixtures like benzene and pyridine. Hosseini et al. [4] expanded on this in 2023, calculating chemical potentials for 15 alkali halide salts, thereby providing a more efficient method to determine osmotic pressures and propose Lennard-Jones parameters for these salts. I have replicated the calculations of osmotic pressure with respect to concentration using both techniques.
Simulating virtual semi-permeable walls using flat-bottomed potentials (figure 1) allows for the calculation of osmotic pressure for a single concentration. The mean force of the system's walls exerted on the ions is directly related to the osmotic pressure as shown in equation 1, where Fwall is the instantaneous force exerted by one wall onto the ions, k is the force constant, N is the number of steps, zi is the position in the z axis of the ions located outside of both the walls, and zwall is the position of the wall. The mean force is averaged between the two half-harmonic walls. The osmotic pressure was calculated as shown in equation 2, where A is the cross-sectional area of the simulation box, and with it, the osmotic coefficients are calculated as a ratio between the observed and ideal results.
Figure 1: 1m NaCl simulation with outlined flat-bottom potential (red line).
A harmonic potential allows for the calculation of osmotic pressure for a nonuniform equilibrium concentration profile for ions (figure 2). The potential curve used in this method uses the form of , where k is the spring constant and z is the z coordinate centered in the box. In equilibrium, slices of a concentration profile are regarded as stationary under the sum of three forces: the gradient of external potential (FU, eq. 3), and the osmotic pressure (FΠ, eq. 4) on the two sides of the slice, which will differ slightly because of the nonuniform concentration. C(z) is the total concentration at height z, and Δz is the length of the slice.
When the system is in equilibrium these forces sum to 0. The change in osmotic pressure across a slice is proportional to the change in the harmonic potential multiplied by the concentration of the slice, as seen in equation 5. One can integrate the force balance equation after plugging in it the harmonic potential, as seen in equation 6, and with it solve for the osmotic pressure at the coordinates of both sides of a slice by averaging over their corresponding concentration.
Figure 2: 3m NaCl simulation with outlines harmonic potential (red line) and concentration gradient (orange line).
- S. Murad and J. G. Powles, “Computer simulation of osmosis and reverse osmosis in solutions,” Chemical Physics Letters, vol. 225, no. 4, pp. 437–440, 1994, doi: 10.1016/0009-2614(94)87108-6.
- Y. Luo and B. Roux, “Simulation of Osmotic Pressure in Concentrated Aqueous Salt Solutions,” J. Phys. Chem. Lett., vol. 1, no. 1, pp. 183–189, 2010, doi: 10.1021/jz900079w.
- C. Gillespie and S. T. Milner, “Using osmotic pressure simulations to test potentials for ions,” Soft Matter, vol. 16, no. 42, pp. 9816–9821, 2020, doi: 10.1039/D0SM00957A.
- A. Hosseini and H. S. Ashbaugh, “Osmotic Force Balance Evaluation of Aqueous Electrolyte Osmotic Pressures and Chemical Potentials,” J. Chem. Theory Comput., vol. 19, no. 23, pp. 8826–8838, 2023, doi: 10.1021/acs.jctc.3c00982.