Prerpints.org logo
Preprint
Article

Gelation Time of Network-Forming Polymer Solutions with Reversible Cross-Link Junctions of Variable Multiplicity

Altmetrics

Downloads

121

Views

25

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

11 April 2023

Posted:

12 April 2023

You are already at the latest version

Alerts
Abstract
The gelation time tg necessary for a solution of functional (associating) molecules to reach its gel point after a temperature jump, or a sudden concentration change, is theoretically calculated as a function of the concentration and temperature for the molecules carrying the number f of functional groups. Quite generally, tg is given by the product of the relaxation time tR and a thermodynamic factor Q. For pairwise cross-linking, and cross-linking with a fixed multiplicity k, precise values of tg are derived. For cross-linking with variable multiplicity, kinetic equation of stepwise association is derived for the reactivity pk of the functional groups, on the basis of which tg and tR are estimated. The thermodynamic factor Q is shown to generate a singularity of logarithmic divergence as the temperature (concentration) approaches the equilibrium gel point, while the relaxation time tR changes continuously across it. The retardation effect on the gelation time due to the reversibility of the cross-linking is explicitly calculated for some specific models of cross-linking. For a micellar cross-linking covering a wide range of the multiplicity, as seen in hydrophobically-modified water-soluble polymers, tR is shown to obey a formula similar to the Aniansson-Wall law.
Keywords: 
Subject: Chemistry and Materials Science  -   Polymers and Plastics

1. Introduction

Thermoreversible gelation in solutions of polymers, as well as of low molecular weight molecules, has been attracting researcher’s interest[1,2,3,4,5,6,7]. Many examples of the phase diagrams with sol–gel transition lines have been reported in the literature, early ones of which are atactic polystyrene in various organic solvents[8,9], poly(n-butylmethacrylate) in 2-propanol[10,11], poly(vinyl chloride) in dimethyl malonate[12] etc. Quite recently, phase diagrams of heat-setting sol–gel transition interfering with LCST phase separation in aqueous solutions of methylcellulose[13], poly(N-isopropylacrylamide)[14,15], and poly(N-isopropylacrylamide-co-n-butylmethacrylate), known as mebiol[16], have been reported. Some reviews and conceptual works have appeared with relation to responsive gels[17,18], hydrogels for biomedical applications[6,7,19], and hydrogen-bonding supramolecular gelators[20,21,22]. However, the kinetic process of cross-linking reaction to reach the gel point, in particular the gelation time, has not been clarified yet.
Here, the gelation time t g is defined by the time necessary for the network-forming solutions of functional molecules to reach their gel points after cross-linking reaction is started. It is related to the relaxation time t R for the solution to go back to their thermal equilibrium state, but must be strictly distinguished from each other because the former goes to infinity on the equilibrium sol–gel transition line while the latter remains continuous across it.
Fundamental investigation of polymer cross-linking reactions as a function of time is very important to obtain better understanding of the gelation mechanism for controlling the production process and final performance of gels. In particular, acceleration of the gelation time in the most efficient way under restricted conditions is critical for the application of gels to biomedical technology, food industry, adhesion industry, etc.
Gelation time is very important also from scientific viewpoint. It includes microscopic information on the kinetics of cross-linking reaction. By the macroscopic measurements of the gel point by rheology, we can speculate microscopic structure and formation process of the network junctions, in particular the rate-controlling steps.
The gelation time critically depends on the thermodynamic conditions under which cross-linking reaction proceeds, such as temperature, concentration of the functional molecules, their mixing ratio if there are more than one component, amount of the cross-linkers if any, pH, etc. The purpose of the present study is to clarify the difference between t g and t R , and find them as functions of the system parameters such as concentration and temperature by focusing on their relations to the rate constants of reaction in the kinetic equation.
Experimental data of the gelation time t g in the literature obtained by the rheological measurements, scattering, calorimetry etc are usually plotted against the temperature and/or polymer concentration to find the activation energy in the rate constant and the power index of the concentration. For instance, Ohkura et al[23,24] reported the data on the gelation time of poly(vynil alcohol) (PVA) in a mixed solvent of water and dimethylsulphoxide. They measured the time by ball-drop method in the temperature quenching experiments for various polymer concentrations. Results were summarized in the form
t g − 1 = A ( T ) ( Δ ϕ ) n
where Δ ϕ is the concentration of PVA measured from the overlap concentration. The critical concentration at which t g becomes infinite was found to be independent of temperature. The exponent n was found to be n = 2 under the assumption that the gel point is given by the overlap concentration. Because the cross-link junctions are formed by micro-crystallites, the temperature factor was found to agree with the nucleation rate
A ( T ) ≃ e − B / T ( Δ T ) 2
where B is a temperature-independent constant.
Investigation was done by Mal et al[25] on the solutions of crystalline polymer poly(viylidene fluoride) (PVF 2 ) in two different solvents. Results were plotted in the form of (1). The critical concentration depends on the temperature, from which they constructed the sol–gel transition line for three samples. The gel point depends on the nature of the solvent; gelation is easier in a poor solvent. The exponent of concentration was found to be n ≃ 0.45 − 0.60 , much smaller than the value for PVA.
Hong et al[26] studied the same polymer PVF 2 in tetra(ethylene glycol) dimethyl ether. They confirmed that the temperature factor takes a crystalline nucleation form, but found a crossover in the power law from n = 2.17 to n = 1.18 as the critical gel point is approached. Based on these data, a new model of gel formation in crystalline polymers was proposed.
Similar study was done by Tobitani et al[27,28] for heat-setting gels of protein, bovine serum albumin (BSA) and beta-lactoglobulin ( β -Lg), in order to construct the equilibrium sol–gel transition curves by the measurement of dynamic-mechanical moduli. Data were again summarized in the form of (1). Because the gelation is accompanied by protein denaturation, the temperature factor A ( T ) took an activation type. They found also the exponent of the concentration much higher than two. To incorporate these results, they proposed a new model of the gelation time by considering multiple reaction of functional groups in the cross-link junctions[29,30,31]. We will refer to their idea again in the section treating multiple cross-links. The activation type temperature dependence of the gelation time has been reported for some other gel-forming polymers, such as polyurethane dispersion[32], silica alkoxides[33,34] etc.
Investigation on the gelation time was reported in the literature for the binary cross-linking in the mixture of functional molecules[35,36,37,38,39]. The main problem of binary gels is to find the optimal mixing ratio for minimizing the gelation time, so that plot of t g against the mixing ratio for different quenching temperatures is fundamental.
From theoretical side, there have been accumulating studies based on the kinetic equations for the cross-linking reaction of functional groups in the classical Flory–Stockmayer pictures of gelation. For instance, Stockmayer[40,41] derived the rate equation of irreversible gelation for the reactivity (conversion) α as a function of time, and found its value at the gel point.
Later, Ziff et al[42,43,44,45] extended Stockmayer’s model on the basis of more general kinetic equation of cluster formation to include the time development of reactivity in the post-gel region. They treated not only unary cross-linking (self-reaction) but also binary cross-linking (copolymerization) in their series of papers[43]. Scaling and universality of gel formation was also discussed[45]. Dongen et al[46] considered reversibility of the cross-linking reactions, and studied retardation effect due to the existence of reverse reaction. However, unfortunately, the effect of temperature and concentration was not explicitly considered in these studies.
In this paper, we focus on the reversibility and variable cross-link multiplicity of stepwise rate equation for the cross-linking reaction. Solving the rate equation for the probability p k ( t ) of cross-linking with multiplicity k, we derive the gelation time t g as a function of the polymer concentration and temperature. In contrast to (1), it includes a temperature–concentration cross-term. Results are expressed in the common form t g = t R Q in terms of the relaxation time t R and a thermodynamic factor Q. The latter shows a singularity near the equilibrium sol–gel transition line. We investigate the retardation effect due to the reverse reaction (dissociation) and the effect of variable multiplicity of the junction. Gelation time for the reversible binary cross-linking will be reported in our forthcoming paper.

2. Theoretical Method

The model solution we consider is a polymer solution in which the number N of reactive (associative) molecules (denoted as R{A f }) with degree of polymerization n are dissolved in the number N 0 of solvent molecules (S). Molecules can be any type, such as high molecular weight linear polymers, star polymers, or low molecular weight polyfunctional molecules, etc. Each molecule carries the number f of functional groups A (functionalityf) which can form interchain cross-links made up of variable number k of A-groups (multiplicityk)[4,29,30].
We are based on the lattice-theoretical picture of polymer solutions[47], and divide the system volume V into cells of size a of the solvent molecule, each of which is assumed to accommodate a statistical repeat unit of the reactive molecules. The volume of a reactive molecule is then given by n, and that of a solvent molecule is n 0 ≡ 1 in the unit of the cell volume. We assume incompressibility of the solution, so that we have Ω = n N + N 0 for the total volume. The volume fraction of each component is then given by ϕ = n N / Ω for the reactive molecule, ϕ 0 = N 0 / Ω for the solvent. In terms of the functional groups, the number concentration of A-groups on the reactive molecules is ψ = f ϕ / n .
We study the time development of the solution after its temperature is changed suddenly from the initial one T i to the final one T ≡ T f (temperature jump), or after the molecules are quickly dissolved to the solvent in the preparation of the solution with concentration ϕ at a constant temperature T (concentration jump) (see Figure 1 (a)). The temperature quench depth is defined by Δ T ≡ T gel − T , while the concentration quench depth is Δ ϕ ≡ ϕ − ϕ g , where T gel and ϕ g are their values on the sol–gel transition line. For simplicity, we consider in this paper only cold-setting gels. For heat-setting gels, other factors, such as dehydration, polymer conformation change, temperature activation of the functional groups, etc must be considered in addition to simple cross-linking. Hence they lie beyond the scope of the present paper.
After a temperature (or concentration) jump, cross-link reactions leading to network formation proceed. Let n k ( t ) be the number of the cross-link junctions of multiplicity k at time t. Because each junction includes the number k of the functional groups A, the probability p k ( t ) for an arbitrarily chosen functional group to belong to the junction of multiplicity k is related to the number of junctions by the equation
ψ p k ( t ) / k = n k ( t )
After a long time, the solution reaches its equilibrium state with equilibrium reactivity p ¯ k ≡ p k ( t → ∞ ) . Figure 1 (b) schematically shows the time development of p ( t ) ≡ p 2 ( t ) , its linearized approximation p ˜ ( t ) , irreversible counterpart p ( IR ) ( t ) , reactivity of the sol part p s ( t ) in the post-gel region. Because of the reverse reaction (dissociation of cross-links), there is a retardation in the gelation time from the irreversible limit.
In our previous work[31], we studied in detail thermoreversible gelation and phase separation in solutions capable of unary (self) cross-linking. We started from the equilibrium condition
ψ p ¯ k / k ( ψ p ¯ 1 ) k = K k ( T )
where K k ( T ) is the equilibrium constant, p 1 ( t ) is the probability for a functional group to remain free from reaction, and p ¯ 1 is its equilibrium value. This equilibrium condition leads to
ψ p ¯ k = k K k z ¯ k
for the reactivity in terms of the number of free groups z ¯ ≡ ψ p ¯ 1 . From the normalization condition of p ¯ k , we find the conservation law
ψ = z ¯ u ( z ¯ )
where
u ( z ¯ ) = ∑ k ≥ 1 k K k z ¯ k − 1

