An Analytical Understanding of Reentrant Condensation of a Polyelectrolyte in the Presence of an Oppositely Charged Surfactant
Huaisong Yong, Holger Merlitz
TL;DR
This paper explains how oppositely charged surfactants affect the condensation and phase transitions of polyelectrolytes, offering insights into liquid-liquid phase separation in biological systems.
Contribution
A novel mean-field theory with analytical solutions explains reentrant condensation and phase transitions in polyelectrolytes with surfactants.
Findings
Strong electrostatic adsorption between ionic monomers and surfactant ions is critical for collapse and re-entry transitions.
A minimum coupling energy is required for phase transitions, explaining the role of surfactant chain length.
The theory rationalizes the 'egg shape'-like phase diagram and applies to biological systems involving DNA/RNA and proteins/peptides.
Abstract
We explore the phase-transition mechanism of the reentrant condensation of a polyelectrolyte in the presence of an oppositely charged surfactant, which is of fundamental importance to the understanding of liquid–liquid phase separation (LLPS) in soft materials and biological systems. We focus on the adsorption and attraction effects of surfactants near/on polymer chains and ignore their own nonessential mixing effects if surfactant molecules are far away from polymer chains. This novel approach allows us to construct a simple mean-field equilibrium theory with closed-form analytical solutions, which can rationalize the essential features of the emergent “egg shape”-like phase diagram. The theory addresses that a strong electrostatic adsorption between the ionic monomers and surfactant ions is critical to understand the peculiar phenomenon that both the collapse and re-entry transitions…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
1
2
3
4
5
6| symbols | physical meaning of the symbol |
|---|---|
|
| the number of monomers in a polyelectrolyte chain |
|
| the surfactant chain size in units of a polyelectrolyte monomer |
|
| the size of a monomer along the direction of the polymer backbone, the size of a solvent molecule, the size of a surfactant counterion, and the size of a polyelectrolyte counterion |
|
| the fraction of charged monomers in a polyelectrolyte chain |
| ∏ | the osmotic pressure acting on the polymer phase by the bulk |
|
| the volume fraction of monomers |
|
| the volume fraction of counterion of a surfactant |
|
| the Boltzmann constant |
|
| the thermodynamic temperature |
| φ | the fraction of ionic monomers occupied preferentially by the surfactant ions |
| λ ≈ 1 | the volume ratio between the ionic head of a surfactant ion and a solvent (or a counterion) |
| χs | the parameter considering the demixing effect between a surfactant and water on polyelectrolyte chains |
| μ | the exchange chemical potential for a surfactant ion on the polymer chains |
| the volume fraction of the surfactant ion next to polyelectrolyte chains in the solvent phase | |
|
| the Bjerrum length of water |
| κ | the inverse Debye screening length |
| ε1 | the excess-adsorption strength of one surfactant ion with respect to the monovalent counterion of the polyelectrolyte on its ionic monomers |
| ε2 | the hydrophobic-aggregation strength among surfactant ions that are adsorbed on ionic monomers |
| εFH,1 | the Flory–Huggins parameter between the solvent and charged monomers |
| εFH,2 | the Flory–Huggins parameter between the solvent and uncharged monomers |
|
| the free energy owing to the mixing of the polymer with solvent and added surfactant |
|
| the adsorption free energy owing to the mixing of surfactant ions and polyelectrolyte counterions on the polymer chains |
|
| the associative attraction between ionized monomers caused by a temporary bridge due to the hydrophobic-aggregation effect of adsorbed surfactant ions on ionic monomers |
|
| the free energy of nonassociative pairwise-like electrostatic interaction due to the |
|
| the energy of nonelectrostatic excluded-volume interactions between the solvent (water) and the charge-neutral part of monomers |
| the total Gibbs free energy
per volume unit under the framework of the isothermal–isobaric ( |
- —Deutsche Forschungsgemeinschaft10.13039/501100001659
- —Southwest Petroleum University10.13039/501100004753
- —Tianfu Emei Scholar Foundation of Sichuan ProvinceNA
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsElectrostatics and Colloid Interactions · Surfactants and Colloidal Systems · Electrowetting and Microfluidic Technologies
Introduction
1
Phase transitions of polyelectrolyte solutions are pivotal to the stability and function of soft materials. The predictability of phase separations such as condensation and aggregation in polyelectrolyte solutions is of fundamental importance to various fields, from materials science ?−? ? to cell biology ?,? and art preservation. ?,? A particular scenario of polyelectrolyte phase separations is the reentrant condensation of a polyelectrolyte in the diluted aqueous solution of an oppositely charged surfactant, as sketched in Figurea, which has attracted significant attention in the past decades. ?−? ? ? ? ? ? ? ? A remarkable feature of the reentrant condensation is that both the collapse and the reentry branches occur at rather low surfactant concentrations, ?−? ? ? where the overall concentration of the added surfactant (such as hexadecyl trimethylammonium bromides) is usually on the order of 5 mmol/L and thus below the bulk critical micelle concentration (CMC) of the surfactant. It is worth mentioning that this reentrant condensation cannot be well-understood merely according to classical polyelectrolyte theories, ?−? ? which predict the absence of any reentrant signature for polyelectrolyte phase transitions when the concentration of added small monovalent salts is on the order of 5 mmol/L. Therefore, new concepts and theoretical formalisms are required to understand this observed peculiar reentrant behavior at low-surfactant concentrations.
*(a) A sketch of the reentrant condensation of polyelectrolytes induced by a diluted oppositely charged surfactant. Typical values of surfactant concentrations for collapse and reentry transitions are on the order of 1 mmol/L and below the CMC of the surfactant. −
Note that the reentrant condensation is not necessarily symmetric with respect to surfactant concentrations. (b) A sketch of the well-studied phase-transition mechanism of electrostatic-adsorption-induced hydrophobic aggregation ,−
for the collapse branch of the reentrant condensation as illustrated in panel (a). The ionic monomers (symbol “⊕” in the figure) are preferentially adsorbed by oppositely charged surfactant ions (symbol “⊖” in the figure), and polyelectrolyte chains further form temporary bridges due to the hydrophobic-aggregation effect of adsorbed surfactant ions (black lines which connect in the figure). The pink lines represent polymer chains, and we do not show counterions of both polyelectrolytes and surfactants for simplification.*
Another unusual feature of the reentrant condensation is that the electrostatic binding between the charged monomer and the oppositely charged surfactant leads to the formation of insoluble polymer–surfactant coacervates once a certain surfactant concentration, known as the critical aggregation concentration (CAC), is exceeded. The CAC can be several orders of magnitude below the CMC of the surfactant,? which implies that the influence of the bulk micellar behaviors of the surfactant is negligible for the reentrant condensation. We note that in the literature, ?,? a surfactant concentration at which the insoluble polymer–surfactant coacervates start to redissolve (i.e., the onset of the reentry branch of the polyelectrolyte condensation in Figurea) is usually termed the second critical aggregation concentration (CAC_2_).
Clarifying the function of oppositely charged surfactants in such a reentrant condensation (as sketched in Figurea) is critical for a deep understanding of biological phase separations if proteins/peptides bound to DNAs/RNAs play an important role, ?−? ? ? ? ? ? ? ? for example, for a better understanding of the DNA length on the co-condensation behaviors of proteins with DNAs. ?,? Furthermore, from the aspect of applied research, the elucidation of the phase-transition mechanism of the reentrant condensation is pivotal for developing new polymer formulations such as for drug delivery? and new soft materials such as porous polymer materials.?
A well-studied phase-transition mechanism ?,?−? ? ? for the collapse branch of the reentrant condensation (as sketched in Figurea) is that the surfactant ion replaces the polyelectrolyte counterions and preferentially adsorbs on the ionic monomers as well as forming electrostatic dipoles. Furthermore, the long hydrophobic tails of the surfactants can undergo hydrophobic aggregation (Figureb). Thus, when the surfactant concentration reaches a critical value, a charge-adsorption-induced premicellar aggregation occurs, which ultimately leads to the formation of insoluble polymer–surfactant coacervates. In contrast, the phase-transition mechanism for the reentry branch, as sketched in Figurea, can be significantly complicated by specific chemistry-related details of the uncharged moieties of the polyelectrolyte. It has been reported that the hydrophobicity of uncharged monomers of the polyelectrolyte affects and decreases the surfactant concentration that is necessary to trigger the reentry transition. ?,?
Related to the nonlinear attraction effect among ionic monomers deriving from the electrostatic adsorption of the surfactant ?,?−? ? ? (as illustrated by Figureb), theoretical rationalizations in the past decades such as chemical reaction-like theories ?,? and lattice/field-like theories ?,?,?,?−? ? ? ? ? ? ? ? ? ? ? have significantly promoted the understanding of phase behaviors of polyelectrolytes in the presence of oppositely charged surfactants, such as recent interesting studies on the complicated micellar structures ?−? ? ? and the nonequilibrium/metastable/kinetically trapped states ?,?,? of polymer–surfactant coacervates. However, there remain two fundamental problems? to be addressed clearly for the reentrant condensation as sketched in Figurea: the first problem is why polyelectrolytes show phase transition only if the surfactant chain length is longer than a minimum length. ?−? ? ? The second problem is why both the collapse and re-entry transitions of polyelectrolytes can occur when the concentration of the surfactant is lower than its bulk critical micelle concentration (CMC).
The primary goal of this work is to try to resolve the above two problems by exploiting the well-studied phase-transition mechanism, as sketched in Figureb. Compared with other published works, ?,?,?−? ? ? ? ? ? ? ? ? ? ? ? ? the novelty of the present work lies in the selection of the contributions to the free energy function. Experiments have shown that a reentrant condensation of a polyelectrolyte can take place when the surfactant concentration is several orders of magnitude below its bulk critical micelle concentration (CMC).? Thus, although the surfactant’s own entropic mixing term is included in the construction of the free energy function in our mean-field equilibrium theory, it would significantly increase the complexity of the subsequent analytical description without actually contributing to the reentrant phenomenon that we want to describe. We therefore omitted this term in the subsequent analytical calculation. Instead, we focus on an explicit contribution from the hydrophobic interaction among the tails of the adsorbed surfactants, which turns out to play a pivotal role in the properties of the resulting phase transition. This novel approach allows us to successfully solve our theory analytically with closed-form solutions and to rationalize the essential features such as the emergent “egg shape” of spinodal and binodal phase diagrams for the reentrant condensation, as illustrated by Figurea.
In the rest of this article, a mean-field theory for the polyelectrolyte solution in terms of the free energy is constructed in Section. Its analytical solution will be considered in detail in Section, where we outline some general results of mean-field theory. Furthermore, a simplified phase diagram of the polyelectrolyte solution will be discussed in this section, as well as the applicability of the mean-field theory. Final concluding remarks are given in Section.
Methods and Theory
2
In this section, we construct the mean-field theory as shown by eq for the polyelectrolyte solution in the presence of an oppositely charged surfactant based upon the Gibbs free energy and clarify its physics foundation. The theory is unfolded as follows.
As sketched in Figureb, we consider flexible polyelectrolytes with monovalent ionic monomers and monovalent counterions in an aqueous solution. We denote N as the degree of polymerization of the polyelectrolyte chain and a as the extent of each monomer along the direction of the polymer backbone. The charged monomers are distributed randomly on the polyelectrolyte chain, and their fraction is denoted by p. In general, the size of an ionized monomer, the size of a counterion of polyelectrolyte, the size of a counterion of the surfactant, and the size of a solvent molecule can be very different. Since we are interested in the general physical understanding of the theory, we restrict ourselves here to the symmetric case, i.e., we consider that these sizes are identical in the spirit of the classical Flory–Huggins lattice model, ?,? which simplifies the analytical arguments substantially. The overall volume fraction of monovalent ionic monomers and neutral monomers is denoted by c, so that the volume fraction of monovalent counterions is pc. We denote the volume fraction of the counterions of the surfactant as c _ x _, and then the volume fraction of the surfactant ion is nc _ x _, where n denotes the surfactant chain length in units of a polyelectrolyte monomer.
Here and in the following, let us consider the free energy per unit of volume for an incompressible system if not otherwise noted. The volume unit is given by the size of the solvent molecules in the spirit of the classical Flory–Huggins lattice model. ?,? We separate the system into two parts: the polymer chains with their enclosed solvents and small ions and the bulk without polymer chains. The contribution to the free energy from the three-dimensional mixing of polymers with the solvent and the added surfactant is given by G sol, see the first line of eq. The bulk is just described by an osmotic pressure (∏) acting on the polymer phase by the bulk. If the volume of the polymer coils changes, then mechanical work against the external pressures is involved in the free energy change. This will become important when we minimize the free energy with respect to the monomer concentration (c). This approach is thus based on the framework of the isothermal–isobaric (NPT) ensemble. Here and in the following, we consider energies in units of k B T, where k B is the Boltzmann constant and T is the thermodynamic temperature.
The energy of nonelectrostatic excluded-volume interactions between the solvent (water) and the charge-neutral part of monomers is given by the classical Flory–Huggins formalism ?,? and denoted as G FH, see the second line of eq. The Flory–Huggins parameters between the solvent and charged monomers are denoted by ε_FH,1_, and ε_FH,2_ are the Flory–Huggins parameters between the solvent and uncharged monomers. In experiments, these Flory–Huggins parameters can be directly estimated by using Hansen solubility parameters? according to molecular structures of both monomers and solvents.
The isothermal adsorption free energy per unit of volume owing to the one-dimensional mixing of surfactant ions and monovalent polyelectrolyte counterions on ionic monomers is given by G ads, see the third line of eq. The fraction of ionic monomers occupied preferentially by the large surfactant ions is denoted by φ. The excess-adsorption strength of one surfactant ion compared to a monovalent counterion of the polyelectrolyte on its ionic monomers is denoted by ε_1_. The exchange chemical potential of a surfactant ion on the polymer chains is denoted by μ, which scales as μ ∼ ln(nc _ x ) if the surfactant concentration (n + 1)c _ x _ in the bulk is sufficiently small. Because there is a strong demixing tendency between the surfactant and solvent molecules, we take care of this demixing effect on polyelectrolyte chains by the parameter χ_s, which can be as large as 2.0 for alkyl trimethylammonium bromides.? We denote the volume ratio between the ionic head of a surfactant and a solvent (or a counterion) by λ, which is usually on the order of unity, though it can be varied significantly in some specific experiments. ?,?
As sketched in Figureb, we consider the associative attraction, G attr, between ionized monomers caused by temporary bridges due to the hydrophobic-aggregation effect of adsorbed surfactant ions on polymer chains, see the fourth line of eq. Here, the bridge is a kind of short-range attractive interaction derived from the formation of electrostatic dipoles between monovalent surfactant ions and monovalent ionic monomers, which leads us to a rather simple statistical construction of G attr. Notice that unbound surfactant molecules can also assemble hydrophobically with a bound surfactant molecule. Thus, the presence of a free surfactant molecule close to polymer chains will inevitably frustrate the effective temporary cross-linking effect. Under a mean-field consideration, we assume that a bound surfactant molecule meets a second bound surfactant molecule with a probability (φ^2^ c). However, the effective bridge can be formed only if free surfactant molecules are not present, which is expressed by the probability (1 – ρ), where ρ is the volume fraction of the surfactant ion close to polyelectrolyte chains in the solvent phase. The coupling energy (ε_2_) stands for the hydrophobic energy gain per surfactant that is adsorbed on a polyelectrolyte chain. Since this hydrophobic energy gain originates in the overlap between hydrophobic tails (see Figureb), it increases with the degree of polymerization of the surfactant (n). Thus, ε_2_ can be viewed as being a parameter which depends on n but is independent of the solution concentration without restricting generality.
It is worth noting that there is a significant enrichment of oppositely charged surfactants in networks/gels in the phase transition of polyelectrolyte networks/gels, which was found by Khokhlov? and Rumyantsev? with their co-workers. We also note experiments? showed that the preferential adsorption of large-sized surfactant ions on ionic monomers inevitably leads to the enrichment of surfactant ions around polyelectrolyte chains, which indicates that there is a relation of nc _ x _ ≪ ρ → φ. To simplify our following analytical calculations, these facts prompt us to assume that the saturation of the bound surfactant on polyelectrolyte chains in semidilute and concentrated polymer solutions is largely governed by the number of the already adsorbed surfactants; in other words, we use the approximation ρ ≈ φ in the construction of G attr. This approximative approach will avoid heavy calculations without compromising on the physical conclusions because the conformations of polymer chains are network-like in semidilute and concentrated polymer solutions. Notice that a physical boundary condition is explicitly embedded in the approximation ρ ≈ φ, that is, both φ → 0 and φ → 1 leading to the vanishing of G attr → 0. The approximation correctly reflects the fact that the polyelectrolyte is miscible within aqueous solutions at both low and high concentrations of oppositely charged surfactants in reentrant condensation, respectively.
We note that previous experiments ?−? ? reported that the accumulation or preferential adsorption of surfactant ions on polyelectrolyte chains can lead to charge inversion of some polyelectrolytes at certain concentrations of oppositely charged surfactants. This implies in physics that the free energy of the polyelectrolyte solution must have a minimum (or minimum) or exhibit symmetry breaking around these surfactant concentrations. Notice that this consideration is reflected by the statistical construction of G attr in an interpolative way with a minimum and a symmetry breaking around φ = 2/3 where the coupling factor φ^2^(1 – φ) shows its maximum, which corresponds to the fact that the symmetry of simple adsorption and desorption states of ionic monomers (φ or 1 – φ states) by surfactant ions is broken and which further corresponds to the maximum coupling between adsorption and attraction effects of surfactant ions on polyelectrolyte chains when the phase transition occurs from a soluble to a collapsed state or vice versa (including both collapse and reentry transitions). This viewpoint will become clear in detail in the analytical analyses in Sections and 3.2. Because of this feature, the fraction of polyelectrolyte chains adsorbed preferentially by surfactant ions (φ) is also the order parameter for the phase transition in our model.
The surfactant concentration at which the reentrant condensation of the polyelectrolyte occurs is usually below about 100 mmol/L. ?−? ? ? This fact indicates that the free energy of nonassociative pairwise-like electrostatic interaction due to the long-range correlation of all ions can be approximated by the classical “double screening theory” ?,? in the low-salt limit as G DS with its truncated form, see the fifth line of eq. Here, the nonassociative attraction correlation of all small ions is approximated by the Debye–Hückel theory? for the isothermal–isobaric (NPT) ensemble. The inverse Debye screening length κ is given by (κa)^2^ = 4π(l B/a)(pc + (n + 1)c _ x _) for polyelectrolytes in diluted surfactant solutions, where l B is the Bjerrum length of water. There is a strong charge neutralization effect ?,? when monovalent surfactant ions adsorb on monovalent ionic monomers because of the formation of electrostatic dipoles. This results in the contribution of the electrostatic repulsion between ionized monomers for G DS being insignificant compared with other terms in free energy because the charge–charge interaction between ionized monomers is highly screened (see a similar argument in chapter six of the excellent monograph by Muthukumar?). Therefore, in the following discussion, we will ignore it to simplify analytical calculations.
Then, the total Gibbs free energy per volume unit is simply considered as G(φ, c, c _ x _) = G sol + G FH + G ads + G attr + G DS. The essential idea of our model is that the short-range hydrophobic-aggregation effect of a polyelectrolyte solution is statistically constructed by the key formalism of G attr in eq. And the long-range electrostatic correlation of all ions is considered by the classical “double screening theory”, ?,? i.e., G DS in eq. We note that this approach for long-range electrostatic interactions does not require that the polyelectrolyte has to be dilute, as indicated by a recent comprehensive study.? We list the symbols used in the theory in Table for the reader’s convenience. Notice that in this study, we deliberately neglect the possible hydrophobic adsorption between the surfactant tail and a hydrophobic polymer backbone, as well as the related electrostatic repulsion effect among adsorbed surfactant ions. This approach will significantly simplify our analytical calculations and guide us to focus on the well-studied phase-transition mechanism, ?,?−? ? ? as illustrated by Figureb, and will finally lead us to qualitative but rather simple answers for the two fundamental questions raised in the Introduction.
1: List of Symbols Used in the Theory
Results and Discussion
3
The Minimum Coupling Energy for the Hydrophobic-Aggregation
Effect in the Phase Transition
3.1
Following a similar computing approach as our previous work? and the work by Sommer,? we can estimate a minimum coupling energy (ε_2_) for the hydrophobic-aggregation effect that is necessary for a phase transition. Based on the construction of G attr in eq, we see that the maximum contribution to the adsorption and attraction effects of surfactant ions on polyelectrolyte chains occurs around the symmetry-breaking point φ = 2/3. As long as this contribution is dominant so that it defines the overall symmetry properties of the entire free energy function, this simplification allows us to mathematically cast our model to a canonical ensemble-like model by setting φ = 2/3 for the case of a very diluted solution of the surfactant ((n + 1)c _ x _ → 0), at which both the chemical potential μ and the osmotic pressure ∏ approach their limiting values μ_0_ and ∏0, respectively:
Here, the canonical ensemble-like free energy is given by G – ∏0.
To show a noticeable hydrophobic effect, the surfactant chain length must be above a minimum length. This implies that there exists a minimum coupling energy for the hydrophobic-aggregation effect (ε_2_) that is necessary for a phase transition. By eq, we can roughly estimate a minimum value for the coupling energy (ε_2_). This can be realized by determining the boundary condition for the spinodal decomposition of polyelectrolyte solution, which is given by d^2^(G – ∏0)/dc ^2^ = 0 in analogy to the framework of the canonical ensemble and follows with refs ? and ?
Here, we define an overall effective interaction parameter 2χ_0_ ≡ 8ε_2_ p ^2^/27 + 2ε_FH,1_(1 + p)p + 2ε_FH,2_(1 – p ^2^). Its critical or minimum value to allow phase transition is given by d(2χ_0_)/dc = 0, which follows the approach by classical monographs. ?,? This approach is mathematically equivalent to d^3^(G – ∏0)/dc ^3^ = 0, which? finds the critical point for the free energy G – ∏0,
There exists no exact explicit analytical solution for c in eq, but an acceptable approximation can be obtained by ignoring the fourth term since it is much smaller than the sum of the other terms within experimental values of the parameter l B/a. For the critical point of the polyelectrolyte collapse, we obtain the approximation for the critical volume fraction of monomers c at
The exact solution of eq is recovered by eq when the parameter p approaches zero. One can improve upon the analytical solution of eq by employing the method of fixed-point iteration ?−? ? with the initial value defined by eq. However, the approximative solution of eq is sufficient to deduce the key features of the overall effective interaction parameter χ_0_ without compromising the physical conclusions. The details behind the computational solution of this approach can be found in Section A of Supporting Information.
By insertion of eq into eq, we get an estimation of the critical or minimum value of χ_0_,
Equation shows that a shorter polymer chain (N) prescribes a higher minimum coupling energy (ε_2_) that is necessary for a phase transition. This correctly reflects the obvious but nontrivial fact that a shorter polymer chain requires a higher enthalpic attraction among polymer chains to overcome/compensate for a larger translational entropy of polymer chains in phase transition. By setting p = 0 in eq, we recover the boundary condition for the spinodal decomposition of an uncharged polymer solution, i.e., . We derive the minimum of ε_2_ for the case of a very long polyelectrolyte chain (N → ∞) by eq,
As discussed before, the strength of the hydrophobic-aggregation effect obviously depends on the surfactant’s chain length. It is also well-known that the strength of hydrophobic attraction is about 1.5 k B T between two methyl/methylene groups.? Therefore, there is a minimum surfactant chain length (n min) ?−? ? ? to induce phase separation of the polyelectrolyte, which is numerically estimated by n min ≈ 2(ε_2_)min/3 according to eqs and ?. As an example of a hydrophilic polyelectrolyte, we choose model parameters p ≈ 0.9, l B/a ≈ 2, N = 6000 (or N = 30), and ε_FH,1_ = ε_FH,2_ ≈ 0.41 for ionized poly(acrylic acid); ?,? then, eq predicts that the minimum surfactant chain length to induce phase separation is about n min ≈ 7 for the mixtures of ionized poly(acrylic acid) and alkyl trimethylammonium bromides. It is noticeable that this prediction is quite close to its experimental value (n min ≈ 8)? for the chosen set of model parameters. For the coacervation of another polyelectrolyte poly(diallyl dimethylammonium chloride) with surfactant sodium alkyl sulfate,? one can estimate a similar minimum surfactant chain length (n min ≈ 8) to initiate phase transition by using similar experimental values of model parameters.?
A feature of the constructions of eqs and ? is that the minimum coupling energy (ε_2_)min for the hydrophobic-aggregation effect estimated by eqs and ? is not related to the model parameters ε_1_, λ, and χ_s_. This is an expected result because one cannot expect a phase transition merely according to the simple exchange effect of surfactant ions on polymer chains based on a formal analogy with the well-understood 1D-Ising model. ?,? Here, if we solely take into account pairwise-like interactions for preferential adsorption between surfactant ions and ionic monomers, linear polymers are 1D substrates for ions.
Another feature of the constructions of eqs and ? is that, as shown in Figure, no matter the values of l B/a ≥ 0 and ε_FH,1_ ≥ 0, there is a unique local maximum of (ε_2_)min when ε_FH,2_ is larger than , but no local maximum of (ε_2_)min exists when 0 ≤ ε_FH,2_ ≤ . This feature is coincidental with the well-known Θ condition for uncharged polymers. Thus, we can see that the minimum coupling energy (ε_2_)min for the hydrophobic-aggregation effect to lead to a phase transition is related to the solvent quality (parameter ε_FH,2_) for the uncharged part of the polyelectrolyte. We note that this analytical result reflects the experimental fact ?,? that the hydrophobicity of uncharged monomers of polyelectrolytes has a noticeable effect on the necessary surfactant concentration to trigger the reentry transition.
According to eq , the minimum coupling energy (ε2) with respect to the fraction of charged monomer (p) for typical values of the parameter εFH,2 with l B/a = 2, εFH,1 = 0.45, and N → ∞.
The Nonmonotonic Effective Flory–Huggins
χ Parameter
3.2
In order to find the equilibrium state of the polyelectrolyte phase with respect to the bulk solvent phase, we follow a similar computing approach as in our previous work.? First, we minimize the free energy G(φ, c, c _ x _) with respect to the adsorption fraction of surfactant ions φ. Then, the replacement of the equilibrium solution for φ(μ, c, c _ x _) into the free energy leads to an effective free energy for the polyelectrolyte, in which the effect of surfactant ions is mapped onto an effective monomer–monomer interaction, which depends on the concentration of surfactant ions. Finally, the minimization of the free energy per monomer, i.e., G(φ, c, c _ x _)/c instead of G(φ, c, c _ x _), with respect to the monomer concentration (c) leads to the equilibrium solution of the model.
Because we are interested in the case of a very diluted solution of surfactant ((n + 1)c _ x _ → 0) in this work, we ignore the influence of surfactant concentration in the bulk solution in the construction of G sol, G FH, and G DS. The physical foundation behind this assumption is that the surfactant concentration for the reentrant condensation of a polyelectrolyte can be several orders of magnitude below the bulk critical micelle concentration (CMC) of the surfactant,? which implies that the influence of mixing behaviors of the surfactant in bulk solution is negligible for the reentrant condensation. This assumption allows us to focus on the adsorption and attraction effects of surfactants near/on polymer chains and ignore their own nonessential mixing effects if surfactant molecules are far away from polymer chains, which can avoid heavy calculations without losing generality of the key physical conclusions.
The construction of G attr in eq implies that the maximum contraction of polymer chains is reached around the maximum-coupling point φ = 2/3, where the hydrophobic-aggregation effect reaches its maximum, which corresponds to the phase transition from a soluble to a collapsed state or vice versa (including both collapse and reentry transitions). This peculiar feature of the model leads us to the introduction of a perturbation (δ) from the maximum-coupling point φ = 2/3 according to
This perturbation approach follows closely our previous work.? With the Taylor expansion of δ-containing terms in the logarithm function up to the accuracy of square terms (δ^2^) under the constraint of |δ^3^| ≪ 1 and ignoring constant terms, we obtain
with
The numerical term “ln 2” in eq is due to the fact that the symmetry of simple adsorption and desorption states of ionic monomers (φ or 1 – φ states) by surfactant ions is broken. This can be seen from the transformation of G ads in eq by the manipulation of eq.
Minimizing the free energy in eq with respect to δ yields
When the surfactant’s bulk chemical potential μ is close to μ_C_ = −((ε_1_ + ln 3 – ln 2)/λ – ln 3 + χ_s_/3 + 1 – 1/λ), the value of δ is close to zero. This case ideally follows the small perturbation expansion introduced by eq and is also the most important situation we are going to focus on in this research. This is discussed in further detail in the following subsections. Resubstituting eq into eq, yields
with the effective Flory–Huggins parameter χ_eff_ as follows:
and the entropic term of free energy per monomer g sol as follows:
The above constructed χ-function is a function of the square of the chemical potential μ, which indicates that a χ_eff_ corresponds to two values of μ and thus captures the reentrant signature of polyelectrolyte condensation at lower and higher surfactant concentrations. According to eq, we know that the reentrant signature of polyelectrolyte condensation is controlled by the hydrophobic-aggregation effect, since only this effect is nonmonotonic with respect to the surfactant concentration.
The General Spinodal and Binodal Phase Diagrams
3.3
In order to analytically discuss the phase transition, the pressure isotherm Π(c, μ; N, l B/a, p, ε_1_, ε_2_, ε_FH,1_, ε_FH,2_) can be calculated from eq by the condition ∂(G/c)/∂c = 0, which leads to
By setting p = 0 in eq for an uncharged polymer solution, with ∂Π/∂c = 0 and the constraint of 0 < c < 1, as it should be, we recover the boundary condition for the spinodal decomposition of the uncharged polymer solution, i.e., .
For the general case of a discontinuous phase transition, the osmotic pressure must display an unstable region of negative compressibility given by ∂Π/∂c < 0. In Figurea, we display two examples showing the two critical points of the pressure isotherms given by eq, i.e., corresponding to the states of ∂Π/∂c = ∂^2^Π/∂c^2^ = 0. In Figureb, we display an example that shows the pressure isotherms for a discontinuous condensation transition given by eq. The coexistence region of the binodal is defined by the Maxwell construction as indicated by the equal area criterion (the area of the envelope “ABC” is equal to that of the envelope “CDE”: S 1 = S 2) with the horizontal isobaric line in the figure (the black dashed line “ACE”). This cannot be obtained analytically in an exact way but must be computed by numerical calculations. However, we are able to calculate the spinodal state analytically at which the solution begins to turn unstable, the existence of which is the necessary condition for a discontinuous transition scenario. The spinodal of the polyelectrolyte solution is given by ∂Π/∂c = 0 and can be written in the following form:
Here, we denote the external adsorption field U as
This defines the spinodal phase diagram in the “μc” space with the eight parameters p, N, l B/a, ε_1_, ε_2_, χ_s_, ε_FH,1_, and ε_FH,2_.
According to eq , the osmotic pressure of the polyelectrolyte solution is plotted as a function of the inverse volume fraction of polyelectrolyte monomers for the parameters p = 0.6, l B/a = 2, ε1 = 6, ε2 = 15, εFH,1 = 0.4, εFH,2 = 0.4, λ = 1, χs = 2, and N → ∞. In panel (a), we show the two critical points (indicated by filled circles with symbols “A” and “B”) of the pressure isotherms by choosing μ = −6.80 and μ = −5.15, which also correspond to μA and μB, respectively, in Figure . In panel (b), we choose μ = −((ε1 + ln 3 – ln 2)/λ – ln 3 + χs/3 + 1 – 1/λ) = −5.97, which also corresponds to μC in Figure . The coexistence pressure of the binodal, obtained by the Maxwell construction, is indicated by the horizontal dashed line in the figure (the black dashed line “ACE”) with the equal area criterion (the area of the envelope “ABC” is equal to that of the envelope “CDE”: S 1 = S 2), and the spinodal points are indicated by filled circles in the figure (symbols “B” and “D”).
In Figure, we display the spinodal phase diagram (closed red curve “ACBD”) in the “μc” space given by eq and its corresponding binodal diagram (closed blue curve “AEBF”) by numerical calculations for a typical case of polyelectrolyte with parameters p = 0.6, l B/a = 2, ε_1_ = 6, ε_2_ = 15, ε_FH,1_ = ε_FH,2_ = 0.4, λ = 1, χ_s_ = 2, and N → ∞. The two solutions at the same value of μ correspond to the two extrema of the pressure isotherm, which imply the coexistence of a condensed polymer phase and a dissolved polymer phase. The region of phase separation is closed topologically. For our theory, we obtain collapse and re-entry transitions that are symmetric; for example, the closed curves in Figure are symmetric with respect to the tie-line “ECDF”. The lower part of μ defines the collapse transition as indicated by the lower half of the “egg-shape” curve, while the higher part of μ defines the reentry transition. In Figure, μ_A_ indicates the surfactant concentration above which the collapse transition can occur, which is close to the critical aggregation concentration (CAC). μ_B_ indicates the surfactant concentration below which the reentry transition can take place, which is close to the second critical aggregation concentration (CAC_2_). μ_C_ indicates the surfactant concentration at which the polyelectrolyte solution shows the maximum extent of phase separation and the effective Flory–Huggins parameter χ_eff_ reaches its maximum (see eq).
According to eq , a typical spinodal phase diagram (closed red curve “ACBD”) in the “μc” space and its corresponding binodal diagram (closed blue curve “AEBF”), obtained by numerical calculations for a polyelectrolyte in the dilute solution of an oppositely charged surfactant for the case of p = 0.6, l B/a = 2, ε1 = 6, ε2 = 15, εFH,1 = εFH,2 = 0.4, λ = 1, χs = 2, and N → ∞. Both the binodal diagram and the spinodal diagram are symmetric with respect to the tie-line “ECDF”. In the figure, μA indicates the lower limit surfactant chemical potential for the collapse transition, μB indicates the upper limit surfactant chemical potential for the reentry transition, and μC indicates the chemical potential at which the polyelectrolyte solution shows the maximum extent of phase separation.
From Figure, we see that the external adsorption field U has a minimum and a maximum at μ_A_ and μ_B_, respectively. Because we obtain symmetric collapse and reentry transitions for our theory, these facts indicate that there is a unique monomer concentration (c AB) corresponding to both μ_A_ and μ_B_. It can be solved in a formal way by
Then, the values of μ_A_ and μ_B_ can be obtained by the insertion of the unique solution of the monomer concentration (c AB) into eq. It is remarkable that the unique monomer concentration (c AB) at both μ_A_ and μ_B_ are also the two critical points of the osmotic pressure in eq, i.e., corresponding to the solution of ∂Π/∂c = ∂^2^Π/∂c^2^ = 0. As shown in Figure, this property can be seen graphically since there are two symmetric tangent points, and each is at (c AB, μ_A_) and (c AB, μ_B_) for the spinodal and its corresponding binodal phase diagrams. Here, we point out that in general, the best way to obtain μ_A_ and μ_B_ is via numerical methods because it is challenging to derive simple analytical approximations. Nevertheless, because we obtain symmetric collapse and reentry transitions for our theory, we can derive the following exact and simple relations:
The values of μ_A_, μ_B_, and μ_C_ are usually comparable with each other in experimental observations; ?,?−? ? ? this can also be easily seen from the example in Figure.
As shown in Figure, when the bulk chemical potential of surfactant ions is close to μ_C_, the tie-line “ECDF” shows that there is an obvious difference in the polymer concentration of the dilute polymer phase as predicated by the spinodal and binodal phase diagrams (symbols “C” and “E” in Figure). A similar noticeable difference is also observed for the polymer concentration of the condensed polymer phase (symbols “D” and “F” in Figure). This basic characteristic of our model implies that there may exist rich nonequilibrium/metastable/kinetically trapped states when the polymer–surfactant coacervates evolve from a spinodal to a binodal stateconsidering the fact that the time scales to reach a true equilibrium state for a polymer solution can span months or even years,? which were indeed observed in experiments. ?,?,?
Another basic characteristic of our model is that it predicts that complete polyelectrolyte compaction can occur at far less than 100% of the degree of charge compensation of ionic monomers. This prediction can be easily seen from the tie-line “ECDF” in Figure, where the surfactant concentration (μ_C_) is insufficient to compensate for the charges of all ionic monomers when the polyelectrolyte solution shows the maximum extent of phase separation. It is remarkable that this prediction was confirmed in detail by experimental studies ?,? on the phase behaviors of both natural and synthetic polyelectrolytes (DNA and poly(acrylic acid)) in the presence of oppositely charged surfactants. According to eqs and ?, we see that the model shows some flexibility in realizing charge compensations and charge regulations of polyelectrolytes by varying the adsorption strength of surfactant ions on ionic monomers (ε_1_), the demixing strength between the surfactant and water (χ_s_), and the size of the ionic head of a surfactant (λ). The experimental studies ?,? and our theory clearly show that it is not necessary to stick to the concept of full charge compensation to explain the collapse transition of a polyelectrolyte in the presence of an oppositely charged surfactant, a concept which was employed by Rumyantsev and co-workers? for the study of polyelectrolyte gels.
Comparison to Additional Experiments
3.4
In this subsection, we compare analytical predictions of our theory with additional experimental results reported in the literature.
In Figure, we display the spinodal phase diagram according to eq for various polyelectrolyte chain lengths (N) under the condition of p = 0.6, l B/a = 2, ε_1_ = 6, ε_2_ = 18, ε_FH,1_ = ε_FH,2_ = 0.4, λ = 1, and χ_s_ = 1. As shown in Figure, the phase transition can be influenced drastically by the chain length (N) of a polyelectrolyte. We see that an increase in polyelectrolyte chain length promotes the coexistence region of collapse transition to a lower surfactant concentration and shifts the coexistence region of reentry transition to a higher surfactant concentration. However, this chain-length dependence diminishes for long polyelectrolyte chains. In contrast to uncharged linear polymer, ?,?,? by setting the external adsorption field U = 0 in eq, our model shows that it is hard to realize a real dilute phase for the reentrant condensation of polyelectrolyte when the charge fraction is not sufficiently low. As indicated in both Figures and ? for the binodal and spinodal diagrams, this is particularly noticeable in the limiting case of an infinite polyelectrolyte chain length. There remains no small monomer concentration (c) in the dilute phase for small values of the external field U when phase separation occurs, i.e., close to the optimally loaded state of the polyelectrolyte with surfactant ions where the effective Flory–Huggins parameter χ_eff_ reaches its maximum (see eq).
Spinodal phase diagrams of a polyelectrolyte in the dilute solution of an oppositely charged surfactant for different polyelectrolyte chain lengths (N) under the condition of p = 0.6, l B/a = 2, ε1 = 6, ε2 = 18, εFH,1 = εFH,2 = 0.4, λ = 1, and χs = 1, as computed from eq .
Our theory also predicts that there is no phase transition if the polyelectrolyte chain length is too short. With the parameter set used for Figure, no phase transition occurs with chain lengths of N < 10. Actually, this is a consequence of the fact that a shorter polyelectrolyte chain requires a higher minimum coupling energy for the hydrophobic-aggregation effect in phase transition; see eq. Furthermore, the impact of translational entropy of the chain increases with shorter chain length, which generally works against the formation of aggregates. It is remarkable that these analytical predictions of our theory were confirmed in detail by previous experimental studies ?,?,?,? on the influence of polyelectrolyte chain length in the phase separation of polyelectrolyte admixed with an oppositely charged surfactant.
Remarkably, similar to the above discussion, our theoretical approach can be leveraged to explain and analyze some co-condensation behaviors of DNAs/RNAs with proteins/peptides (such as FUS and HP1 proteins). ?−? ? ? ? ? ? ? ? Here, a relatively short protein–peptide chain acts analogously to a surfactant. For example, the phase separation of DNA induced by protein binding is dependent on the DNA length, ?,? which can be explained by our theory on the effect of polyelectrolyte chain length (see discussion for Figure). Here, only one specific binding site of a protein anchors on the DNA chains, whereas other nonspecific sections of the protein can interact with other proteins via physical cross-linking, a physical picture which is similar to the phase-transition mechanism shown in Figureb.
In Figurea, we display the spinodal phase diagram according to eq for different strengths of the hydrophobic-aggregation parameter (ε_2_) under the condition of p = 0.6, l B/a = 2, ε_1_ = 6, ε_FH,1_ = ε_FH,2_ = 0.4, λ = 1, χ_s_ = 1, and N → ∞. We see that an increase in the strength of the hydrophobic-aggregation effect (ε_2_), which is related to an increase in the surfactant chain length (n), shifts the coexistence region of the collapse transition to a lower concentration of the surfactant but shifts the coexistence region of the reentry transition to a higher concentration of the surfactant, which is in agreement with experimental results of highly hydrophilic polyelectrolytes in aqueous solutions of oppositely charged surfactants. ?,?−? ? ? However, when the surfactant chain length exceeds a maximum length, it is expected that the chain-length effect of the surfactant will saturate due to the lack of an efficient stack of steric surfactant tails in their premicellar aggregation (see Figureb for a sketch of the premicellar aggregation). But in this situation, the demixing effect between the surfactant and water (χ_s_) still plays a role on polyelectrolyte chains. In Figureb, we display the spinodal phase diagram according to eq for different demixing effects between the surfactant and water (χ_s_) with a parameter set of p = 0.6, l B/a = 2, ε_1_ = 6, ε_2_ = 15, ε_FH,1_ = ε_FH,2_ = 0.4, λ = 1, and N → ∞. We see that an increase in the demixing effect between the surfactant and water (χ_s_) shifts the coexistence region of both collapse and reentry transitions to lower concentrations of the surfactant. This analytical result indeed qualitatively explains the phase-behavior difference of a highly hydrophilic polyelectrolyte in aqueous solutions of hexadecyl- and dodecyl-trimethylammonium bromides (both surfactants have a very long steric alkyl tail).?
(a) According to eq , spinodal phase diagrams of a polyelectrolyte in the dilute solution of an oppositely charged surfactant for different strengths of hydrophobic-aggregation parameter (ε2) under the condition of p = 0.6, l B/a = 2, ε1 = 6, εFH,1 = εFH,2 = 0.4, λ = 1, χs = 1, and N → ∞. (b) According to eq , spinodal phase diagrams of a polyelectrolyte in the dilute solution of an oppositely charged surfactant for different demixing effects between the surfactant and water (χs) under the condition of p = 0.6, l B/a = 2, ε1 = 6, ε2 = 15, εFH,1 = εFH,2 = 0.4, λ = 1, and N → ∞.
Interestingly, by an increase in charge fraction (p) of a polyelectrolyte chain while keeping other model parameters unchanged, the corresponding spinodal phase diagrams (similar to Figurea; also see eq) indicate that the coexistence region of the collapse transition shifts to a lower concentration of the surfactant, but the coexistence region of the reentry transition shifts to a higher concentration of the surfactant. It is noticeable that this analytical result is corroborated in detail by experimental studies of a recent dissertation? on the phase behaviors of poly(acrylic acid) and sodium polyacrylate in aqueous solutions of hexadecyl trimethylammonium bromide.
An Estimation of the Critical Aggregation
Concentration (CAC)
3.5
We note that during the process of charge neutralization between a polyelectrolyte and an oppositely charged surfactant, a strong electrostatic dipole between a monovalent ionic monomer and a monovalent surfactant ion can form because of a noticeable reduction of the dielectric constant around the polyelectrolyte chains, arising from the long alkyl surfactant tails. ?,? This indicates a significant shift of ε_1_ on the left-hand side of the spinodal construction, i.e., eq, which is often not small for a polyelectrolyte (on the order of about 5 k B T for the strength of an ionic bond? in water at low salt concentrations and on the order of about 10 k B T for the strength of the electric dipole interaction?). This fact implies that eq melts down when the chemical potential μ is far away from zero, which is approximately given by the condition that the external adsorption field U is close to zero, i.e., eq. This way, we can estimate the critical aggregation concentration (CAC) of the surfactant in the surfactant–polyelectrolyte binding isotherm for the occurrence of a reentrant condensation, which is given by
We can extract some interesting physics from eq. First, the larger size of the ionic head of a surfactant ion (λ) can promote the CAC to a higher surfactant concentration. It is remarkable that this analytical prediction was observed in previous experimental studies ?,?,? on the reentrant condensation of some polyelectrolytes. Second, if the values of ε_1_ and χ_s_ are sufficiently large, as, for example, with a rather moderate choice of parameters ε_1_ = 6, χ_s_ = 1, and λ ≈ 1, the surfactant concentration which permits a polyelectrolyte reentrant condensation to occur is in the order of about (n + 1)c _ x _ ∼ CAC ≈2.0 × 10^–3^ (volume fraction in solution). It is remarkable that this CAC value is lower than the bulk critical micelle concentration (CMC) of some ionic surfactants reported in literature.? We note that this can be the case for both the collapse and reentry branches of the reentrant condensation, as shown from Figures–?. According to this simple estimation, we already see that the ionic surfactant concentrations for polyelectrolyte reentrant condensation can be lower than its bulk critical micelle concentration (CMC) even if we do not evaluate the CMC of an ionic surfactant by its molecular structure parameters (a, n, λ, and ε_2_ in our model; see Section B of Supporting Information for details). Third, because λ ≈ 1 and ε_1_ ≫ χ_s_/3 for most common cases, by eq, we see that the primary factor to induce both collapse and reentry transitions of polyelectrolytes at low surfactant concentrations is the strong electrostatic adsorption between ionic monomers and surfactant ions (ε_1_). Notice that these predictions concur with all experimental observations reported in the literature ?−? ? ? ? ? ? ? ? on the reentrant condensation of polyelectrolytes in dilute aqueous solutions of oppositely charged surfactants.
Conclusions and Remarks
4
In summary, in this work, we explore the reentrant condensation of polyelectrolytes in the presence of an oppositely charged surfactant, a phenomenon whose phase-transition mechanism is of fundamental importance to the understanding of liquid–liquid phase separation (LLPS) in soft materials and biological systems. Compared with previous theoretical formalisms, ?,?,?−? ? ? ? ? ? ? ? ? ? ? ? ? the novelty of the present work is that we focus on the adsorption and attraction effects of surfactants near/on polymer chains and ignore their own nonessential mixing effects if surfactant molecules are far away from polymer chains. This novel approach allows us to construct a simple mean-field equilibrium theory and to solve it analytically with closed-form solutions, which can rationalize the essential features (such as the emergent “egg shape” of spinodal and binodal phase diagrams) of the reentrant condensation of a polyelectrolyte induced by diluted oppositely charged surfactants. This approach also allows us to clearly address the fact that a strong electrostatic adsorption between the ionic monomers and surfactant ions is critical to understand the peculiar phenomenon that both the collapse and reentry transitions of polyelectrolytes can even occur when the concentration of the surfactant is below its bulk critical micelle concentration (CMC). The analytical solution of the theory also indicates that a minimum coupling energy for the nonlinear hydrophobic-aggregation effect of adsorbed surfactants is essential for a phase transition to occur, which explains why polyelectrolytes show that phase transition only if the surfactant chain length is above a certain threshold value. ?−? ? ?
In contrast to uncharged linear polymers, ?,?,? a distinctive finding of the present theory in this work is that it is hard to realize a real dilute polymer phase when the fraction of charged monomers is not sufficiently low. This is clear from our theoretical calculations for the limiting case of an infinite chain length. Our theory justifies that the monomer charge is an important factor in the regulation of polymer liquid–liquid phase separation. In addition, we found that it is not necessary to stick to the concept of charge compensation to explain the reentrant condensation of polyelectrolytes, which has been confirmed in detail in previous experimental studies. ?,? Our theory showed that it can flexibly realize charge compensation and charge regulation of polyelectrolytes in phase transitions by changing the adsorption strength of surfactant ions on ionic monomers (ε_1_), the demixing strength between the surfactant and water (χ_s_), and the size of the ionic head of a surfactant (λ).
Nevertheless, it is necessary to point out that our approximations neglect contributions that would describe the surfactant’s bulk phase transitions, including micellar formation. The surfactant’s entropic mixing terms are included in the free energy function, but these contributions would hamper the simplicity of the analytical solutions of the model. We have thus, for the sake of clarity and transparency of the formalism, neglected these terms in our analytical approximations and rather focused on systems in which polyelectrolyte condensation occurs at concentrations well below the critical micellar concentration of a surfactant. Another consequence of the restriction to low surfactant concentrations is that the surfactant–polymer interaction cannot be weak, i.e., the adsorption strength of surfactant ions on ionic monomers is not small (ε_1_ ∼ 5 k B T). For the case of weak interaction (ε_1_ ∼ 1 k B T), a high ionic surfactant concentration would be required to induce a reentry transition, which can be well-explained by classical polyelectrolyte theories. ?−? ? We also point out that a limitation of the current theory is its neglect of the charge distribution on the polyelectrolyte chains. This approach simplifies our analytical calculations but underestimates the overall entropic effect, which predicts symmetric collapse and reentry transitions, which are, however, not always the experimental observation. A starting point of theoretically accounting for the charge distribution of polyelectrolyte chains can be the calculation framework summarized by Avni, Andelman, and Podgornik.?
Another limitation of the current theory is that it has neglected the interaction between the polyelectrolyte backbone and the hydrophobic tail of the surfactant. This implies that our simplified formalism is primarily confined to a hydrophilic polyelectrolyte backbone, i.e., in the case of and . This approach simplifies our analytical calculations but underestimates the overall hydrophobic effect in the polyelectrolyte phase transition, which may limit its predictions on the phase behaviors of certain amphiphilic polyelectrolytes in the presence of surfactants. For example, the hydrophobic effect can affect and decrease considerably the surfactant concentration that is necessary to trigger the reentry transition, as pointed out by previous experiments. ?,? Under certain conditions, the hydrophobic effect may lead to a liquid-like polyelectrolyte phase in coexistence with a solid-like polyelectrolyte phase (liquid–solid phase separation or LSPS; see two classical review papers on this topic ?,? ), which is beyond the scope of our simplified theoretical formalism. In principle, for the case of polyelectrolyte backbones with strong hydrophobic properties, our theoretical approach can be refined and extended to accommodate additional hydrophobic adsorption between the surfactant tail and the polymer backbone as well as the related electrostatic repulsion effect among adsorbed surfactant ions. However, a detailed investigation of these above aspects to refine our theory is worthy of future consideration and lies beyond the goal of the present study.
Supplementary Material
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Li J.Parakhonskiy B. V.Skirtach A. G.A decade of developing applications exploiting the properties of polyelectrolyte multilayer capsules Chem. Commun.202359780783510.1039/D 2CC 04806 J 36472384 · doi ↗ · pubmed ↗
- 2Guzmán E.Rubio R. G.Ortega F.A closer physico-chemical look to the Layer-by-Layer electrostatic self-assembly of polyelectrolyte multilayers Adv. Colloid Interface Sci.202028210219710.1016/j.cis.2020.10219732579951 · doi ↗ · pubmed ↗
- 3Ivanov A. S.Pershina L. V.Nikolaev K. G.Skorb E. V.Recent Progress of Layer-by-layer Assembly, Free-Standing Film and Hydrogel Based on Polyelectrolytes Macromol. Biosci.202121210011710.1002/mabi.20210011734272830 · doi ↗ · pubmed ↗
- 4Holehouse A. S.Kragelund B. B.The molecular basis for cellular function of intrinsically disordered protein regions Nat. Rev. Mol. Cell Biol.202425318721110.1038/s 41580-023-00673-037957331 PMC 11459374 · doi ↗ · pubmed ↗
- 5Achazi K.Haag R.Ballauff M.Dernedde J.Kizhakkedathu J. N.Maysinger D.Multhaup G.Understanding the Interaction of Polyelectrolyte Architectures with Proteins and Biosystems Angew. Chem., Int. Ed.20216083882390410.1002/anie.202006457 PMC 789419232589355 · doi ↗ · pubmed ↗
- 6Caruso M. R.D’Agostino G.Milioto S.Cavallaro G.Lazzara G.A review on biopolymer-based treatments for consolidation and surface protection of cultural heritage materials J. Mater. Sci.20235832129541297510.1007/s 10853-023-08833-5 · doi ↗
- 7Gîfu I. C.IanchişR.Nistor C. L.Petcu C.Fierascu I.Fierascu R. C.Polyelectrolyte CoatingsA Viable Approach for Cultural Heritage Protection Materials 2023167287310.3390/ma 1607287337049167 PMC 10096418 · doi ↗ · pubmed ↗
- 8Puente-Santamaría A.Ortega F.Maestro A.Rubio R. G.Guzmán E.Non-equilibrium states in polyelectrolyte-surfactant systems at fluid interfaces: A critical review Curr. Opin. Colloid Interface Sci.20247110180410.1016/j.cocis.2024.101804 · doi ↗