Predictive Approach to the Phase Behavior of Polymer-Water-Surfactant-Electrolyte Systems Using a Pseudosolvent Concept

1. Introduction

Interactions between surfactants and polymers is an important phenomenon relevant to several practical applications in paints and coatings, personal care products, chemical, pharmaceutical, mineral processing and petroleum industries [1,2,3,4,5,6]. The simultaneous presence of polymer and surfactant molecules alter the rheological properties of solutions, adsorption characteristics at solid-liquid interfaces, stability of colloidal dispersions, the solubilization capacities in water for sparingly soluble molecules, and liquid-liquid interfacial tensions, and thereby impact either beneficially or adversely the practical applications. An important application that has occupied the attention of researchers over the last five decades is the enhanced recovery of petroleum from existing reservoirs by chemical flooding [2,3,4,5]. In the chemical flooding process for enhanced oil recovery (EOR), a surfactant slug (often composed of anionic petroleum sulfonates) is injected into the ground to displace the formation oil. It is followed by a polymer (typically anionic or nonionic)-thickened mobility buffer and drive water. The surfactant slug is designed to generate ultra-low interfacial tensions against the formation oil. The polymer buffer is designed to possess adequate mobility control characteristics. The simultaneous presence of polymer and surfactant molecules in the reservoir results in interactions between these molecules as well as with the minerals in the reservoir rocks. These interactions among polymer, surfactant and electrolytes can drastically alter the properties of the surfactant slug and/or of the polymer mobility buffer from their designed characteristics, specifically due to changes in their phase behavior, leading to significant reduction in the anticipated oil recovery.

Many experimental studies on the phase behavior of solutions containing anionic or nonionic polymer molecules and anionic surfactant molecules have been carried out since early 1970s, but surprisingly no theoretical treatment has been developed to date. Note that we are not concerned here with oppositely charged polymers and surfactants but polymers and surfactants, typically both being anionic. Most of the early phase behavior studies came from the petroleum industry in the course of the EOR research in laboratory as well as in large scale oilfield reservoirs. The earliest experimental work is that of Trushenski et al. [7,8], who studied the phase behavior of solutions containing mixtures of Mahogony AA sulfonate (a petroleum derived surfactant), isopropyl alcohol, water, sodium chloride and either of the two polymers, Pusher 700 (a hydrolyzed polyacrylamide) or Kelco's Xanflood (a Xanthan biopolymer). The interest was on identifying the domain of single-phase existence in the multicomponent system and the composition boundaries at which a phase separation occurs. They found that the phase behavior of the systems differed very little for the two polymers. The phase behavior was also not affected by the solution concentration of the polymer in the range 500 to 1500 ppm. The single-phase region in the composition space of the phase diagram expanded when the temperature increased from 110°F to 190°F. Further, the addition of small quantities of oil was found to stabilize an otherwise unstable (phase separating) polymer-surfactant-electrolyte system. Szabo [9] examined systems constituted from polymers including Xanflood, hydroxyethyl and methyl celluloses, hydrolyzed and unhydrolyzed polyacrylamide, and polyacrylamido-methyl propane-sulfonate and four different petroleum sulfonate type surfactants. The phase behavior of solutions were monitored in systems containing 5 weight percent sulfonate and 2 weight percent NaCl. In all cases, phase separation into two liquid phases occurred on contacting the different components. The top phase was found to contain most of the polymer while the bottom phase contained most of the sulfonate. As the polymer concentration was increased, the volume fraction of the surfactant rich bottom phase decreased. The polymer concentration in the top phase and the sulfonate concentration in the bottom phase both increased. However, the polymer concentration in the bottom phase and the sulfonate concentration in the top phase remained constant. No systematic differences were seen with the variation in the type of polymer molecules.

Pope et al. [10] studied the effect of salt concentration on the phase behavior of polymer-surfactant solutions. The amount of NaCl required to induce phase separation in polymer-surfactant solutions was measured for combinations of different polymers and surfactants and at various concentrations of them in solution. The electrolyte concentration at which phase separation occurs, designated as the critical electrolyte concentration (CEC), was found to be independent of the polymer type and polymer concentration (in the range 100 to 1000 ppm). No variation of CEC with the molecular weight of the polymer (in the range 4 x 10⁵ to 5 x 10⁶) was observed for polyethylene oxide.

Gupta [11] examined the effect of calcium salts on the phase behavior of Pusher 700 - Mahogony AA sulfonate solution, containing 1600 ppm of the polymer. At 50 ppm concentration of calcium the solution was stable, whereas at 300 ppm the solution separated into a polymer-rich top phase and a surfactant-rich bottom phase. In general, the influence of the divalent electrolyte was more severe on the phase stability compared to that of the monovalent electrolyte.

Kalpakci [12] investigated the phase behavior of several anionic and nonionic surfactant solutions in the presence of polymers. He found that the phase stability of the solutions decreased with increasing concentration of the surfactant and of the electrolyte as well as with increasing equivalent weight of the surfactant. The polymer size in the molecular weight range 7 x 10⁴ to 2 x 10⁶ and polymer concentration between 100 and 1500 ppm were found to have very little effect on the stability of the anionic sulfonate-polymer solutions. It was found that the polymer-surfactant solutions can be stabilized by the addition of low molecular weight additives such as alcohols or low molecular weight sulfonates a well as by the addition of small amounts of oil. Further, unstable polymer-surfactant solutions were stabilized if the solution was subjected to ultrasonication resulting in the breakdown of the molecular size of the polymer.

These experimental results available in the early EOR literature show that phase separation is a common phenomenon occurring in polymer-surfactant solutions in the presence of electrolytes, even at low concentrations of various components. The phase separation patterns are not significantly influenced by the molecular type or the molecular size of the polymer for the high molecular weight polymers studied and by the variations in polymer concentrations over the limited range of 100 to 1500 ppm typically investigated for EOR applications. The phase behavior is, however, critically affected by the equivalent weight of the surfactant and by the type and amount of the electrolyte added. The stability of the polymer-surfactant solution is enhanced if the polymer molecular weight is significantly reduced as in the ultrasonication experiment or by the addition of certain types of additives including small amounts of oil. Beyond these generalizations, little is known concerning the phase behavior of polymer-surfactant solutions from a predictive point of view, even after forty years. It is well recognized in recent oil field simulations, that reliable information on the polymer-surfactant phase behavior is needed for making realistic predictions of oil recovery efficiencies [3,4,5].

The principal goal of this work is to develop a predictive approach to the phase behavior of aqueous solutions containing polymer, surfactant, and electrolyte molecules. Specifically, we seek to develop a thermodynamic treatment of phase behavior and compare the predictions of the treatment with experimental results obtained on reasonably well-defined chemical components in solution. Towards this goal, we postulate a novel approach to treat the four-component system of polymer, surfactant, electrolyte, and water by viewing it as a pseudo binary system. The essence of this approach is as follows. One component of the pseudo binary system is the polymer molecule, treated as the solute molecule. The other component, termed here pseudosolvent, is made up of water, the surfactant, and the electrolyte. The thermodynamics of this solute-pseudosolvent binary system is quantitatively described using the formalism of any suitable polymer solution theory available in the literature. The criterion of phase separation stipulated by the polymer solution theory is then applied to the pseudo binary system in order to find the composition boundary of the pseudo-components at which the solution undergoes phase separation. This composition boundary is then reinterpreted in terms of the actual concentrations of the surfactant and the electrolyte constituting the pseudosolvent.

In Section 2, the experimental methods and the materials used, and the observed phase behavior are briefly described for some polymer and surfactant systems of interest to the enhanced oil recovery application. In Section 3, the conceptual outline of the pseudosolvent model is presented. The pseudosolvent model is quantitatively developed in Section 4 and the estimation of model parameters is detailed. Phase diagrams are computed and compared to the experimentally measured ones in Section 5. Most importantly, for ease of practical applications, a simple approach to constructing phase diagrams using one experimentally measured surfactant-electrolyte composition lying at the phase boundary is proposed, without requiring any information on the polymer. The last section summarizes the principal conclusions from this work.

2. Experimental Study of Phase Behavior

2.1. Materials

The water-soluble polymers used in this study are listed in Table 1. The polyacrylamide and the biopolymer have been explored both in laboratory and field applications for enhanced oil recovery [2,7,8,9,10,11,12]. Flocon is an aqueous solution of the microbially produced heteropolysaccharide xanthan gum, with a molecular weight of about 2 x 10⁶. Polyethylene oxide (PEO) has been used in several studies of polymer-surfactant interactions and is a well-defined polymer for modeling studies [10,12].

Commercially available low molecular weight anionic surfactants as well as deoiled petroleum sulfonates were employed in this study. Table 2 lists the surfactants, their chemical structures, the names of suppliers.

The petroleum sulfonates supplied by Witco Chemical (TRS 40 and TRS 10-410) included a distribution of molecules as well as some amount of oil. These petroleum sulfonates were deoiled and fractionated into two parts, one with a low average equivalent weight (LEW) and the other with a high average equivalent weight (HEW). The HEW samples were used for the phase behavior measurements [13]. Other surfactants were used as received. For the petroleum sulfonates, the chemical formula is given assuming the presence of a single monosulfonate having the measured equivalent weight.

The chemical structures of the polymer repeat units and the surfactants used in this study are presented in Figure 1.

2.2. Phase Behavior Measurements

Aqueous solutions of polymers with specified amounts of salt were prepared using either a magnetic stirrer or a propeller-type mixer. Care was taken to ensure that maximum dispersion of polymer molecules was achieved without simultaneously causing any mechanical degradation. The viscosity of the polymer solution was measured prior to phase behavior studies to make sure that no polymer degradation had occurred. The solution of biopolymer (Flocon) was prefiltered through a Whatman No. 3 filter paper to remove any impurities and bacteria from the Flocon material before the phase behavior studies [13].

Desired amounts of surfactants were weighed in glass containers (bottles) to which previously prepared polymer-electrolyte aqueous solution was added. The sample bottles were mixed with a magnetic stirrer. The solutions were kept at room temperature (24 ± 2°C) for phase behavior observations. The concentration boundaries which separate the one-phase region from the two-phase region for several polymer-surfactant systems were established based on two weeks of solution equilibration.

Figures 2 shows the phase boundaries for different surfactants in solutions of PEO (4M), Pusher 700 (5M), and Flocon (2M). The points represent the experimentally observed compositions at which phase separation occurred. Note that the salt and surfactant concentrations are expressed as weight percent, that is weight (in g) of the component in 100 g of the solution. These can be converted to molar concentrations assuming that the total volume of the solution is practically that of water, given that the concentrations of the surfactant and the salt are not very large.

Figure 2 shows that for a given polymer molecule, the phase stability of the polymer-surfactant solution decreases (namely, phase separation occurs) with increasing equivalent weight of the surfactant, increasing salt concentration, and increasing surfactant concentration. For a given surfactant molecule, the phase stabilities for the polymers decrease in the following order: PEO (4M) ≥ Pusher 700 ≥ Flocon, but the difference among these high molecular weight polymers is not very significant.

3. Development of Pseudosolvent Model for Predicting Phase Behavior

3.1. Pseudosolvent and its Characteristic Energy ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$

An evaluation of polymer solution theories in the literature makes it apparent that an adequate predictive theory for polymer-water system is presently not available [15,16]. The classical Flory-Huggins theory, given the original meaning of the Flory interaction parameter ${\mathsf{\chi }}_{12}$ [14], cannot predict the characteristic lower critical solution temperature (LCST) behavior exhibited by many aqueous polymer solutions, including the PEO-water system. While using the classical Flory theory, usually, the interaction parameter is empirically fitted as a function of the temperature and the volume fraction of the polymer in solution to describe (not predict) the experimentally observed phase behavior of polymer-water systems [15,16,17]. The development of models for aqueous polymer solutions with true predictive capabilities is an area of continuing interest [18,19]. When the aqueous polymer solution also includes surfactant and electrolyte molecules, a formal molecular level modeling of the four-component system becomes even more complex. Since the primary focus of this work is on predicting the phase behavior of these four component systems, a simplified approach tailored to treat phase behavior, is proposed here.

In this approach, the pure solvent (water) is replaced by a pseudosolvent made up of water, electrolyte, and surfactant (present as both singly dispersed surfactant molecules and micellar aggregates). The polymer molecule is the solute. In this manner, the four-component system is treated in terms of an equivalent pseudo two-component system. We then need a polymer solution model as a framework to describe this pseudobinary system and for this purpose we choose the lattice fluid model of polymer solutions developed by Sanchez and Lacombe [20,21,22,23,24]. This model is chosen because it allows representation of the LCST behavior and is thus capable of at least qualitatively or semi-quantitatively describe aqueous polymer solutions. It should be noted that the LCST behavior in lattice fluid model emerges from the consideration of the presence of vacant sites in the lattice and the resulting entropic effects. This is in contrast to the LCST behavior in aqueous polymer solutions that arises from hydrogen bonding interactions. Indeed the lattice fluid model has been extended to account for hydrogen bonding [25], but for our purposes of just needing a polymer solution model framework, the basic lattice fluid model is adequate. The lattice fluid model for pure fluids and polymer solutions is briefly summarized in the Appendix with all key equations used in this study.

In the lattice fluid model, a pure component is defined by three equation of state parameters, a characteristic energy ${\mathsf{\epsilon }}^{*}$, a characteristic size of the component r and a characteristic lattice site of hard-core volume ${\mathrm{v}}^{*}$. One can alternately use the characteristic pressure P*, characteristic temperature T* and the characteristic density ${\mathsf{\rho }}^{*}$ as the defining parameters (see Equations (A-1) to (A-4) in the Appendix). For small molecules, the equation of state parameters for the pure components are usually estimated from experimental saturation vapor pressure data. The characteristic parameters for water determined in this way [20] are, ${\mathrm{P}}^{*}$= 26,520 atm, ${\mathrm{T}}^{*}$= 623^oK and ${\mathsf{\rho }}^{*}$ = 1.105 g/cm³, or equivalently, ${\mathsf{\epsilon }}^{*}$= 1.238 kcal/mol, ${\mathrm{v}}^{*}$= 1.927 cm³/mol and r = M/(${\mathsf{\rho }}^{*}$${\mathrm{v}}^{*}$) = 8.45, the molecular weight of water M being 18. For high molecular weight polymers, vapor pressures are negligibly small and therefore, the equation of state parameters are usually determined from experimental density data. Alternately, when only limited density data are available, the equation of state parameters can be determined from a single experimental density ρ, thermal expansion coefficient α and isothermal compressibility β, at the same temperature and pressure, using Equations (A-6) and (A-7) in the Appendix. For polyethylene oxide of molecular weight M, the equation of state parameters estimated by this single point method [13] are, ${\mathrm{P}}^{*}$= 3,963 atm, ${\mathrm{T}}^{*}$= 597^oK and ${\mathsf{\rho }}^{*}$= 1.202 g/cm³, or equivalently, ${\mathsf{\epsilon }}^{*}$= 1.194 kcal/mol, ${\mathrm{v}}^{*}$= 12.36 cm³/mol and r = M/(${\mathsf{\rho }}^{*}$${\mathrm{v}}^{*}$). Note that the characteristic size parameter r is obviously dependent on the molecular weight M of the polymer.

For the pseudosolvent composed of water, the surfactant, and the electrolyte as the actual components, we assume that two of its equation of state parameters, ${\mathrm{v}}_1^{*}$ and ${\mathrm{r}}_1$, (subscript 1 denotes the solvent and 2 denotes the polymer) are identical to those of water. The third parameter, characteristic energy ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ (the subscript 1P is introduced to distinguish it from 1 used for water) is taken to be different from that of water and is assumed to account for all the surfactant-electrolyte composition-dependent characteristics of the pseudosolvent. Thus the critical postulate underlying our pseudosolvent model is the assumption that a single characteristic energy parameter, ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ contains all the information related to the type and amounts of surfactant and electrolyte molecules present in water. It should be mentioned that this is an ad hoc extra-thermodynamic postulate and is justified only a posteriori by the resulting simplicity of the thermodynamic approach and the usefulness of the results.

3.2. Polymer – Pseudosolvent Binary Parameters ξ and δ

For the binary polymer-pseudosolvent solution two additional parameters characteristic of the mixture properties are introduced, as is common in many theories of solutions to account for non-idealities of mixture behavior. The binary energy parameter ξ represents the deviation of the interaction energy from the geometric mean rule and the binary volume parameter δ represents the deviation of the close packed mixture volume from the arithmetic additivity rule. The introduction of binary parameters ξ and δ adds considerable flexibility to the lattice fluid model for mixtures, as it does in all solution models in the literature. The binary parameters for the polymer-pseudosolvent are assumed to be the same as for polymer-water and determined from experimental thermodynamic data on polymer-water solutions. The estimates for the binary parameters obtained from different solution properties such as activity of the solvent, heat of mixing or volume change on mixing are generally not identical, reflecting the fundamental inadequacies in the current state of the theory of the liquid state. Since we are interested in the polymer solution phase behavior, we use the water activity data to estimate the binary polymer-water parameters. As shown in Figure 3, the experimental activity data [26] for a PEO-water system is fitted well by the theoretical water activity calculated from Equation (A-13) in the Appendix, for the binary parameter values ξ = 1.023 and δ = - 0.25. These binary parameters are taken to be independent of temperature.

3.3 Polymer – Pseudosolvent Phase Diagram in the Parametric Space of ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$

Given the equation of state parameters ${\mathsf{\epsilon }}_2^{*}$, ${\mathrm{v}}_2^{*}$, ${\mathrm{r}}_2$ for the polymer, the two equation of state parameters, ${\mathrm{v}}_1^{*}$, ${\mathrm{r}}_1$ for the pseudosolvent (assumed to be the same as for water) and the binary parameters, ξ, δ for the polymer-pseudosolvent (assumed to be the same as for polymer-water), one can calculate the phase behavior of the polymer solution for various pseudosolvents characterized by differing values for ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$. The phase stability criteria are summarized in Equations (A-20) to (A-23) in the Appendix. Since the equation of state parameters as well as the binary polymer-water binary interaction parameters are known for polyethylene oxide, we were able to construct the polymer-pseudosolvent phase diagram shown in Figure 4 for three different molecular weights of PEO (14 K, 300 K and 4 M) in the ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ vs polymer concentration space. The phase diagram shows that while the characteristic parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the pseudosolvent is below a certain value, the polymer-pseudosolvent solution is stable; and when ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ is above a certain value, the solution is unstable. The characteristic energy of the pseudosolvent ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ at which phase separation occurs is a function of the polymer molecular weight, the concentration of the polymer in solution as well as the temperature. Similar phase diagrams can be constructed for other polymers of interest, knowing their equation of state parameters and their binary interaction parameters with water.