2.1. Pairwise Cross-Linking

Let us first consider the simplest case of pairwise cross-linking reaction
2 A 1 ⇆ A 2
for which k takes only two values: k = 1 (free), and k = 2 (bound). Cross-linking by covalent bond is included as the irreversible limit of such reaction. For such a simple second order reaction, we can find the exact solution of the rate equation. Therefore, the model provides a good starting point for the study of more complex cross-linking.
Let us write as p 2 ≡ p , p 1 = 1 − p . Probability p is the conventional reactivity of the cross-linking reaction. Time development kinetic equation for the number n k of cross-link junctions can be written as
d n 2 ( t ) d t = α n 1 ( t ) 2 − β n 2 ( t )
and hence we have
d p ( t ) d t = 2 ψ α ( 1 − p ) 2 − β p ( t )
where α and β are the rate constant of the forward reaction and backward reaction. (Throughout this paper, we avoid conventional symbols k f and k b etc to prevent confusion with multiplicity k.) By using the scaled time τ ≡ 2 ψ α t , we have a simple second order equation
d p d τ = g ( p )
where
g ( p ) ≡ ( 1 − p ) 2 − p 2 K 2 ψ
Here, a new constant
K 2 = α / β
is the equilibrium constant written in terms of the rate constants.
The solution of this equation with the initial condition p ( 0 ) = 0 is given by
p ( τ ) = η ( + ) η ( − ) ( 1 − e − γ τ ) η ( + ) − η ( − ) e − γ τ
where
η ( ± ) = 1 4 K 2 ψ 1 + 4 K 2 ψ ± 1 + 8 K 2 ψ
and
γ ≡ η ( + ) − η ( − ) = 1 2 K 2 ψ 1 + 8 K 2 ψ
The equilibrium reactivity is given by p ¯ = η ( − ) . Therefore γ has the meaning of relaxation time of the reversible reaction, i.e. the time necessary for the system to reach its chemical equilibrium. Detailed derivation of (14) is given in Appendix A.
Result is drawn in Figure 2 (a) for the functionality f = 3 , 4 , 5 . In Figure 2 (b), the weight- and number-average molecular weight of the three-dimensional cross-linked polymers ( M ¯ w and M ¯ n in the unit of the molecular weight of the primary molecules) are plotted against time. They are explicitly given by[40]
M ¯ w ( τ ) = 1 + p ( τ ) 1 − f ′ p ( τ ) , M ¯ n ( τ ) = 1 1 − f p ( τ ) / 2
In the post-gel region where the gel point is passed, the reactivity of the sol part and that of the gel part become different. We have calculated the former on the basis of Flory’s picture. But, there are other possibilities[42,43,44,45]. In this paper we focus on the process of approaching the gel point, and prevent discussion on the post-gel regime.
In Figure 2 (b), the weight-average molecular weight (broken lines) and the number-average molecular weight (solid lines) are plotted against the scaled time τ for the combined variable (scaled concentration) x ≡ K 2 ψ = 2.0 . The number-average remains finite at the gel point.
Because the relaxation time t R is defined by the speed to reach the equilibrium state, it is given by 1 / γ . In terms of the bare time, we have
1 t R = 2 ψ α γ = α K 2 1 + 8 K 2 ψ = β 1 + 8 K 2 ψ
In the limit of irreversible reaction where β → 0 with finite α , γ τ is small, and η ( ± ) → 1 , so that we go back to the Stockmayer’s result[40,41](see p ( IR ) ( t ) in Figure 1 (b))
p ( t ) = p ( IR ) = 2 ψ α t 1 + 2 ψ α t
So far, the reaction has been assumed to occur independently with equal probability for any functional group (assumption of equal reactivity). Now, we employ an additional assumption such that all cross-linked three-dimensional molecules take tree forms. Cycle formation is not considered. Then, the conventional tree statistics[40,41,47,48,49,50] gives the gel point
p g = 1 / ( f − 1 )
for the reactivity. Substituting this value for p ( τ ) in (14), we find for the gelation time
t g ( x ) = t R ( x ) Q ( x )
where
Q ( x ) ≡ ln η ( − ) ( x ) [ f ′ η ( + ) ( x ) − 1 ] η ( + ) ( x ) [ f ′ η ( − ) ( x ) − 1 ]
is a function of a combined variable x ≡ K 2 ( T ) ψ of the temperature and concentration. Since this factor Q ( x ) drastically depends upon the quench depth Δ T , or Δ ϕ , in the experiments near the sol–gel transition point, we refer to it as thermodynamic factor. In particular, the gel-point condition (20) is equivalent to f ′ η ( − ) ( x ) − 1 = 0 , so that Q ( x ) goes logarithmically to infinity as x approaches the critical value x g = f ′ / 2 f ′ ′ 2 at the gel point. (In what follows we use abbreviated notations f ′ ≡ f − 1 , f ′ ′ ≡ f − 2 , ⋯ for simplicity.) Near the gel point, we have x ≡ x g ( 1 + ϵ ) ( ϵ < < 1 ) . A lengthy calculation leads to the form
Q ( x ) ≃ ln f f ′ ′ f ′ 2 2 1 ϵ ≃ − ln ϵ + finite const .
Therefore, in the limit of shallow quenching, we have to wait for a logarithmically long time for the solution to turn into a gel.
As for the relaxation time, we have to notice that t R can be expressed in terms of the equilibrium reactivity. The solution of the conservation law (6) for a pairwise cross-linking is given by
z ¯ = 1 4 K 2 − 1 + 1 + 8 K 2 ψ
We therefore have p ¯ 2 / p ¯ 1 = 2 K 2 z ¯ , and hence
1 t R = β 1 + 2 p ¯ 2 p ¯ 1
We will have some extensions of this formula in the following sections for the study of multiple cross-links.
We next consider retardation of the gelation time due to the reversible reaction. We start from the gelation time (21) in the form
ψ α t g = x 1 + 8 x Q ( x )
and take β → 0 limit while α is kept finite. Because x takes a large value, we expand the right hand side of this equation in powers of 1 / x . After a quite a lengthy calculation, we find
t g = 1 2 α f ′ ′ ψ 1 + R f , 2 x + O ( 1 x 2 )
where the front factor
t g ( IR ) ≡ 1 2 α f ′ ′ ψ
is the gelation time of irreversible reaction (covalent bonding). It can directly be obtained from Stockmayer’s form (19) by fixing p ( t ) = 1 / f ′ . The first correction due to reversible reaction is O ( 1 / x ) . Its coefficient is found to be
R f , 2 = 3 f − 4 6 f ′ ′ 2
Let us refer to it as retardation coefficient. In the next section, we shall derive the retardation coefficient R f , k for the cross-linking with arbitrary multiplicity k.
To compare with experimental data, we plot in Figure 3 the reciprocal of the gelation time
1 β t g = 1 + 8 x Q ( x ) ≡ Γ ( x )
as a function of the combined variable x for varied functionality f.
Near the sol–gel transition point x = x g , it goes logarithmically to zero as
Γ ( x ) ≃ 1 + 8 x g ln ( 1 / Δ x )
For larger values of the functionality f, all plots remain qualitatively the same although quantitatively very different. They approach much closer to the vertical axis.
Experimental data[23,24] on PVA in a mixed solvent of water and DMSO suggest that the gelation time near the transition point obeys a power low of concentration with the power index 2. The paper attributed this observation to the binary collision of the polymers at the overlap concentration, which is necessary for the formation of microcrystalline junctions for PVA solutions. In the present models of cross-linking reaction among functional groups, the gel point is located at much lower concentration than the overlap concentration.
For sufficiently high x, it approaches the irreversible limit
Γ ( x ) ≃ 2 f ′ ′ x
In this limit, the gelation time is separable into a product of the temperature factor λ ( T ) and the concentration factor ψ = f ϕ / n . The linear behavior in the high concentration region has not been experimentally reported.
So far, the gelation time has been derived from the exact solution of the rate equation. Only assumption is that the gel point is assumed to be given by the conventional tree approximation. We now examine the validity of linearized equation of (11) before moving to more complex cross-links with variable multiplicity. As time goes on, the polymer solution approaches its equilibrium state. Assuming that it is sufficiently close to equilibrium, let the reactivity be
p ( t ) = p ¯ ( 1 + ξ ( t ) )
and consider only linear terms of the deviation ξ ( t ) in (11). We have
d ξ d τ = g ′ ( p ¯ ) ξ
where
g ′ ( p ¯ ) = − 2 ( 1 − p ¯ ) + 1 2 K 2 ψ
is the derivative of g ( p ¯ ) . Substituting the equilibrium value of p leads to
g ′ ( p ¯ ) = − 1 2 K 2 ψ 1 + 2 p ¯ 2 p ¯ 1
Hence the relaxation time in the linear approximation agrees with the exact one. If we assume p ( 0 ) = 0 , the reactivity in the linear approximation takes the form (see Figure 1 (b))
p ˜ ( t ) = p ¯ ( 1 − e − t / t R )
The gel-point condition p ( t g ) = 1 / f ′ then leads to (21), but with Q ( x ) of slightly different form
Q ( x ) = ln f ′ η ( + ) ( x ) f ′ η ( − ) ( x ) − 1 ≃ − ln ϵ + finite const .
for the thermodynamic factor leading to the same logarithmic divergence. Unfortunately, the linear approximation does not give a correct irreversible limit, so that retardation effect is impossible to study.

