We chose the force field developed for silicalite-1 by Desbiens and co-workers to model interactions for zeosils in this study. This force field is tailored to reproduce intrusion pressure of TIP4P water around ∼ 100 MPa, as observed in experiments. Water molecules are modeled with the rigid TIP4P model and the electrolytes are represented as charged beads that differs by their Lennard-Jones parameters and charges. We considered the following monovalent electrolytes which are: LiCl, NaCl, KCl, RbCl, CsCl, KF, KBr and KI. We also explored the following divalent electrolytes to probe influence of charge ion: MgCl We used Lorentz-Berthelot mixing rules for cross-terms in Lennard-Jones potential, and Ewald summation to account efficiently for electrostatic interactions. We used a cutoff of 8.0 Å for both the Lennard-Jones potential and the separation between real space and reciprocal space in the Ewald summation.