Figure 4 shows that for high molecular weights polymer, PEO (4M) and PEO (300K), the value of ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ at the phase boundary is almost constant in the polymer concentration range of 100 ppm to 1500 ppm (equivalently 0.01 wt% to 0.15 wt%). In contrast, for the low molecular weight polymer PEO (14K), the value of ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ is significantly higher than that for the high molecular weight polymer, especially at low polymer concentrations.

3.4. Mapping Characteristic Energy ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the Pseudosolvent to its Composition

The characteristic energy ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the pseudosolvent needs to be translated to the actual composition of the pseudosolvent (namely the chemical structures and concentrations of the surfactant and the electrolyte). This is done by considering the free energy difference ΔG between the pseudosolvent and water. In the formalism of the lattice fluid model, given the equation of state parameters for a pure fluid, the free energy G of a pure fluid can be calculated using Equation (A-1) in the Appendix by introducing in it the fluid density calculated from the equation of state, Equation (A-4). Consequently, the free energy difference ΔG between the pseudosolvent and water is calculated as follows:

$$\Delta \mathrm{G}\hspace{0.17em}=\hspace{0.17em}\mathrm{G}\hspace{0.17em}{(\mathsf{\epsilon }}_{1\mathrm{P}}^{*},\hspace{0.17em}{\mathrm{v}}_1^{*},\hspace{0.17em}{\mathrm{r}}_1)\hspace{0.17em}-\hspace{0.17em}\mathrm{G}\hspace{0.17em}{(\mathsf{\epsilon }}_1^{*},\hspace{0.17em}{\mathrm{v}}_1^{*},\hspace{0.17em}{\mathrm{r}}_1)$$

(1)

The calculated relation between ΔG and ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$is presented in Figure 5. Note that this relation is independent of the polymer and is thus a universal relation for pseudosolvents in the framework of the lattice fluid model and can be used to explore the phase behavior involving any polymer.

The characteristic energy parameter for water ${\mathsf{\epsilon }}_1^{*}$ is 1.238 at which ΔG = 0. For ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ < ${\mathsf{\epsilon }}_1^{*}$, ΔG > 0 and for ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ > ${\mathsf{\epsilon }}_1^{*}$, ΔG < 0. The characteristic energy parameter is a measure of the cohesive energy of the solvent. The higher the energy parameter, the larger the cohesive energy for the solvent, greater the probability of polymer exclusion and therefore the larger the likelihood of phase separation.

4. Calculation of ΔG from Thermodynamic Data on Surfactants and Electrolytes

The thermodynamic quantity of interest for predicting the phase behavior of polymer-surfactant-electrolyte-water system, in the framework of the pseudosolvent concept, is the free energy difference ΔG between the pseudosolvent and pure water. We need to develop a quantitative relation between ΔG (which is taken to represent through its magnitude all the essential characteristics of the pseudosolvent) and the actual chemical composition of the pseudosolvent. The free energy difference ΔG can be related to the actual chemical composition of the pseudosolvent using thermodynamic models and data (entirely independent of the polymers) available in the literature for surfactant and electrolyte solutions, as detailed in this section. ΔG can be calculated for different surfactant and electrolytes as:

$$\Delta \mathrm{G}\hspace{0.17em}=\hspace{0.17em}\mathrm{G}\hspace{0.17em}(\mathrm{water}\hspace{0.17em}+\hspace{0.17em}\mathrm{surfactant}\hspace{0.17em}+\hspace{0.17em}\mathrm{electrolyte})\hspace{0.17em}-\hspace{0.17em}\mathrm{G}\hspace{0.17em}\left(\mathrm{water}\right)$$

(2)

In this manner, the actual composition of the pseudosolvent can be related to the free energy difference ΔG and thereby to the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$of the pseudosolvent.

As mentioned earlier, the pseudosolvent consists of surfactant molecules, any electrolyte added and water. The surfactant molecules are present both in their singly dispersed state as well as in the form of micellar aggregates. The pseudosolvent can be thought as a mixture of water (component 1) in its pure water reference state with the other components, electrolyte counterion (component 2), electrolyte coion (component 3), surfactant counterion (component 4) and surfactant coion (component 5) in their infinitely dilute reference states. Since components 2 and 3 come from the salt and components 4 and 5 come from the surfactant, their mole fractions x or molar concentrations C are related as x₂ = x₃ and x₄ = x₅, or equivalently, C₂ = C₃ and C₄ = C₅. Expressions for estimating the magnitude of ΔG are developed in this section.

4.1. Water + Electrolyte Systems

For the aqueous electrolyte solutions, we adopt the treatment of Kawaguchi et al [27] formulated based on the Analytic Solution of Group (ASOG) model [28] to calculate the activity of the components. In this treatment, the electrolyte solution is considered to consist of free water molecules (W), hydrated cations (C) and hydrated anions (A). The electrolyte is assumed to be dissociated completely to produce the cations and the anions with which water molecules are bound. The structure of the hydrated ion is determined by the number of hydration water molecules around it. A hydration water molecule is distinguished from a free water molecule because the water of hydration is tightly bound such that its motion and intramolecular states are different from those of the free water. The hydrated cation (n_C water molecules and a cation) and the hydrated anion (n_A water molecules and an anion) are treated as the two solute components present in the solvent water. Using this visualization, the activities of the components of the electrolyte solution can be viewed as composed of (i) the ideal entropy of mixing, (ii) an excess entropic contribution arising from the presence of the hydrated anions, the hydrated cations, and the free water molecules, described using the Flory-Huggins expression [14] and (iii) an enthalpic contribution arising from the interactions between the free water and the hydration water molecules that surround the ions. The long-range ion-ion electrostatic interactions were found to be relatively small and therefore neglected in the present study. Therefore, the activity of component i can be written as

$${\mathrm{a}}_{\mathrm{i}} \hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{\mathrm{i}}\hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}\hspace{0.17em}\hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{\mathrm{i}}\hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{FH}}\hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{G}}$$

(3)

where ${\mathrm{x}}_{\mathrm{i}}$ is the mole fraction of component i (i refers to free water, hydrated cation and the hydrated anion), ${\mathsf{\gamma }}_{\mathrm{i}}$ is the activity coefficient and ${\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{FH}}$ and ${\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{G}}$are the excess entropic and enthalpic contributions due to the non-idealities.

The expression for the excess entopic contribution to the activity coefficient is given by the ASOG model [27] as follows:

$$\ln \hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{FH}}\hspace{0.17em}=\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\nu }}_{\mathrm{i}}/{\displaystyle \sum _{\mathrm{j}}{\mathrm{x}}_{\mathrm{j}}^{\mathrm{h}}{\mathsf{\nu }}_{\mathrm{j}})\hspace{0.17em}+\hspace{0.17em}1\hspace{0.17em}-\hspace{0.17em}}{(\mathsf{\nu }}_{\mathrm{i}}/{\displaystyle \sum _{\mathrm{j}}{\mathrm{x}}_{\mathrm{j}}^{\mathrm{h}}{\mathsf{\nu }}_{\mathrm{j}})}$$

(4)

In Equation (4), ${\mathrm{x}}_{\mathrm{i}}^{\mathrm{h}}$ is the mole fraction of component i in the mixture, with the superscript h denoting that the water of hydration is treated as an integral part of the species. ${\mathsf{\nu }}_{\mathrm{i}}$ is the number of atoms (other than hydrogen) in the component i and equals 1, ( ${\mathrm{n}}_{\mathrm{C}}$+ 1), and (${\mathrm{n}}_{\mathrm{A}}$+ 1), for free water, hydrated cation and hydrated anion respectively. Accounting for the hydrated water surrounding the ions, the mole fractions of free water (${\mathrm{x}}_{\mathrm{W}}^{\mathrm{h}}$), hydrated electrolyte cation (${\mathrm{x}}_{\mathrm{C}}^{\mathrm{h}}$) and hydrated electrolyte anion (${\mathrm{x}}_{\mathrm{A}}^{\mathrm{h}}$) are given by:

$$\begin{array}{c}{\mathrm{x}}_{\mathrm{W}}^{\mathrm{h}}\hspace{0.17em}=\hspace{0.17em}\frac{{\mathrm{x}}_{\mathrm{W}}\hspace{0.17em}-\hspace{0.17em}{\mathrm{n}}_{\mathrm{C}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}-\hspace{0.17em}{\mathrm{n}}_{\mathrm{A}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}}{{\mathrm{x}}_{\mathrm{W}}\hspace{0.17em}-\hspace{0.17em}\hspace{0.17em}{\mathrm{n}}_{\mathrm{C}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}-\hspace{0.17em}{\mathrm{n}}_{\mathrm{A}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}},{\mathrm{x}}_{\mathrm{C}}^{\mathrm{h}}\hspace{0.17em}=\hspace{0.17em}\frac{\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}}{{\mathrm{x}}_{\mathrm{W}}\hspace{0.17em}-\hspace{0.17em}\hspace{0.17em}{\mathrm{n}}_{\mathrm{C}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}-\hspace{0.17em}{\mathrm{n}}_{\mathrm{A}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}},\\ {\mathrm{x}}_{\mathrm{A}}^{\mathrm{h}}\hspace{0.17em}=\hspace{0.17em}\frac{\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}}{{\mathrm{x}}_{\mathrm{W}}\hspace{0.17em}-\hspace{0.17em}\hspace{0.17em}{\mathrm{n}}_{\mathrm{C}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}-\hspace{0.17em}{\mathrm{n}}_{\mathrm{A}}\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{C}}\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_{\mathrm{A}}}\end{array}$$

(5)

where ${\mathrm{x}}_{\mathrm{W}}$, ${\mathrm{x}}_{\mathrm{C}}$ and ${\mathrm{x}}_{\mathrm{A}}$ are the actual mole fractions of the water, the electrolyte cation, and the electrolyte anion in solution. If we consider only a water - electrolyte solution, in the notations introduced for all components of the pseudosolvent, ${\mathrm{x}}_{\mathrm{W}}$ = ${\mathrm{x}}_1$, ${\mathrm{x}}_{\mathrm{C}}$ = ${\mathrm{x}}_2$, and ${\mathrm{x}}_{\mathrm{A}}$ = ${\mathrm{x}}_3$. Equation (5) refers to a single cation and a single anion, but if there are multiple distinct ions, they can be accounted through additional terms for each distinct ion.

The enthalpic contribution to the activity coefficient in the ASOG model [27] is expressed as:

$$\begin{array}{l}\ln \hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{G}}\hspace{0.17em}=\hspace{0.17em}{\displaystyle \sum _{\mathrm{k}}{\mathsf{\nu }}_{\mathrm{ki}}\hspace{0.17em}(}\ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{k}}-\hspace{0.17em}\ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{k}}^{\left(\mathrm{i}\right)})\\ \ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{k}}\hspace{0.17em}=\hspace{0.17em}-\ln \hspace{0.17em}({\displaystyle \sum _{\mathrm{e}}{\mathrm{q}}_{\mathrm{e}}^{\mathrm{h}}\hspace{0.17em}{\mathrm{a}}_{\mathrm{k}/\mathrm{e}})\hspace{0.17em}+\hspace{0.17em}1\hspace{0.17em}-\hspace{0.17em}}{\displaystyle \sum _{\mathrm{e}}{(\mathrm{q}}_{\mathrm{e}}^{\mathrm{h}}\hspace{0.17em}{\mathrm{a}}_{\mathrm{e}/\mathrm{k}}/{\displaystyle \sum _{\mathrm{m}}{\mathrm{q}}_{\mathrm{m}}^{\mathrm{h}}\hspace{0.17em}{\mathrm{a}}_{\mathrm{e}/\mathrm{m}})}}\\ {\mathrm{q}}_{\mathrm{k}}\hspace{0.17em}=\hspace{0.17em}{\displaystyle \sum _{\mathrm{k}}{\mathrm{x}}_{\mathrm{i}}^{\mathrm{h}}\hspace{0.17em}{\mathsf{\nu }}_{\mathrm{ki}}\hspace{0.17em}}/{\displaystyle \sum _{\mathrm{i}}{\mathrm{x}}_{\mathrm{i}}^{\mathrm{h}}\hspace{0.17em}{\displaystyle \sum _{\mathrm{k}}{\mathsf{\nu }}_{\mathrm{ki}}}\hspace{0.17em}}\end{array}$$

(6)

where ${\mathsf{\nu }}_{\mathrm{ki}}$ is the number of interacting groups of kind k in component i, ${\mathsf{\Gamma }}_{\mathrm{k}}$ is the activity coefficient of the interacting group k in the mixture, ${\mathsf{\Gamma }}_{\mathrm{k}}^{\left(\mathrm{i}\right)}$ is the standard state activity coefficient of group k. and ${\mathrm{q}}_{\mathrm{k}}$ is the group fraction of group k. The expression for the activity coefficient for group k, ${\mathsf{\Gamma }}_{\mathrm{k}}$, follows the Wilson model [29].

There are only two interacting groups in this system, namely free water (denoted by subscript k = OH) and the hydrated water (k = OH*) since the interaction between water molecules and ions is entirely accounted for, solely through this hydration. Each component has only one kind of group in it, k = OH in free water and k = OH* in hydrated cation and in hydrated anion. The group activity coefficients are given [27] by

$$\begin{array}{l}\ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{OH}}\hspace{0.17em}=\hspace{0.17em}- \ln \hspace{0.17em}{(\mathrm{q}}_{\mathrm{OH}}{+\mathrm{q}}_{\mathrm{OH}* }{\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*})+1 - (\frac{{\mathrm{q}}_{\mathrm{OH}}}{{\mathrm{q}}_{\mathrm{OH}}{+\mathrm{q}}_{\mathrm{OH}*}{ \mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}}+\frac{{\mathrm{q}}_{\mathrm{OH}* }{\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}}{{\mathrm{q}}_{\mathrm{OH} }{\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}{+\mathrm{q}}_{\mathrm{OH}*}})\\ \ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{OH}*}\hspace{0.17em}=\hspace{0.17em}- \ln \hspace{0.17em}{(\mathrm{q}}_{\mathrm{OH}*}{+\mathrm{q}}_{\mathrm{OH} }{\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}})+1 - (\frac{{\mathrm{q}}_{\mathrm{OH}*}}{{\mathrm{q}}_{\mathrm{OH}*}{+\mathrm{q}}_{\mathrm{OH}}{ \mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}}+\frac{{\mathrm{q}}_{\mathrm{OH} }{\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}}{{\mathrm{q}}_{\mathrm{OH}* }{\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}{+\mathrm{q}}_{\mathrm{OH}}})\end{array}$$

(7)

In Equation (7), ${\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}$ and ${\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}$ are the Wilson interaction energy parameters and ${\mathrm{q}}_{\mathrm{k}}$ is the group mole fraction of group k. For water, the standard state is pure water and for the anion and the cation, the standard state is the infinite dilution condition. Therefore from Equation (7),

$$\ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{OH}}^{\left(\mathrm{i}\right)}\hspace{0.17em}=\hspace{0.17em}0 ,\ln \hspace{0.17em}{\mathsf{\Gamma }}_{\mathrm{OH}*}^{\left(\mathrm{i}\right)}\hspace{0.17em}\hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}-\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}\hspace{0.17em}+1-\hspace{0.17em}\hspace{0.17em}{\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}\hspace{0.17em}$$

(8)

Kawaguchi [27] has reported the hydration numbers of 11 ions and the interaction energy parameters (${\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}$ and ${\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}$) between the free water denoted by OH and the hydration water denoted by OH*. These parameter estimates were obtained by correlating the activities of water in electrolyte solutions very accurately, taking ${\mathsf{\nu }}_{\mathrm{ki}}$ in Equation (6) as 1.6 for free water, ${\mathrm{n}}_{\mathrm{C}}$ for the hydrated cation and ${\mathrm{n}}_{\mathrm{A}}$ for the hydrated anion [27,30]. For 14 electrolyte solutions with 1:1 and 2:1 electrolytes, the activities of water calculated by ASOG model using these optimized parameters were found to deviate on the average by less than 1% from the experimental data [31] , for salt concentrations of up to 5 molality [27]. Therefore, we directly make use of the parameter estimates reported in Kawaguchi’s work. The interaction energy parameters ${\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}*}$ and ${\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}}$ where OH denotes a free water and OH* a hydration water were determined to be 1.82 and 1.78, respectively [27,30]. Note that the interaction parameters between identical groups will be unity, that is ${\mathrm{a}}_{\mathrm{OH}/\mathrm{OH}}$ = ${\mathrm{a}}_{\mathrm{OH}*/\mathrm{OH}*}$ = 1. The values of the hydration numbers ${\mathrm{n}}_{\mathrm{C}}$ and ${\mathrm{n}}_{\mathrm{A}}$ optimized by Kawaguchi by fitting the water activity data are presented in Table 3. The value of ${\mathrm{n}}_{\mathrm{C}}$ decreases with the increase of ionic radius. This tendency qualitatively coincides with that anticipated from the literature [32].

In the Kawaguchi treatment of electrolyte solutions, the hydrated anion and cation are taken as the distinct chemical species to develop expressions for the activity coefficients. But in the pseudosolvent model, we treat the electrolyte ions in the unhydrated state as the components. Therefore, while writing an expression for the free energy difference ΔG between the pseudosolvent composed of water + electrolyte and water, we should introduce a correction term ΔG_ref to account for this difference in the reference states. On this basis, the free energy difference ΔG between the pseudosolvent composed of water + electrolyte and pure water can be written as

$$\begin{array}{l}\Delta \mathrm{G}=\mathrm{RT}\hspace{0.17em}\hspace{0.17em}[{\mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_1+\hspace{0.17em}{\mathrm{x}}_2\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_2+\hspace{0.17em}{\mathrm{x}}_3\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_3\hspace{0.17em}]\hspace{0.17em}\hspace{0.17em}+\\ \mathrm{RT}\hspace{0.17em}\hspace{0.17em}[{\mathrm{x}}_1\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_1^{\mathrm{FH}}{\mathsf{\gamma }}_1^{\mathrm{G}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_2\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_2^{\mathrm{FH}}{\mathsf{\gamma }}_2^{\mathrm{G}})+{\mathrm{x}}_3\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_3^{\mathrm{FH}}{\mathsf{\gamma }}_3^{\mathrm{G}})\hspace{0.17em}]+{\mathrm{x}}_2{ \Delta \mathrm{G}}_{\mathrm{ref}}\end{array}$$

(9)

The correction term ΔG_ref accounts for the change in energy when unhydrated ions become hydrated with associated water molecules. This process involves breaking of water-water hydrogen bonds and replacement by ion-dipole interactions between the ion and water. The corresponding energy changes cannot be unambiguously estimated and we have assumed a value of 5000 cal/mole based on typical magnitudes for hydrogen bonding energies (which can range between 2 to 10 kcal/mole) and ion-dipole interactions reported in the literature. We present computed results that show the sensitivity of the calculated free energy difference ΔG between the pseudosolvent and water, to variations in this assumed value. All other parameters for calculating ΔG in Equation (9) have been established and validated by Kawaguchi [27] based on accurate fitting of water activity data for numerous electrolyte solutions over a wide range of concentrations.

4.2. Water + Surfactant Systems

The surfactant solution is composed of water (component 1), surfactant counterion (component 4) and surfactant coion (component 5). The formation of micelles in the solution can be formally represented by the law of mass action as

$$\mathrm{m}\hspace{0.17em}{\mathrm{S}}^{-}\hspace{0.17em}+\hspace{0.17em}\mathrm{n}\hspace{0.17em}{\mathrm{A}}^{+}\hspace{0.17em}\hspace{0.17em}\rightleftarrows \hspace{0.17em}\hspace{0.17em}{\left(\mathrm{SA}\right)}^{-\hspace{0.17em}(\mathrm{m}\hspace{0.17em}-\hspace{0.17em}\mathrm{n})}$$

(10)

where S^- represents the free surfactant ion (representing an anionic surfactant and can be similarly written for cationic surfactants), A⁺ represents the free counterion, m is the aggregation number of the micelle, n is the number of counterions bound to the micelle surface. Below the critical micelle concentration (CMC), the surfactant molecules are present mostly in the singly dispersed form, and above the CMC the added surfactant molecules appear as micellar aggregates (denoted by the subscript mic). The singly dispersed surfactant is completely dissociated into surfactant counterion (denoted as 4S) and surfactant coion (denoted as 5S). The activity of the surfactant in solution is equal to the activity of the singly dispersed surfactant because of the monomer-micelle equilibrium.

The free energy difference ΔG between the surfactant solution (pseudosolvent) and pure water can be written similar to Equation (9), in the form

$$\begin{array}{l}\Delta \mathrm{G}=\hspace{0.17em}{\mathrm{RT} [ \mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}({\mathrm{x}}_1)+\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{4\hspace{0.17em}}\ln \hspace{0.17em}({\mathrm{x}}_{4\mathrm{S}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_5\hspace{0.17em}\ln \hspace{0.17em}({\mathrm{x}}_{5\mathrm{S}}{)+\mathrm{x}}_{\mathrm{mic}}\hspace{0.17em}\ln \hspace{0.17em}({\mathrm{x}}_{\mathrm{mic}}\left)\right]\hspace{0.17em}\\ +\hspace{0.17em}{\mathrm{RT} [ \mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\gamma }}_1^{\mathrm{S}})+\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{4\hspace{0.17em}}\ln \hspace{0.17em}{(\mathsf{\gamma }}_{4\mathrm{S}}^{\mathrm{S}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_5\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\gamma }}_{5\mathrm{S}}^{\mathrm{S}}{)+\mathrm{x}}_{\mathrm{mic}}\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\gamma }}_{\mathrm{mic}}^{\mathrm{S}}\left)\right]\end{array}$$

(11)

where the superscript S is added to indicated that these activity coefficients appear in surfactant solutions. In Equation (11), ${\mathrm{x}}_{\mathrm{mic}}$ and ${\mathsf{\gamma }}_{\mathrm{mic}}^{\mathrm{S}}$ refer to the mole fraction and activity coefficient for micelles. The monomer-micelle equilibrium, Equation (10), implies that the micelle activity can be related to the monomer activity in the form

$$\ln \hspace{0.17em}{(\mathrm{x}}_{\mathrm{mic}}{ \mathsf{\gamma }}_{\mathrm{mic}}^{\mathrm{S}})=\mathrm{m} \ln \hspace{0.17em}{(\mathrm{x}}_{5\mathrm{S}}{ \mathsf{\gamma }}_{5\mathrm{S}}^{\mathrm{S}})+\mathrm{n} \ln \hspace{0.17em}{(\mathrm{x}}_{4\mathrm{S}}{ \mathsf{\gamma }}_{4\mathrm{S}}^{\mathrm{S}})\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}$$

(12)

Further, the concentration of singly dispersed surfactant and micelles can be related to the total concentration via the mass balance

$${\mathrm{x}}_4\hspace{0.17em}={\mathrm{x}}_{4\mathrm{S}}{+\mathrm{n} \mathrm{x}}_{\mathrm{mic}} ,{\mathrm{x}}_5\hspace{0.17em}={\mathrm{x}}_{5\mathrm{S}}{+\mathrm{m} \mathrm{x}}_{\mathrm{mic}}\hspace{0.17em}$$

(13)

Introducing Equations (12) and (13) into Equation (11), we get

$$\begin{array}{l}\Delta \mathrm{G}=\hspace{0.17em}{\mathrm{RT} [ \mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}({\mathrm{x}}_1)+\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{4\hspace{0.17em}}\ln \hspace{0.17em}({\mathrm{x}}_{4\mathrm{S}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_5\hspace{0.17em}\ln \hspace{0.17em}({\mathrm{x}}_{5\mathrm{S}}) ]\hspace{0.17em}\\ +\hspace{0.17em}{\mathrm{RT} [ \mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\gamma }}_1^{\mathrm{S}})+\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{4\hspace{0.17em}}\ln \hspace{0.17em}{(\mathsf{\gamma }}_{4\mathrm{S}}^{\mathrm{S}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_5\hspace{0.17em}\ln \hspace{0.17em}{(\mathsf{\gamma }}_{5\mathrm{S}}^{\mathrm{S}}) ]\end{array}$$

(14)

The amount of the free surfactant ions$\hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}$ can be estimated from the critical micelle concentration, CMC. The mole fraction of surfactant in the solution near the CMC is generally quite small and hence the total moles in 1 liter of solution can be approximated as 55.5 moles based on water alone. Therefore, with the CMC expressed as a molar concentration (mole/liter), we can estimate

$${\mathrm{x}}_{5\mathrm{S}}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}\mathrm{CMC}/55.5$$

(15)

The CMC of homologous straight-chain ionic surfactants in aqueous medium, in the absence of any added salt (C₂ = 0), displays a dependence on the number of carbon atoms N in the hydrophobic chain in the form

$${\ln \mathrm{CMC}|}_{\mathrm{C}2=0}=\hspace{0.17em}\hspace{0.17em}\mathrm{A} - \mathrm{B} \mathrm{N}$$

(16)

where A and B are constants [33].

The number of counterions binding on the micelles depends on the concentration of counterions in the solution. The binding of counterions on the micelles is assumed to follow the simple Langmuir adsorption isotherm. Correspondingly, the degree of dissociation of the micelle, denoted by α, can be written as

$$\begin{array}{l}\mathsf{\alpha }=\frac{\mathrm{m} - \mathrm{n}}{\mathrm{m}}=\frac{{\mathsf{\alpha }}^{*}}{{1+\mathrm{K}}_{\mathrm{ad}} \mathrm{I}},\\ \mathrm{I}=\frac{{\displaystyle \sum {\mathrm{C}}_{\mathrm{i}}{\mathrm{Z}}_{\mathrm{i}}^2}}{2}=\frac{{(\mathrm{C}}_{4\mathrm{S}}{+\mathrm{C}}_{5\mathrm{S}}{+(\mathsf{\delta } \mathsf{\alpha } \mathrm{m})}^2{ (\mathrm{C}}_5{-\mathrm{C}}_{5\mathrm{S}})/\mathrm{m} )}{2}\hspace{0.17em}\\ \mathrm{I}=\hspace{0.17em}\hspace{0.17em}\frac{55.5 {(\mathrm{x}}_{4\mathrm{S} }{+\mathrm{x}}_{5\mathrm{S}}{+(\mathsf{\delta } \mathsf{\alpha } \mathrm{m})}^2{ (\mathrm{x}}_5{-\mathrm{x}}_{5\mathrm{S}})/\mathrm{m})}{2}\end{array}$$

(17)

In equation (17), α^* is the degree of dissociation of micelles in infinitely dilute solution, $\hspace{0.17em}{\mathrm{K}}_{\mathrm{ad}}$ is the Langmuir adsorption equilibrium constant (liter/mole) for counterion binding to micelles, I is the ionic strength of the solution, ${\mathrm{x}}_{4\mathrm{S}},\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}\hspace{0.17em}$ are the mole fractions and ${\mathrm{C}}_{4\mathrm{S}},\hspace{0.17em}\hspace{0.17em}{\mathrm{C}}_{5\mathrm{S}}\hspace{0.17em}$ are the molar concentrations of the counterion and coion of the free surfactant, respectively and ${(\mathrm{C}}_5-\hspace{0.17em}{\mathrm{C}}_{5\mathrm{S}})/\mathrm{m}\hspace{0.17em}$ is the molar concentration of micelles of aggregation number m. Below the CMC, ${\mathrm{x}}_{4\mathrm{S}}=\hspace{0.17em}{\mathrm{x}}_4,\hspace{0.17em}\hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}{=\mathrm{x}}_5$. Above the CMC, ${\mathrm{x}}_{5\mathrm{S}}\hspace{0.17em}$= CMC/55.5 and ${\mathrm{x}}_{4\mathrm{S}}\hspace{0.17em}=\hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}{+\mathsf{\alpha } (\mathrm{x}}_5{ - \mathrm{x}}_{5\mathrm{S}})$.

In the absence of any added electrolyte, the ionic strength is determined by the concentration of the singly dispersed surfactant which is dissociated into surfactant coion and surfactant counterion (the first two terms in the expression for I) and the micelles which contain (αm) multiple charges (the last term in the expression for I). It is generally assumed that the charges on the micelles are partially shielded and the shielding factor δ < 1 is introduced to reduce the effect of the micellar charge on the ionic strength of the solution. It has been found that using δ = 1 is not in agreement with the experimentally determined values of activity and osmotic coefficients. For several ionic surfactants δ is estimated to be around 0.5 [34] and that value is assumed in this work.

The activity coefficient of the surfactant coion and the surfactant counterion are difficult to determine independently and the experimental measurements usually focus on an average activity coefficient defined as ${\mathsf{\gamma }}_{\pm S}^S\hspace{0.17em}={{(\mathsf{\gamma }}_{4S}^S\hspace{0.17em}{\mathsf{\gamma }}_{5S}^S)}^{1/2}$[34,35]. There is very little data in the literature on experimentally determined average activity coefficients for surfactants. The average activity coefficient of the surfactant in the presence of ionic species is attributed mainly to the "salting out" or "salting in" of the hydrophobic group of the surfactant in the aqueous solvent. For nonpolar solutes in aqueous electrolyte solutions of ionic strength I, McDevit and Long [36] developed a theoretical expression for the activity coefficient in the form, $\ln \mathsf{\gamma } \hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}\hspace{0.17em}\mathrm{I}$. This linear dependence is identical to the well-known empirical Setschenow relation [37], with$\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}$ designated as the Setschenow constant. Although the equation of McDevit and Long was developed for nonpolar solutes, it has been applied to various polar and polar organic compounds [38] and nonionic surfactants [39] as well. To retain simplicity, we will assume that this relation can be applied to calculate the average activity coefficient of the surfactant.

$${\ln \mathsf{\gamma }}_{\text{±}\mathrm{S}}^{\mathrm{S}} \hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}\hspace{0.17em}\mathrm{I}\hspace{0.17em}$$

(18)

The activity coefficient of water$\hspace{0.17em}{\mathsf{\gamma }}_1^S$ can then be derived through the Gibbs-Duhem equation. When the molar concentration of the surfactant$\hspace{0.17em}{\mathrm{C}}_5$ is less than the CMC, no micellization occurs. When $\hspace{0.17em}{\mathrm{C}}_5$ is greater than the CMC, $\hspace{0.17em}{\mathsf{\gamma }}_1^S$ includes contributions arising from the presence of micelles. One thus gets the surfactant contribution to the activity coefficient of water to be

$$\begin{array}{l}{\ln \mathsf{\gamma }}_1^{\mathrm{S}} \hspace{0.17em}=\hspace{0.17em}-\hspace{0.17em}\frac{{\mathrm{k}}_{\mathrm{S}}\hspace{0.17em}{\mathrm{C}}_5^2}{{\mathrm{C}}_1}{ \mathrm{when} \mathrm{C}}_5 \text{<} \mathrm{CMC}\\ {\ln \mathsf{\gamma }}_1^{\mathrm{S}} \hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}-\frac{{\mathrm{k}}_{\mathrm{S}}\hspace{0.17em}{\mathrm{CMC}}^2}{{\mathrm{C}}_1}-\frac{{\mathrm{k}}_{\mathrm{S}} \mathsf{\alpha }}{{2\mathrm{C}}_1}{(\mathrm{C}}_5^2{- \mathrm{CMC}}^2)+\frac{1}{{\mathrm{C}}_1}{(\mathrm{C}}_5-\mathrm{CMC})+\\ \hspace{0.17em}\frac{(1-\mathsf{\alpha }) \mathrm{CMC}}{{\mathrm{C}}_1 \mathsf{\alpha }}\ln (\frac{{\mathrm{C}}_5 \mathsf{\alpha }+(1-\mathsf{\alpha }) \mathrm{CMC}}{\mathrm{CMC}})\hspace{0.17em}{\mathrm{when} \mathrm{C}}_5 \text{>} \mathrm{CMC}\end{array}$$

(19)

The free energy difference ΔG between the water + surfactant solution (pseudosolvent) and pure water can now be calculated from Equation (14) by introducing Equations (15) to (19) in it.

4.3. Water + Surfactant + Electrolyte Systems

For water-surfactant-electrolyte solutions, the expression for the free energy difference ΔG between the pseudosolvent and pure water can be written by extending Equation (9) for the free energy difference of the water + electrolyte system and Equation (14) for the free energy difference of the water + surfactant system. On this basis, we can write

$$\begin{array}{l}\Delta \mathrm{G}=\mathrm{RT}\hspace{0.17em}\hspace{0.17em}{[\mathrm{x}}_1\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_1+\hspace{0.17em}{\mathrm{x}}_2\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_2+\hspace{0.17em}{\mathrm{x}}_3\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_3\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_{4\hspace{0.17em}}\ln \hspace{0.17em}{\mathrm{x}}_{4\mathrm{S}}\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_5\hspace{0.17em}\ln \hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}]\hspace{0.17em}\\ +\mathrm{RT}\hspace{0.17em}\hspace{0.17em}{[\mathrm{x}}_1\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_1^{\mathrm{FH}}{\mathsf{\gamma }}_1^{\mathrm{G}}{\mathsf{\gamma }}_1^{\mathrm{S}})\hspace{0.17em}+\hspace{0.17em}{\mathrm{x}}_2\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_2^{\mathrm{FH}}{\mathsf{\gamma }}_2^{\mathrm{G}}{)+\mathrm{x}}_3\hspace{0.17em}\ln (\hspace{0.17em}{\mathsf{\gamma }}_3^{\mathrm{FH}}{\mathsf{\gamma }}_3^{\mathrm{G}}{)+(\mathrm{x}}_{4\hspace{0.17em}}{+\mathrm{x}}_5) \ln \hspace{0.17em}{\mathsf{\gamma }}_{\text{±}\mathrm{S}}^{\mathrm{S}}\hspace{0.17em}]\\ +\hspace{0.17em}{\mathrm{x}}_2\hspace{0.17em}{\Delta \mathrm{G}}_{\mathrm{ref}}\end{array}$$

(20)

As mentioned earlier, subscripts 2 and 3 refer to the counterion and the coion of the electrolyte. Subscripts 4 and 5 refer to the counterion and the coion of the surfactant. The$\hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{FH}}$contribution is calculated from Equation (4), the$\hspace{0.17em}{\mathsf{\gamma }}_{\mathrm{i}}^{\mathrm{G}}$contribution from equation (6), the$\hspace{0.17em}{\mathsf{\gamma }}_{\text{±}\mathrm{S}}^{\mathrm{S}}$ contribution from Equation (18) and the $\hspace{0.17em}{\mathsf{\gamma }}_1^{\mathrm{S}}$ contribution from Equation (19). The ionic strength I now includes the added electrolyte.

$$\begin{array}{l}\mathrm{I}=\frac{{(\mathrm{C}}_2{+\mathrm{C}}_3{+\mathrm{C}}_{4\mathrm{S}}{+\mathrm{C}}_{5\mathrm{S}}{+(\mathsf{\delta } \mathsf{\alpha } \mathrm{m})}^2{ (\mathrm{C}}_5{-\mathrm{C}}_{5\mathrm{S}})/\mathrm{m} )}{2}\hspace{0.17em}\\ \mathrm{I}=\hspace{0.17em}\hspace{0.17em}\frac{55.5 {(\mathrm{x}}_2{+\mathrm{x}}_3+\hspace{0.17em}{\mathrm{x}}_{4\mathrm{S} }{+\mathrm{x}}_{5\mathrm{S}}{+(\mathsf{\delta } \mathsf{\alpha } \mathrm{m})}^2{ (\mathrm{x}}_5{-\mathrm{x}}_{5\mathrm{S}})/\mathrm{m})}{2}\end{array}$$

(21)

The presence of electrolyte in aqueous solution causes a decrease in the CMC, which will influence the concentrations$\hspace{0.17em}{\mathrm{x}}_{4\mathrm{S}}$and$\hspace{0.17em}{\mathrm{x}}_{5\mathrm{S}}$. The depression of CMC in ionic surfactant solutions is due mainly to the decrease in the thickness of the ionic double layer surrounding the ionic head groups in the presence of the additional electrolyte and the consequent decrease in the electrostatic repulsions between them at the micelle surface. Experimental data [40] suggest that for ionic surfactants, the effect of the concentration of electrolyte on the CMC is given by,

$${\ln \mathrm{CMC}|}_{\mathrm{C}2} - {\ln \mathrm{CMC}|}_{\mathrm{C}2=0}\hspace{0.17em}=\hspace{0.17em}\hspace{0.17em}-{\mathrm{a} [\ln ( \mathrm{C}}_2+{\mathrm{CMC}|}_{\mathrm{C}2=0}) - \ln ({\mathrm{CMC}|}_{\mathrm{C}2=0})]$$

(22)

where a is a constant for a given ionic head group at a particular temperature and is independent of the hydrophobic tail length, and$\hspace{0.17em}{\mathrm{C}}_2$ is the concentration of the added electrolyte in mole per liter. Introducing Equation (16) into equation (22) we can calculate the CMC for a surfactant with N carbon atoms in the hydrophobic tail at an added electrolyte concentration of$\hspace{0.17em}{\mathrm{C}}_2$ from,

$${\ln \mathrm{CMC}|}_{\mathrm{C}2}=(1+\mathrm{a}) (\mathrm{A}-\mathrm{BN})\hspace{0.17em}-{\mathrm{a} \ln ( \mathrm{C}}_2+\mathrm{A}-\mathrm{BN})$$

(23)

4.4. Estimation of Thermodynamic Parameters from Literature

For calculating the free energy difference ΔG between the pseudosolvent (composed of water-surfactant-electrolyte) and the pure water using Equation (20), the values of the following thermodynamic parameters for the surfactant-electrolyte solutions are needed.

(a)

A and B in Equation (16) for the hydrophobic tail length dependence of the CMC

(b)

a in Equation (22) for the ionic strength dependence of the CMC

(c)

$\hspace{0.17em}{\mathrm{K}}_{\mathrm{ad}}$ and α* in Equation (17) for counterion binding at micellar surface

(d)

$\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}$ in Equation (18) for the activity coefficient of surfactant.

The parameters A and B in Equation (16) for calculating the critical micelle concentration as a function of surfactant chain length, in the absence of any added salt, are 3.45 and 0.69 for alkyl sulfates, and 3.68 and 0.67 for alkyl sulfonates and alkylbenzene sulfonates [40]. The parameter a in Equation (22) to calculate the effect of salt on the CMC for sodium alkyl sulfate and sulfonates is 0.458 [40, 41]. All these parameter values are based on the CMC expressed in molar concentration (moles/liter). The CMC of surfactants with different hydrocarbon chain lengths and at various added salt concentrations can be readily calculated from Equation (23) incorporating the values for A, B, b, N and$\hspace{0.17em}{\mathrm{C}}_2$.

The micelle aggregation number m appearing in Equations (17) and (23) is taken to be 50 for all surfactants considered here to avoid adding additional parameters. Changing the value of m does not significantly affect the computed free energy differences or the phase diagrams. The adsorption constant$\hspace{0.17em}{\mathrm{K}}_{\mathrm{ad}}$ in Equation (17) for the counter-ion adsorption at micellar surface is taken as 1.0 (mole/liter)^-1 and α* in Equation (17) is taken to be 0.5. Both the values of$\hspace{0.17em}{\mathrm{K}}_{\mathrm{ad}}$ and α* are chosen to agree with the experimental observation that the fraction of counter-ions dissociation on the micelle surface for alkyl sulfonates and alkyl sulfates is in the range of 0.3 to 0.6 for the surfactant and electrolyte concentrations studied [33,41,42].

For alkyl sulfonate and alkyl sulfate surfactants, experimental average activity coefficient data are not available for the homologous series of surfactants. Experimental data for short chain alkyl carboxylates [35] suggest that the mean molar activity coefficients for different carboxylate homologs approach a constant value of -0.25 at their critical micellar concentrations. If this behavior is displayed by the alkyl sulfates and sulfonates, the Setschenow constant$\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}$ for differing hydrophobic chain lengths can be estimated from the following equation,

$${\ln \mathsf{\gamma }}_{\text{±}\mathrm{S}}^{\mathrm{S}}{ (\mathrm{at} \mathrm{CMC})=\mathrm{k}}_{\mathrm{S}} \mathrm{CMC}=\mathrm{constant}$$

(24)

For the present study, the value of constant in Equation (24) for alkyl sulfate, alkyl sulfonate and alkyl benzene sulfonate surfactants is chosen as -0.5 since it is in the range of the experimental average activity coefficient for sodium dodecyl sulfate [34,42,43]. Incorporating this constant value of -0.5 and the cmc calculated as described above, the Setschenow constant k_S is determined using Equation (24). The validity of Equation (24) for other surfactants (in other words, whether$\hspace{0.17em}{\mathsf{\gamma }}_{\text{±}\mathrm{S}}^{\mathrm{S}}$ at CMC is equal to a constant value for surfactant homologs) needs experimental verification. However, at the present time, due to lack of experimental activity coefficient data for various surfactant homologs, the validity of Equation (24) is assumed and a constant value of -0.5 is assigned to calculate the values of$\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}$ for alkyl sulfate, alkyl sulfonate and alkyl benzene sulfonate surfactant homologs.

Since the parameters $\hspace{0.17em}{\mathrm{K}}_{\mathrm{ad}}$, α* and $\hspace{0.17em}{\mathrm{k}}_{\mathrm{S}}$ have been assigned values that are reasonable, but without any direct supporting experimental data for the specific molecular systems, we perform parametric sensitivity analysis in Section 5.4 on how changes in these parameter values affect the calculation of the free energy difference ΔG, and thereby the phase diagram.

5. Phase Behavior Predictions and Construction of Phase Diagrams

5.1. Calculation of the Free Energy Difference Between the Pseudosolvent and Water

The free energy difference ΔG defined in Equation (20) has been calculated for sodium pentyl sulfonate (R₅SO₃Na), sodium decyl sulfate (R₁₀SO₄Na), sodium dodecyl sulfate (R₁₂SO₄Na), sodium tetradecyl sulfate (R₁₄SO₄Na), TRS 40 HEW (RB₁₇SO₃Na), TRS 10-410 HEW (RB₂₄SO₃Na), and TRS 18 HEW (RB₃₄SO₃Na) at different salt and surfactant concentrations using the parameter values provided above. For illustrative purposes, we have plotted in Figure 6, the calculated values of ΔG for the pseudosolvent containing sodium pentyl sulfonate, sodium dodecyl sulfate and TRS 40 as surfactants and sodium chloride as the electrolyte, in the concentration range of interest. Similar plots have been constructed for each of the surfactant mentioned above and, in all cases, using NaCl as the electrolyte.

Figures 5 discussed earlier in Section 3.3 shows the relation between the characteristic energy parameter, ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the pseudosolvent, and the free energy difference ΔG between the pseudosolvent and water, calculated from Equation (A-1), in the framework of the lattice fluid theory. Figure 5 shows that as ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ increases, ΔG decreases. The implication of this for the calculated results shown on Figure 6 is that the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the pseudosolvent increases as the concentration of surfactant and/or salt is increased. Since the increase in ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ results in phase separation in polymer solutions, we have the general conclusion that an increase in the concentration of surfactant and/or salt will lead to phase separation in aqueous solutions containing polymers. For any concentration of salt and surfactant, the value of ΔG can be found from Figure 6, then the corresponding value of ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ can be read off from Figure 5, and thus one can obtain the equation of state parameters (${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$, ${\mathrm{v}}_1^{*}$, r) of the pseudosolvent.

5.2. Construction of Ternary Phase Diagram from Theory Based on iso-ΔG values

To construct the ternary phase diagram for any polymer - pseudosolvent system, we start from Figure 4 constructed for the polymer of interest (in this case, PEO). For a given polymer molecular weight and concentration, one can find the value of the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ when phase separation would occur from Figure 4. The corresponding free energy difference ΔG is then determined from Figure 5, which is a universal relation independent of the polymer. All water-surfactant-electrolyte compositions with free energy differing from water by the same magnitude of ΔG will lie at the polymer-pseudosolvent phase boundary, allowing us to construct the phase diagram in the actual compositional space of polymer + surfactant + electrolyte + water system. These iso-ΔG surfactant + electrolyte + water compositions are found from Figure 6.

Ternary phase diagrams can thus be constructed in this manner by combining the results for polymer + water system such as in Figure 4 obtained from the lattice fluid theory, the universal relation between the characteristic energy of the pseudosolvent and the free energy difference obtained in the framework of the lattice fluid model shown in Figure 5 and the relation between and the actual composition of the solution obtained from the polymer independent models of surfactant + electrolyte systems, as shown on Figure 6. The ternary phase diagram constructed in the manner are shown on Figure 7 for polyethylene oxide + sodium dodecyl sulfate + NaCl + water system at two concentrations of the electrolyte, NaCl.

One can observe the influence of polymer molecular weight, polymer concentration, surfactant concentration, and electrolyte concentration from the ternary diagram. Both Figures 7a and 7b show that increasing surfactant concentration leads to phase separation. Comparing the two figures, we see that increasing the salt concentration reduces the single phase compositional domain. Similarly, increasing the equivalent weight of the surfactant (not shown in this figure) also decreases the single phase compositional domain.

In addition to predicting the effect of salt concentration, surfactant concentration, and the equivalent weight of surfactant on the phase behavior of the polymer-surfactant solutions in agreement with known experimental behavior, the pseudosolvent model also correctly describes the effect of polymer molecular weight and polymer concentration. Both Figures 7a and 7b show that for high molecular weight PEO samples, the phase behavior is not changed appreciably by the change in polymer concentration while for the lower molecular weight PEO, there is a change with polymer concentration. One can observe that as the polymer concentration is changed, the surfactant concentration at the phase boundary remains practically unchanged in both figures at the given salt concentrations. Both figures show that the compositional space of the single phase domain is not much altered by the polymer molecular weight, as long as the polymers have high molecular weights.

Both Trushenski et al. [7,8] and Pope et al. [10] found that the influence of polymer concentration (in the range 100 to 1500 ppm) and polymer molecular weight (for molecular weights larger than 400K)on the phase behavior of the polymer-surfactant-electrolyte systems were negligible. Kalpakci [6] made the same observation and further noted that as the high molecular weight polymer was degraded using an ultrasonic mixer and/or orifice mixer, the phase stability of the system improved and the aqueous solution remained stable at higher salt and surfactant concentrations.

Our calculated results in Figure 7 are in agreement with these experimental observations. For the high molecular weight polymers, PEO (4M) and PEO (300K), for the polymer concentration in the range of 100 ppm to 1500 ppm (equivalently 0.01 wt% to 0.15 wt%), the phase behavior is the same, namely the phase boundaries have the same concentration of the surfactant and electrolyte. In contrast Figure 7 shows that for the lower molecular weight polymer, PEO (14K), the higher the concentrations of surfactant and salt, at which phase separation would occur. This model prediction is in agreement with the experimental findings of Kalpacki that decreasing polymer molecular weight during ultrasonication resulted in enhanced phase stability at higher surfactant and salt concentration.

5.4. Construction of Phase Diagram Using a Single Experimental Data

In the absence of polymer equation of state parameters or the binary parameters of the polymer-water system, it will not be possible to calculate the phase diagram similar to Figure 4. In the absence of knowledge of ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ values, it is not possible to construct the ternary phase diagram such as in Figure 7. However, for these situations, we suggest a simple approach to constructing the phase diagram similar to Figure 2. If we know from experiments one composition lying at the phase boundary (any one experimental point in Figure 2), we can calculate the corresponding ΔG of the surfactant + electrolyte + water system from Figure 6 (which can be constructed for that surfactant + electrolyte system). Various combinations of surfactant + electrolyte systems can give rise to identical magnitudes of ΔG (or to identical values of the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$) are expected to have similar phase behavior in the framework of the pseudosolvent model. All surfactant + electrolyte + water systems having iso-ΔG values can be identified from Figure 6 in order to construct the phase diagram without requiring any information about the polymer. Illustrative phase diagrams calculated in this manner, using one experimental phase composition data lying at the phase boundary from Figure 2 are shown on Figure 8 for the three polymers considered in this work.

For example, one of the solution compositions on the phase boundary for the Pusher 700-R₁₂SO₄Na system contains 0.76 wt % surfactant and 2.0 wt % NaCl. It can be found from Figure 6 that the free energy difference ΔG for this pseudosolvent is -33.6 cal/mol. This implies, in the framework of the pseudosolvent model, that all compositions of the pseudosolvent with the free energy difference of -33.6 cal/mol, which can be found from Figure 6, must lie on the phase boundary for the Pusher 700 - surfactant- electrolyte system. The iso-ΔG also occurs at 0% NaCl and 2.82 wt% surfactant and 4% NaCl and 0.35 wt% surfactant for sodium dodecyl sulfate, providing other compositions lying at the phase boundary. This procedure has been used for the three polymer systems, to establish their predicted phase boundaries as shown by the continuous lines in Figure 8. The predicted and experimental phase boundaries (shown as points) are found to be in reasonable agreement.

The most significant feature of the pseudosolvent model is thus the possibility of constructing phase diagrams using knowledge of just one experimental phase composition lying at the phase boundary, without requiring any information about the polymer. Using this one experimental point available for any one surfactant and/or electrolyte, one can theoretically establish the phase boundary for the polymer with any other surfactant and electrolyte systems.

5.4. Assessment of Parametric Sensitivity to Predicted Results

To carry out predictive computations we have used estimated parameter values where direct experimental data for specific molecules are not available. The choice of values for the degree of dissociation on micelle surface at infinite dilution α* and the shielding parameter δ are reasonable based on available data for multiple surfactants. In contrast, the counterion adsorption equilibrium constant K_ad, the Setschenow constant to describe salt effect on surfactant activity k_S and the correction term ΔG_ref to account for the consideration of hydrated ions as the actual species in place of the unhydrated ions have all been assigned values based on limited or no direct experimental data. Therefore, we have performed illustrative calculations to assess the sensitivity of the free energy calculations to variations in these parameter values. The results are shown on Figure 9.

Changes in the values of K_ad and k_S only change the free energy values to a small extent and therefore the phase diagrams will change only marginally. The correction term ΔG_ref does not affect the dependence of the free energy curve on the surfactant concentration but displaces the curves up or down, that is, to smaller or larger ΔG magnitudes. This will have the effect of shifting the phase boundary curves in Figure 8 up or down. For quantitative comparisons of predictions against experiments, it would be possible to update the present calculations whenever the currently unavailable surfactant parameter values are determined experimentally

The free energy calculations show that the predominant contributions come from the ideal free energy of mixing and the non-ideal free energy contribution arising from the salt effect on the surfactant. All other terms are practically unimportant for the concentrations of salt and surfactant considered in this study, even though they have been rigorously included here. Therefore improved estimation of the Setschenow constant k_S will assure that true predictions can be obtained.

6. Conclusions

Phase separation in four-component systems made up of polymer + surfactant + electrolyte + water, is treated by introducing the concept of a pseudosolvent. This concept transforms the problem to that of a pseudo-binary solution and thus avoids the unmanageability of the thermodynamic description of phase stability criteria in four-component systems. The pseudosolvent is taken to be made up of water, electrolyte, singly dispersed surfactant molecules, and micellar aggregates. The polymer molecule is regarded as the solute. The lattice fluid theory is then used to describe the phase behavior of the polymer + pseudosolvent system. All important characteristics of the pseudosolvent are assumed to be described by a single parameter, namely the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$. This is an extra thermodynamic postulate underlying the pseudosolvent treatment. This characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ of the pseudosolvent is related to the free energy difference ΔG between the pseudosolvent and water. In order to calculate the free energy difference ΔG between the pseudosolvent and water, electrolyte solution theory and surfactant solution theory available in the literature have been used.

The free energy calculations show that the main contributions to the free energy difference ΔG between the pseudosolvent and water come from the ideal free energy of mixing and the non-ideal free energy contribution arising from the salt effect on the surfactant. All other terms are practically unimportant for the concentrations of salt and surfactant considered in this study. These terms which have been rigorously modeled in this work could become important at higher concentrations of electrolytes when studies focus on somewhat smaller polymer molecular weights.

Ternary phase diagrams have been constructed by combining the value of the characteristic energy parameter ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ when phase separation would occur from Figure 4, the polymer-independent universal relation between ${\mathsf{\epsilon }}_{1\mathrm{P}}^{*}$ and the free energy difference ΔG between the pseudosolvent and water shown on Figure 5, and the dependence of ΔG on the surfactant-electrolyte composition shown on Figure 6. The ternary phase diagrams correctly show agreement with the observations in the literature that the phase behavior of polymer-surfactant solutions show only slight dependence on molecular weight of the polymer for high molecular weight polymers, they are not very sensitive to polymer concentrations for high molecular weight polymers at low polymer concentration range, the phase stability increases significantly as the polymer molecular weight is drastically decreased and the compositional domain of the single phase region decreases with increasing polymer molecular weight, increasing surfactant equivalent weight, increasing surfactant concentration and increasing added salt concentration.

More interestingly, we present an approach predicting the phase behavior without any information about the polymer, as long as a single experimental composition data for a system lying on the phase boundary is known. Using the appropriate parameter values in the pseudosolvent model, the free energy difference ΔG between the pseudosolvent and water is calculated. All other chemical compositions required to complete the phase boundary are then determined from surfactant-electrolyte compositions having iso- ΔG values. The comparison between the phase boundaries calculated in this manner from the pseudosolvent model and the experimental phase boundaries show satisfactory agreement.