2.2. Fixed-Multiplicity Model

Because most of physical gels have multiple cross-link junctions, let us next consider the effect of cross-link multiplicity. We first study an extreme case where simultaneous formation of k junctions takes place
k A 1 ⇄ A k
before we move onto more complex case of stepwise association. A functional group is either free ( k = 1 ) or reacted ( k ≥ 2 ). Let us refer to it as fixed-multiplicity model (see Figure 4).
We have a rate equation for the fixed-multiplicity reversible reaction
d n k ( t ) d t = α n 1 ( t ) k − β n k ( t )
Or, equivalently
d p k ( t ) d t = α k ψ k ′ p 1 ( t ) k − β p k ( t )
( k ′ ≡ k − 1 ). Hence we have
d p d τ = g ( p )
where p ≡ p k , p 1 = 1 − p . Time is scaled as
τ ≡ α k ψ k ′ t
and
g ( p ) ≡ ( 1 − p ) k − β ′ p
The coefficient β of the reverse reaction is changed to
β ′ ≡ β α k ψ k ′ = 1 k K k ψ k ′
where
K k = α β
is the equilibrium constant.
In the study of heat-setting β -Lactoglobulin, Tobitani and Ross-Murphy[27,28] solved this equation in the irreversible limit of β → 0 . The rate equation can easily be integrated in this limit. We find
τ = 1 k ′ 1 ( 1 − p ) k ′ − 1 ( 1 − p ( 0 ) ) k ′
The relaxation is not exponential, but obeys a power law. Let us assume that all functional groups are free in the initial state ( p ( 0 ) = 0 ).
In our previous study of gelation with cross-links of variable multipliicity[30,31], we showed that the gel point is in general given by the condition
1 − f ′ ( 〈 k 〉 w − 1 ) = 0
where
〈 k 〉 w ≡ ∑ k ≥ 1 k p k
is the average multiplicity of the junctions. The average ∑ k ( k − 1 ) p k = 〈 k 〉 w − 1 was referred to as branching number[30,31]. In the case of fixed-multiplicity model, the gel point is found by the condition
p k ( t g ) = 1 f ′ k ′ ≡ p g
Substituting this value into (47), we find
1 t g ( IR ) = α k ψ k ′ k ′ ( 1 − p g ) k ′ 1 − ( 1 − p g ) k ′ ≡ A ( T ) ϕ k ′
Therefore, it turns out that temperature and concentration is separable for irreversible gelation. The temperature factor is
A ( T ) = α ( T ) f n k ′ k k ′ [ f ′ k ′ / ( f ′ k ′ − 1 ) ] k ′ − 1
It goes back to (28) for the irreversible pairwise cross-linking of k = 2 .
Let us study reversible rate equation. Because rigorous integration of (42) is difficult, let us employ the linear approximation. The linearized kinetic equation leads to the solution
p ˜ ( τ ) = p ¯ ( 1 − e − γ τ )
where
γ ≡ − g ′ ( p ¯ ) = k ( 1 − p ¯ ) k ′ + β ′
and hence the relaxation time is
1 t R = α k ψ k ′ γ = β 1 + k 2 K k z ¯ k ′ = β 1 + k p ¯ k ( z ¯ ) p ¯ 1 ( z ¯ )
by using equilibrium values of the reactivities. They can be written as functions of the solution z ¯ of the conservation equation (6) with
u ( z ¯ ) ≡ 1 + k K k z ¯ k ′
The relaxation time similar form to (55) was first proposed by Kresheck et al[51] (referred to as KHDS), and later by Muller[52] for the self-assembled micelle formation in solutions of surfactant molecules studied by temperature-jump experiments. The KHDS form was derived under the assumption that, in the stepwise association of the surfactant molecules, the last step is slowest compared to the intermediate steps, and hence it is the rate-controlled step.
Substituting the gel-point condition (48) into the linear solution (53), we find the gelation time takes the same form as (21) but with the different thermodynamic factor
Q ( z ¯ ) = ln f ′ k ′ p ¯ k ( z ¯ ) f ′ k ′ p ¯ k ( z ¯ ) − 1
The equilibrium gel point is given by the condition 1 − f ′ k ′ p ¯ k ( z ¯ ) = 0 and hence
( K k ( T ) ψ k ′ ) g = ( f ′ k ′ ) k ′ k ( f ′ k ′ − 1 ) k
We find again a logarithmic divergence of the thermodynamic factor near the sol–gel transition line.
To find specific results, let us assume that the equilibrium constant takes a form[30]
K k ( T ) = λ ( T ) k ′
where λ ( T ) is a binding constant per one functional group. The conservation law becomes
x ≡ λ ψ = z ¯ ( 1 + k z ¯ k ′ )
where z ¯ ≡ x ( 1 − p ¯ k ) is the scaled concentration of the free groups. Solving this equation for z ¯ , and substituting the solution to
τ g ≡ α k ψ k ′ t g = k x k ′ Q ( z ¯ ) 1 + k p ¯ k ( z ¯ ) / p ¯ 1 ( z ¯ )
we find the gelation time as functions of x. Let us plot it in our standard form
Γ ( x ) ≡ 1 β t g = ( 1 + k 2 z ¯ k ′ ) ln f ′ k ′ k z ¯ k f ′ k ′ k z ¯ x k − 1
where z ¯ = z ¯ ( x ) is the solution of the conservation law (60).
Figure 4 (b) plots 1 / β t g , 1 / β t R , and Q as functions of x for varied multiplicity k for tetra-functional molecules f = 4 as a typical example. Multiplicity k = 2 goes back to the pairwise association studied in the previous section. All curves diverge logarithmically as x approaches the critical value x g = f ′ k ′ / k 1 / k ′ ( f ′ k ′ − 1 ) k / k ′ from above. For f = 4 , we know that the gel point is not a monotonically decreasing function of k, but takes a minimum at k = 3 (see Figure 7.5 (a) in the reference[4]).
The retardation coefficient R f , k for f-functional primary molecules with cross-link multiplicity k can be found in a similar method of expanding t g in powers of the reverse rate constant β . We find
t g = t g ( IR ) 1 + R f , k x k ′ + O 1 x 2 k ′
with
R f , k = ( 1 − p g ) 2 k − 1 + ( 2 k − 1 ) p g − 1 2 k ( 2 k − 1 ) [ 1 − ( 1 − p g ) k ′ ] ( 1 − p g ) k
where t g ( IR ) is the irreversible gelation time (51) by Tobitani–Ross-Murphy[27]. Detailed calculation is given in Appendix B.

2.3. Stepwise Association

In most of physical gels, we expect that cross-linking proceeds via step-by-step association of the free functional groups:
A 1 + A k − 1 ⇆ A k ( k = 2 , 3 , ⋯ s )
d n k ( t ) d t = J k − J k + 1
= − { β k + α k n 1 ( t ) } n k ( t ) + β k + 1 n k + 1 ( t ) + α k − 1 n 1 ( t ) n k − 1 ( t )
In some physical gels, a particular value of the multiplicity k is most stable, while in others cross-links are monotonically destabilized with increasing k. The gelation time and physical properties of the networks thus depend on the stepwise association constants. Let us therefore study thermoreversible gelation by such step-by-step cross-linking for a given set of stepwise constants as a model reversible gelation (see Figure 5).
The rate equation is described by
d n k ( t ) d t = J k − J k + 1
= − { β k + α k n 1 ( t ) } n k ( t ) + β k + 1 n k + 1 ( t ) + α k − 1 n 1 ( t ) n k − 1 ( t )
for k ≥ 2 , and
d n 1 ( t ) d t = − 2 J 2 − ∑ k = 3 s J k
= − 2 α 1 n 1 ( t ) 2 + 2 β 2 n 2 ( t ) + ∑ k = 3 s { β k n k ( t ) − α k − 1 n 1 ( t ) n k − 1 ( t ) }
for k = 1 . Here,
J k ( t ) ≡ α k − 1 n 1 ( t ) n k − 1 ( t ) − β k n k ( t )
for k ≥ 2 is the flux between k − 1 and k state ( J 1 ≡ 0 ), i.e., the number of k − 1 junctions changing to k junctions per unit time as a result of forward and backward reaction with the rate constant α k and β k (see Figure 5).
Similar rate equations were proposed by Aniansson and Wall[53,54,55] for the study of micelle formation in surfactant solutions. In our gelation problem, associating groups are attached to the polymer main chains (indicated by arrows), while in micellization problem they move freely. But within the assumption of equal reactivity, the basic kinetics governing the association process can be regarded as fundamentally the same in both cases.
Let us first transform these equations by using reactivity p k ( t ) in order to find the equilibrium solution easily. Using equation (3), we find
d p k ( t ) d t = − { β k + α k ψ p 1 ( t ) } p k ( t ) + k k + 1 β k + 1 p k + 1 ( t ) + k k − 1 α k − 1 ψ p 1 ( t ) p k − 1 ( t )
for k ≥ 2 , and
d p 1 ( t ) d t = − 2 α 1 ψ p 1 ( t ) 2 + β 2 p 2 ( t ) + ∑ k = 3 s β k k p k ( t ) − α k − 1 k − 1 ψ p 1 ( t ) p k − 1 ( t )
for k = 1 .
Let us next confirm that the equilibrium distribution (5) satisfies these kinetic equations. On substitution of (5) into (69), we find that the relation
( β k K k − α k − 1 K k − 1 ) − ( α k K k − β k + 1 K k + 1 ) z ¯ = 0
must be fulfilled. We therefore propose the detailed balance condition such that the equlibrium constants satisfy the relation
K k = α k − 1 β k K k − 1
By repeated use of this relation, we have a well known relation
K k = α k − 1 α k − 2 β k β k − 1 K k − 2 = ⋯ = λ k λ k − 1 ⋯ λ 2
between the equilibrium constants K k and the stepwise association constants
λ k ≡ α k − 1 β k
for k − 1 and k state. The kinetic equations (69) and (70) together with the detailed balance conditions (73) provide a complete set to find the solution for the gelation time. The gel-point condition is given by (48).

