Next Article in Journal / Special Issue
Small Oligonucleotides Detection in Three-Dimensional Polymer Network of DNA-PEG Hydrogels
Previous Article in Journal / Special Issue
Hyaluronic Acid Functionalization with Jeffamine® M2005: A Comparison of the Thermo-Responsiveness Properties of the Hydrogel Obtained through Two Different Synthesis Routes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermoreversible Gelation with Two-Component Mixed Cross-Link Junctions of Variable Multiplicity in Ternary Polymer Solutions

Department of Polymer Chemistry, Graduate School of Engineering, Kyoto University, Katsura, Kyoto 615-8510, Japan
Submission received: 7 June 2021 / Revised: 7 July 2021 / Accepted: 8 July 2021 / Published: 11 July 2021
(This article belongs to the Special Issue Gels Horizons: From Science to Smart Materials)

Abstract

:
Theoretical scheme is developed to study thermoreversible gelation interfering with liquid–liquid phase separation in mixtures of reactive f-functional molecules R{A f } and g-functional ones R{B g } dissolved in a common solvent. Formed polymer networks are assumed to include multiple cross-link junctions containing arbitrary numbers k 1 and k 2 of functional groups A and B of each species. Sol-gel transition lines and spinodal lines are drawn on the ternary phase plane for some important models of multiple cross-link junctions with specified microscopic structure. It is shown that, if the cross-link structure satisfies a certain simple condition, there appears a special molar ratio of the two functional groups at which gelation takes place with a lowest concentration of the solute molecules, as has been often observed in the experiments. This optimal gelation concentration depends on f and g (functionality) of the solute molecules and the numbers k 1 and k 2 (multiplicity) of the functional groups in a cross-link junction. For cross-links which allow variable multiplicity, special attention is paid on the perfectly immiscible cross-links leading to interpenetrating polymer networks, and also on perfectly miscible cross-links leading to reentrant sol-gel-sol transition. Results are compared with recent observations on ion-binding polymer solutions, polymer solutions forming recognizable biomolecular complexes, polymer/surfactant mixtures, hydrogen-bonding polymers, and hydrophobically-modified amphiphilic water-soluble polymers.

1. Introduction: Binary Gels with Mixed Cross-Links

A majority of physical gels made up of natural polymers contain more than two component functional molecules and/or groups in their cross-link junctions [1,2,3,4,5,6]. Synthetic gels are also often made by mixing chemical cross-linkers, ions, other reactive molecules into polymer solutions [2,4,6]. More recently, dual polymer networks with coexisting chemical and physical cross-links of different structures have been attracting researchers interest [7,8,9,10]. To find the molecular distribution of cross-linked clusters (three-dimensional polymers), their average molecular weight, and the gel point, in terms of their system parameters such as concentration, temperature, reaction time, etc. in such multi-component polymer gels with complex cross-links, we have to know the relation between microscopic structure of the cross-link junctions and macroscopic change in thermodynamic state of the system. For a special case of chemical co-reaction between two species of polydisperse functional molecules, Stockmayer [11] derived, on the basis of the classical tree statistics, the molecular distribution and the weight-average molecular weight of the generated three-dimensional polymers in terms of the degree of reaction p , q for each functional group. From the result he presented the gel point condition for the formation of binary polymer networks. In this paper, we eliminate the assumption of pairwise reaction, and generalize his result to cross-links of arbitrary multiplicity in order to apply the theoretical model to real physical gels. Here the multiplicity is defined by a set of numbers of functional groups contained in a cross-link junction.
One of the important examples of such mixed multiple junctions is ion-binding cross-linking in which variable number of ligands attached on the polymer chains form ion complexes [12] around an added metal ion as a nucleus. Ions attract functional groups (ligands) on the polymer chains to form coordination complexes, which serve as network junctions [12,13]. When the concentration of added ions reaches a critical value, polymer solutions turn into gels. Typical examples of such ionically cross-linked gels are metalic cations (Ca 2 + , Cu 2 + , Zn 2 + , Fe 3 + , ⋯) and basic anions (B (OH) 4 , Sb (OH) 6 , ⋯) complexed by synthetic polymers (poly(vinyl alcohol), polyacrylamide, poly(acrylic acid), ⋯) [14,15,16,17,18,19], or natural polymers (polysaccharides, proteins, ⋯) [20,21,22,23,24,25] bearing functional groups. The characteristic properties, such as reversible gelation, phase separation, viscoelasticity, concentration fluctuations, etc. of these aqueous polymer solutions, originate in the detailed structure of the coordination complexes and variation of their stepwise equilibrium constants. For instance, gels are often observed to turn back to sols under excess presence of ions (reentrant sol-gel-sol transition) [26]. Similar reentrant transitions were also observed in polymer solutions with biomolecularly recognizable cross-link junctions [27].
Other important examples are amphiphilic co-networks [6] such as seen in aqueous solutions of hydrophobic polymers mixed with low-molecular weight surfactants. Interaction between polymers and surfactants has been a subject of great interest [28,29]. The problem was initially laid in studies of proteins associated with natural lipids and with synthetic surfactants. More recently, interaction of water-soluble synthetic polymers such as poly(ethylene oxide) (PEO) with ionic and non-ionic surfactants [30,31,32,33] have attracted interest of researchers because of scientific and technological implications. When polymers carry small fraction of hydrophobic groups, profound influence of added surfactants on the rheological properties has been reported. For example, the plateau modulus of hydrophobically modified PEO solutions exhibits a peak when sodium dodecyl sulfate (SDS) is added to solutions of low polymer concentrations [34]. It was aslo reported that the viscosity of hydrophobically modified methylcellulose solutions mixed with SDS exhibits similar peaks at low polymer concentrations [35,36]. These non-monotonic rheologies suggest that there is an optimal concentration for the surfactant binding depending on the miscibility of different species of hydrophobes within the same cross-link junction.
In view of these important multi-component gelling polymer solutions, general equation for finding the gel point in terms of their compositions and temperature is desirable. Historically, the gel point conditions of single component trifunctional molecules, and mixture of trifunctional and bifunctional molecules, were studied by Flory [37,38,39,40]. It was later generarized by Stockmayer to reactive molecules of arbitrary functionality for monodisperse [41] and polydisperse [42] cases. These studies (referred to as FS) were based on the fundamental assumptions such that (i) cross-link reactions between functional groups occur pairwisely by chemical (covalent) bonds with equal reactivity, and (ii) intramolecular reactions leading to cyclic structures are prevented in the pre-gel regime. FS was later extended by Fukui and Yamabe [43] to include multiple reactions of the one component functional groups in which cross-link junctions contain arbitrary number of functional groups. Gelation with multiple cross-linking was later treated for thermoreversible gels with physical cross-link junctions by Tanaka and Stockmayer [44,45] (referred to as TS).
On the other hand, gel point condition and molecular weight distribution in the systems of binary multifunctional molecules with co-polymerization between different functional groups were first studied by Stockmayer [11] for the mixtures of polyfunctional A and B molecules. Systematic methods for studying their average properties were later developed by Macosko and Miller [46]. Thermoreversible gelation with binary pairwise cross-links was studied by Tanaka et al. [45,47,48,49,50], and later extended to thermoreversible gels with multiple cross-link junctions [51,52].
It is well known that the decrease in the mixing entropy with growing tree-dimensional molecules drives the solution into a tendency to liquid–liquid phase separation. The relationship between gelation and phase separation has been relatively well understood [1,2,45], but progress continues to be made with some open interesting questions.
In this paper we describe detailed models and their applications to thermoreversible gelation of binary reactive molecules in a common solvent. We present equations for the gel point condition and the phase transition boundaries (spinodal lines) in terms of the microscopic parameters related to the structure of network junctions. Hence we can analyze the structure and strength of the network junctions from macroscopic measurements of phase transitions.

2. Stoichiometric Definition of the Model Solution

The model polymer solution we study here consists of two species of reactive molecules, referred to as R{A f }(A molecule) and R{B g } (B molecule), in a common solvent, mostly water. Molecules can be any type, such as high molecular weight linear polymers, star polymers, or low molecular weight polyfunctional molecules (cross-linkers), etc. An A molecule carries the number f of functional groups A (A groups), and a B molecule carries the number g of groups B (B groups). We can assume g f without losing generality, and refer to R{A f } as primary component, R{B g } as secondary one of the ternary solution. Functional groups are capable of forming complexes (mixed cross-link junctions) with binding numbers controlled by the thermodynamic conditions with equilibrium constants. Specific examples are polymers bearing hydrophobic short chains with different species or different sizes at their chain ends (mixtures of telechelic associating polymers). Hydrophobes form mixed micelles at cross-link junctions. Another examples are star polymers bearing functional groups (ligands) at their arm ends forming complexes with metal cations, etc. In the present study, solvent is assumed to be inactive.
Let n A be the number of statistical repeat units on an A molecule, and n B on a B molecule. The molecular weights of them are then M A = M 0 ( A ) n A and M B = M 0 ( B ) n B , where M 0 ( A ) and M 0 ( B ) are the molecular weights of their statistical repeat units.
To treat concentration, we take the volume of a solvent molecule as the unit of volume, and assume for simplicity that the volume of the statistical repeat units on a functional molecule is the same as that of the solvent molecule, as conventionally assumed in a simple polymer solution theory. The total volume of the solution is then given by Ω = n A N A + n B N B + N 0 . Here, N α is the number of molecules of the type α in the solution. The volume fraction of each component is ϕ A = n A N A / Ω for R{A f }, ϕ B = n B N B / Ω for R{B g }, and ϕ 0 = N 0 / Ω for the solvent. The number concentration of A groups and B groups are then given by ψ A = f ϕ A / n A and ψ B = g ϕ B / n B .

3. Equilibrium Condition for the Formation of Cross-Link Junctions

Let us consider formation of the network junctions J ( k 1 , k 2 ) containing the number k 1 of A groups and k 2 of B groups in the reversible chemical reaction (see Figure 1)
k 1 J ( 1 , 0 ) + k 2 J ( 0 , 1 ) J ( k 1 , k 2 )
Let p k 1 , k 2 be the probability for an arbitrarily chosen A group to belong to a J ( k 1 , k 2 ) junction, and let q k 1 , k 2 be that for a B group. They are the counterparts of the conventional reactivity of the functional groups. Obviously, trivial equality
ψ A p k 1 , k 2 / k 1 = ψ B q k 1 , k 2 / k 2
holds for the number of junctions specifid by ( k 1 , k 2 ) .
For chemical cross-links, the reactivity obeys kinetic equation of the reaction, and hence p k 1 , k 2 is a function of the reaction time. For physical cross-links, we may impose the equilibrium condition for their formation, and in an ideal case the reactivity must fulfil the equation
ψ A p k 1 , k 2 / k 1 ( ψ A p 1 , 0 ) k 1 ( ψ B q 0 , 1 ) k 2 = K k 1 , k 2
for A groups in the reaction (1), where K k 1 , k 2 is the equilibrium constant of the junctions of type ( k 1 , k 2 ) . Hence we have
ψ A p k 1 , k 2 = k 1 K k 1 , k 2 ( ψ A p 1 , 0 ) k 1 ( ψ B q 0 , 1 ) k 2
Therefore, the reaction probability is given by
p k 1 , k 2 = p 1 , 0 k 1 K k 1 , k 2 z A k 1 1 z B k 2
where
z A ψ A p 1 , 0 z B ψ B q 0 , 1
are the concentration of the free functional groups that remain unreacted in the solution. Summing up for all possible k 1 1 and k 2 0 , and using the normalization condition for the probability p k 1 , k 2 , we find
ψ A = z A u A ( z A , z B )
where the function defined by
u A ( z A , z B ) k 1 1 , k 2 0 k 1 K k 1 , k 2 z A k 1 1 z B k 2
characterizes the microscopic structure of a junction containing at least one A group.
Similarly, from the identity (2), we find
q k 1 , k 2 = q 0 , 1 k 2 K k 1 , k 2 z A k 1 z B k 2 1
and hence
ψ B = z B u B ( z A , z B )
for B groups with
u B ( z A , z B ) k 1 0 , k 2 1 k 2 K k 1 , k 2 z A k 1 z B k 2 1
The structure functions u A , u B have physical meanings of the reciprocal unreactivity u A ( z A , z B ) = 1 / p 1 , 0 , u B ( z A , z B ) = 1 / q 0 , 1 .
Coupled Equations (7) and (10) show that the total numbers of A groups and B groups do not change in the process of cross-link formation, but only their connectivity is changed. Therefore in what follows we call them materials conservation law. To study solution properties, these coupled equations must be solved for the two unknown variables z A , z B as functions of the concentration ψ A , ψ B given in the preparation stage of the experiments.
Let us next introduce the reactivity p and q of the functional groups of each species by the definitions p 1 p 1 , 0 and q 1 q 0 , 1 . They have the meaning of conventional reactivities, but extended to multiple reactions, such that for an arbitralily chosen functional group to be in any kind of cross-links. Obviously we have the relation
p = k 1 1 , k 2 1 p k 1 , k 2
q = k 1 1 , k 2 1 q k 1 , k 2
In terms of these reactivities, we can write the materials conservation law as
p = k 1 1 , k 2 1 ( k 1 K k 1 , k 2 ) ψ A k 1 ψ B k 2 ( 1 p ) k 1 ( 1 q ) k 2
q = k 1 1 , k 2 1 ( k 2 K k 1 , k 2 ) ψ A k 1 ψ B k 2 ( 1 p ) k 1 ( 1 q ) k 2
where abbreviated notations k 1 k 1 1 and k 2 k 2 1 have been introduced for simplicity.

4. The Average Branching Number and Response to the Concentration Fluctuations

Let us consider the average structure of cross-link junctions. Using Equations (5) and (9) for the reactivities p k 1 , k 2 , q k 1 , k 2 , the average number of A groups in the cross-link junctions is calculated as
μ ¯ A , A k 1 , k 2 k 1 p k 1 , k 2 = p 1 , 0 k 1 , k 2 k 1 2 K k 1 , k 2 z A k 1 z B k 2 = p 1 , 0 1 + z A z A u A ( z A , z B ) = 1 + κ A , A ( z A , z B )
where the element of κ ^ -matrix is defined by the logarithmic derivative
κ α , β ln u α ln z β
The relation p 1 , 0 u A ( z A , z B ) = 1 has been used.
Similarly, we have
μ ¯ B , B k 1 , k 2 k 2 q k 1 , k 2 = 1 + κ B , B ( z A , z B )
for the average number of B groups in the junctions. Other elements of the μ ^ -matrix are defined by
μ ¯ A , B k 1 , k 2 k 2 p k 1 , k 2 = κ A , B ( z A , z B )
μ ¯ B , A k 1 , k 2 k 1 q k 1 , k 2 = κ B , A ( z A , z B )
although these do not have direct geometrical implication of the junction structure. The determinant of the μ ^ -matrix
D μ ( 1 + κ A , A ) ( 1 + κ B , B ) κ A , B κ B , A
is related to the determinant of κ ^ -matrix
D κ κ A , A κ B , B κ A , B κ B , A
by the equation
D μ = 1 + κ A , A + κ B , B + D κ
Let us next consider infinitesimal variation of the concentrations of unreacted groups in response to the variations d ψ A , d ψ B of the total concentrations. This depends on the average structure of the cross-link junction, and can be found by taking logarithmic derivatives of the materials conservation law (7) and (10). We have
d ln ψ A d ln ψ B = 1 + κ A , A κ A , B κ B , A 1 + κ B , B d ln z A d ln z B
Solving for d ln z A and d ln z B , we find
d ln z A d ln z B = 1 D μ 1 + κ B , B κ A , B κ B , A 1 + κ A , A d ln ψ A d ln ψ B

5. The Weight Average Molecular Weight of the Cross-Linked Clusters and the Gel Point Condition

To study solution properties, in particular gelation and phase separation, let us consider the distribution of cross-linked clusters (three-dimensional polymers) in the solution. Let N l , m be the number of clusters made up of the number l of R{A f } (A molecules) and m of R{B g } (B molecules). Their number concentration is given by
ν l , m N l , m / Ω
The number-average molecular weight of the clusters is defined by
M ¯ n l , m ( M A l + M B m ) N l , m / l , m N l , m
and similarly the weight-average molecular weight is by
M ¯ w l , m ( M A l + M B m ) 2 N l , m / l , m ( M A l + M B m ) N l , m
In our previous paper [51,52], we derived both averages for the mixtures capable of forming any types of multiple cross-links on the basis of classical tree statistics [11,37,38,39,40,41,42,43]. The assumption of pairwise reaction is eliminated to include multiple cross-links made up of any numbers of A and B groups. Under the simplifying assumption that the molecular weight M 0 ( A ) of a statistical repeat unit of A molecules and M 0 ( B ) of B molecules are the same ( M 0 ), the result (Equation (3.22) in the literature [52]) is
ϕ M ¯ w M 0 = n A ϕ A + n B ϕ B + 1 D n A 2 ψ A [ κ A , A ( g 1 ) D κ ] + n B 2 ψ B [ κ B , B ( f 1 ) D κ ] + n A n B D ψ A κ A , B + ψ B κ B , A
where ϕ ϕ A + ϕ B is the total solute volume fraction, κ α , β is the κ ^ -matrix defined above, and D κ is its determinant. The denominator in M ¯ w is defined by
D ( z A , z B ) 1 f κ A , A g κ B , B + f g D κ
(different from D κ , D μ ). This function corresponds to the determinant of the matrix introduced by Gordon [53] in his cascade theory of gelation for mixtures of multi-component reactive molecules. Abbreviated notations f f 1 and g g 1 have been used since they will frequently appear in the following. This result is a generalization of the old Stockmayer’s average molecular weight [11] for the pairwise reaction to multiple ones. At the gel point, the weight average molecular weight goes to infinity, and hence we have
D ( z A , z B ) = 0
for a gel to appear (gel point condition). At the point satisfying this gel point condition, the solution changes from a sol state to a gel state, and vice versa. Such a reversible phase transition is called thermoreversible gelation, or sol-gel transition (referred to as SGT). Materials conservation law (7) and (10), together with the gel point condition (28), leads to the relation between ψ A and ψ B , and therefore gives the SGT line on the ternary phase plane when parameters z A and z B are eliminated.

6. The Average Number of Molecules in the Cross-Linked Clusters and the Spinodal Condition for the Solution

6.1. Chemical Potentials Described in Terms of the Number-Average Degree of Polymerization

To study thermodynamic properties (in particular phase separation) of a polymer solution, we have to derive the chemical potentials of all components in the solution. In our present ternary reactive solutions, we can apply our general theoretical framework [45] for the sudy of associating polymer solutions. We omit its details here by decribing only the outline of derivation.
We start from the free energy of the model solution. It consists of the sum of the two parts
Δ F = Δ cross F + Δ mix F
The first term (per unit volume) is the free energy of cluster formation by cross-linking
Δ cross F / Ω = l , m Δ A l , m ν l , m
where Δ A l , m is the free energy to form an ( l , m ) cluster from the number l of A molecules and m of B molecules in the separated free state.
The second term is the free energy change
β Δ mix F / Ω = l , m ν l , m ln ϕ l , m + ϕ 0 ln ϕ 0 + g ( { ϕ } )
upon mixing the formed clusters with the solvent. In this equation, ϕ l , m ( n A l + n B m ) ν l , m is the volume fraction of the clusters, and β 1 / k B T is the inverse temperature. The interaction term
g ( { ϕ } ) ( χ A , 0 ϕ A + χ B , 0 ϕ B ) ϕ 0 + χ A , B ϕ A ϕ B
gives enthalpy change due to the contact between different species of molecules. The interaction parameters χ α , β ( T ) are Flory’s χ -parameters in the conventional polymer solution theory [40].
By the standard thermodynamic procedure of taking partial derivatives of the free energy, we find the chemical potentials Δ μ l , m for the clusters of ( l , m ) type. We then impose equilibrium condition
Δ μ l , m = l Δ μ 1 , 0 + m Δ μ 0 , 1
for the cluster formation, and find the overall independent values Δ μ A Δ μ 1 , 0 and Δ μ B Δ μ 0 , 1 for each species. Together with the chemical potential for the solvent, we have
β Δ μ A n A = 1 + ln x n A ν S + g A ( { ϕ } )
β Δ μ B n B = 1 + ln y n B ν S + g B ( { ϕ } )
β Δ μ 0 = 1 + ln ϕ 0 ν S + g 0 ( { ϕ } )
Here, x f ν 1 , 0 and y g ν 0 , 1 are the concentration of the functional groups carried by the unassociated free molecules of each species, and
ν S = ν + ϕ 0
is the total translational degree of freedom of solute and solvent molecules, where
ν l , m ν l , m
is the total number of connected clusters in a unit volume of the solution. The sum on the right hand side includes unrected molecules ( 1 , 0 ) and ( 0 , 1 ) . The gel part is excluded because the upper limit goes to infinity for a gel.
The interaction terms are explicitly given by
g A ( { ϕ } ) = χ A , 0 ϕ 0 ( 1 ϕ A , 0 ) χ B , 0 ϕ 0 ϕ B + χ A , B ϕ B ( 1 ϕ A )
g B ( { ϕ } ) = χ B , 0 ϕ 0 ( 1 ϕ B , 0 ) χ A , 0 ϕ 0 ϕ A + χ A , B ϕ A ( 1 ϕ B )
g 0 ( { ϕ } ) = ( χ A , 0 ϕ A + χ B , 0 ϕ B ) ( 1 ϕ 0 ) χ A , B ϕ A ϕ B
The number concentrations x , y of unassociated free molecules can be related to the number concentration z A , z B of unassocated functional groups through the equations
x f ν 1 , 0 = ψ A p 1 , 0 f = ψ A u A ( z A , z B ) f = z A u A ( z A , z B ) f
y g ν 0 , 1 = ψ B q 0 , 1 g = ψ B u B ( z A , z B ) g = z B u B ( z A , z B ) g
Owing to the materials conservation law, x , y can be regarded as functions of the concentration { ϕ } .

6.2. Average Degree of Polymerization and Number Concentration of the Cross-Linked Clusters

Let us next consider the total degree ν of translational motion of clusters, i.e., the concentration of their centers of mass. Since a cluster of the type ( l , m ) includes the number l of A molecules and m of B molecules, the identity
l , m ( l + m ) ν l , m = ϕ A n A + ϕ B n B
holds for the total number of molecules in the pre-gel regime. Hence we have
ν = ϕ A n A + ϕ B n B P ¯ n 1
where
P ¯ n l , m ( l + m ) ν l , m / l , m ν l , m
is the number-average degree of polymerization of cross-linked clusters. Because the number-average molecular weight is defined by Equation (24), P ¯ n can be found by formally substituting M A = M B = 1 in the calculation of M ¯ n . We have reported the latter in our preceding paper [52]. It is
1 M ¯ n = f m A M A 1 f + 1 μ ¯ n ( A ) 1 + g m B M B 1 g + 1 μ ¯ n ( B ) 1
where m A M A N A / ( M A N A + M B N B ) , m B M B N B / ( M A N A + M B N B ) are the mass fractions, and
1 μ ¯ n ( A ) k 1 1 , k 2 0 p k 1 , k 2 k 1 + k 2
1 μ ¯ n ( B ) k 1 0 , k 2 1 q k 1 , k 2 k 1 + k 2
are the reciprocal number-average multiplicity of the cross-link junctions. Setting M A = M B = 1 in these equations, we find
ν = ψ A 1 f + 1 μ ¯ n ( A ) 1 + ψ B 1 g + 1 μ ¯ n ( B ) 1

6.3. Spinodal Condition

The Gibbs condition for a liquid–liquid phase equilibrium (binodal condition) requires a balance of these three chemical potentials in two distinct liquid phases with different concentrations. Because of the many difficulties in solving the coupled equations for binodals, we instead consider the phase boundary of stability limit (spinodal condition) by constructing the Gibbs matrix G ^ from these chemical potentials. Inside of the spinodal line, the solution becomes unstable against infinitesimal concentration fluctuations. The spinodal condition is simply given by a single equation | G ^ | = 0 , and easy to solve. In general, a cloud-point curve observed in an experiment is expected to be identical to the binodal, and may be different from the spinodal, but usually they lie close to each other.
Due to the Gibbs–Dühem relation, two of the three chemical potentials are independent. We choose the solvent as the reference component, and consider A and B as independent components. To avoid cumbersome treatment involving ν S , we make subtraction and consider the difference measured from the reference chemical potential of the solvent. Thus, we have
β Δ μ A β ( Δ μ A Δ μ 0 ) = 1 n A ln x ln ϕ 0 + g A g 0 + const .
β Δ μ B β ( Δ μ B Δ μ 0 ) = 1 n B ln y ln ϕ 0 + g B g 0 + const .
The change of ln x and ln y in response to the infinitesimal variation d ψ A , d ψ B is found by simple differentiation as
d ln x d ln y = 1 f κ A , A f κ A , B g κ B , A 1 g κ B . B d ln z A d ln z B
By using the transition matrix (22) from z to ψ , we find
d ln x d ln y = 1 D μ 1 f κ A , A + κ B , B f D κ f κ A , B g κ B , A 1 g κ B . B + κ A , A g D κ d ln ψ A d ln ψ B
where D κ is the determinant (19) of κ ^ -matrix.
Substituting this relation into the chemical potentials, and together with the equations
d ( g A g 0 ) = 2 χ A , 0 d ϕ A Δ χ d ϕ B
d ( g B g 0 ) = Δ χ d ϕ A 2 χ B , 0 d ϕ B
for the mixing enthalpy, with Δ χ χ A , 0 + χ B , 0 χ A , B , we finally find
d ( β Δ μ A ) d ( β Δ μ B ) = 1 + κ B , B f ( κ A , A + D κ ) n A ϕ A D μ + 1 ϕ 0 2 χ A , 0 f κ A , B n A ϕ B D μ + 1 ϕ 0 Δ χ g κ B , A n B ϕ A D μ + 1 ϕ 0 Δ χ 1 + κ A , A g ( κ B , B + D κ ) n B ϕ B D μ + 1 ϕ 0 2 χ B , 0 d ϕ A d ϕ B
The first matrix factor on the right hand side is the Gibbs matrix G ^ .
After a lengthy calculation of its determinant, we find the spinodal condition as
| G ^ | = K A K B f g κ A , B κ B , A n A n B ϕ A ϕ B D μ 2 + 1 ϕ 0 D μ K A n A ϕ A + K B n B ϕ B + f κ A , B n A ϕ B + g κ B , A n B ϕ A 2 χ A , B D μ   2 D μ χ B , 0 K A n A ϕ A + χ A , 0 K B n B ϕ B + Δ χ 2 f κ A , B n A ϕ B + g κ B , A n B ϕ A χ ˜ = 0
where
K A = ( 1 + κ B , B ) ( 1 f κ A , A ) + f κ A , B κ B , A
K B = ( 1 + κ A , A ) ( 1 g κ B , B ) + g κ A , B κ B , A
These are functions of z A and z B , and hence functions of reactivities p and q. In what follows we therefore write the spinodal condition as
S ( p , q ) | G ^ ( z A , z B ) | = 0
The ternary effective interaction parameter in the last term of (50) is defined by
χ ˜ χ A , B 2 + χ A , 0 2 + χ B , 0 2 2 χ A , B χ A , 0 2 χ A , B χ B , 0 2 χ A , 0 χ B , 0
From the materials conservation law, we can find p and q as functions of the volume fraction ϕ A and ϕ B . Then the spinodal lines can be drawn on the ternary phase plane.
We can easily confirm that (50) reduces to the conventional Flory–Huggins equation (XIII-1d of [40]) for ternary spinodal condition
S ( p , q ) = 1 n A n B ϕ A ϕ B + 1 n A n 0 ϕ A ϕ 0 + 1 n B n 0 ϕ B ϕ 0 2 χ B , 0 n A ϕ A + χ A , 0 n B ϕ B + χ A , B n 0 ϕ 0 χ ˜ = 0
by setting all κ elements to be zero.
In the special case of no solvent, we take the limit ϕ 0 0 in (50), and find
K A n A ϕ A + K B n B ϕ B 2 χ A , B = 0
with
K A K A + g n A n B κ B , A / D μ
K B K B + f n B n A κ A , B / D μ
This is the standard form of spinodal condition for binary associating solutions we have studied over the past decades [45].

6.4. Treatment of the Temperature

Although exploring the effect of temperature is not the main purpose of this study, we briefly summarize how to treat temperature within the present theoretical framework. Temperature comes in through the conventional interaction parameters χ α , β and the new parameters K k 1 , k 2 (equilibrium constants for cross-link formation). We regard A molecules as the primary component and B molecules as the secondary component because the latter is often an additional component mixed into the solution of the primary one. Therefore, we assume the conventional Shultz–Flory form
χ A , 0 ( T ) = 1 / 2 ψ ( 1 Θ / T )
for the interaction between A molecules and solvent. Here, ψ is a material parameter of order unity, and Θ is the theta temperature of the primary solution. Because Θ does not include the effect of cross-linking, it is referred to as unperturbed theta temperature. The equilibrium constants can be in principle described in terms of the binding constants
λ ( T ) = λ 0 exp ( ϵ A / k B T ) , μ ( T ) μ 0 exp ( ϵ B / k B T )
of functional group of each species in the cross-link junctions (expamples given below), where λ 0 , μ 0 are the entropy contribution and ϵ A , ϵ B are the enthalpy of binding.

7. Fixed Multiplicity Model

Let us first consider the simplest model junction whose numbers k 1 and k 2 are fixed (fixed multiplicity model). Functional groups in the solution are therefore either in an unreacted state (1,0), (0,1), or belong to a junction of ( k 1 , k 2 ) type. For special case of pairwise reaction ( k 1 , k 2 ) = ( 1 , 1 ) , this model reduces to the conventional chemical (covalent) bonds for which molecular distribution function of the cross-linked polymers and the gel point condition were studied by Stockmayer [11].
For the fixed multiplicity model, the number of cross-link junctions formed in the solution is ψ A p / k 1 as seen from A groups, and ψ B q / k 2 as seen from B groups. Since they are the same, we have a relation
q = r R p
where
R ψ B / ψ A
is the relative concentration of the number of B groups per an A group, and
r k 2 / k 1
is the multiplicity ratio. The relative concentration R is often used for the gelation of polymers (A) induced by added agencies (B), such as metallic cations, basic anions, surfactants, and other cross-linkers. The stoichiometric concentration is given by the condition R = r , and hence at this concentration the reactivities of the functional groups take the same value p = q . In what follows we use the scaled relative concentration
ξ R / r
in order to see the physical meanings of the concentrations easily. For the low concentration region 0 < ξ 1 , the reactivities change in the range 0 < p ξ , 0 < q 1 , while for the high concentration region 1 < ξ , they change in the region 0 < p 1 , 0 < q 1 / ξ . The stoichiometric concentration is given by ξ = 1 .

7.1. Sol-Gel Transition (SGT) and Spinodals for the Fixed Multiplicity Model

For the fixed multiplicity model, the materials conservation law takes the form
p = ( k 1 K ) ψ A k 1 ψ B k 2 ( 1 p ) k 1 ( 1 q ) k 2
q = ( k 2 K ) ψ A k 1 ψ B k 2 ( 1 p ) k 1 ( 1 q ) k 2
where only one equilibrium constant K K k 1 , k 2 is necessary. Eliminating q by using the stoichiometric relation (59), we find the equation
F ( p ) p ( k 1 K ) ψ A k r k 2 ( 1 p ) k 1 ( ξ p ) k 2 = 0
which should be solved for the unknown variable p as a function of the controlled parameters ψ A and ξ fixed in the preparation stage of the solution. The new multiplicity sum
k k 1 + k 2 1
has been introduced.
We next find the gel point condition. Simple calculation leads to
κ ^ = k 1 p , k 2 p k 1 q , k 2 q
and hence
D κ = k p q
and
D μ = 1 + k 1 p + k 2 q k p q
From these results we find
D ( p , q ) = 1 f k 1 p g k 2 q f g k p q = 0
for the gel point.
The K functions for the study of spinodal condition are then found to be
K A = 1 f k 1 p + k 2 q + f k p q
K B = 1 g k 2 q + k 1 p + g k p q
The spinodal condition is also found in the form of S ( p , q ) = 0 by using these specific forms.
Since the gel point condition (69) is a second order algebraic equations for p, the solution can be found easily as
p = p c ( ξ ) = f k 1 ξ + g k 2 2 f g k 1 + 1 + 4 f g k ξ ( f k 1 ξ + g k 2 ) 2
and q = q c ( ξ ) = p c ( ξ ) / ξ . The concentration of A groups at the gel point is then given by
ψ A , c = p c ( ξ ) ( k 1 K ) r k 2 ( 1 p c ( ξ ) ) k 1 ( ξ p c ( ξ ) ) k 2 1 / k
This result of an SGT line on the ( ψ A , ξ ) plane can be mapped onto the ternary phase plane ( ϕ A , ϕ B , ϕ 0 ) of the volume fraction in the following two ways in accordance with the methods by which the solutions are prepared in the experiments.
The first method of preparation is mixing a solution of R{A f }/S and a solution of R{B g }/S with a common concentration (volume fraction) ϕ with the volume ratio 1:u. The result leads to
ϕ A = ϕ ( 1 u ) , ϕ B = ϕ u , ϕ 0 = 1 ϕ
where ϕ ϕ A + ϕ B is the total solute volume fraction. We then have
R = α u / ( 1 u ) , ξ = β u / ( 1 u )
where α g n A / f n B depends only on the functionality and molecular weight of the solute molecules, and β α / r = g n A k 2 / f n B k 1 . (Do not confuse with 1 / k B T .) Since ξ depends only on u, we find ϕ c = ϕ c ( u ) from (72).
The second method of preparation is to add B molecules R{B g }, or its solution R{B g }/S of a known concentration, to a pure R{A f }/S solution, as in titration experiments. In particular for adding B molecules, we have
ϕ A = ϕ ( 1 u ) , ϕ B = u , ϕ 0 = ( 1 ϕ ) ( 1 u )
and hence
R = α u / ϕ ( 1 u ) , ξ = β u / ϕ ( 1 u )

7.2. Optimal Gelation Concentration

We have developed a very general strategy to study the relationship between microscopic structure of the cross-link junctions and macroscopic SGT in ternary polymer solutions. As an example of its practical applications, we try to find a minimum concentration of the functional A molecules, mostly functional polymers, necessary for driving the solutions into gel when cross-linking agency is added. If a minimum exists, it provides the most efficient concentration to obtain gelling solution from a viewpoint of materials designing.
To study this problem of optimizing gelation, let us examine the gel point condition (69). For an infinitesimal change d ψ A and d ψ B , the reactivities changes along SGT line under the condition
d D = f ( g k q + k 1 ) p d ln p + g ( f k p + k 2 ) q d ln q = 0
On the other hand, from the materials conservation law (63a) and (63b), we find
1 + k 1 p 1 p , k 2 q 1 q k 1 p 1 p , 1 + k 2 q 1 q d ln p d ln q = k 1 , k 2 k 1 , k 2 d ln ψ A d ln ψ B
Solving for d ln p and d ln q , we find
d ln p d ln q = ( 1 p ) ( 1 q ) 1 + k 1 p + k 2 q k p q k 1 k q 1 q , k 2 1 q k 1 1 p , k 2 k p 1 p d ln ψ A d ln ψ B
Upon substitution into minimum SGT condition (77), we find
H ( p , q ) H A ( p , q ) d ln ψ A + H B ( p , q ) d ln ψ B = 0
where
H A ( p , q ) f p ( g k q + k 1 ) k 1 k q 1 q + g q ( f k p + k 2 ) k 1 1 p
H B ( p , q ) f p ( g k q + k 1 ) k 2 1 q + g q ( f k p + k 2 ) k 1 k p 1 p
To find the concentration of B molecules which minimizes the concentration of A molecules along the SGT line, we fix d ln ψ A = 0 in this equation, and have
H B ( p , q ) = 0
The solution ( p * , q * , ξ * ) of the coupled Equations (64), (69) and (82) gives the optimal SGT.

7.3. Specific Models of the Cross-Link Junctions with Fixed Multiplicity

7.3.1. Pairwise Cross-Links

Let us start from the simplest case of pairwise reaction ( k 1 , k 2 ) = ( 1 , 1 ) . We have r = 1 , k = 1 , ξ = R . The elements of κ ^ -matrix are κ A , A = κ B , B = 0 , κ A , B = p , κ B , A = q , and therefore lead to D κ = p q and K A = 1 + f p q , K B = 1 + g p q . Obviously, the gel point condition reduces to the monodisperse case of Stockmayer’s result [11]
D ( p , q ) = 1 f g p q = 0
The optimal condition (82) then gives
p * = f g + 1 2 f g
From two other equations we find
ψ A * = 4 f g K ( f g 1 ) 2 , R * = ( f g + 1 ) 2 4 f g
Gel formation is limited to a certain concentration region between R min and R max , where R min = 1 / f g and R max = f g . It is natural that we need a minimum amount of B component (lower bound R min ) to drive the solution into gel, but a certain condition is necessary to see that the gel turns back to sol for excess B component at high B concentration (called reentrant sol-gel-sol transition, referred to as SGST). Under excess B molecules, A molecules are covered by B molecules and separated from each other rather than forming connected clusters because reactions among B groups ( ( 0 , 2 ) etc.) are prohibited (see Figure 2).
In the following study, we attempt to find, in terms of the multiplicity of the cross-links, the necessary condition for the appearance of an optimal SGT point followed by a reentrant sol phase. The optimal SGT is practically important because at this point the solution can turn into a gel with a minimum amount of A molecules.
The spinodal condition for the fixed multiplicity model is given by (50) with K A , K B of Equations (70a) and (70b). We find that it takes the form
1 + ( f + g ) p q n A n B ϕ A ϕ B ( 1 p q ) 2 + 1 ϕ 0 ( 1 p q ) 1 + ( f p + n A g / n B ) q n A ϕ A + 1 + ( g q + n B f / n A ) p n B ϕ B 2 χ A , B ( 1 p q ) 1 1 p q 2 χ A , 0 1 + f p q n A ϕ A + 2 χ B , 0 1 + g p q n B ϕ B + Δ χ f p n A ϕ B + g q n B ϕ A χ ˜ = 0
for pairwise association.
In the particular case of no solvent, this reduces to
1 + ( f p + n A g / n B ) q n A ϕ A ( 1 p q ) + 1 + ( g q + n B f / n A ) p n B ϕ B ( 1 p q ) 2 χ A , B = 0
If we regard B molecules as monofunctional molecules such as hydrogen-bonding water, low molecular weight surfactant, etc., we can set g = 0 in this equation, and find our previous results [45,54,55] for the spinodal
1 n A ϕ A + 1 + [ f n B / n A + ( f n B / n A ) 2 ] p n B ϕ B ( 1 p q ) 2 χ A , B = 0
with q = ( f n B ϕ A / g n A ϕ B ) p .
Figure 3a,b show how SGST and spinodals change on the ternary phase plane with the association constant K and the interaction parameters χ . In contrast to the conventional triangular phase plane, we use throughout this paper two new independent concentration variables, the volume fraction of ϕ A of the A molecules (vertical axis) and the concentration of the B functional groups ξ R / r = β ϕ B / ϕ A relative to that of A functional groups (horizontal axis). The stoichiometric concentration of the functional groups is given by ξ = 1 . The origin of the coordinate corresponds to the pure solvent, the point ( 0 , 1 ) to pure A component, and ( , 0 ) to pure B component. The phase region is limited from above by the condition ϕ 0 = 1 ( ϕ A + ϕ B ) = 0 (thin line). The main reason for using this phase plane is that in the experiments B molecules are often the third components added to the already prepared solutions of A molecules. However, we can use it for any ternary solutions instead of the conventional triangular plane without losing generality.
In genaral, the association constant K changes with solvent quality, so that it depends on the interaction parameters. The change is mainly caused by the exclusion of solvent molecules from the junction zones upon forming cross-linking bonds, leading to the change in the number of molecular contact. In the present numerical study, however, we have fixed it for simplicity independently of the χ parameters. Further scrutiny on this point from physical viewpoint will provide a clue to clarify the dependence of thermoreversible gelation on the solvent quality.
Figure 3a shows how the spinodal lines shift with increase of χ A , B , the interaction parameter between A- and B molecules, while the association constant is fixed, for a symmetric case of f = g = 3 and n A = n B = 5 . In this test calculation, we have fixed K = 30 and changed χ A , B from 2.0 to 5.0. The sol region (below the red solid line) and the gel region (above it) are separated by a nonmonotonic line with a minimum, so that the sol-gel transition is a reentrant SGST. The minimum point is specifically given by the coordinates (25/16,8/81) = (1.56,0.099) according to Equation (85). With increase in χ A , B the phase-separated two-phase region (the region above the broken line) expands because repelsive force between the two components becomes stronger.
Figure 3b shows the same as (a), but K is changed from 30 to 100 while χ A , B is fixed at 5.0. We can see how the gel region expands with K, with flatter bottom part. The spinodals also shift with K because the repulsive force is compensated by increasing association [49].
Figure 4 shows the same as Figure 3a but here the interaction parameters between components are different. For a solvent which is poor for A but good for B, we have fixed as χ A , 0 = 2.0 , χ A , B = χ B , 0 = 0 (blue broken line), and for an opposite case, we have fixed as χ B , 0 = 2.0 , χ A , B = χ A , 0 = 0 (green broken line). In the former case, the solution starts from a two-phase state, whose phase-saparated region shrinks with increase in the concentration of B molecules, and eventually merges into one phase. The latter case shows an opposite tendency.
The turning points of the spinodal curves in the phase diagram, and in all phase diagrams presented in the following, lie close to the crititical point of phase separation. However, they do not exactly correspond to the critical points because our ternary co-ordinates are different from the conventional triangular ones.

7.3.2. Mononuclear A Cross-Links

In order to study the reentrant SGST in more detail, we next consider cross-links containing only one A functional group (functional group carried by the molecules of higher functionality). It is an example of mononuclear cross-links since it has only one A group (referred to as nucleus) in a junction. We have r = k 2 , k 1 = 0 , k = k 2 and K A = 1 + ( k 2 + f k 2 p ) q , K B = 1 g ( k 2 k 2 p ) q , D μ = 1 + ( k 2 k 2 p ) q . The gel point condition leads to
p c ( ξ ) = k 2 2 f k 2 1 + 1 + 4 f k 2 ξ g k 2 2
for the critical reactivity. The SGST line on the concentration plane takes the form
ψ A , c ( ξ ) = p c ( ξ ) ( k 2 K ) ( 1 p c ( ξ ) ) ( ξ p c ( ξ ) ) k 2 1 / k 2
The optimal SGST condition becomes
H B ( p ) = 2 f g ( k 2 p ) 2 + [ g ( f k 2 + f k 2 ) + 1 ] ( k 2 p ) + k 2 ( g k 2 1 ) = 0
where f f 2 , and hence
k 2 p * = g ( f k 2 + f k 2 ) + 1 2 f g 1 + 1 + 4 f g k 2 ( g k 2 1 ) [ g ( f k 2 + f k 2 ) + 1 ] 2
The lower bound of the gel region can be found by putting q c ( ξ ) = 1 . We find ξ min = 0 except a special case of k 2 = 1 (pairwise cross-link studied above). Since the number concentration of A groups is limited in the region below the maximum value K f / n A corresponding to the volume fraction ϕ A = 1 , ξ min is practically limited by this finite bound. Gelation becomes reentrant one (see Figure 5. The reentrant part lies beyond the region shown in the graph). The upper bound of the gel region can be found by setting p c ( ξ ) = 1 , from which condition we have
ξ max = g ( f k 2 + k 2 )
Figure 5a shows SGST line (red solid line) and spinodal lines for three cases of the interaction parameters as above (broken lines) for a fixed association constants of a symmetric case of f = g = 3 and n A = n B = 5 . In this test calculation, we have fixed K = 20 and changed the quality of the solvent. The sol region (below the red solid line) and the gel region (above it) are separated by a line with a minimum (lying beyond the region shown in the graph). Figure 5b shows the reactivity p and q as functions of the scaled B concentration for three fixed volume fractions of A molecules ϕ A = 0.05 , 0.10 , 0.15 . They all start from 0. The reactivity of B groups shows a maximum near the stoichiometric concentration.

7.3.3. Mononuclear B Cross-Links

We next consider model cross-links containing only one B group (functional group carried by the molecules of lower functionality), referred to as mononuclear B cross-links. We have k 2 = 0 , k = k 1 and
K A = 1 f ( k 1 + k 1 q ) p , K B = 1 + ( k 1 + g k 1 q ) p , D μ = 1 + ( k 1 k 1 q ) p
The gel point condition becomes
D ( p ) = ξ f k 1 ξ p f g k 1 p 2 = 0
so that the reactivity of A groups along the SGT line is given by
p c ( ξ ) = k 1 ξ 2 g k 1 1 + 1 + 4 g k 1 f k 1 2 ξ
The concentration of A groups along this SGT is
ψ A , c ( ξ ) = p c ( ξ ) K ( 1 p c ( ξ ) ) k 1 ( ξ p c ( ξ ) ) 1 / k 1
Therefore, the gel region is limited from below by
ξ min = 1 / f ( g k 1 + k 1 )
but there is no upper bound except the special case of k 1 = 1 (pairwise cross-link studied above). The SGT decreases monotonically with the concentration of B groups (see Figure 6). Mixing B molecules monotonically promotes gelation.
For special case of k 1 = 2 , the cross-link takes the -ABA- form. Let us call it sandwich cross-link. Poly(vinyl alcohol) and polysaccharides cross-linked by borax ions fall on this categoly. In fact, phase diagrams similar to this one were already reported for garactomannan/borax solutions [21]. There is a uniform clear gel region at low concentration of B molecules before the solution becomes turbid. However, there is a difference such that monocomplexes of the type -AB may coexist in such ion-binding cross-links [26] while in the present fixed multiplicity model they are excluded. SGT line may have a minimum outside of the borax concentration region measured in the experiments [21] although it is not evident in the reported phase diagram.
In Figure 6, incidentally f , n A and g , n B are the same for these examples. Therefore, if vertical coordinate ϕ A is replaced by ϕ B and horizontal coordinate ξ = β ϕ B / ϕ A is replaced by ξ = β ϕ A / ϕ B , Figure 5 (mononuclear A) and Figure 6 (mononuclear B) give the same phase diagram. However, because the concentration region covered is completely different, Figure 5 shows A-solvent edge and Figure 6 shows B-solvent edge of the conventional triangular diagram.

7.3.4. Monofunctional B Molecules

In some important gelling polymer solutions, the added cross-linking molecules are monofunctional molecules ( g = 1 ). Metal ions forming coordination complexes with polymers carrying various ligands such as acrylic groups, hydroxy groups etc. fall onto this category. Surfactant molecules mixed with hydrophobically-modified water-soluble polymers (associating polymers) are also regarded as a typical example of such gelling solutions, because a surfactant molecule usually carries only one hydrophobe. We therefore here summarize the nature of SGT for the special case of g = 1 . The functionality f and the junction multiplicity k 2 can be arbitrary, but k 1 is assumed to be k 1 2 in order to have gels. We have
K A = 1 f k 1 p + k 2 q + f k p q , K B = 1 + k 1 p , D μ = 1 + k 1 p + k 2 q k p q
Because g = 0 , the gel point condition gives a constant reactivity
p c = 1 / f k 1
along the SGT line, and hence the concentration of A groups along it is
ψ A , c ( ξ ) = f k 1 ( k 1 K ) r k 2 ( f k 1 1 ) k 1 ( f k 1 ξ 1 ) k 2 1 / k
There is only a lower bound ξ min = 1 / f k 1 . The SGT line monotonically decreases. The gelation becomes easier in proportion to the concentration of added cross-linkers (see Figure 7).
In our previous study on ion-binding polymer solutions [26], we showed the existence of an optimal gelation point where a polymer solution turns into gel at the minimum polymer concentration when metal ions, regarded as g = 1 , are added. The appearance of such an optimal gel point was attributed to the effect of monocomplexes, i.e., complexes including only one ligand. From the viewpoint of ions, a monocomplex is formed by the adsorption of an ion onto one of the ligands carried by polymer chains without forming a cross-link. In the present case of g = 1 , it corresponds to ( 1 , 1 ) cross-link junction. If ( 1 , 1 ) is allowed in addition to ( 2 , 1 ) , most of A groups are covered by B groups at their high concentration region, so that connected clusters are separated into smaller ones. In the fixed multiplicity model studied in this section, however, formation of ( 1 , 1 ) junction is excluded because only one type of cross-links with fixed multiplicity ( k 1 , k 2 ) for k 1 2 is allowed. In other words, if cross-link junctions of ( 1 , k 2 ) are allowed to coexist with those of k 1 2 , an optimal gel point may appear. Details of such a variable multiplicity problem will be studied in the following sections.
Figure 7 shows the phase diagram of polyfunctional molecules ( f = n A = 10 ) cross-linked by monofunctional small B molecules ( n B = 1 , g = 1 ) by a sandwich type ( 2 , 1 ) multiplicity. The interaction parameter χ A , B is varied (blue broken line) in a good solvent. SGT is monotonic with the concentration of the cross-linkers. Clear gel appears first, and then phase separation takes place. The clear gel region shrinks with increase in χ A , B . This type of phase diagrams are also frequently observed in the experiments.

8. Variable Multiplicity Model

Let us next study gelation with simultaneous formation of cross-links of different type. The multiplicity index ( k 1 , k 2 ) are allowed to vary under the thermodynamic equilibrium condition determined by the intermolecular interactions among the functional groups. In what follow, we consider several important cases frequently observed in experiments.

8.1. Completely Immiscible Cross-Link Junctions

In this case, there is no association between functional groups of different species. Only cross-links of the pure type ( k 1 , 0 ) and ( 0 , k 2 ) are allowed. The polymer solution turns into an interpenetrating polymer networks (IPNs) when concentrations of both A and B groups exceed the critical values. Let us assume that multiplicities k 1 , k 2 are variable. Then, the IPN considered here is an extension of the conventional chemical IPN to more general physical gels with multiple junctions.
We have
u A ( z A ) = k 1 1 k 1 K k 1 , 0 z A k 1
u B ( z B ) = k 2 1 k 2 K 0 , k 2 z B k 2
Hence,
κ A , A = ln u A ( z A ) ln z A
κ B , B = ln u B ( z B ) ln z B
and κ A , B = κ B , A = 0 . For spinodal conditions, K functions are given by
K A = ( 1 + κ B , B ) ( 1 f κ A , A )
K B = ( 1 + κ A , A ) ( 1 g κ B , B )
and
D μ = ( 1 + κ A , A ) ( 1 + κ B , B )
The gel point condition is factorized as
D ( z A , z B ) = { ψ A f z A 2 u A ( z A ) } { ψ B g z B 2 u B ( z B ) } = 0
where abbreviated notation u ( z ) d u ( z ) / d z has been used. Accordingly, the concentration plane is separated into four regions of SOL, A-GEL, B-GEL and IPN(A/B-GEL).
To find SGT of each component, let us introduce association constants λ ( T ) and μ ( T ) for A- and B groups as in our previous study on one component thermoreversible gels [44,47]. They depend on the temperature T. We then have
K k 1 , 0 = γ k 1 ( A ) λ ( T ) k 1 , K 0 , k 2 = γ k 2 ( B ) μ ( T ) k 2
with dimensionless numerical constants γ k 1 ( A ) and γ k 2 ( B ) . The gel condition becomes
D ( z A , z B ) = { a f z A 2 ϕ 2 ( A ) ( z A ) } { b g z B 2 ϕ 2 ( B ) ( z B ) } = 0
with polynomials
ϕ 2 ( A ) ( z A ) k 1 2 k 1 k 1 γ k 1 ( A ) z A k 1
ϕ 2 ( B ) ( z B ) k 2 2 k 2 k 2 γ k 2 ( B ) z B k 2
and scaled concentrations
a = λ ( T ) ψ A
b = μ ( T ) ψ B
and the scaled variables z A , z B instead of λ z A , μ z B . If the multiplicities k 1 and k 2 are fixed, the gel point condition, together with the materials conservation law, gives two solutions
a = f k 1 ( γ ( A ) ) 1 / k 1 ( f k 1 1 ) k 1 / k 1
b = g k 2 ( γ ( B ) ) 1 / k 2 ( g k 2 1 ) k 2 / k 2
corresponding to the SGT line of A-GEL and B-GEL. For pairwise case, we already reported on the thermoreversible IPN in detail [48,49]. An example IPN with multiple cross-linking is presented in Figure 8.
Let us briefly summarize the structure of binary gels after passing the gel pint. For binary gels, there are four regions in the post-gel regime: A-A connected region (referred to as A A ) which has infinite paths on the gel connected by a homo A-A reaction, B-B connected region (B B ) with infinite paths connected by a homo B-B reaction, A-B connected region (A B ) with infinite paths connected by co-reaction only, and A-A/B-B connected region (A B ) with infinite paths connected by both A-A reaction and B-B reaction: bi-continuous gel). The simplest example of A B is presented in Figure 8 for immiscible cross-links. The network is bi-continuous but separated with each other, and therefore it is called interpenetrating polymer network (IPN). More detailed analysis of the binary post-gel properties will soon be reported.

8.2. Completely Miscible Cross-Link Junctions

We next consider the opposite case where A- and B groups are similar and completely miscible with each other within a cross-link junction. Typical examples of such cases are associating polymers carrying many hydrophobes (large f) mixed with surfactant molecules ( g = 1 ) [35,56], and also mixed with telechelic polymers ( g = 2 ) (polymers carrying hydrophobes at both chain ends)[57]. If hydrophobes on both species are the same, they can completely mix with each other within miceller cross-link junctions formed by their aggregation [56].
For such a completely miscible junction, we can assume that k 1 , k 2 may vary under a total multiplicity k 1 + k 2 = k which can also vary in a range fixed by the physical nature of interaction. We therefore have
K k 1 , k 2 = γ k k ! k 1 ! k 2 ! λ ( T ) k 1 μ ( T ) k 2 δ ( k 1 + k 2 k )
for the equilibrium constant for k 1 1 , where λ ( T ) , μ ( T ) are the binding (association) constants of A- and B group, and γ k is a numerical constant. In what follows, we use the symbols z A , z B instead of λ z A , μ z B for simplicity. The binding polynomial for A groups takes the form
u A ( z A , z B ) = 1 + k 2 γ k k 1 = 1 k k ! k 1 ! k 2 ! k 1 z A k 1 z B k k 1 = 1 + k 2 k γ k ( z A + z B ) k = k 1 k γ k ( z A + z B ) k
Similarly, we have
u B ( z A , z B ) = k 1 k γ k ( z A + z B ) k
for B groups. Since u A = u B , we write it simply as u ( z A , z B ) . The materials conservation law is given by
a = z A u ( z A , z B )
b = z B u ( z A , z B )
in terms of a λ ( T ) ψ A , b μ ( T ) ψ B . Taking the sum, these coupled equations reduce to the one
ψ = z u ( z )
where ψ a + b , and z z A + z B . This relation is apparently the same as the materials conservation law for one component gels with multiple junctions studied in TS [44]. We have only to replace the concentration of the functional group by the weighted sum of the concentration of each component, and z variable by the sum of z for each component. By using the solution z of this equation, we can find z A = ( 1 w ) z and z B = w z , where
w b / ( a + b )
is the weighted number fraction of B groups among the total functional groups.
Now, from a simple calculation of κ ^ -matrix, we find the condition
D ( z ) 1 [ f ( 1 w ) + g w ] κ ( z ) = 0
for SGT with one-component κ -function
κ ( z ) d ln u ( z ) d ln z
Because the prefactor f ( 1 w ) + g w can be regarded as weight-average functionality f w f w 1 , where
f w ( ξ ) f ( 1 w ) + g w = f + g ξ 1 + ξ
the problem has entirely reduced to the gelation of one component functional molecules. In what follows, we refer to this fact as mixing law.
In order to study the existence of an optimal gelation point, we derive an additional condition d D ( z ) = 0 . Together with the relation d ln ψ = [ 1 + κ ( z ) ] d z / z obtained from the materials conservation law, we find
H B ( z ) κ 2 ( z ) ( 1 + ξ ) ( f g ) κ ( z ) [ 1 + κ ( z ) ] = 0
for the appearance of the minimum in ψ A along SGT line, where a new κ 2 -function is defined by the logarithmic derivative of κ ( z )
κ 2 ( z ) d ln κ ( z ) d ln z
By using u ( z ) only, we can write this condition explicitly as
H B ( z ) = g u ( z ) u ( z ) 2 + g z u ( z ) 3 + u ( z ) 2 u ( z ) = 0
where abbreviated notations g g 3 , u ( z ) d 2 u ( z ) / d z 2 have been used.
For the spinodal condition, we have
K A = 1 + [ w f ( 1 w ) ] κ ( z )
K B = 1 + [ 1 w g w ] κ ( z )
and
D μ = 1 + κ ( z )

8.2.1. Pairwise Cross-Links

We first study the simplest case of pairwise cross-linking. We have only ( 2 , 0 ) (A 2 cross-link), ( 1 , 1 ) (AB cross-link), and ( 0 , 2 ) (B 2 cross-link). In general we have randomly mixed polymer networks if either of the functionality is larger than (or equal to) 3. This case was studied in detail in our previous paper [55]. If B molecule is monofunctional ( g = 1 ), it serves as a hinderer for gelation, i.e., gelation of A molecules is retarded, or prevented when B molecules are added to the solution (monotonically increasing ψ A , c ). If B molecule is bifunctional ( g = 2 ), generated networks have long AA subchains between neighboring cross-links [40]. In both cases, there is no optimal gel point.

8.2.2. Fixed Multiplicity Model

We therefore move onto the multiple junctions. Let us first study a fixed multiplicity case of k ( 3 ) , so that we have
u ( z ) = 1 + γ k z k , ( k k 1 , γ γ k )
From the gel point condition, we find
γ k z c ( ξ ) k = 1 / [ f w ( ξ ) k 1 ]
and hence
λ ψ A , c = 1 ( γ k ) 1 / k f w ( ξ ) k ( 1 + ξ ) [ f w ( ξ ) k 1 ] k / k
Simple analysis of the condition (123) for optimal gel point leads to the conclusion that for monofunctional B molecules ( g = 1 ) it has a solution
z * = k γ k 2 1 / k
( k k 2 ) and therefore
ξ * μ λ R * = 1 2 f k 1
but for multifunctional B molecules ( g 2 ) there is no solution. One of the most interesting examples of such monofunctional B systems is the solutions of associating polymers mixed with surfactant molecules. There have been accumulated experimental data indicating the existence of a maximum in viscosity, and also in plateau modulus, as functions of the surfactant concentration [34,35]. Such maxima occur at the optimal concentration ξ * of surfactants for a given concentration of associating polymers.

8.2.3. Mini-Max Model

Let us generalize the model to the multiplicity k allowed to vary in a certain interval between k min and k max . We therefore have
u ( z ) = 1 + d ϕ ( z ) d z
where
ϕ ( z ) k = k min k max z k = z k min 1 z β 1 z
where β k max k min + 1 . Gelation of such a mini-max junction was already studied for one-component networks [44]. We don’t repeat it here but instead focus on the appearance of the optimal gel point. Obviously g 2 is necessary for the condition (123) to have a solution, but g = 2 leads to a monotonic SGT line as studied above. Therefore B molecule must be a monofunctional one ( g = 1 ). A straightforward but lengthy calculation of H B ( z ) shows that the minimum multiplicity k min must be larger than 3. If pairwise junction k min = 2 is allowed, B molecules attach on the functional groups of A molecules to cover the surface, and prevent gelation. As a result, the gel concentration of A molecules monotonically increases. Speciality of k min = 2 was already pointed out in our previous paper [56]. Upper bound in k has only weak effect so that the separation of cross-linked clusters into smaller pieces by the coverage of B molecules (monocomplexes) is not prevented by k max .
Figure 9 shows an example of the phase diagram with k min = 3 , k max = 5 for a mixture of trifunctiona A molecules ( f = 3 ) and bifunctional B molecules ( g = 2 ). The association constant is fixed at λ = μ = 2.0 , while spinodal lines of two cases of the solvent χ A , B = 5.0 , χ A , 0 = χ B , 0 = 0 (green broken line) and χ A , 0 = 0.9 , χ A , B = χ B , 0 = 0 (blue broken line) are drawn.
Figure 10a shows phase diagram of telechelic associating polymers ( f = 2 , n A = 100 ) mixed with surfactant molecules ( g = 1 , n B = 1 ) with k min = 5 , k max = 20 . The association constant is fixed at λ = μ = 80 , while three spinodal lines for the solvent χ A , 0 = 0.60 , 0.603 , 0.7 (poor for A molecules) with χ A , B = χ B , 0 = 0 (blue broken lines) are drawn. Below χ A , 0 = 0.603 (weakly poor solvent for the telechelic polymers), the two-phase region has a closed shape with a lower boud of B concentration. Because such closed-loop two-phase regions take place by adding cross-linking molecules which are perfectly miscible to the solvent, they are caused not by the enthalpy change but by the decrease in the mixing entropy of molecules due to association.

8.3. Partially Miscible Cross-Link Junctions

In mixtures of associating polymers and surfactants, the nature of the hydrophobes they carry may sometimes very different; with different species, different sizes, different shapes, etc. They are therefore only partially miscible within the miceller junctions. Another example of junctions with partial miscibility is formation of outer-sphere complexes in ion-binding cross-links. In this example, a metal ion forms a complex with fixed number of ligands ( ( k 1 , 1 ) type primary complex), and such primary complexes form higher order secondary complexes of the ( k 1 , 1 ) k , where k = 1 , 2 , . We can study thermoreversible gelation in such extended types of cross-link junctions with complex structure within the theoretical framework presented here. Some of the results on mixed micelles will be reported in our forthcoming paper.

9. Conclusions and Discussion

We have presented a very general theoretical framework for the study of binary thermoreversible gels in common solvents with specified multiple structures of mixed cross-links. All results are presented from a unified point of view in terms of the ternary phase diagrams of a new style by using scaled concentrations. The phase plane is different from that of the conventional triangular diagram; one axis of the plane is the concentration of the functional groups of the primary solute component (mainly polymers), and the other is the concentration of the functional groups of the secondary component (cross-linkers of various types) measured relative to that of the primary component. The method proposed is particularly convenient for the treatment of low-concentration region of the secondary component added to the solution.
The model solutions proposed in this study has obvious advantages in finding the microscopic parameters regarding the network junctions from macroscopic measurements of phase transitions. Hence, it can be regarded as a generalization of the conventional Eldridge–Ferry method [58,59] to binary multiple cross-links. From the theoretical and numerical results, together with qualitative comparison with experiments, the following conclusions can be drawn:
  • There exists an optimal concentration of B (secondary) molecules for gelation if monocomplexes (surface coverage by B molecules) are formed at its high concentration region. The condition for optimization depends on the functionalities of both components, the multiplicity of the cross-links, and the association constants of complexes.
  • In general, there are lower and upper bound in the concentration of B molecules between which the gel phase appears (reentrant sol-gel-sol transition). Both bounds depend on the functionalities of both components, the multiplicity of the cross-links, and the association constants of complexes.
  • Gelation interferes with liquid–liquid phase separation. The relative position of the gel region and the two-phase region systematically shifts with the association constants of the cross-links and Flory’s interaction parameters. In one case, gelation takes place before phase separation (from clear gel to turbid gel), and in another vice versa (only turbid gel) with increase in the concentration of B molecules. The crossing point between the sol-gel transition line and the phase separation line is a higher order critical point.
  • Addition of B molecules to the phase separated solutions of A molecules in a poor solvent may induce homogenization of the solution above a critical concentration and form a clear gel if the solvent is sufficiently good for B molecules. Oppositely, an already formed clear gel may turn into turbid one if the solvent is poor for B molecules.
Our theoretical framework may directly be applicable to some important thermoreversible gels, such as ion-binding PVA, galactomannan, etc. with borate, Ca 2 + , Fe 3 + ions etc. The primary component in these studies are linear chains. Some of the results were already reported in our previous paper [26]. It may also be applicable to the problem of polymer-surfactant interaction. Solutions of associating polymers mixed with low molecular weight surfactants have been a focus of the researchers interest over the past decades. However, relationship between mixing property of the hydrophobes carried by them and their thermoreversible gelation has not been fully explored so far. Our model solution provides a good starting point to study the problem.
We have assumed that the equilibrium constants ( K , λ , μ ) for association are independent of Flory’s interaction parameters χ ’s at our present stage of study. In the experiments they are often observed to depend on each other, because SGT of the same polymer takes place at different polymer concentrations in different solvents. In mixed solvents of different quality, in particular, there is a particular solvent composition at which gelation becomes easiest [60,61]. Study of the solvent effect will therefore require some more precise analysis of the structure of cross-link junctions in contact with solvent molecules.
Further important problem which remains open is the effect of temperature on the phase behavior. The main interest lies in the occurrence of UCST and LCST by association of solvent molecules [45,54]. A particularly important example is the association of water molecules onto polymer chains. In such water-soluble polymer solutions, the primary component is hydrogen-bonding polymer, the secondary component is an active solvent (ethanol, methanol, etc. or organic solvents), and the common solvent is water. Phase separation of such water-soluble polymers in mixed solvents have recently been studied in many researchers, but their gelation property remains as an open question.
Finally, we have to remark on the post-gel region. There are several treatments of the post-gel region in the phase diagrams in the theoretical study of polycondensation by tree statistics; one assumes tree structure for a gel network as for the sol, but the other permits cycle formation within the network. The former was proposed by Stockmayer [41], and the latter by Flory [37,38,39]. Later Ziff and Stell [62] examined another possibility from a kinetic point of view. There appear additional terms in the chemical potentials of each component depending on the treatment of sol-gel interaction in the post-gel region. For thermoreversible gels of one component with pairwise cross-linking in a solvent, these two conventional models were successfully applied without violating the thermodynamic requirements [50]. Therefore, the same treatments have to be also applied to the binary gels with multiple junctions. However, we have used in this study the same forms of the chemical potentials for the post-gel region as well as for the pre-gel region to avoid mathematical complexity. Detailed attempt to improve this point lies, however, outside of the scope of this study.
Generalization of the present theoretical method to multicomponent mixtures of functional molecules is straightforward. We have attempted its application to ternary polymer solutions in which all A,B,C functional groups are allowed to involve in a cross-link junction. It is hoped to pursue these topics in later publications.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guenet, J.M. Thermoreversible Gelation of Polymers and Biopolymers; Academic Press, Harcourt Brace Jovanovich Pub.: London, UK, 1992. [Google Scholar]
  2. Te Nijenhuis, K. Thermoreversible Networks. Adv. Polym. Sci. 1997, 130, 1–261. [Google Scholar]
  3. Zhang, J.; Hu, Y.; Li, Y. Gel Chemistry: Interactions, Structures and Properties; Springer: Heidelberg, Germany, 2018. [Google Scholar]
  4. Thakur, V.K.; Thakur, M.K. Gels Horizons: From Science to Smart Materials. In Polymer Gels: Science and Fundamentals; Springer: Heidelberg, Germany, 2018. [Google Scholar]
  5. Thakur, V.K.; Thakur, M.K. Gels Horizons: From Science to Smart Materials. In Hydrogels: Recent Advances; Springer: Heidelberg, Germany, 2018. [Google Scholar]
  6. Patrickios, C.S. (Ed.) Amphiphilic Polymer Co-Networks: Synthesis, Properties, Modelling and Applications; Royal Society of Chemistry: Cambridge, UK, 2020. [Google Scholar]
  7. Fan, H.; Wang, J.; Jin, Z. Tough, Swelling-Resistant, Self-Healing, and Adhesive Dual-Cross-Linked Hydrogels Based on Polymer-Tannic Acid Multiple Hydrogen Bonds. Macromolecules 2018, 51, 1696–1705. [Google Scholar] [CrossRef]
  8. Zhang, C.; Yang, Z.; Duong, N.T.; Li, X.; Nishiyama, Y.; Wu, Q.; Zhang, R.; Sun, P. Using Dynamic Bonds to Enhance the Mechanical Performance: From Microscopic Molecular Interactions to Macroscopic Properties. Macromolecules 2019, 52, 5014–5025. [Google Scholar] [CrossRef]
  9. Yu, H.C.; Li, C.Y.; Du, M.; Song, Y.; Wu, Z.L.; Zheng, Q. Improved toughness and stability of κ-carrageenan/polyacrylamide double-network hydrogels by dual cross-linking of the first network. Macromolecules 2019, 52, 629–638. [Google Scholar] [CrossRef]
  10. Cristiani, T.R.; Filippidi, E.; Behrens, R.L.; Valentine, M.T.; Eisenbach, C.D. Tailoring the toughness of elastomers by incorporating Ionic cross-linking. Macromolecules 2020, 53, 4099–4109. [Google Scholar] [CrossRef]
  11. Stockmayer, W.H. Molecular distribution in condensation polymers. J. Polym. Sci. 1952, 4, 69–71. [Google Scholar] [CrossRef]
  12. Piepenbrock, M.-O.M.; Lloyd, G.O.; Clarke, N.; Steed, J.W. Metal- and Anion-Binding Supramolecular Gels. Chem. Rev. 2010, 110, 1960–2004. [Google Scholar] [CrossRef] [PubMed]
  13. Christensen, T.; Gooden, D.M.; Kung, J.E.; Toone, E.J. Additivity and the Physical Basis of Multivalency Effects: A Thermodynamic Investigation of the Calcium EDTA Interaction. J. Am. Chem. Soc. 2003, 125, 7357–7366. [Google Scholar] [CrossRef] [PubMed]
  14. Shibayama, M.; Yoshizawa, H.; Kurokawa, H.; Fujiwara, H.; Nomura, S. Sol-gel transition of poly(vinyl alcohol)-borate complex. Polymer 1988, 29, 2066–2071. [Google Scholar] [CrossRef]
  15. Kurokawa, H.; Shibayama, M.; Ishimaru, T.; Nomura, S.; Wu, W.-L. Phase behaviour and sol-gel transition of poly (vinyl alcohol)-borate complex in aqueous solution. Polymer 1992, 33, 2182–2188. [Google Scholar] [CrossRef]
  16. Keita, G.; Ricard, A.; Audebert, R.; Pezron, E.; Leibler, L. The poly(vinyl alcohol)-borate system: Influence of polyelectrolyte effects on phase diagrams. Polymer 1995, 36, 49–54. [Google Scholar] [CrossRef]
  17. Shibayama, M.; Ikkai, F.; Moriwaki, R.; Nomura, S. Complexation of Poly(vinyl alcohol)-Congo Red Aqueous Solutions. 1. Viscosity Behavior and Gelation Mechanism. Macromolecules 1994, 27, 1738–1743. [Google Scholar] [CrossRef]
  18. Shibayama, M.; Ikkai, F.; Nomura, S. Complexation of Poly(vinyl alcohol)-Congo Red Aqueous Solutions. 2. SANS and SAXS Studies on Sol-Gel Transition. Macromolecules 1994, 27, 6383–6388. [Google Scholar] [CrossRef]
  19. Kohri, M.; Yanagimoto, K.; Kohaku, K.; Shiomoto, S.; Kobayashi, M.; Imai, A.; Shiba, F.; Taniguchi, T.; Kishikawa, K. Magnetically Responsive Polymer Network Constructed by Poly(acrylic acid) and Holmium. Macromolecules 2018, 51, 6740–6745. [Google Scholar] [CrossRef]
  20. Pezron, E.; Ricard, A.; Lafuma, F.; Audebert, R. Reversible Gel Formation Induced by Ion Comlexation. 1. Borax-Galactomannan Interactions. Macromolecules 1988, 21, 1121–1125. [Google Scholar] [CrossRef]
  21. Pezron, E.; Leibler, L.; Ricard, A.; Audebert, R. Reversible Gel Formation Induced by Ion Comlexation. 2. Phase Diagram. Macromolecules 1988, 21, 1126–1131. [Google Scholar] [CrossRef]
  22. Pezron, E.; Leibler, L.; Ricard, A.; Lafuma, F.; Audebert, R. Complex Formation in Polymer-Ion Solutions. 1. Polymer Concentration Effects. Macromolecules 1989, 22, 1169–1174. [Google Scholar] [CrossRef]
  23. Pezron, E.; Leibler, L.; Lafuma, F. Complex Formation in Polymer-Ion Solutions. 2. Polyelectrolyte Effects. Macromolecules 1989, 22, 2656–2662. [Google Scholar] [CrossRef]
  24. Lips, A.; Clark, A.H.; Cutler, N.; Durand, D. Measurement of cooperativity of binding of calcium to neutral sodium pectate. Food Hydrocoll. 1991, 5, 87–99. [Google Scholar] [CrossRef]
  25. Zhang, W.; Piculell, L.; Nilsson, S. Effects of Specific Anion Binding on the Helix-Coil Transition of Lower Charged Carrageenans. NMR Data and Conformational Equilibria Analyzed within the Poisson-Boltzman Cell Model. Macromolecules 1992, 25, 6165–6172. [Google Scholar] [CrossRef]
  26. Tanaka, F.; Nakagawa, Y.; Ohta, S.; Ito, T. Thermoreversible gelation with ion-binding cross-links of varible multiplicity. J. Chem. Phys. 2019, 150, 174904. [Google Scholar] [CrossRef]
  27. Norioka, C.; Okita, K.; Mukada, M.; Kawamura, A.; Miyata, T. Biomolecularly stimuli-responsive tetra-poly (ethylene glycol) that undergoes sol-gel transition in response to a target biomolecule. Polym. Chem. 2017, 8, 6378–6385. [Google Scholar] [CrossRef]
  28. Dubin, D. (Ed.) Microdomains in Polymer Solutions; Plenum Press: New York, NY, USA, 1985. [Google Scholar]
  29. Goddard, E.D.; Ananthapadmanabhan, K.P. (Eds.) Interactions of Surfactants with Polymers and Proteins; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar]
  30. Cabane, B. Structure of Some Polymer-Detergent Aggregates in Water. J. Phys. Chem. 1977, 81, 1639–1645. [Google Scholar] [CrossRef]
  31. Brown, G.; Chakrabarti, A. Structure formation in self-associationg polymer and surfactant systems. J. Chem. Phys. 1992, 96, 3251–3254. [Google Scholar] [CrossRef]
  32. Feitosa, E.; Brown, W.; Hansson, P. Interactions between the Non-Ionic Surfactant C12E5 and Poly(ethylene oxide) Studied Using Dynamic Light Scattering and Fluorescsnce Quenching. Macromolecules 1996, 29, 2169–2178. [Google Scholar] [CrossRef]
  33. Feitosa, E.; Brown, W.; Vasilescu, M.; Swason-Vethamuthu, M. Effect of Temperature on the Interaction between the Nonionic Surfactant C12E5 and Poly(ethylene oxide) Investigated by Dynamic Light Scattering and Fluorescence Methods. Macromolecules 1996, 29, 6837–6846. [Google Scholar] [CrossRef]
  34. Annable, T.; Buscall, R.; Ettelaie, R.; Shepherd, P.; Whittlestone, D. Influence of Surfactants on the Rheology of Associating Polymers in Solution. Langmuir 1994, 10, 1060–1070. [Google Scholar] [CrossRef]
  35. Nilsson, S. Interactions between Water-Soluble Cellulose Derivatives and Surfactants. 1. The HPMC/SDS/Water System. Macromolecules 1995, 28, 7837–7844. [Google Scholar] [CrossRef]
  36. Piculell, L.; Thuresson, K.; Ericsson, O. Surfactant Binding and Micellization in Polymer Solutions and Gels: Binding Isotherms and their Consequences. Faraday Disc. 1995, 101, 307–318. [Google Scholar] [CrossRef]
  37. Flory, P.J. Molecular Size Distribution in Three Dimensional Polymers I. Gelation. J. Am. Chem. Soc. 1941, 63, 3083–3090. [Google Scholar] [CrossRef]
  38. Flory, P.J. Molecular Size Distribution in Three Dimensional Polymers II. Trifunctional Branching Units. J. Am. Chem. Soc. 1941, 63, 3091–3096. [Google Scholar] [CrossRef]
  39. Flory, P.J. Molecular Size Distribution in Three Dimensional Polymers III. Tetrafunctional Branching Units. J. Am. Chem. Soc. 1941, 63, 3096–3100. [Google Scholar] [CrossRef]
  40. Flory, P.J. Principles of Polymer Chemistry; Cornell University Press: Ithaca, NY, USA, 1953; Chapter XI. [Google Scholar]
  41. Stockmayer, W.H. Theory of Molecular Size Distribution and Gel Formation in Branched-Chain Polymers. J. Chem. Phys. 1943, 11, 45–55. [Google Scholar] [CrossRef]
  42. Stockmayer, W.H. Theory of Molecular Size Distribution and Gel Formation in Branched Polymers II. General Cross Linking. J. Chem. Phys. 1944, 12, 125–131. [Google Scholar] [CrossRef]
  43. Fukui, K.; Yamabe, T. A General Theory of Gel Formation with Multifunctional Interunit Junctions. Bull. Chem. Soc. Jpn. 1967, 40, 2052–2063. [Google Scholar] [CrossRef] [Green Version]
  44. Tanaka, F.; Stockmayer, W.H. Thermoreversible Gelation with Junctions of Variable Multiplicity. Macromolecules 1994, 27, 3943–3954. [Google Scholar] [CrossRef]
  45. Tanaka, F. Polymer Physics—Applications to Molecular Association and Thermoreversible Gelation; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  46. Macosko, C.W.; Miller, D.R. A New Derivation of Average Molecular Weights of Nonlinear Polymers. Macromolecules 1976, 9, 199–205. [Google Scholar] [CrossRef] [Green Version]
  47. Tanaka, F. Thermodynamic Theory of Network-Forming Polymer Solutions. 1. Macromolecules 1990, 23, 3784–3789. [Google Scholar] [CrossRef]
  48. Tanaka, F. Thermodynamic Theory of Network-Forming Polymer Solutions. 2. Equilibrium Gelation by Conterminous Cross-Linking. Macromolecules 1990, 23, 3790–3795. [Google Scholar] [CrossRef]
  49. Tanaka, F.; Ishida, M. Phase Formation of Two-Component Physical Gels. Phys. A 1994, 204, 660–672. [Google Scholar] [CrossRef]
  50. Tanaka, F.; Ishida, M. Thermoreversible Gelation with Two-Component Networks. Macromolecules 1999, 32, 1271–1283. [Google Scholar] [CrossRef]
  51. Tanaka, F. Gel Formation with Multiple Interunit Junctions. I—Molecules Carrying Different Functional Groups. J. Polym. Sci. Part B Polym. Phys. 2003, 41, 2405–2412. [Google Scholar] [CrossRef]
  52. Tanaka, F. Gel Formation with Multiple Interunit Junctions. II—Mixture of Different Functional Molecules. J. Polym. Sci. Part B Polym. Phys. 2003, 41, 2413–2421. [Google Scholar] [CrossRef]
  53. Gordon, M. Good’s Theory of Cascade Processes applied to the Statistics of Polymer Distribution. Proc. R. Soc. 1962, A268, 240–257. [Google Scholar]
  54. Matsuyama, A.; Tanaka, F. Theory of Solvation-Induced Reentrant Phase Separation in Polymer Solutions. Phys. Rev. Lett. 1990, 65, 341–344. [Google Scholar] [CrossRef] [Green Version]
  55. Ishida, M.; Tanaka, F. Theoretical Study of Hydrogen-Bonded Supramolecular Liquid Crystals. Macromolecules 2002, 35, 7460–7472. [Google Scholar]
  56. Tanaka, F. Polymer-Surfactant Interaction in Thermoreversible Gels. Macromolecules 1998, 31, 384–393. [Google Scholar] [CrossRef]
  57. Renou, F.; Benyahia, L.; Nicoli, T. Structure and Rheology of Mixed Polymeric Micelles Formed by Hydrophobically End-Capped Poly(ethylene oxide). Macromolecules 2008, 41, 6523–6530. [Google Scholar] [CrossRef]
  58. Eldridge, J.E.; Ferry, J.D. Studies on the Cross-Linking Process in Gelatin Gels. III. Dependence of Melting Point on Concentration and Molecular Weight. J. Phys. Chem. 1954, 58, 992–995. [Google Scholar] [CrossRef]
  59. Tanaka, F.; Nishinari, K. Junction Multiplicity in Thermoreversible Gelation. Macromolecules 1996, 29, 3625–3628. [Google Scholar] [CrossRef]
  60. Ohkura, M.; Kanaya, T.; Kaji, K. Gels of poly(vinyl alcohol) from dimethyl sulphoxide/water solutions. Polymer 1992, 33, 3686–3690. [Google Scholar] [CrossRef]
  61. Ohkura, M.; Kanaya, T.; Kaji, K. Gelation rates of poly(vinyl alcohol) solution. Polymer 1992, 33, 5044–5048. [Google Scholar] [CrossRef]
  62. Ziff, R.M.; Stell, G. Kinetics of Polymer Gelation. J. Chem. Phys. 1980, 73, 3492–3499. [Google Scholar] [CrossRef]
Figure 1. Connected tree-like clusters consisting of tetrafunctional molecules R{A 4 } and trifunctional molecules R{B 3 } in a solvent. Reversible chemical reaction in forming cross-link junctions contains arbitrary number k 1 of functional groups A and k 2 of groups B.
Figure 1. Connected tree-like clusters consisting of tetrafunctional molecules R{A 4 } and trifunctional molecules R{B 3 } in a solvent. Reversible chemical reaction in forming cross-link junctions contains arbitrary number k 1 of functional groups A and k 2 of groups B.
Gels 07 00089 g001
Figure 2. Reentrant sol region at high concentration of B molecules. Most of the clusters consist of small numbers of A molecules whose surfaces are covered by B molecules.
Figure 2. Reentrant sol region at high concentration of B molecules. Most of the clusters consist of small numbers of A molecules whose surfaces are covered by B molecules.
Gels 07 00089 g002
Figure 3. Sol-gel-sol transition lines (red solid lines) and spinodal lines (blue broken lines) for pairwise reaction k 1 = k 2 = 1 drawn on the ternary phase plane. The horizontal axis is the scaled concentration ξ of the secondary functional groups B, and the vertical axis is the volume fraction ϕ A of the primary functional groups A. The upper bound of the graph (thin solid line) shows the limiting solution without solvent. The functionalities are fixed at f = g = 3 . The number of repeat units are fixed at n A = n B = 5 . (a) The interaction parameter χ A , B between A and B molecules is changed from 2.0 to 5.0 while the association constant is fixed at K = 30 . The two phase region is indicated by 2 Φ . (b) The same as (a) but the association constant is changed from 30 to 100 while the interaction parameter is fixed at χ A , B = 5.0 .
Figure 3. Sol-gel-sol transition lines (red solid lines) and spinodal lines (blue broken lines) for pairwise reaction k 1 = k 2 = 1 drawn on the ternary phase plane. The horizontal axis is the scaled concentration ξ of the secondary functional groups B, and the vertical axis is the volume fraction ϕ A of the primary functional groups A. The upper bound of the graph (thin solid line) shows the limiting solution without solvent. The functionalities are fixed at f = g = 3 . The number of repeat units are fixed at n A = n B = 5 . (a) The interaction parameter χ A , B between A and B molecules is changed from 2.0 to 5.0 while the association constant is fixed at K = 30 . The two phase region is indicated by 2 Φ . (b) The same as (a) but the association constant is changed from 30 to 100 while the interaction parameter is fixed at χ A , B = 5.0 .
Gels 07 00089 g003
Figure 4. The same as Figure 3a but the interaction parameters are changed: (blue broken line) χ A , 0 = 2.0 , χ A , B = χ B , 0 = 0 , (green broken line) χ B , 0 = 2.0 , χ A , B = χ A , 0 = 0 . In the former case, the solution starts from a two-phase state, whose region shrinks with increase in the concentration of B molecules (good solvent for B), and eventually merges into one phase. Turbid gel turns into clear one with increase in B concentration. The latter case shows an opposite tendency; from clear gel to turbid gel.
Figure 4. The same as Figure 3a but the interaction parameters are changed: (blue broken line) χ A , 0 = 2.0 , χ A , B = χ B , 0 = 0 , (green broken line) χ B , 0 = 2.0 , χ A , B = χ A , 0 = 0 . In the former case, the solution starts from a two-phase state, whose region shrinks with increase in the concentration of B molecules (good solvent for B), and eventually merges into one phase. Turbid gel turns into clear one with increase in B concentration. The latter case shows an opposite tendency; from clear gel to turbid gel.
Gels 07 00089 g004
Figure 5. (a) Phase diagram of a mononuclear A cross-link junction; f = g = 3 , ( k 1 , k 2 ) = ( 1 , 2 ) . The gel region is bound by ξ min = 0 from below and ξ max = g ( f k 2 + k 2 ) from above. SGST line (red line), spinodals for χ A , B = 2.0 (red broken line), for χ A , 0 = 2.0 (blue broken line), and for χ B , 0 = 2.0 (green broken line) are shown. For the repulsive interaction χ A , B = 2.0 in a good solvent χ A , 0 = χ B , 0 = 0 , phase separation takes place before the solution gels when the concentration of added B molecules is increased. (b) Reactivity p (solid lines) and q (broken lines) shown as functions of the scaled B concentration for three fixed A concentrations. They all start from 0. The reactivity of B groups shows a maximum near the stoichiometric concentration.
Figure 5. (a) Phase diagram of a mononuclear A cross-link junction; f = g = 3 , ( k 1 , k 2 ) = ( 1 , 2 ) . The gel region is bound by ξ min = 0 from below and ξ max = g ( f k 2 + k 2 ) from above. SGST line (red line), spinodals for χ A , B = 2.0 (red broken line), for χ A , 0 = 2.0 (blue broken line), and for χ B , 0 = 2.0 (green broken line) are shown. For the repulsive interaction χ A , B = 2.0 in a good solvent χ A , 0 = χ B , 0 = 0 , phase separation takes place before the solution gels when the concentration of added B molecules is increased. (b) Reactivity p (solid lines) and q (broken lines) shown as functions of the scaled B concentration for three fixed A concentrations. They all start from 0. The reactivity of B groups shows a maximum near the stoichiometric concentration.
Gels 07 00089 g005
Figure 6. (a) Phase diagram for the mononuclear B cross-link junction. f = g = 3 , ( k 1 , k 2 ) = ( 2 , 1 ) (sandwich cross-link). The gel region is bound at ξ min = 1 / f ( g k 1 + k 1 ) from below. There is no optimal concentration. SGT line (red line), spinodals for χ A , B = 2.0 (red broken line), for χ A , 0 = 2.0 (blue broken line), and for χ B , 0 = 2.0 (green broken line) are shown. Adding trifunctional B molecules promotes gelation monotonically with B concentration. For the repulsive interaction χ A , B = 2.0 in a good solvent χ A , 0 = χ B , 0 = 0 , gelation takes place before the solution phase separates when the concentration of added B molecules is increased. Therefore, there is a uniform clear gel region at low concentration of B molecules. (b) Reactivity p (solid lines) and q (broken lines) shown as functions of the scaled B concentration for three fixed A concentrations. They are monotonic functions of B concentration, but different from those of ( 1 , 2 ) type cross-links in that q starts from finite values.
Figure 6. (a) Phase diagram for the mononuclear B cross-link junction. f = g = 3 , ( k 1 , k 2 ) = ( 2 , 1 ) (sandwich cross-link). The gel region is bound at ξ min = 1 / f ( g k 1 + k 1 ) from below. There is no optimal concentration. SGT line (red line), spinodals for χ A , B = 2.0 (red broken line), for χ A , 0 = 2.0 (blue broken line), and for χ B , 0 = 2.0 (green broken line) are shown. Adding trifunctional B molecules promotes gelation monotonically with B concentration. For the repulsive interaction χ A , B = 2.0 in a good solvent χ A , 0 = χ B , 0 = 0 , gelation takes place before the solution phase separates when the concentration of added B molecules is increased. Therefore, there is a uniform clear gel region at low concentration of B molecules. (b) Reactivity p (solid lines) and q (broken lines) shown as functions of the scaled B concentration for three fixed A concentrations. They are monotonic functions of B concentration, but different from those of ( 1 , 2 ) type cross-links in that q starts from finite values.
Gels 07 00089 g006
Figure 7. Phase diagram of monofunctional low molecular weight B molecules ( n B = 1 , g = 1 ) with χ A , B varied (blue broken line) in a good solvent for f = n A = 10 , ( k 1 , k 2 ) = ( 2 , 1 ) (sandwich cross-link of longer molecules). The lower bound of the gel region is ξ min = f k 1 . Many A functional groups remain active on the surface of clusters under excess of B molecules. Clear gel region shrinks with increase in χ A , B , and eventually disappears.
Figure 7. Phase diagram of monofunctional low molecular weight B molecules ( n B = 1 , g = 1 ) with χ A , B varied (blue broken line) in a good solvent for f = n A = 10 , ( k 1 , k 2 ) = ( 2 , 1 ) (sandwich cross-link of longer molecules). The lower bound of the gel region is ξ min = f k 1 . Many A functional groups remain active on the surface of clusters under excess of B molecules. Clear gel region shrinks with increase in χ A , B , and eventually disappears.
Gels 07 00089 g007
Figure 8. A typical example of phase diagrams for a symmetric IPN shown on the triangular volume fraction plane at a given temperature with fixed multiplicities k 1 = 3 , k 2 = 3 for trifunctional A- and B molecules f = g = 3 with λ = μ = 3.0 and γ ( A ) = γ ( B ) = 1 . SOL region, A-GEL region, B-GEL region and IPN region (A/B-GEL) are indicated. Spinodal lines for AB interaction (blue broken line), for poor solvent of A (black broken lines) and for poor solvent of B (green broken lines) are also shown.
Figure 8. A typical example of phase diagrams for a symmetric IPN shown on the triangular volume fraction plane at a given temperature with fixed multiplicities k 1 = 3 , k 2 = 3 for trifunctional A- and B molecules f = g = 3 with λ = μ = 3.0 and γ ( A ) = γ ( B ) = 1 . SOL region, A-GEL region, B-GEL region and IPN region (A/B-GEL) are indicated. Spinodal lines for AB interaction (blue broken line), for poor solvent of A (black broken lines) and for poor solvent of B (green broken lines) are also shown.
Gels 07 00089 g008
Figure 9. Phase diagram of a mixture of trifunctional A molecules and bifunctional B molecules in a solvent, whose functional groups are completely miscible in a range of multiplicity between 3 and 5. Solution of A molecules has SGT without B molecules, whose gel concentration monotonically decreases by mixing B molecules. Spinodal line for an interaction parameter χ A , B = 5.0 between A- and B molecules in a good solvent (green broken line) and in a solvent poor to A molecules χ A , 0 = 0.9 (blue broken line) are also drawn. For the latter, there is a critical B concentration above which the solution becomes uniform.
Figure 9. Phase diagram of a mixture of trifunctional A molecules and bifunctional B molecules in a solvent, whose functional groups are completely miscible in a range of multiplicity between 3 and 5. Solution of A molecules has SGT without B molecules, whose gel concentration monotonically decreases by mixing B molecules. Spinodal line for an interaction parameter χ A , B = 5.0 between A- and B molecules in a good solvent (green broken line) and in a solvent poor to A molecules χ A , 0 = 0.9 (blue broken line) are also drawn. For the latter, there is a critical B concentration above which the solution becomes uniform.
Gels 07 00089 g009
Figure 10. (a) Phase diagram of telechelic associating polymers ( f = 2 , n A = 100 ) mixed with surfactant molecules ( g = 1 , n B = 1 ). Hydrophobes carried by both molecules are assumed to be completely miscible within a formed micelles covering an allowed multiplicity range from k min = 5 to k max = 20 . The polymer-solvent interaction parameter χ A , 0 is changed for the calculation of spinodals (blue broken lines). The sol-gel transition (red line) has a minimum at a certain concentration of the surfactant near the stoichiometric one, so that it is an example of the reentrant sol-gel-sol transition. (b) The minimum multiplicity is allowed to change from 2 to 8 while the maximum multiplicity is fixed at 20. If k min = 2 is allowed, the minimum concentration does not appear.
Figure 10. (a) Phase diagram of telechelic associating polymers ( f = 2 , n A = 100 ) mixed with surfactant molecules ( g = 1 , n B = 1 ). Hydrophobes carried by both molecules are assumed to be completely miscible within a formed micelles covering an allowed multiplicity range from k min = 5 to k max = 20 . The polymer-solvent interaction parameter χ A , 0 is changed for the calculation of spinodals (blue broken lines). The sol-gel transition (red line) has a minimum at a certain concentration of the surfactant near the stoichiometric one, so that it is an example of the reentrant sol-gel-sol transition. (b) The minimum multiplicity is allowed to change from 2 to 8 while the maximum multiplicity is fixed at 20. If k min = 2 is allowed, the minimum concentration does not appear.
Gels 07 00089 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tanaka, F. Thermoreversible Gelation with Two-Component Mixed Cross-Link Junctions of Variable Multiplicity in Ternary Polymer Solutions. Gels 2021, 7, 89. https://doi.org/10.3390/gels7030089

AMA Style

Tanaka F. Thermoreversible Gelation with Two-Component Mixed Cross-Link Junctions of Variable Multiplicity in Ternary Polymer Solutions. Gels. 2021; 7(3):89. https://doi.org/10.3390/gels7030089

Chicago/Turabian Style

Tanaka, Fumihiko. 2021. "Thermoreversible Gelation with Two-Component Mixed Cross-Link Junctions of Variable Multiplicity in Ternary Polymer Solutions" Gels 7, no. 3: 89. https://doi.org/10.3390/gels7030089

APA Style

Tanaka, F. (2021). Thermoreversible Gelation with Two-Component Mixed Cross-Link Junctions of Variable Multiplicity in Ternary Polymer Solutions. Gels, 7(3), 89. https://doi.org/10.3390/gels7030089

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop