We estimate the probability of ion intrusion by using the raw values of the alchemical work and chemical potential reported in Tables and. At bulk molality of 2.41 mol/kg for monovalent electrolytes CA the probability of ion intrusion is respectively about ∼ 10 -6 , ∼ 10 -13 and ∼ 10 -18 -10 -20 for faujasite, chabazite and silicalite-1. The probability of ion intrusion for divalent-cation electrolytes CA 2 at bulk molality of 2.52 mol/kg is lower than for monovalent electrolytes with probabilities of ∼ 10 -18 , ∼ 10 -40 and ∼ 10 -60 respectively, because it is less favorable to solvate simultaneously three ions in the zeosils than an ion pair in case of monovalent electrolytes. For bulk electrolyte at high molalities (i.e. at 10.40 mol/kg for monovalent electrolytes), the probability of electrolyte intrusion increases for all zeosils, but remains insignificant for silicalite (∼10 -16 -10 -17 ) and chabazite (∼10 -10 -10 -12 ). On the contrary, the probability of electrolyte intrusion in faujasite becomes significant (∼10 -4 -∼10 -1 ), indicating that electrolyte may penetrate the zeosil at such concentration. Finally, in the presence of divalent electrolytes, insertion probability in the three zeosils also increases but remains low so that no electrolyte may be intruded as well. In Figure , we presented the distribution of the alchemical work for inserting a first NaCl unit in silicalite-1, chabazite and faujasite. The histogram queue for silicalite-1 and chabazite reveals that several intrusion of NaCl electrolyte are associated with insertion probabilities P acc of ∼ 10 -1 for bulk molality of 10.09 mol/kg. However, after eventual intrusion of a NaCl pair, it is unlikely that the pair remains in the framework because alchemical work associated with the reverse move for transforming the NaCl into water molecules (denoted as W -N aCl in our former work ) is the opposite of W N aCl resulting in a very high probability (P del ∝ e -β(∆µ CnAm -W CnAm ) ) to remove the ion pair.