3. Results

3.1. Three-State Model

First we study the three-state model of cross-linking. In this model we have free ( k = 1 ), double ( k = 2 ), and triple ( k = 3 ) cross-link junctions. Some biopolymers form cross-links of either double helix or triple helix[2,56,57]. Depending on the environmental condition, there is a competition between them. We can study such competition in forming cross-links by using three-state rate equation.
The kinetic equations are
d p 1 d t = − 2 α 1 p 1 2 + ( β 2 − 1 2 α 2 ψ p 1 ) p 2 + β 3 3 p 3
d p 2 d t = 2 α 1 ψ p 1 2 − ( β 2 + α 2 ψ p 1 ) p 2 + 2 β 3 3 p 3
d p 3 d t = 3 α 2 2 ψ p 1 p 2 − β 3 p 3
The conservation law of the total number of functional groups is described by the normalization condition p 1 + p 2 + p 3 = 1 . The equilibrium solutions are given by ψ p ¯ 1 = z ¯ , ψ p ¯ 2 = 2 K 2 z ¯ 2 , ψ p ¯ 3 = 3 K 3 z ¯ 3 . Detailed balance conditions lead to the relation K 2 = α 1 / β 2 ≡ λ 2 , and K 3 = α 1 α 2 / β 2 β 3 ≡ λ 2 λ 3 .
Let us first eliminate triple association p 3 by using the normalization condition. And then we transform the kinetic equation in terms of the deviation ξ 1 , ξ 2 from the equilibrium state defined by p k ( t ) ≡ p ¯ k ( 1 + ξ k ( t ) ) . After a lengthy calculation, we find
d d t ξ 1 ξ 2 = − a 1 , 1 a 1 , 2 a 2 , 1 a 2 , 2 ξ 1 ξ 2 + f 1 ( ξ 1 , ξ 2 ) f 2 ( ξ 1 , ξ 2 )
where
a 1 , 1 ≡ β 3 3 + 4 α 1 z ¯ + α 2 K 2 z ¯ 2
a 1 , 2 ≡ K 2 z ¯ 2 β 3 3 − 2 β 2 + α 2 z ¯
a 2 , 1 ≡ 1 K 2 z ¯ β 3 3 − 2 α 1 z ¯ + α 2 K 2 z ¯ 2
a 2 , 2 ≡ β 2 + β 3 2 3 + λ 3 z ¯
for the linear terms, and
f 1 ( ξ 1 , ξ 2 ) ≡ − ( 2 α 1 ξ 1 + α 2 K 2 z ¯ ξ 2 ) z ¯ ξ 1
f 2 ( ξ 1 , ξ 2 ) ≡ ( α 1 ξ 1 − α 2 K 2 z ¯ ξ 2 ) ξ 1 / K 2
for the nonlinear (quadratic) terms.
We can numerically solve these equations for a given set of rate constants. But here we confine to the linear approximation in order to obtain physical picture on the factors controlling the gelation time. Nonlinear terms are used only for the stability analysis of the linear solution.
Let us find the eigen-values γ of the matrix A ^ ≡ ( a i , j ) . The equation to find them is
γ 2 − A ( z ¯ ) γ + B ( z ¯ ) = 0
where
A ( z ¯ ) ≡ Tr ( A ^ ) = β 2 + β 3 + ( 4 α 1 + α 2 ) z ¯ + α 2 K 2 z ¯ 2
B ( z ¯ ) ≡ Det ( A ^ ) = β 2 β 3 + 4 α 1 β 3 z ¯ + 9 α 1 α 2 z ¯ 2
Note that B ( z ¯ ) can be written as
B ( z ¯ ) = β 2 β 3 ( z ¯ u ( z ¯ ) ) ′ = β 2 β 3 1 + 2 p ¯ 2 p ¯ 1 + 3 p ¯ 3 p ¯ 1
so that it is regarded as an extension of KHDS decay rate (55) to the three-state model. (Prime ′ indicates differentiation.) The linear relaxation has two modes
γ ( ± ) ( z ¯ ) ≡ A ( z ¯ ) 2 1 ± 1 − 4 B ( z ¯ ) A ( z ¯ ) 2
The gel point is given by the condition
1 − f ′ { p 2 ( t ) + 2 p 3 ( t ) } = 0
Using the deviation ξ k ( t ) , this condition is transformed to
ξ 1 ( t ) + K 2 z ¯ ξ 2 ( t ) = K 2 + 3 K 3 z ¯ − 1 / 2 f ′ z ¯

3.1.1. Rate-Determining Step

We first look into the cases where the two rate constants of the backward reaction are largely different. If β 2 < < β 3 , the cross-linking is controlled by the first step association (pair formation), and vice versa. In both cases we have the situation B / A 2 < < 1 , so that the relaxation time is approximately given by
1 t R ≃ β 2 β 3 ( z ¯ u ( z ¯ ) ) ′ β 2 b 2 ( z ¯ ) + β 3 b 3 ( z ¯ ) = β ¯ ( z ¯ ) 1 + 2 p ¯ 2 p ¯ 1 + 3 p ¯ 3 p ¯ 1
where β ¯ is the average of β defined by
1 β ¯ ( z ¯ ) ≡ b 2 ( z ¯ ) β 2 + b 3 ( z ¯ ) β 3
The relaxation time apparently takes an extended KHDS form. Here, we have eliminated α k in favor of β k and λ k from A ( z ¯ ) , and write
A ( z ¯ ) = β 2 b 2 ( z ¯ ) + β 3 b 3 ( z ¯ )
by using
b 2 ( z ¯ ) ≡ 1 + 4 λ 2 z ¯
b 3 ( z ¯ ) ≡ 1 + λ 3 z ¯ + λ 2 λ 3 z ¯
If the first association (pairing) is the rate-determining step, the relaxation time is approximately
1 t R ≃ β 2 b 3 ( z ¯ ) ( z ¯ u ( z ¯ ) ) ′
and goes back to the pairwise association in the limit of λ 3 → 0 . If the second associaition (triple junction) is the rate-determining step, the relaxation time is approximately given by
1 t R ≃ β 3 b 2 ( z ¯ ) ( z ¯ u ( z ¯ ) ) ′
From these results, we may expect the extended KHDS form for the multiplicity more than three
1 t R ≃ β 1 + 1 p ¯ 1 ∑ k ≥ 2 k p ¯ k = β 〈 k 〉 w p ¯ 1
with some average β , because the average multiplicity is related to the average branching number of the cross-link junctions.

3.1.2. Quasi-Stationary State Approximation

Let us next consider the special case where the first step reaction is so fast that a stationary state is easily established for the binary cross-linking. Since we have d ξ 2 / d t ≃ 0 , the kinetic equation reduces to
d ξ 1 d t ≃ − B ( z ¯ ) a 2 , 2 ξ 1
The rexation time is then given by
γ ( z ¯ ) = β 2 β 3 a 2 , 2 ( z ¯ u ( z ¯ ) ) ′
The front factor can be regarded as an average backward rate constant β ¯ defined by
1 β ¯ ( z ¯ ) ≡ 1 β 3 + 2 / 3 + λ 3 z ¯ β 2
Hence we find again the same form
1 t R = β ¯ ( z ¯ ) 1 + 2 p ¯ 2 p ¯ 1 + 3 p ¯ 3 p ¯ 1
for the relaxation time, with slightly different average of β .

3.1.3. Slow-Mode Approximation

The general solution for the linear equation is given by
ξ 1 ξ 2 = C 1 e 1 e − γ ( − ) t + C 2 e 2 e − γ ( + ) t
where e 1 is the normalized eigen-vector for the eigen-value γ ( − ) (slow mode) and e 2 for the fast mode. The constants C 1 , C 2 can be fixed by the initial conditions p 1 ( 0 ) = 1 , p 2 ( 0 ) = 0 . In the slow mode approximation, we take only the slow mode into consideration, and neglect the contribution from the fast mode. Then the relaxation time is given by 1 / t R = γ ( − ) , and the gelation time is t g = t R Q ( z ¯ ) with a thermodynamic factor
Q ( z ¯ ) = ln 2 f ′ z ¯ D 1 ( z ¯ ) f ′ ( 〈 k 〉 w − 1 ) − 1
from (84), where D 1 ( z ¯ ) ≡ C 1 ( e 1 1 + K 2 z ¯ e 2 1 ) is a function of z ¯ which is non-singular across the gel point. Therefore, we have again logarithmic divergence of the gel time.

3.1.4. Some Numerical Results of the Relaxation Time

In order to find specific results, let us assume that the association constant of the second step is μ times larger than the first one:
λ 2 ≡ λ , λ 3 = μ λ
Then the forward rate constants are given by
α 1 = β 2 λ , α 2 = μ λ β 3
The conservation law takes the form
x ≡ λ ψ = z ¯ ( 1 + z ¯ + μ z ¯ 2 )
where λ z ¯ is written as z ¯ for simplicity. Now we have
A ( z ¯ ) = β 2 ( 1 + 4 z ¯ ) + β 3 ( 1 + μ z ¯ + μ z ¯ 2 )
B ( z ¯ ) = β 2 β 3 ( 1 + 4 z ¯ + 9 μ z ¯ 2 )
In particular, for a special case of β 2 = β 3 ≡ β the relaxation time is given by
1 β t R = 1 2 A ( z ¯ ) 1 − 1 − 4 B ( z ¯ ) A ( z ¯ ) 2
with
A ( z ¯ ) = 2 + ( μ + 4 ) z ¯ + μ z ¯ 2
B ( z ¯ ) = 1 + 4 z ¯ + 9 μ z ¯ 2
Figure 6 (a) shows a schematic free energy surface of association drawn against the reaction coordinates for the three-state model. With decrease in the ratio μ , the surface curve of the triple junction changes from A to D. For the limit of μ = 0 (curve D), the model reduces to the pairwise one. Figure 6 (b) plots The reciprocal gelation time 1 / β t g , the reciprocal relaxation time 1 / β t R , and the thermodynamic factor Q ( x ) plotted against a scaled concentration x ≡ λ ( T ) ψ for three ratios μ = 0 , 1 , 2 of the stepwise constants for the tetra-functional molecules ( f = 4 ). Due to the existence of k = 2 state, relaxation is slower than the fixed-multiplicity model of k = 3 .

3.2. Micellar Cross-Linking

Let us next consider cross-link junctions whose multiplicity is stabilized around a large number k 0 due to the physical reason of associative force, such as geometrical packing, saturation of the interaction force, etc. Typical examples of thermoreversible gels cross-linked by such micellar junctions are hydrophobically-modified water-soluble polymers (associating polymers)[4], for instance, aqueous solutions of polymers with short hydrophobic chains attached at their both chain ends (telechelic polymers), such as hydrophobic ethoxylated urethane (called HEUR)[58,59,60], hydrophobic poly(N-isopropylacrylamide)[61,62,63,64]. Hydrophobes form spherical micelles and serve as the network junctions. Triblock coploymers of the type ABA, ABC, etc, for instance, poly(propylene oxide)-poly(ethylene oxide)-poly(propylene oxide) triblock copolymers[65,66,67], are another important examples whose cross-links are micelles of various shape, spherical, cylindrical, planar etc, consisting of the block segments.
The size distribution of the cross-links changes as reaction proceeds step by step, from unimer (free functional group) dominant one to the final micelle (large junction) dominant one. Such time-development is similar to that observed for the micelle formation in surfactant solutions[51,52,53,54]. The functional groups treated here are, however, attached to the polymer main chains. Therefore, their association takes place under the constraint of chain movement. But here we assume the equal reactivity, and start from the stepwise kinetic equations (69) and (70). Their linearized forms are
d ξ k d t = − ( β k + α k z ¯ ) + β k ξ k − 1 + α k z ¯ ξ k + 1 + ( β k − α k z ¯ ) ξ 1 ( k ≥ 2 )
d ξ 1 d t = − 2 α 1 z ¯ ( 2 ξ 1 − ξ 2 ) − ∑ k ≥ 3 β k K k z ¯ k ( ξ k − 1 − ξ k + ξ 1 )
The gel-point condition (48) turns into
f ′ ∑ k ≥ 2 ( k − 1 ) p ¯ k ξ k ( t g ) = 1 − f ′ ∑ k ≥ 2 ( k − 1 ) p ¯ k
Following the idea by Aniansson and Wall[53](referred to as AW), we take an analogy between micelle formation and a heat transfer from one metal to another through a connected thin metal wire. Let us choose judiciously a multiplicity k 1 below which the population distribution is dominant in the initial stage, and k 2 ( < k 0 ) above which the population is dominant in the final stage. Stationary state is approximately retained between them for k 1 + 1 ≤ k ≤ k 2 in the most time during the approach to equilibrium (assumption of quasi-stationary flow). We then have
J k ≃ 0 ( 1 ≤ k ≤ k 1 ) independent of k ≡ J ( k 1 + 1 ≤ k ≤ k 2 ) 0 ( k 2 + 1 ≤ k )
In Appendix C, we show that the flux J is proportional to the probability of free groups as
R J ≃ m 0 + M 1 M 0 ξ 1
where m j and M j are the j-th moments for small aggregates and large micelles defined by
m j ≡ ∑ k = 1 k 1 k j p ¯ k
M j ≡ ∑ k = k 2 + 1 ∞ k j p ¯ k
for j = 0 , 1 , 2 ⋯ . Here the resistance R in the heat flow analogy is given by
R ≡ ∑ k = k 1 + 1 k 2 1 β k K k z ¯ k = 1 ψ ∑ k = k 1 + 1 k 2 k β k p ¯ k
We also find that the probability of large micelles for k ≥ k 2 + 1 is approximately given by
ξ k ≃ k ξ 1 − R J
The kinetic equation for k = 1 then takes the form
d ξ 1 d t = − 2 J 2 − ∑ k = 3 ∞ J k ≃ − m 1 + M 1 R M 0 ( k 2 − k 1 ) ξ 1
Hence, the relaxation time is
1 t R = ( k 2 − k 1 ) m 1 R M 0 1 + M 1 m 1
which is similar to AW formula. The solution of (111) is given by
ξ 1 ( t ) = C e − t / t R
with C ≡ ( 1 − p ¯ 1 ) / p ¯ 1 .
Substituting (110) into the gel-point condition (105), we have
f ′ ∑ k = 2 k 1 ( k − 1 ) k p ¯ k + ∑ k ≥ k 2 + 1 ( k − 1 ) k p ¯ k ξ 1 ( t g ) − f ′ R J ∑ k ≥ k 2 + 1 p ¯ k = 1 − f ′ ∑ k ≥ 2 ( k − 1 ) p ¯ k
and hence we finally find the form
t g = t R Q ( z ¯ )
for the gelation time with the thermodynamic factor
Q ( z ¯ ) = ln f ′ C M f ′ ∑ k ≥ 2 ( k − 1 ) p ¯ k − 1
where
M ≡ M 1 M 0 ( m 1 + M 1 ) − m 2 − M 2
To see the relaxation time more specifically, we consider the special case in which the multiplicity of the cross-links is limited from above such that only k = 1 , 2 , ⋯ s are allowed. Then, the last step from s − 1 to s is expected to be slower than other steps, and hence we can choose k 1 = 2 and k 2 = s − 1 . Figure 7 shows the quasi-stationary approximation schematically. Circles connected by dotted lines indicate the probability deviation ξ k ( t ) at time t. Since the moments are reduced to m 1 = m 2 = p ¯ 1 and M 0 = p ¯ s , M 1 = s p ¯ s , we have
1 t R = ( s − 2 ) p ¯ 1 R p ¯ s 1 + s p ¯ s p ¯ 1
This is again KHDS form with an effective rate constant
β eff = ( s − 2 ) p ¯ 1 R p ¯ s
written in terms of the resistance constant
R ≡ ∑ k = 2 s − 1 k β k p ¯ k = ∑ k = 2 s − 1 1 β k K k z ¯ k
Let us compare this with the relaxation time of the fixed multiplicity model (55) whose association takes place simultaneously. In the stepwise association, cross-link junctions grow one by one from smaller to larger until they are saturated. All intermediate states must reach equilibrium before the last one (largest junction), and hence they serve as a resistance for the solution to go back to equilibrium. As a result, the dissociation rate β is replaced by an effective one β eff . Relaxation is much slower, in particular, for high concentration region.
Figure 8 shows the gelation time t g , the relaxation time t R , and the thermodynamic factor Q plotted against the scaled concentration x ≡ λ ( T ) ϕ for the functionality f = 4 with the varied upper bound s = 3 , 5 , 7 , 9 , 11 of the junction multiplicity. For simplicity, uniform β k is assumed. The case of s = 3 goes back to the three-state model studied above. With increase of the upper limit s, relaxation becomes slower because the resistance R due to the stepwise time development increases by the existence of many intermediate states. This tendency is opposite to the simultaneous formation of the fixed-multiplicity junctions by one step, for which relaxation becomes faster with k. In the high concentration region, the reciprocal relaxation time behaves with a power index 2 / s as 1 / t R ∼ x 2 / s . We can study these results in more detail by comparing the experimental data of associating polymer solutions measured by using temperature-jump technique.

4. Discussion

Throughout this paper, we have confined our study to the pre-gel region t ≤ t g before the solution reaches the gel point. In this region, the conservation law (6) of the functional groups, equilibrium reactivities p ¯ k (5), and the gel-point condition (48) all hold without ambiguity as they are given in the text. From the equilibrium constants, we can derive the detailed balance condition for the rate constants in the kinetic equation. However, in the post-gel region t g ≤ t where gel networks exist, there is a possibility such that the reaction within the gel part may be different from that of the sol part.
There are several treatments of the reaction in the post-gel region 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[40], and the latter by Flory[48,49,50]. Later Ziff and Stell[42,43,44,45,46] examined another possibility from a kinetic point of view. We have, however, avoided this problem because, for the estimation of the gelation time, we need only information on the reactivities before the gel point. After the gel point is passed, the solution, in particular its sol part, may change very differently depending on the mechanism of cross-linking reaction.

5. Conclusions

We have presented a very general theoretical framework for the study of the gelation time of thermoreversible gels with specified multiple structures of cross-link junctions. It is based on the kinetic equation for the stepwise association of functional groups. All results are presented from a unified point of view in terms of the gelation-time diagram — plot of 1 / t g , 1 / t R and Q against the scaled concentration variable x ≡ λ ( T ) ψ .
From the theoretical modeling, the following conclusions can be drawn:
(1)
The gelation time t g , the relaxation time t R , and the thermodynamic factor Q are all functions of a single variable x ≡ λ ( T ) ψ (scaled concentration), where λ ( T ) is the stepwise association constant at the final temperature T at which cross-linking reaction proceeds. Therefore, temperature and concentration are not separable, but give the same effect if they are properly scaled under a fixed value of x.
(2)
These three factors obey a fundamental relation t g ( x ) = t R ( x ) Q ( x ) . The thermodynamic factor Q ( x ) is logarithmically singular at the equilibrium gel point x = x g , while the relaxation time t R ( x ) is continuous across the gel point. They are calculated for some important models of cross-link junctions, such as pairwise cross-linking, three-state model, cross-linking with fixed high multiplicity, and micellar cross-linking.
(3)
The gelation time t g ( x ) of reversible cross-linking approaches that of irreversible one in the asymptotic region of large x (either high concentration ϕ or high values of the association constant λ ( T ) ). The power index of 1 / t g lies at somewhere between k − 1 (simultaneous cross-linking) and 2 / s (stepwise cross-linking).
(4)
For a large micellar cross-link junctions, the gelation time is derived on the basis of the quasi-stationary approximation (Aniansson–Wall formula) for the relaxation time. Combination with the singular part of the thermodynamic factor estimated by our preceding equilibrium gelation theory provides an accurate estimation of the gelation time.
The model solutions proposed in this study have obvious advantages in finding the microscopic parameters regarding the cross-linking reaction, such as stepwise rate constants and cross-link multiplicity, from macroscopic measurements on the gelation time and the relaxation time.
Our theoretical framework may directly be applicable to some important thermoreversible gels for which the equilibrium sol–gel transition lines are established. We hope detailed experimental data on the gelation time for the systems treated here will be reported in the near future.

Appendix A. Solution of the Kinetic Equation

Because of the second order nature of the reaction, we have a factorized form
g ( p ) = ( p − η ( + ) ) ( p − η ( − ) )
with equation (15). Equation (11) can easily be integrated. The solution is given by
p ( τ ) = [ η ( + ) − p ( 0 ) ] η ( − ) − [ η ( − ) − p ( 0 ) ] η ( + ) e − γ τ [ η ( + ) − p ( 0 ) ] − [ η ( − ) − p ( 0 ) ] e − γ τ
where γ ≡ η ( + ) − η ( − ) is the relaxation rate, and p ( 0 ) is the initial value of the reactivity. We assume for simplicity that the initial temperature is sufficiently high, so that all functional groups are free p ( 0 ) = 0 . Hence we have the result (14).

Appendix B. Retardation Coefficient

To find the retardation coefficient and investigate the behavior of the gelation time in high concentration region, we consider the integral of the equation (42)
τ = ∫ 0 p d p ( 1 − p ) k − β ′ p
where β ′ = 1 / k x k ′ . Expanding the integrand in powers of β ′ , and fixing the reactivity at its gel-point value p g = 1 / f ′ k ′ , we have
τ g = ∑ m = 0 ∞ I m ( p g ) β ′ m
where integrals are defined by
I m ( p g ) ≡ ∫ 0 p g p m d p ( 1 − p ) k ( m + 1 )
Then, the gelation time takes the form
β t g = I 0 ( p g ) k x k ′ 1 + R f , k x k ′ + O 1 x 2 k ′
and hence we have
R f , k = I 1 ( p g ) / k I 0 ( p g )
for the retardation coefficient. We have already found it in (29) for the pairwise cross-linking k = 2 . For a general k, it is explicitly given by equation (64). The first order correction to gelation time is proportional to 1 / x k − 1 .

Appendix C. Quasi-Stationary Approximation

Because the flux J k is given by
J k ≃ β k K k z ¯ k ( ξ k − 1 − ξ k + ξ 1 ) for k ≥ 2
in the linear approximation, we find for k ≥ k 2 + 1
∑ k = 2 k ( ξ k − 1 − ξ k + ξ 1 ) = ∑ k = 2 k J k β k K k z ¯ k ≃ J ∑ k = k 1 + 1 k 2 1 β k K k z ¯ k
The left hand side is simply − ξ k + k ξ 1 by cancellation, so that we have
ξ k ≃ k ξ 1 − R J
where R defined by (109) is the analogy of the resistance in the heat flow. Similarly, for k ≤ k 1 , we find
ξ k ≃ k ξ 1
To find J, we substitute these relations into the materials conservation law ∑ p k = 1 . We find
∑ k = 1 k 1 p ¯ k ( k ξ 1 ) + ∑ k = k 1 + 1 ∞ p ¯ k ( k ξ 1 − R J ) ≃ 0
and hence we have (107).

References

  1. Guenet, J. M. Thermoreversible Gelation of Polymers and Biopolymers., 2nd ed.; Academic Press, Harcourt Brace Jovanovich Publishers: London, 1992. [Google Scholar]
  2. te Nijenhuis, K. Thermoreversible Networks. Adv. Polym. Sci. 1997, 130, 1–252. [Google Scholar]
  3. Winter, H. H.; Mours, M. Rheology of Polymers near Liquid–Solid Transitions. Adv. Polym. Sci. 1997, 134, 165–234. [Google Scholar]
  4. Tanaka, F. Polymer Physics—Applications to Molecular Association and Thermoreversible Gelation. Cambridge University Press: Cambridge, 2011.
  5. Zhang, J.; Hu, Y.; Li, Y. Gel Chemistry: Interactions, Structures and Properties. Springer: Singapore, 2018.
  6. Thakur, V. K.; Thakur, M. K. Polymer Gels: Science and Fundamentals., 1st ed.; Springer: Singapore, 2018. [Google Scholar]
  7. Thakur, V. K.; Thakur, M. K. Hydrogels: Recent Advances., 1st ed.; Springer: Singapore, 2018. [Google Scholar]
  8. Tan, H. M.; Moet, A.; Hiltnet, A.; Baer, E. Thermoreversible Gelation of Atactic Polystyrene Solutions. Macromolecules 1983, 16, 28–34. [Google Scholar] [CrossRef]
  9. Boyer, R. F.; Baer, E.; Hiltner, A. Concerning Gelation Effects in Atactic Polystyrene Solutions. Macromolecules 1985, 18, 427–434. [Google Scholar] [CrossRef]
  10. Jelich, L. M.; Nunes, S. P.; Paul, N. E.; Wolf, B. A. On the Cooccurrence of Demixing and Thermoreversible Gelation of Polymer Solutions. 1. Experimental Observations. Macromolecules 1987, 20, 1943–1947. [Google Scholar] [CrossRef]
  11. Schneider, T.; Wolf, B. A.; Kasten, H.; Kremer, F. Thermoreversible Gelation and Vitrification of Highly Concentrated Polymer Solutions under Poor Thermodynamic Conditions. Macromolecules 1991, 24, 5387–5392. [Google Scholar] [CrossRef]
  12. Mutin, P. H.; Guenet, J. M. Physical Gels from PVC: Aging and Solvent Effects on Thermal Behavior, Swelling, and Compression Modulus. Macromolecules 1989, 22, 843–848. [Google Scholar] [CrossRef]
  13. Arvidson, S. A.; Lott, J. R.; McAllister, J. W.; Zhang, J.; Bates, F. S.; Lodge, T. P. Interplay of Phase Separation and Thermoreversible Gelation in Aqueous Methylcellulose Solutions. Macromolecules 2013, 46, 300–309. [Google Scholar] [CrossRef]
  14. Nakano, S.; Ogiso, T.; Kita, R.; Shinyashiki, N.; Yagihara, S.; Yoneyama, M.; Katsumoto, Y. Thermoreversible gelation of isotactic-rich poly(?it N-isopropylacrylamide) in water. J. Chem. Phys. 2011, 135, 114903. [Google Scholar] [CrossRef]
  15. Wang, C.; Hashimoto, T.; Chuang, Y.-C.; Tanaka, T.; Chang, Y.-P.; Yang, T.-W.; Huang, M.-T. Physical Gelation of Aqueous Solutions of Atactic Poly(N-isopropylacrylamide). Macromolecules 2022, 55, 9152–9167. [Google Scholar] [CrossRef]
  16. Piazza, R.; Campello, M.; Buzzaccaro, S.; Sciortino, F. Phase Behavior and Microscopic Dynamics of a Thermosensitive Gel-Forming Polymer. Macromolecules 2021, 54, 3897–3906. [Google Scholar] [CrossRef]
  17. Taylor, M. L.; Paul Tomlins, P.; Sahota, T. S. Thermoresponsive Gels. Gels 2017, 3, 4–1. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, K.; Kun Xue, K.; Loh, X. J. Thermo-Responsive Hydrogels: From Recent Progress to Biomedical Applications. Gels 2021, 7, 77–1. [Google Scholar] [CrossRef] [PubMed]
  19. Patrickios, C. S. Amphiphilic Polymer Co-Networks: Synthesis, Properties, Modelling and Applications. Royal Society of Chemistry: London, 2020.
  20. Weiss, R. G.; Terech, P. Molecular Gels: Materials with Self-Assembled Fibrillar Networks. Springer: Dordrecht, 2006; p 978.
  21. Weiss, R. G. Controlling Variables in Molecular Gel Science: How Can We Improve the State of the Art? Gels 2018, 4, 25–1. [Google Scholar] [CrossRef]
  22. Morris, J.; Bietsch, J.; Bashaw, K.; Wang, G. Recently Developed Carbohydrate Based Gelators and Their Applications. Gels 2021, 7, 24–1. [Google Scholar] [CrossRef]
  23. Ohkura, M.; Kanaya, T.; Kaji, K. Gels of poly(vinyl alcohol) from dimethyl sulphoxide / water solutions. Polymer 1992, 33, 3686–3690. [Google Scholar] [CrossRef]
  24. Ohkura, M.; Kanaya, T.; Kaji, K. Gelation rates of poly(vinyl alcohol) solution. Polymer 1992, 33, 5044–5048. [Google Scholar] [CrossRef]
  25. Mal, S.; Maiti, P.; Nandi, A. K. On the Gelation Rates of Thermoreversible Poly(vinylidene fluoride) gels. Macromolecules 1995, 28, 2371–2376. [Google Scholar] [CrossRef]
  26. Hong, P.-D.; Chou, C.-M. Phase separation and gelation behaviors in poly(vinylidene fluoride)/tetra(ethylene glycol) dimethyl ether solutions. Polymer 2000, 41, 8311–8320. [Google Scholar] [CrossRef]
  27. Tobitani, A.; Ross-Murphy, S. B. Heat-Induced Gelation of Globular Proteins.1. Model for the Effects of Time and Temperature on the Gelation Time of BSA Gels. Macromolecules 1997, 30, 4845–4854. [Google Scholar] [CrossRef]
  28. Tobitani, A.; Ross-Murphy, S. B. Heat-Induced Gelation of Globular Proteins.2. Effect of Environmental Factors on Single-Component and Mixed-Protein Gels. Macromolecules 1997, 30, 4855–4862. [Google Scholar] [CrossRef]
  29. 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]
  30. Tanaka, F.; Stockmayer, W. H. Thermoreversible Gelation with Junctions of Variable Multiplicity. Macromolecules 1994, 27, 3943–3954. [Google Scholar] [CrossRef]
  31. Tanaka, F. Thermoreversible Gelation Interfering with Phase Separation in Multicomponent Mixtures of Associating Polymers. Macromolecules 2022, 55, 5233–5248. [Google Scholar] [CrossRef]
  32. Madbouly, S. A.; Otaigbe, J. U. Kinetic Analysis of Fractal Gel Formation in Waterborne Polyurethane Dispersions Undergoing High Deformation Flows. Macromolecules 2006, 39, 4144–4151. [Google Scholar] [CrossRef]
  33. Ponton, A.; Griesmar, P.; Barboux-Doeuff, S.; Sanchez, C. Rheological investigation of the sol–gel transition: effect of hydrolysis variation in silicon oxide and titanium oxide based on matrices. J. Mat. Chem. 2001, 11, 3125–3129. [Google Scholar] [CrossRef]
  34. Ponton, A.; Warlus, S.; Griesmar, P. Rheological Study of Sol-Gel Transition in Silica Alkoxides. J. Coll. Interf. Sci. 2002, 249, 209–216. [Google Scholar] [CrossRef] [PubMed]
  35. Zhou, J.-F.; Johari, G. P. Gelation Time during Polymerization by Ultrasonic Shear Waves Propagation. Macromolecules 1997, 30, 8085–8087. [Google Scholar] [CrossRef]
  36. Normand, V.; Muller, S.; Ravey, J. C.; Parker, A. Gelation Kinetics of Gelatin: A Master Curve and Network Modeling. Macromolecules 2000, 33, 1063–1071. [Google Scholar] [CrossRef]
  37. Sperinde, J. J.; Griffith, L. G. Control and Prediction of Gelation Kinetics in Enzymatically Cross-Linked Poly(ethylene glycol) Hydrogels. Macromolecules 2000, 33, 5746–5480. [Google Scholar] [CrossRef]
  38. Saalwachter, K.; Gottlieb, M.; Liu, R.; Oppermann, W. Gelation as Studied by Proton Multiple-Quantum NMR. Macromolecules 2007, 40, 1555–1561. [Google Scholar] [CrossRef]
  39. Kurakazu, M.; Takuya Katashima, T.; Chijiishi, M.; Nishi, K.; Akagi, Y.; Matsunaga, T.; Shibayama, M.; Chung, U.; Sakai, T. Evaluation of Gelation Kinetics of Tetra-PEG Gel. Macromolecules 2010, 43, 3935–3940. [Google Scholar] [CrossRef]
  40. 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]
  41. 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]
  42. Ziff, R. M.; Stell, G. Kinetics of Polymer Gelation. J. Chem. Phys. 1980, 73, 3492–3499. [Google Scholar] [CrossRef]
  43. Ziff, R. M. Kinetics of Polymerization. J. Stat. Phys. 1980, 23, 241–263. [Google Scholar] [CrossRef]
  44. Ziff, R. M.; Hendriks, E. M.; Ernst, M. H. Critical Properties for Gelation: A Kinetic Approarch. Phys. Rev. Lett. 1982, 49, 593–595. [Google Scholar] [CrossRef]
  45. Ziff, R. M.; Ernst, M. H.; Hendriks, E. M. Kinetics of Gelation and Universality. J. Phys. A: Math. Gen. 1983, 16, 2293–2320. [Google Scholar] [CrossRef]
  46. Dongen, P. G. J. v.; Ernst, M. H. Kinetics of Reversible Polymerization. J. Stat. Phys. 1984, 37, 301–324. [Google Scholar] [CrossRef]
  47. Flory, P. J. , Principles of Polymer Chemistry. Cornell University Press: Ithaca, New York, 1953.
  48. Flory, P. J. Molecular Size Distribution in Three Dimensional Polymers I. Gelation. J. Am. Chem. Soc. 1941, 63, 3083–3090. [Google Scholar] [CrossRef]
  49. 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]
  50. 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]
  51. Kresheck, G. C.; Hamori, E.; Davenport, G.; Scheraga, H. A. Determination of the Dissociation Rate of Dodecylpyridinium Iodide Micelles by a Temperature-Jump Technique18. J. Am. Chem. Soc. 1966, 88, 246–253. [Google Scholar] [CrossRef]
  52. Muller, N. Kinetics of Micelle Dissociation by Temperature-Jump Techniques. A Reinterpretation. J. Phys. Chem. 1972, 70, 3017–3020. [Google Scholar] [CrossRef]
  53. Aniansson, E. A.; Wall, S. N. On the Kinetics of Step-Wise Micelle Association. J. Phys. Chem. 1974, 78, 1024–1030. [Google Scholar] [CrossRef]
  54. Aniansson, E. A.; Wall, S. N.; Almgren, M. Theory of the Kinetics of Micellar Equilibria and Quantitative Interpretation of Chemical Relaxation Studies of Micellar Solutions of Ionic Surfactants. J. Phys. Chem. 1976, 80, 905–922. [Google Scholar] [CrossRef]
  55. Patist, A.; Oh, S. G.; Leung, R.; Shah, D. O. Kinetics of micellization: its significance to technological processes. Coll. Surf. A: Physicochem. Eng. Aspects 2001, 176, 3–16. [Google Scholar] [CrossRef]
  56. Clark, A. H.; Ross-Murphy, S. B. Structural and Mechanical Properties of Biopolymer Gels. Adv. Polym. Sci. 1987, 83, 57–192. [Google Scholar]
  57. Djabourov, M.; Nishinari, K.; Ross-Murphy, S. B. , Physical Gels from Biological and Synthetic Polymers. Cambridge University Press: New York, 2013.
  58. Annable, T.; Buscall, R.; Ettelaie, R.; Whittlestone, D. The Rheology of Solutions of Associating Polymers: Comparison of Experimental Behavior with Transient Network Theory. J. Rheol. 1993, 37, 695–726. [Google Scholar] [CrossRef]
  59. 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]
  60. Yekta, A.; Xu, B.; Duhamel, J.; Adiwidjaja, H.; Winnik, M. A. Fluorescence Studies of Associating Polymers in Water: Determination of the Chain End Aggregation Number and a Model for the Association Process. Macromolecules 1995, 28, 956–966. [Google Scholar] [CrossRef]
  61. Kujawa, P.; Watanabe, H.; Tanaka, F.; Winnik, F. M. Amphiphhilic Telechelic Poly(?it N-isopropylacrylamide) in Water: From Micelles to Gels. Eur. Phys. J. E 2005, 17, 129–137. [Google Scholar] [CrossRef]
  62. Kujawa, P.; Segui, F.; Shaban, S.; Diab, C.; Okada, Y.; Tanaka, F.; Winnik, F. M. Impact of End-Group Association and Main-Chain Hydration on the Thermosensitive Properties of Hydrophobically Modified Telechelic Poly(?it N-isopropylacrylamide) in Water. Macromolecules 2006, 39, 341–348. [Google Scholar] [CrossRef]
  63. Kujawa, P.; Tanaka, F.; Winnik, F. M. Temperature-Dependent Properties of Telechelic Hydrophobically Modified Poly(?it N-isopropylacrylamides) in Water: Evidence from Light Scattering and Fluorescence Spectroscopy for the Formation of Stable Mesoglobules at Elevated Temperatures. Macromolecules 2006, 39, 3048–3055. [Google Scholar] [CrossRef]
  64. Kujawa, P.; Aseyev, V.; Tenhu, H.; Winnik, F. M. Temperature-Sensitive Properties of Poly(?it N-isopropylacrylamides) Mesoglobules Formed in Dilute Aqueous Solutions Heated above Their Demixing Point in Water: Evidence from Light Scattering and Fluorescence Spectroscopy for the Formation of Stable Mesoglobules at Elevated Temperatures. Macromolecules 2006, 39, 7686–7693. [Google Scholar]
  65. Quellet, C.; Eicke, H.-F.; Xu, G.; Hauger, Y. Transient Networks in ABA Block Copolymer-Microemulsion Systems. Macromolecules 1990, 23, 3347–3352. [Google Scholar] [CrossRef]
  66. Mortensen, K.; Brown, W.; Norden, B. Inverse Melting Transition and Evidence of Three-Dimentional Cubatic Structure in a Block-Copolymer Miceller System. Phys. Rev. Lett. 1992, 68, 2340–2343. [Google Scholar] [CrossRef]
  67. Odenwald, M.; Eicke, H.-F.; Meier, W. Transient Networks in ABA Block Copolymers and Microemulsions: A Rheological Study. Macromolecules 1995, 28, 5069–5074. [Google Scholar] [CrossRef]
Preprints 70841 g001
Figure 1. (a) Schematic diagram for a temperature (or concentration) jump experiment. A polymer solution is kept in equilibrium at sufficiently high temperature T i in the sol state P. At time t = 0 , it is quenched to a low temperature T f ( ≡ T ) in the gel region Q, and then time development of cross-linking reaction is observed. The quench depth is defined by Δ T ≡ T gel − T . For a concentration jump experiment, polymers are quickly mixed with solvent at a constant temperature T to reach the state Q. The concentration depth is defined by Δ ϕ ≡ ϕ − ϕ g . (b) Reactivity of functional groups for pairwise cross-linking schematically shown as a function of the time after the reaction is started: p ( t ) exact solution, p ( IR ) ( t ) its irreversible limit, p ˜ ( t ) linear approximation. The gelation time t g is found by the gel-point condition for the reactivity p ( t g ) = 1 / ( f − 1 ) . Due to the backward reaction, there is a retardation time t g − t g ( IR ) for reversible cross-linking. Dotted lines show the reactivity of the sol part after the gel point is passed. They are considered in this paper on the basis of Flory’s picture.
Figure 1. (a) Schematic diagram for a temperature (or concentration) jump experiment. A polymer solution is kept in equilibrium at sufficiently high temperature T i in the sol state P. At time t = 0 , it is quenched to a low temperature T f ( ≡ T ) in the gel region Q, and then time development of cross-linking reaction is observed. The quench depth is defined by Δ T ≡ T gel − T . For a concentration jump experiment, polymers are quickly mixed with solvent at a constant temperature T to reach the state Q. The concentration depth is defined by Δ ϕ ≡ ϕ − ϕ g . (b) Reactivity of functional groups for pairwise cross-linking schematically shown as a function of the time after the reaction is started: p ( t ) exact solution, p ( IR ) ( t ) its irreversible limit, p ˜ ( t ) linear approximation. The gelation time t g is found by the gel-point condition for the reactivity p ( t g ) = 1 / ( f − 1 ) . Due to the backward reaction, there is a retardation time t g − t g ( IR ) for reversible cross-linking. Dotted lines show the reactivity of the sol part after the gel point is passed. They are considered in this paper on the basis of Flory’s picture.
Preprints 70841 g001
Preprints 70841 g002
Figure 2. (a) Reactivity p ( τ ) for the functionality f = 3 , 4 , 5 plotted against the dimensionless scaled time τ ≡ 2 ψ α t for various scaled concentration x ≡ K 2 ( T ) ψ = 2 , 4 , 6 (solid lines). The relaxation time is independent of f as a function of τ , so that p ( τ ) depends only on x. Broken lines are the reactivity of the sol part p s ( τ ) in the post-gel region. (b) The weight-average molecular weight (broken lines) and the number-average molecular weight (solid lines) plotted against the scaled time τ for the concentration x = 2.0 . The number-averages are finite at the gel point.
Figure 2. (a) Reactivity p ( τ ) for the functionality f = 3 , 4 , 5 plotted against the dimensionless scaled time τ ≡ 2 ψ α t for various scaled concentration x ≡ K 2 ( T ) ψ = 2 , 4 , 6 (solid lines). The relaxation time is independent of f as a function of τ , so that p ( τ ) depends only on x. Broken lines are the reactivity of the sol part p s ( τ ) in the post-gel region. (b) The weight-average molecular weight (broken lines) and the number-average molecular weight (solid lines) plotted against the scaled time τ for the concentration x = 2.0 . The number-averages are finite at the gel point.
Preprints 70841 g002
Preprints 70841 g003
Figure 3. The reciprocal gelation time 1 / β t g (solid lines), reciprocal relaxation time 1 / β t R (broken lines), and thermodynamic factor Q (broken dotted lines) plotted against a combined variable x ≡ K 2 ( T ) ψ for various functionality f. 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point x g = f ′ / 2 f ′ ′ 2 , while it is proportional to x at high concentration region. t R is independent of f, and remains finite at the gel point, while Q goes to infinity.
Figure 3. The reciprocal gelation time 1 / β t g (solid lines), reciprocal relaxation time 1 / β t R (broken lines), and thermodynamic factor Q (broken dotted lines) plotted against a combined variable x ≡ K 2 ( T ) ψ for various functionality f. 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point x g = f ′ / 2 f ′ ′ 2 , while it is proportional to x at high concentration region. t R is independent of f, and remains finite at the gel point, while Q goes to infinity.
Preprints 70841 g003
Preprints 70841 g004
Figure 4. (a) Schematic picture of the fixed-multiplicity cross-linking model with multipliicity k. (b) The reciprocal gelation time 1 / β t g , the reciprocal relaxation time 1 / β t R , and the thermodynamic factor Q ( x ) plotted against a combined variable x ≡ λ ( T ) ψ for various multiplicity k for the tetra-functional molecules ( f = 4 ). 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point given by x g = f ′ k ′ / k 1 / k ′ ( f ′ k ′ − 1 ) k / k ′ , while it approaches the irreversible limit ∼ x k ′ at high concentration region. Note that the gel point concentration changes non-monotonically as a function of the multiplicity.
Figure 4. (a) Schematic picture of the fixed-multiplicity cross-linking model with multipliicity k. (b) The reciprocal gelation time 1 / β t g , the reciprocal relaxation time 1 / β t R , and the thermodynamic factor Q ( x ) plotted against a combined variable x ≡ λ ( T ) ψ for various multiplicity k for the tetra-functional molecules ( f = 4 ). 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point given by x g = f ′ k ′ / k 1 / k ′ ( f ′ k ′ − 1 ) k / k ′ , while it approaches the irreversible limit ∼ x k ′ at high concentration region. Note that the gel point concentration changes non-monotonically as a function of the multiplicity.
Preprints 70841 g004
Preprints 70841 g005
Figure 5. Schematic picture to construct the rate equation for a stepwise association during cross-linking. The reactivity p k is the probability for an arbitrarily chosen functional group to belong to a cross-link junction of the multiplicity k. The flux J k is the number concentration of the junctions growing from k − 1 to k in a unit time.
Figure 5. Schematic picture to construct the rate equation for a stepwise association during cross-linking. The reactivity p k is the probability for an arbitrarily chosen functional group to belong to a cross-link junction of the multiplicity k. The flux J k is the number concentration of the junctions growing from k − 1 to k in a unit time.
Preprints 70841 g005
Preprints 70841 g006
Figure 6. (a) Schematic picture of the three-state cross-linking model with multiplicity k = 1 , 2 , 3 . (b) The reciprocal gelation time 1 / β t g (solid lines), the reciprocal relaxation time 1 / β t R (broken lines), and the thermodynamic factor Q ( x ) (broken dotted lines) plotted against a scaled concentration variable x ≡ λ ( T ) ψ for various ratio μ = 0 , 1 , 2 of the stepwise constants for the tetra-functional molecules ( f = 4 ). The reciprocal gelation time 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point, while it is eventually proportional to x 2 / 3 at high concentration region. With decrease of μ , stability of the triple cross-link junction is weakened as shown from A to D in Figure (a). For μ = 0 , the model reduces to the pairwise cross-linking.
Figure 6. (a) Schematic picture of the three-state cross-linking model with multiplicity k = 1 , 2 , 3 . (b) The reciprocal gelation time 1 / β t g (solid lines), the reciprocal relaxation time 1 / β t R (broken lines), and the thermodynamic factor Q ( x ) (broken dotted lines) plotted against a scaled concentration variable x ≡ λ ( T ) ψ for various ratio μ = 0 , 1 , 2 of the stepwise constants for the tetra-functional molecules ( f = 4 ). The reciprocal gelation time 1 / β t g goes logarithmically to zero near the equilibrium sol–gel transition point, while it is eventually proportional to x 2 / 3 at high concentration region. With decrease of μ , stability of the triple cross-link junction is weakened as shown from A to D in Figure (a). For μ = 0 , the model reduces to the pairwise cross-linking.
Preprints 70841 g006
Preprints 70841 g007
Figure 7. Quesi-stationary state of the stepwise cross-linking for the multiplicity with upper bound. The intermediate cross-link junctions ( k = 2 , ⋯ , s − 1 ) are close to equilibrium with ξ k ≃ 0 , and hence the fluxes J k ( k = 2 , ⋯ , s − 1 ) are kept nearly constant.
Figure 7. Quesi-stationary state of the stepwise cross-linking for the multiplicity with upper bound. The intermediate cross-link junctions ( k = 2 , ⋯ , s − 1 ) are close to equilibrium with ξ k ≃ 0 , and hence the fluxes J k ( k = 2 , ⋯ , s − 1 ) are kept nearly constant.
Preprints 70841 g007
Preprints 70841 g008
Figure 8. The reciprocal gelation time 1 / β t g (solid lines), relaxation time 1 / β t R (dotted lines) and the thermodynamic factor Q (broken dotted lines) plotted against the scaled concentration x ≡ λ ( T ) ϕ for the functionality f = 4 with varied upper limit s = 3 , 5 , 7 , 9 , 11 of the junction multiplicity. With increasing s, the gelation time becomes slower with smaller slope at high concentration region.
Figure 8. The reciprocal gelation time 1 / β t g (solid lines), relaxation time 1 / β t R (dotted lines) and the thermodynamic factor Q (broken dotted lines) plotted against the scaled concentration x ≡ λ ( T ) ϕ for the functionality f = 4 with varied upper limit s = 3 , 5 , 7 , 9 , 11 of the junction multiplicity. With increasing s, the gelation time becomes slower with smaller slope at high concentration region.
Preprints 70841 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated