Impact of Ionic Strength and Charge Density on Donnan Potential in the NaCl-Cation Exchange Membrane System

: This work aims to theoretically investigate the effect of both the ﬁxed charge density of ion exchange membranes and the ionic strength of the treated aqueous NaCl solution on the generated Donnan potential at thermodynamic equilibrium conditions. The direct objective of our work is to calculate the equilibrium concentration of the Cl − co-ion inside a swelled cation-exchange membrane equilibrated with a water/NaCl system. Two activity coefﬁcient models are employed, i


Introduction
Ion-exchange membranes (IEMs) are vital in electrically driven membrane processes.They are used in a number of membrane-based separation processes, such as electrodialysis and electrochemical desalination [1][2][3].In IEMs, stationary charged functional groups are attached to the polymer backbone of the polymeric membrane material.Ions of the same charge (i.e., co-ions) are completely or partially excluded from passing through the membrane.In contrast, ions of the opposite charge (i.e., counter-ions) are able to pass.Therefore, anion-exchange membranes have positively charged groups fixed to their polymeric skeleton and exclude positive ions (cations) but are freely permeable to negatively charged ions (anions).Similarly, cation exchange membranes with fixed negatively charged groups exclude anions but are freely permeable to cations.Electrical neutrality must be satisfied in all of these phenomena.
Industrial applications of ion-exchange membranes emerged more than half a century ago.Several IEM-based technologies have reached a satisfactory level of maturity that has advanced their utility to the large-scale commercial phase.For example, electrodialysis is now an essential technology in the brackish water desalination industry [4].IEM found practical applications in many fields, such as the demineralization of industrial waste and whey and sugar liquor [5] and the purification of amino acid solutions [6].Additionally, IEMs are utilized in several recent technologies, including electrodialysis reversal [7], bipolar membrane electrodialysis [8], electro-deionization [9], electrolysis [10], and diffusion dialysis [11].Yoshinobu Tanaka [12] authored an excellent book that addresses both the fundamentals and applications of IEMs.
When an IEM is equilibrated with a saline aqueous phase (including Na + , Mg +2 , Cl − , etc.), an unequal distribution of ions arises between the hydrated (swelled) membrane and the aqueous phase, resulting in a phenomenon known as the Donnan electrical potential at the membrane/solution interface [13].The Donnan potential results in a partial (or ideally a complete) exclusion of co-ions, providing the basis for the permselectivity of IEMs.
Experimental measurement of the Donnan potential is expected to give important insight into the IEM's effective charge density, as the Donnan potential is strongly linked to the IEM charge density.
Despite the recent attempts to experimentally measure individual IEM-interphase Donnan's potentials, it is generally acknowledged that such measurements are extremely difficult [14].Approximate methods based on measuring the overall potential difference in Donnan dialysis experiments were used in the past to predict the individual membraneinterphase Donnan potential [15].Therefore, well-established theoretical methods to estimate the Donnan potential are usually employed [15][16][17][18].
In modeling the performance of IEMs, the treated aqueous phase (e.g., seawater) is usually assumed to be either an ideal solution, where values of the ionic activity coefficients become equal to one, or a real solution with equal ionic activity coefficients in both, the membrane and solution phases.In both cases, ion activity coefficients in the membrane and solution phases are eliminated from the analysis; the activities of the ions in both phases are replaced by their corresponding concentrations.
This work investigates the effect of the ionic strength of the treated solution and the fixed charge density of the IEM on the Donnan potential and IEM permselectivity at thermodynamic equilibrium conditions.Specifically, the co-ion concentration inside the IEM will be computed for different values of the treated solution compositions and the IEM fixed charge values.This study will employ the Debye-Huckel (DH) and Meissner corresponding state models to examine the Donnan electrochemical potential at the membrane/solution interface, which is assumed to be in thermodynamic equilibrium.The DH model, used here as a reference model, is rigorous and becomes exact at the infinite dilution limit of the electrolyte concentration.As such, it provides an anchor asymptotic framework for all other models of non-ideal systems.On the other hand, the Meissner corresponding state model, although an old model, is of high simplicity and applicability for highly concentrated aqueous solutions with molalities exceeding saturation [19].Despite this, the model was not adequately addressed in the technical literature.In terms of the ionic mean activity coefficients (γ ± ), the model showed excellent predictive capability [19].Therefore, this "somewhat ignored" model will be used to compute the Donnan potential associated with an IEM equilibrated with an aqueous system of monovalent electrolyte, i.e., sodium chloride (NaCl).Specifically, the electrolyte mean activity coefficient in both phases, the aqueous feed and the swelled membrane, will be evaluated using the Meissner model, thus providing uniformity in calculating liquid phase nonideality.Experimental data for Donnan's potential of IEMs available in the literature will be used to verify the model's theoretical prediction.

The Donnan Equilibrium Theory
Consider a cationic exchange membrane in equilibrium with a saline aqueous phase (with ions like Na + , K + , Mg +2 , Ca +2 , Cl − , OH − , CO 3 −2 , etc.).The chemical potential of each distributed ion (µ i ) in the swelled membrane phase (m) and the aqueous phase (s) is given as: where µ o i is a reference chemical potential of species (i) at the system temperature (T) and a specified reference pressure (p o ), a i is the activity of species (i), and z i is the valence number of ion (i).Additionally, ψ is the electrical potential, F is Faraday constant (96,500 C/mole), p is the pressure, and V i is partial molar volume of species (i).Superscripts (s) and (m) refer to the solution (feed) phase and the membrane phase, respectively.The last term in each of Equations ( 1) and ( 2) is usually neglected based on its very small contribution [13], and it will be neglected initially in our analysis here.However, the effect of their inclusion will be discussed separately in the next section.
At equilibrium, equality of the chemical potential of each ion requires that: Substituting for µ s i and µ m i using Equations ( 1) and ( 2) (neglecting the last term in each) gives: The Donnan electric potential (E Don ) is (ψ m − ψ s ).Therefore, using Equation ( 4), we obtain the following: C i stands for the concentration of ions (i), and γ i is the activity coefficient of species (i).For a general salt C ν C A ν A (e.g., CaCl 2 ) that ionizes completely: where ν C and ν A stand for the stoichiometric numbers of the cation (C) and the anion (A) in the electrolyte C ν C A ν A .Writing Equation (4) for the cation and the anion and eliminating F (ψ m − ψ s )/RT between the two equations gives (with a i = C i γ i ): The equations representing electrical neutrality for the membrane and the external solution are: In Equation ( 8), q stands for the fixed charge density of the membrane in eq/liter.Equation ( 7), together with the electrical neutrality conditions for the membrane (Equation ( 8)) and the aqueous phase (Equation ( 9)), can be solved for the anion and cation concentrations inside the membrane and the Donnan potential as per Equation (5).
Substituting for C m C from Equation (8) and for C S C from Equation (9) into Equation (7) gives: In the case of a (1-1) electrolyte (e.g., NaCl), z C = +1 and z A = −1, and for the case of a (2-1) electrolyte (e.g., CaCl 2 ), z C = +2 and z A = −1, Equation (10) gives: ( 2 − 1) electrolyte ( 12) The mean activity coefficients in the membrane and aqueous phases γ s ± are functions of the phase composition and temperature.These are calculated using specific liquid-phase models.In this work, the Debye-Huckel and Meissner models are used.

Effect of Osmotic Pressure on Donnan's Potential
As shown in Section 2.1 above, the last term in Equations ( 1) and ( 2) is usually neglected based on its very small contribution [13].Neglecting this term is essentially equivalent to neglecting the coupling between ion concentrations and the osmotic pressure.Julian et al. [20] pointed out that neglecting osmotic pressure effects and considering only ionic equilibrium may lead to an incorrect Donnan potential with an error dependent on the osmotic pressure.Gang Chen [21] also addressed the same issue and indicated that the existing coupling between ion concentrations and the osmotic pressure leads to a membrane potential that consists of not only the classical Donnan potential term but also an additional term due to the osmotic pressure.This implies the existence of a membrane potential even when the impermeable species are not charged.
In order to take into account the effect of the osmotic pressure on Donnan potential, the steps implemented above, shown in Equations ( 1) through (12), are employed again, but the last term in each of Equations ( 1) and ( 2) is not neglected.The procedure followed is the same as that mentioned in [20,21].The equivalents of Equations ( 5), (11) and (12) when the pressure term is included are, respectively: where π stands for the osmotic pressure difference across the membrane (p m − p s ), which arises because of the electrolyte concentration difference between the swelled membrane phase and the feed phase.It is to be noted that the osmotic pressure π depends on the electrolyte concentration in both phases, including C m A , which is the unknown variable in Equations ( 15) and ( 16).The functional dependence of π on C m A and C s A is obtained by implementing the equilibrium criterion for the solvent (water) by equating the chemical potential of water in the swelled membrane phase and the feed phase (µ m water = µ s swater ).
Substituting for π from Equation ( 17) into Equations ( 15) and ( 16) for f(π) and g(π) gives: In the Results and Discussion section, we will show that under the prevailing conditions of the membrane fixed charge density and the external solution concentration, there is practically no effect for f(π) and g(π), i.e., both have an approximate value of unity.

The Debye-Huckel Theory
According to this theory, the nonideality of an electrolyte solution is expressed as the sum of the separate contributions of coulombic forces from solvent-ion and ion-ion interactions.These contributions fundamentally differ from the dispersive van der Waals forces acting between neutral molecules.The DH theory provides an exact model for the mean activity coefficient (γ ± ) at the electrolyte infinite dilution limit.Thus, it stands as a rigorous anchor point for most engineering models of electrolyte solutions.Debye-Huckel theory incorporates several key simplifications [21]: (1) the solvent is a structureless continuum with constant permittivity (ε) and dielectric strength (D s ); (2) the ions are spherically symmetric; and (3) the charge density around each specified ion is obeying a continuous Poisson-Boltzmann distribution.
The derivation of a mathematical model for the mean electrolyte activity coefficient γ ± , which is based on the DH theory, can be phenomenological and based on the concepts of classical electrostatics [21].It can also be derived based on the concepts of statistical mechanics [22].The phenomenological approach gives the following mathematical expression for the mean activity coefficient γ ± of an electrolyte in a solution: where I is the ionic strength of the solution, m i is the molality of species (i) (mole i/kg solvent), which becomes approximately c i (mole i/L solvent) for dilute solutions, ρ s is the solvent density, and D s is the dielectric constant for the solvent.

The Meissner Corresponding State Model
The DH model gets progressively worse as ionic strength increases beyond 0.01.The non-ideal behavior of electrolyte solutions, expressed in terms of the electrolyte mean activity coefficient γ ± , generally shows a very complex functional relationship with concentrations and temperature.Ideally, a robust theoretical model that predicts this complex functional relationship becomes a necessity.One example of such robust models, which was not adequately emphasized in the technical literature, is the Meissner corresponding state model.
Meissner and Tester [23] recognized a possible similar approach to generalize the non-ideal mean ionic activity coefficient (γ ± ) behavior of strong electrolytes.They were inspired by the great success of the law of the corresponding states in correlating the PVTN behavior of gases and liquids.The developed generalized model ably captured the very complex nature of the functional dependence of γ ± on electrolyte concentration and electrolyte chemical nature.For example, some electrolytes like HCl and LiCl show strong positive deviations (γ ± > 1), while others like CuSO 4 show negative deviations (γ ± < 1).The model maintains consistency with the DH at the infinite dilution limit of the electrolyte, i.e., γ ± = 1 as m = 0. Specifically, the reduced activity coefficient ( was found to provide a limiting result that was independent of electrolyte type.Meissner and Tester [23] explored the idea of using Γ o ij outside the DH region to correlate pure single electrolyte γ ± data for over 120 different electrolytes.It was found that log 10 Γ o ij versus I formed a family of non-intersecting curves quite analogous to the compressibility factor chart (Z = f(T r , P r ; ω)).What is remarkable about this model is its extrapolative nature to significantly higher concentrations, even up to ionic strengths of 20 or more [19].Further refinements by Meissner and others produced the final generalized correlation of the model, as follows [23]: The parameter q o ij was reported for many electrolytes at 25 • C and was extrapolated to other temperatures using the following formula: One of the most important features of Meissner's generalized correlation is that with only one data point for γ ± at a particular electrolyte concentration, extrapolation (or interpolation) to any other concentration value can be made.This is sometimes valid for compositions that even exceed the saturation value (solubility) of that salt, which is a valuable feature for multicomponent systems [21].However, one of the downsides of the Meissner model is its limited capability for extrapolation to temperatures higher than the base temperature of 298 K.This is a direct result of the fact that only one adjustable parameter is employed.The model is also limited in accurately capturing the effects of other cosolvents.Nonetheless, when considering water desalination as the intended use, these constraints do not pose significant limitations.

Numerical Solution Strategy
For a given aqueous phase with a known electrolyte concentration, (γ s ± ) is calculated using the selected liquid phase electrolyte model for the mean ionic activity coefficient, i.e., the DH and Meissner models.Together with the known charge density of the IEM (q) as an input parameter, Equation (11) or Equation (12) are solved iteratively for C m A , and hence Equation ( 5) is used to solve for E Don .The concentration of the cation inside the membrane (C m C ) is obtained directly from Equation (8).The flow chart for the numerical solution employed in this work and implemented using Excel spreadsheets (Microsoft ® Excel ® 2016) is shown in Figure 1.

Model Verification
Verifying Donnan's potential associated with IEMs in a certain electrolyte system is highly prioritized in developing novel perm-selective membranes.However, due to the lack of appropriate interfacial experimental techniques, direct measurement of Donnan's potential is complicated and difficult to attain.Most recently, Gokturk et al. [14] developed an experimental methodology for directly measuring Donnan potential for commercial IEMs equilibrated with electrolyte solutions and presented experimental results for the cases of NaCl and MgCl2 aqueous salt solutions.Their measurements are based on the binding energy shift in membrane-related core levels detected by a tender ambient pressure X-ray.The authors compared their experimental measurements with model predictions.They applied Manning's counter-ion condensation theory to predict the ion activity coefficients inside a commercial IEM and the Pitzer model to calculate the ion activity coefficients in the external salt solutions [24,25].This approach requires the input of basic membrane properties, but no adjustable parameters are needed.Reasonable agreement was observed between their experimental measurements and model predictions.

Model Verification
Verifying Donnan's potential associated with IEMs in a certain electrolyte system is highly prioritized in developing novel perm-selective membranes.However, due to the lack of appropriate interfacial experimental techniques, direct measurement of Donnan's potential is complicated and difficult to attain.Most recently, Gokturk et al. [14] developed an experimental methodology for directly measuring Donnan potential for commercial IEMs equilibrated with electrolyte solutions and presented experimental results for the cases of NaCl and MgCl 2 aqueous salt solutions.Their measurements are based on the binding energy shift in membrane-related core levels detected by a tender ambient pressure X-ray.The authors compared their experimental measurements with model predictions.They applied Manning's counter-ion condensation theory to predict the ion activity coefficients inside a commercial IEM and the Pitzer model to calculate the ion activity coefficients in the external salt solutions [24,25].This approach requires the input of basic membrane properties, but no adjustable parameters are needed.Reasonable agreement was observed between their experimental measurements and model predictions.
In this work, the Meissner model, a generalized corresponding state activity coefficient model that proved very reliable as a predictive model for treating highly concentrated solutions, is used to predict activity coefficients.Figure 2 compares the experimental measurements of Gokturk et al. [14] for both NaCl and MgCl 2 with model predictions using the Meissner model.Reasonable agreement is observed, which verifies to some extent both the applicability of the experimental approach of Gokturk and the suitability of the Meissner model for handling highly concentrated electrolyte aqueous solutions.Model prediction for both electrolytes overestimates Donnan's potential.In view of Equation ( 5), this is equivalent to saying that the Meissner model either underestimates the ionic activity inside the membrane relative to the actual value or overestimates the ionic activity in the external aqueous feed solution.However, the model was proven to give excellent predictions for aqueous phases of electrolytes up to and even exceeding a molality of six [21].The predicted concentration of the ion inside the membrane is probably underestimated by the model, which leads to lower than actual ionic activity inside the membrane.It is also observed that the Donnan potential is lower for higher-valence counter-ions (e.g., Mg +2 ), leading to reduced co-ion exclusion, which is in agreement with Gokturk et al. [14].
Water 2023, 15, x FOR PEER REVIEW 8 of 14 In this work, the Meissner model, a generalized corresponding state activity coefficient model that proved very reliable as a predictive model for treating highly concentrated solutions, is used to predict activity coefficients.Figure 2 compares the experimental measurements of Gokturk et al. [14] for both NaCl and MgCl2 with model predictions using the Meissner model.Reasonable agreement is observed, which verifies to some extent both the applicability of the experimental approach of Gokturk and the suitability of the Meissner model for handling highly concentrated electrolyte aqueous solutions.Model prediction for both electrolytes overestimates Donnanʹs potential.In view of Equation (5), this is equivalent to saying that the Meissner model either underestimates the ionic activity inside the membrane relative to the actual value or overestimates the ionic activity in the external aqueous feed solution.However, the model was proven to give excellent predictions for aqueous phases of electrolytes up to and even exceeding a molality of six [21].The predicted concentration of the ion inside the membrane is probably underestimated by the model, which leads to lower than actual ionic activity inside the membrane.It is also observed that the Donnan potential is lower for higher-valence counter-ions (e.g., Mg +2 ), leading to reduced co-ion exclusion, which is in agreement with Gokturk et al. [14].

Effects of Osmotic Pressure, Ionic Adsorption, and Ionic Pairing on Donnan Potential
As mentioned in Section 2.2, researchers have reported that coupling between ion concentrations and the osmotic pressure () can give rise to a membrane potential that consists of not only the classical Donnan potential term but also an additional term due to the osmotic pressure [20,21].According to these investigators, this coupling can lead to the existence of a membrane potential even when the impermeable species are not charged [21].For a (1-1) electrolyte like NaCl, Equations ( 14) and (17) show that the osmotic pressure effect is given by the dimensionless function (), which is simply the ratio of the activity of water in the feed and the membrane phases ( / ) raised to the exponent ( /2 ).The activity of water in each phase is a function of the water mole fraction in that phase, which in turn depends on the equilibrium composition of the electrolyte.Since the water mole fraction in both phases, the swelled membrane phase (m) and the feed phase (s), is not drastically different from each other, the activity ratio is not expected to significantly differ from unity.Our computation shows that for an extreme case of a membrane fixed charge density of 15 eq/L,  (the feed phase mole fraction of water) is less than 5% different from  (the membrane phase mole fraction of water).

Effects of Osmotic Pressure, Ionic Adsorption, and Ionic Pairing on Donnan Potential
As mentioned in Section 2.2, researchers have reported that coupling between ion concentrations and the osmotic pressure (π) can give rise to a membrane potential that consists of not only the classical Donnan potential term but also an additional term due to the osmotic pressure [20,21].According to these investigators, this coupling can lead to the existence of a membrane potential even when the impermeable species are not charged [21].For a (1-1) electrolyte like NaCl, Equations ( 14) and (17) show that the osmotic pressure effect is given by the dimensionless function f(π), which is simply the ratio of the activity of water in the feed and the membrane phases (a s water /a m water ) raised to the exponent (V solute /2V water ).The activity of water in each phase is a function of the water mole fraction in that phase, which in turn depends on the equilibrium composition of the electrolyte.Since the water mole fraction in both phases, the swelled membrane phase (m) and the feed phase (s), is not drastically different from each other, the activity ratio is not expected to significantly differ from unity.Our computation shows that for an extreme case of a membrane fixed charge density of 15 eq/L, x s w (the feed phase mole fraction of water) is less than 5% different from x m w (the membrane phase mole fraction of water).This implies that, within the parameters examined in this study, the impact of osmotic pressure is expected to be negligible.To verify this, the Donnan potential was recalculated with the osmotic pressure taken into account for two extreme cases of the membrane fixed charge density, i.e., 1 eq/L and 15 eq/L.The difference between Donnan potential corrected for osmotic pressure and Donnan potential without osmotic pressure correction is shown in Figure 3.The required partial molar volume of water V water is taken as 18 cm 3 /mol, and the partial molar volume of the (NaCl) electrolyte V solute is calculated using the experimental data recently reported by Khons et al. [26].No effect of the osmotic pressure is observed for a fixed charge density of 1 eq/L.Even for the extreme case of 15 eq/L, the maximum difference in Donnan potential is 1 mV, less than 5% of the conventional Donnan potential.The observed maximum in Figure 3 for the case of 15 eq/L is a result of the opposing values of the two terms appearing in Equation ( 14), i.e., the activity term and the osmotic pressure term.
Water 2023, 15, x FOR PEER REVIEW 9 of 14 This implies that, within the parameters examined in this study, the impact of osmotic pressure is expected to be negligible.To verify this, the Donnan potential was recalculated with the osmotic pressure taken into account for two extreme cases of the membrane fixed charge density, i.e., 1 eq/L and 15 eq/L.The difference between Donnan potential corrected for osmotic pressure and Donnan potential without osmotic pressure correction is shown in Figure 3.The required partial molar volume of water  is taken as 18 cm 3 /mol, and the partial molar volume of the (NaCl) electrolyte  is calculated using the experimental data recently reported by Khons et al. [26].No effect of the osmotic pressure is observed for a fixed charge density of 1 eq/L.Even for the extreme case of 15 eq/L, the maximum difference in Donnan potential is 1 mV, less than 5% of the conventional Donnan potential.The observed maximum in Figure 3 for the case of 15 eq/L is a result of the opposing values of the two terms appearing in Equation ( 14), i.e., the activity term and the osmotic pressure term.The effects of adsorbed ions on the Donnan potential, which develops at the IEM interphase, are particularly interesting in ion exchange membranes [27,28].It is worth mentioning that ion adsorption on localized membrane sites gives rise to a non-ideal or modified Donnan phenomenon.This phenomenon arises because of the electrostatic interaction between adsorbed ions and the dependence of the concentration of adsorbed ions on external solution conditions such as electrolyte concentration and pH [29].The modified Donnan potential is usually computed by augmenting an adsorption isotherm that incorporates the interaction between adsorbed ions into the Donnan formalism [30].Such theoretical studies are of special importance as they delineate the conditions under which interaction between adsorbed ions is important enough to be considered.The interaction of water molecules with charged surfaces and the influence of the surface on the structure of the extended hydrogen bond network is a subject of ongoing extensive studies [31,32].
Another phenomenon that has a similar effect on Donnan equilibrium is that of the counterion-membrane-fixed ion pairing in charged membranes with low dielectric constant regions.Here, the ideal Donnan equilibrium is extended to account for the ion pairing effect, such as the Fuoss approach for contact ion pairs in electrolyte solutions with a correction term for the entropy change in the counter-ion undergoing ion pairing [33].The effects of adsorbed ions on the Donnan potential, which develops at the IEM interphase, are particularly interesting in ion exchange membranes [27,28].It is worth mentioning that ion adsorption on localized membrane sites gives rise to a non-ideal or modified Donnan phenomenon.This phenomenon arises because of the electrostatic interaction between adsorbed ions and the dependence of the concentration of adsorbed ions on external solution conditions such as electrolyte concentration and pH [29].The modified Donnan potential is usually computed by augmenting an adsorption isotherm that incorporates the interaction between adsorbed ions into the Donnan formalism [30].Such theoretical studies are of special importance as they delineate the conditions under which interaction between adsorbed ions is important enough to be considered.The interaction of water molecules with charged surfaces and the influence of the surface on the structure of the extended hydrogen bond network is a subject of ongoing extensive studies [31,32].
Another phenomenon that has a similar effect on Donnan equilibrium is that of the counterion-membrane-fixed ion pairing in charged membranes with low dielectric constant regions.Here, the ideal Donnan equilibrium is extended to account for the ion pairing effect, such as the Fuoss approach for contact ion pairs in electrolyte solutions with a correction term for the entropy change in the counter-ion undergoing ion pairing [33].

Effects of Membrane Fixed Charge Density and Electrolyte Concentration
In our analysis, the effects of osmotic pressure, ion adsorption, and ion pairing are not included.The osmotic pressure shown above does not have a significant effect under the conditions of this study.Adsorption and ion pairing interactions with Donnan potential could be an interesting topic that can be addressed in a separate study; including them in this manuscript could potentially lead to a substantial increase in its overall length.It is worth mentioning that the saturation solubility of NaCl in water at 298 K is about 6 moles of salt per kg of water and that it is not sensitive to moderate temperature changes [34].Despite this, it is not uncommon to find in the technical literature research studies on NaCl solutions with compositions exceeding the saturation solubility [35].Such supersaturated aqueous salt solutions, where solutes in excess of the expected equilibrium solubility are encountered, are a very well-known phenomenon, especially in the crystallizationrelated technical literature.However, such systems are thermodynamically metastable [19].Another phenomenon that can contribute to supersaturation is ion association.In a simulation study by Soniat et al. [36], it was found that for NaCl, NaI, and KCl, larger clusters are present but with a small amount of pairing.
The effect of the IEM fixed charge density (q) on the generated Donnan potential, using the classical DH and Meissner models, is shown in Figure 4.The results in this figure reveal excellent agreement between the DH and Meissner models for electrolyte concentrations up to a molality of 0.1 (approximately 6000 ppm).Beyond this concentration, however, the DH model overestimates the predicted (absolute) Donnan potential relative to the Meissner model.The difference between the predictions of the two models grows with the electrolyte feed concentration.The overestimate of the predicted Donnan potential by the DH model is attributed to the basic assumptions inherent in the DH model.Also, at higher feed concentrations of the electrolyte, the mismatch between the two models grows with the level of the IEM fixed charge.For example, at an electrolyte concentration of 0.2 M (11,688 ppm), the DH model overestimates Donnan potential by approximately the same amount of 5 mV relative to the Meissner model for the three levels of fixed charge density investigated in this study, i.e., 1, 2, and 5 equivalent/L.However, the corresponding overestimates of Donnan potential for an electrolyte concentration of 5 M (292,200 ppm) are 13, 29, and 50 mV for the fixed charge densities of 1, 2, and 5 equivalent/L, respectively.This can result from a higher degree of deviation from ideal solution behavior (i.e., the equal activity coefficient of the ions inside the membrane and in the aqueous feed phase), an inherent assumption in the DH model.Figure 4 also reveals that a higher (absolute) Donnan potential is obtained by a higher IEM fixed charge.This is expected given the higher exclusions of co-ions at the high IEM fixed charge [12].
In our analysis, the effects of osmotic pressure, ion adsorption, and ion pairing are not included.The osmotic pressure shown above does not have a significant effect under the conditions of this study.Adsorption and ion pairing interactions with Donnan potential could be an interesting topic that can be addressed in a separate study; including them in this manuscript could potentially lead to a substantial increase in its overall length.It is worth mentioning that the saturation solubility of NaCl in water at 298 K is about 6 moles of salt per kg of water and that it is not sensitive to moderate temperature changes [34].Despite this, it is not uncommon to find in the technical literature research studies on NaCl solutions with compositions exceeding the saturation solubility [35].Such supersaturated aqueous salt solutions, where solutes in excess of the expected equilibrium solubility are encountered, are a very well-known phenomenon, especially in the crystallization-related technical literature.However, such systems are thermodynamically metastable [19].Another phenomenon that can contribute to supersaturation is ion association.In a simulation study by Soniat et al. [36], it was found that for NaCl, NaI, and KCl, larger clusters are present but with a small amount of pairing.
The effect of the IEM fixed charge density (q) on the generated Donnan potential, using the classical DH and Meissner models, is shown in Figure 4.The results in this figure reveal excellent agreement between the DH and Meissner models for electrolyte concentrations up to a molality of 0.1 (approximately 6000 ppm).Beyond this concentration, however, the DH model overestimates the predicted (absolute) Donnan potential relative to the Meissner model.The difference between the predictions of the two models grows with the electrolyte feed concentration.The overestimate of the predicted Donnan potential by the DH model is attributed to the basic assumptions inherent in the DH model.Also, at higher feed concentrations of the electrolyte, the mismatch between the two models grows with the level of the IEM fixed charge.For example, at an electrolyte concentration of 0.2 M (11,688 ppm), the DH model overestimates Donnan potential by approximately the same amount of 5 mV relative to the Meissner model for the three levels of fixed charge density investigated in this study, i.e., 1, 2, and 5 equivalent/L.However, the corresponding overestimates of Donnan potential for an electrolyte concentration of 5 M (292,200 ppm) are 13, 29, and 50 mV for the fixed charge densities of 1, 2, and 5 equivalent/L, respectively.This can result from a higher degree of deviation from ideal solution behavior (i.e., the equal activity coefficient of the ions inside the membrane and in the aqueous feed phase), an inherent assumption in the DH model.Figure 4 also reveals that a higher (absolute) Donnan potential is obtained by a higher IEM fixed charge.This is expected given the higher exclusions of co-ions at the high IEM fixed charge [12].predictions beyond that concentration.The DH model predicts higher exclusion, i.e., a lower concentration of the co-ion inside the IEM. Figure 5 presents the effect of the electrolyte concentration in the feed on the co-ion concentration inside the IEM for three levels of the fixed charge concentration as predicted by the DH and the Meissner model.The two models agree in their predictions for electrolyte feed concentrations up to 5000 ppm.However, the two models seriously disagree in their predictions beyond that concentration.The DH model predicts higher exclusion, i.e., a lower concentration of the co-ion inside the IEM.Molality of coion inside the IEM Molality of NaCl in the feed aquous phase q = 1 eq/l; Meissner q = 2 eq/l; Meissner q = 1 eq/l; Debye-Huckel q = 2 eq/l; Debye-Huckel q = 5 eq/l; Meissner q = 5 eq/l; Debye-Huckel

Conclusions
Ion exchange membranes play a vital role in electrically driven membrane processes.The effect of the ionic strength of the treated feed solutions and the fixed charge density of the ion exchange membrane on the associated Donnan potential is investigated using two predictive models, i.e., the Debye-Huckel and the Meissner models.The latter model

Figure 1 .
Figure 1.The flow chart implemented in calculating the Donnan potential in the NaCl-IEM system.

Figure 1 .
Figure 1.The flow chart implemented in calculating the Donnan potential in the NaCl-IEM system.

Figure 2 .
Figure 2. Comparison between experimental Donnan potential in the NaCl−IEM system and Meissner model prediction.

Figure 2 .
Figure 2. Comparison between experimental Donnan potential in the NaCl−IEM system and Meissner model prediction.

Figure 3 .
Figure 3.Effect of electrolyte molality and membrane fixed charge density on the difference in Donnan potential with and without osmotic pressure correction in the NaCl−IEM system.

Figure 3 .
Figure 3.Effect of electrolyte molality and membrane fixed charge density on the difference in Donnan potential with and without osmotic pressure correction in the NaCl−IEM system.

Figure 4 .
Figure 4. Predicted Donnan potential for the NaCl−IEM system at 298 K for three levels of the fixed charge.

Figure 5
Figure5presents the effect of the electrolyte concentration in the feed on the co-ion concentration inside the IEM for three levels of the fixed charge concentration as predicted by the DH and the Meissner model.The two models agree in their predictions for electrolyte feed concentrations up to 5000 ppm.However, the two models seriously disagree in their

Figure 4 .
Figure 4. Predicted Donnan potential for the NaCl−IEM system at 298 K for three levels of the fixed charge.

Figure 5 .
Figure 5. Predicted co-ion concentration inside the IEM at 298 K for three levels of the fixed charge.

Figure 6
Figure 6 shows the effect of the IEM fixed charge density on the membrane anion-tocation concentration ratio ( ) / m m Cl Na C C − + , as predicted by both the DH and Meissner models.This ratio is indicative of the IEM perm-selectivity.The DH model agreement with the Meissner model breaks after approximately an electrolyte concentration of 1 M.The DH model predicts a total exclusion (zero concentration) of the co-ion for the fixed charge of 5 equivalents per liter and a minimum exclusion (maximum concentration) at a bulk electrolyte concentration of 3 M for the case of a fixed charge of 1 equivalent per liter.On the other hand, the Meissner model correctly predicts the 100% asymptotic limiting value of ( ) / m m Cl Na C C − +as the bulk concentration of the electrolyte increases.

Figure 5 .
Figure 5. Predicted co-ion concentration inside the IEM at 298 K for three levels of the fixed charge.

Figure 6 14 Figure 6 .
Figure 6 shows the effect of the IEM fixed charge density on the membrane anionto-cation concentration ratio C m Cl − /C m Na + , as predicted by both the DH and Meissner models.This ratio is indicative of the IEM perm-selectivity.The DH model agreement with the Meissner model breaks after approximately an electrolyte concentration of 1 M.The DH model predicts a total exclusion (zero concentration) of the co-ion for the fixed charge of 5 equivalents per liter and a minimum exclusion (maximum concentration) at a bulk electrolyte concentration of 3 M for the case of a fixed charge of 1 equivalent per liter.On the other hand, the Meissner model correctly predicts the 100% asymptotic limiting value of C m Cl − /C m Na + as the bulk concentration of the electrolyte increases.Water 2023, 15, x FOR PEER REVIEW 12 of 14

Figure 6 .
Figure 6.Effect of the IEM fixed charge on C m Cl − /C m Na + as predicted by the DH and Meissner models.