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

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_{\mathrm{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_{\mathrm{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_{\mathrm{g}}$ and $t_{\mathrm{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_{\mathrm{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_{\mathrm{g}}^{-1}=A\left(T\right){\left(\mathsf{\Delta }\varphi \right)}^n$$

(1)

where $\mathsf{\Delta }\varphi$ is the concentration of PVA measured from the overlap concentration. The critical concentration at which $t_{\mathrm{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\left(T\right)\simeq e^{-B/T{\left(\mathsf{\Delta }T\right)}^2}$$

(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\simeq 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 ($\beta$-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\left(T\right)$ 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_{\mathrm{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) $\alpha$ 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\left(t\right)$ of cross-linking with multiplicity k, we derive the gelation time $t_{\mathrm{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_{\mathrm{g}}=t_{\mathrm{R}}Q$ in terms of the relaxation time $t_{\mathrm{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\mathsf{\equiv }1$ in the unit of the cell volume. We assume incompressibility of the solution, so that we have $\Omega =nN+N_0$ for the total volume. The volume fraction of each component is then given by $\varphi =nN/\Omega$ for the reactive molecule, ${\varphi }_0=N_0/\Omega$ for the solvent. In terms of the functional groups, the number concentration of A-groups on the reactive molecules is $\psi =f\varphi /n$.

We study the time development of the solution after its temperature is changed suddenly from the initial one $T_{\mathrm{i}}$ to the final one $T\mathsf{\equiv }{T}_{\mathrm{f}}$ (temperature jump), or after the molecules are quickly dissolved to the solvent in the preparation of the solution with concentration $\varphi$ at a constant temperature T (concentration jump) (see Figure 1 (a)). The temperature quench depth is defined by $\mathsf{\Delta }T\mathsf{\equiv }{T}_{\mathrm{gel}}-T$, while the concentration quench depth is $\mathsf{\Delta }\varphi \mathsf{\equiv }\varphi -{\varphi }_{\mathrm{g}}$, where $T_{\mathrm{gel}}$ and ${\varphi }_{\mathrm{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\left(t\right)$ 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\left(t\right)$ for an arbitrarily chosen functional group to belong to the junction of multiplicity k is related to the number of junctions by the equation

$$\psi p_k\left(t\right)/k=n_k\left(t\right)$$

(3)

After a long time, the solution reaches its equilibrium state with equilibrium reactivity ${\overline{p}}_k\mathsf{\equiv }{p}_k(t\to \infty )$. Figure 1 (b) schematically shows the time development of $p\left(t\right)\mathsf{\equiv }{p}_2\left(t\right)$, its linearized approximation $\tilde{p}\left(t\right)$, irreversible counterpart $p^{\left(\mathrm{IR}\right)}\left(t\right)$, reactivity of the sol part $p_{\mathrm{s}}\left(t\right)$ 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

$$\frac{\psi {\overline{p}}_k/k}{{\left(\psi {\overline{p}}_1\right)}^k}=K_k\left(T\right)$$

(4)

where $K_k\left(T\right)$ is the equilibrium constant, $p_1\left(t\right)$ is the probability for a functional group to remain free from reaction, and ${\overline{p}}_1$ is its equilibrium value. This equilibrium condition leads to

$$\psi {\overline{p}}_k=k{K}_k{\overline{z}}^k$$

(5)

for the reactivity in terms of the number of free groups $\overline{z}\mathsf{\equiv }\psi {\overline{p}}_1$. From the normalization condition of ${\overline{p}}_k$, we find the conservation law

$$\psi =\overline{z}u\left(\overline{z}\right)$$

(6)

where

$$u\left(\overline{z}\right)=\sum _{k\ge 1}k{K}_k{\overline{z}}^{k-1}$$

(7)

2.1. Pairwise Cross-Linking

Let us first consider the simplest case of pairwise cross-linking reaction

$$2{\mathrm{A}}_1\leftrightarrows {\mathrm{A}}_2$$

(8)

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\mathsf{\equiv }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

$$\frac{d{n}_2\left(t\right)}{dt}=\alpha n_1{\left(t\right)}^2-\beta n_2\left(t\right)$$

(9)

and hence we have

$$\frac{dp\left(t\right)}{dt}=2\psi \alpha {(1-p)}^2-\beta p\left(t\right)$$

(10)

where $\alpha$ and $\beta$ are the rate constant of the forward reaction and backward reaction. (Throughout this paper, we avoid conventional symbols $k_{\mathrm{f}}$ and $k_{\mathrm{b}}$ etc to prevent confusion with multiplicity k.) By using the scaled time $\tau \mathsf{\equiv }2\psi \alpha t$, we have a simple second order equation

$$\frac{dp}{d\tau }=g\left(p\right)$$

(11)

where

$$g\left(p\right)\mathsf{\equiv }{(1-p)}^2-\frac{p}{2K_2\psi }$$

(12)

Here, a new constant

$$K_2=\alpha /\beta$$

(13)

is the equilibrium constant written in terms of the rate constants.

The solution of this equation with the initial condition $p\left(0\right)=0$ is given by

$$p\left(\tau \right)=\frac{{\eta }^{(+)}{\eta }^{(-)}(1-e^{-\gamma \tau })}{{\eta }^{(+)}-{\eta }^{(-)}{e}^{-\gamma \tau }}$$

(14)

where

$${\eta }^{(\pm )}=\frac{1}{4K_2\psi }\left\{1+4K_2\psi \pm \sqrt{1+8K_2\psi }\right\}$$

(15)

and

$$\gamma \mathsf{\equiv }{\eta }^{(+)}-{\eta }^{(-)}=\frac{1}{2K_2\psi }\sqrt{1+8K_2\psi }$$

(16)

The equilibrium reactivity is given by $\overline{p}={\eta }^{(-)}$. Therefore $\gamma$ 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 (${\overline{M}}_{\mathrm{w}}$ and ${\overline{M}}_{\mathrm{n}}$ in the unit of the molecular weight of the primary molecules) are plotted against time. They are explicitly given by[40]

$${\overline{M}}_{\mathrm{w}}\left(\tau \right)=\frac{1+p\left(\tau \right)}{1-f^{\prime }p\left(\tau \right)},{\overline{M}}_{\mathrm{n}}\left(\tau \right)=\frac{1}{1-fp\left(\tau \right)/2}$$

(17)

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 $\tau$ for the combined variable (scaled concentration) $x\mathsf{\equiv }{K}_2\psi =2.0$. The number-average remains finite at the gel point.

Because the relaxation time $t_{\mathrm{R}}$ is defined by the speed to reach the equilibrium state, it is given by $1/\gamma$. In terms of the bare time, we have

$$\frac{1}{t_{\mathrm{R}}}=2\psi \alpha \gamma =\frac{\alpha }{K_2}\sqrt{1+8K_2\psi }=\beta \sqrt{1+8K_2\psi }$$

(18)

In the limit of irreversible reaction where $\beta \to 0$ with finite $\alpha$, $\gamma \tau$ is small, and ${\eta }^{(\pm )}\to 1$, so that we go back to the Stockmayer’s result[40,41](see $p^{\left(\mathrm{IR}\right)}\left(t\right)$ in Figure 1 (b))

$$p\left(t\right)=p^{\left(\mathrm{IR}\right)}=\frac{2\psi \alpha t}{1+2\psi \alpha t}$$

(19)

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_{\mathrm{g}}=1/(f-1)$$

(20)

for the reactivity. Substituting this value for $p\left(\tau \right)$ in (14), we find for the gelation time

$$t_{\mathrm{g}}\left(x\right)=t_{\mathrm{R}}\left(x\right)Q\left(x\right)$$

(21)

where

$$Q\left(x\right)\mathsf{\equiv }ln\left\{\frac{{\eta }^{(-)}\left(x\right)[f^{\prime }{\eta }^{(+)}\left(x\right)-1]}{{\eta }^{(+)}\left(x\right)[f^{\prime }{\eta }^{(-)}\left(x\right)-1]}\right\}$$

(22)

is a function of a combined variable $x\mathsf{\equiv }{K}_2\left(T\right)\psi$ of the temperature and concentration. Since this factor $Q\left(x\right)$ drastically depends upon the quench depth $\mathsf{\Delta }T$, or $\mathsf{\Delta }\varphi$, 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^{\prime }{\eta }^{(-)}\left(x\right)-1=0$, so that $Q\left(x\right)$ goes logarithmically to infinity as x approaches the critical value $x_{\mathrm{g}}=f^{\prime }/2{f^{\prime \prime }}^2$ at the gel point. (In what follows we use abbreviated notations $f^{\prime }\mathsf{\equiv }f-1,f^{\prime \prime }\mathsf{\equiv }f-2,\cdots$ for simplicity.) Near the gel point, we have $x\mathsf{\equiv }{x}_{\mathrm{g}}(1+ϵ)(ϵ<<1)$. A lengthy calculation leads to the form

$$Q\left(x\right)\simeq ln\left\{{\left(\frac{f{f}^{\prime \prime }}{{f^{\prime }}^2}\right)}^2\frac{1}{ϵ}\right\}\simeq -lnϵ+\mathrm{finite}\mathrm{const}.$$

(23)

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_{\mathrm{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

$$\overline{z}=\frac{1}{4K_2}\left\{-1+\sqrt{1+8K_2\psi }\right\}$$

(24)

We therefore have ${\overline{p}}_2/{\overline{p}}_1=2K_2\overline{z}$, and hence

$$\frac{1}{t_{\mathrm{R}}}=\beta \left(1+2\frac{{\overline{p}}_2}{{\overline{p}}_1}\right)$$

(25)

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

$$\psi \alpha t_{\mathrm{g}}=\frac{x}{\sqrt{1+8x}}Q\left(x\right)$$

(26)

and take $\beta \to 0$ limit while $\alpha$ 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_{\mathrm{g}}=\frac{1}{2\alpha f^{\prime \prime }\psi }\left(1+\frac{R_{f,2}}{x}+O\left(\frac{1}{x^2}\right)\right)$$

(27)

where the front factor

$$t_{\mathrm{g}}^{(\mathrm{IR})}\mathsf{\equiv }\frac{1}{2\alpha f^{\prime \prime }\psi }$$

(28)

is the gelation time of irreversible reaction (covalent bonding). It can directly be obtained from Stockmayer’s form (19) by fixing $p\left(t\right)=1/f^{\prime }$. The first correction due to reversible reaction is $O(1/x)$. Its coefficient is found to be

$$R_{f,2}=\frac{3f-4}{6{f^{\prime \prime }}^2}$$

(29)

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

$$\frac{1}{\beta t_{\mathrm{g}}}=\frac{\sqrt{1+8x}}{Q\left(x\right)}\mathsf{\equiv }\Gamma \left(x\right)$$

(30)

as a function of the combined variable x for varied functionality f.

Near the sol–gel transition point $x=x_{\mathrm{g}}$, it goes logarithmically to zero as

$$\Gamma \left(x\right)\simeq \frac{\sqrt{1+8x_{\mathrm{g}}}}{ln(1/\mathsf{\Delta }x)}$$

(31)

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

$$\Gamma \left(x\right)\simeq 2f^{\prime \prime }x$$

(32)

In this limit, the gelation time is separable into a product of the temperature factor $\lambda \left(T\right)$ and the concentration factor $\psi =f\varphi /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\left(t\right)=\overline{p}(1+\xi \left(t\right))$$

(33)

and consider only linear terms of the deviation $\xi \left(t\right)$ in (11). We have

$$\frac{d\xi }{d\tau }=g^{\prime }\left(\overline{p}\right)\xi$$

(34)

where

$$g^{\prime }\left(\overline{p}\right)=-\left\{2(1-\overline{p})+\frac{1}{2K_2\psi }\right\}$$

(35)

is the derivative of $g\left(\overline{p}\right)$. Substituting the equilibrium value of p leads to

$$g^{\prime }\left(\overline{p}\right)=-\frac{1}{2K_2\psi }\left(1+2\frac{{\overline{p}}_2}{{\overline{p}}_1}\right)$$

(36)

Hence the relaxation time in the linear approximation agrees with the exact one. If we assume $p\left(0\right)=0$ , the reactivity in the linear approximation takes the form (see Figure 1 (b))

$$\tilde{p}\left(t\right)=\overline{p}(1-e^{-t/t_{\mathrm{R}}})$$

(37)

The gel-point condition $p\left(t_{\mathrm{g}}\right)=1/f^{\prime }$ then leads to (21), but with $Q\left(x\right)$ of slightly different form

$$Q\left(x\right)=ln\left\{\frac{f^{\prime }{\eta }^{(+)}\left(x\right)}{f^{\prime }{\eta }^{(-)}\left(x\right)-1}\right\}\simeq -lnϵ+\mathrm{finite}\mathrm{const}.$$

(38)

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{\mathrm{A}}_1\rightleftarrows {\mathrm{A}}_k$$

(39)

before we move onto more complex case of stepwise association. A functional group is either free ($k=1$) or reacted ($k\ge 2$). Let us refer to it as fixed-multiplicity model (see Figure 4).

We have a rate equation for the fixed-multiplicity reversible reaction

$$\frac{d{n}_k\left(t\right)}{dt}=\alpha n_1{\left(t\right)}^k-\beta n_k\left(t\right)$$

(40)

Or, equivalently

$$\frac{d{p}_k\left(t\right)}{dt}=\alpha k{\psi }^{k^{\prime }}{p}_1{\left(t\right)}^k-\beta p_k\left(t\right)$$

(41)

($k^{\prime }\mathsf{\equiv }k-1$). Hence we have

$$\frac{dp}{d\tau }=g\left(p\right)$$

(42)

where $p\mathsf{\equiv }{p}_k,p_1=1-p$. Time is scaled as

$$\tau \mathsf{\equiv }\alpha k{\psi }^{k^{\prime }}t$$

(43)

and

$$g\left(p\right)\mathsf{\equiv }{(1-p)}^k-{\beta }^{\prime }p$$

(44)

The coefficient $\beta$ of the reverse reaction is changed to

$${\beta }^{\prime }\mathsf{\equiv }\frac{\beta }{\alpha k{\psi }^{k^{\prime }}}=\frac{1}{k{K}_k{\psi }^{k^{\prime }}}$$

(45)

where

$$K_k=\frac{\alpha }{\beta }$$

(46)

is the equilibrium constant.

In the study of heat-setting $\beta$-Lactoglobulin, Tobitani and Ross-Murphy[27,28] solved this equation in the irreversible limit of $\beta \to 0$. The rate equation can easily be integrated in this limit. We find

$$\tau =\frac{1}{k^{\prime }}\left\{\frac{1}{{(1-p)}^{k^{\prime }}}-\frac{1}{{(1-p\left(0\right))}^{k^{\prime }}}\right\}$$

(47)

The relaxation is not exponential, but obeys a power law. Let us assume that all functional groups are free in the initial state ($p\left(0\right)=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^{\prime }({\langle k\rangle }_{\mathrm{w}}-1)=0$$

(48)

where

$${\langle k\rangle }_{\mathrm{w}}\mathsf{\equiv }\sum _{k\ge 1}k{p}_k$$

(49)

is the average multiplicity of the junctions. The average ${\sum }_k(k-1)p_k={\langle k\rangle }_{\mathrm{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\left(t_{\mathrm{g}}\right)=\frac{1}{f^{\prime }{k}^{\prime }}\mathsf{\equiv }{p}_{\mathrm{g}}$$

(50)

Substituting this value into (47), we find

$$\frac{1}{t_{\mathrm{g}}^{\left(\mathrm{IR}\right)}}=\alpha k{\psi }^{k^{\prime }}\frac{k^{\prime }{(1-p_{\mathrm{g}})}^{k^{\prime }}}{1-{(1-p_{\mathrm{g}})}^{k^{\prime }}}\mathsf{\equiv }A\left(T\right){\varphi }^{k^{\prime }}$$

(51)

Therefore, it turns out that temperature and concentration is separable for irreversible gelation. The temperature factor is

$$A\left(T\right)=\alpha \left(T\right){\left(\frac{f}{n}\right)}^{k^{\prime }}\frac{k{k}^{\prime }}{{[f^{\prime }{k}^{\prime }/(f^{\prime }{k}^{\prime }-1)]}^{k^{\prime }}-1}$$

(52)

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

$$\tilde{p}\left(\tau \right)=\overline{p}(1-e^{-\gamma \tau })$$

(53)

where

$$\gamma \mathsf{\equiv }-g^{\prime }\left(\overline{p}\right)=k{(1-\overline{p})}^{k^{\prime }}+{\beta }^{\prime }$$

(54)

and hence the relaxation time is

$$\frac{1}{t_{\mathrm{R}}}=\alpha k{\psi }^{k^{\prime }}\gamma =\beta \left(1+k^2{K}_k{\overline{z}}^{k^{\prime }}\right)=\beta \left(1+k\frac{{\overline{p}}_k\left(\overline{z}\right)}{{\overline{p}}_1\left(\overline{z}\right)}\right)$$

(55)

by using equilibrium values of the reactivities. They can be written as functions of the solution $\overline{z}$ of the conservation equation (6) with

$$u\left(\overline{z}\right)\mathsf{\equiv }1+k{K}_k{\overline{z}}^{k^{\prime }}$$

(56)

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\left(\overline{z}\right)=ln\left(\frac{f^{\prime }{k}^{\prime }{\overline{p}}_k\left(\overline{z}\right)}{f^{\prime }{k}^{\prime }{\overline{p}}_k\left(\overline{z}\right)-1}\right)$$

(57)

The equilibrium gel point is given by the condition $1-f^{\prime }{k}^{\prime }{\overline{p}}_k\left(\overline{z}\right)=0$ and hence

$${\left(K_k\left(T\right){\psi }^{k^{\prime }}\right)}_{\mathrm{g}}=\frac{{\left(f^{\prime }{k}^{\prime }\right)}^{k^{\prime }}}{k{(f^{\prime }{k}^{\prime }-1)}^k}$$

(58)

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\left(T\right)=\lambda {\left(T\right)}^{k^{\prime }}$$

(59)

where $\lambda \left(T\right)$ is a binding constant per one functional group. The conservation law becomes

$$x\mathsf{\equiv }\lambda \psi =\overline{z}(1+k{\overline{z}}^{k^{\prime }})$$

(60)

where $\overline{z}\mathsf{\equiv }x(1-{\overline{p}}_k)$ is the scaled concentration of the free groups. Solving this equation for $\overline{z}$, and substituting the solution to

$${\tau }_{\mathrm{g}}\mathsf{\equiv }\alpha k{\psi }^{k^{\prime }}{t}_{\mathrm{g}}=\frac{k{x}^{k^{\prime }}Q\left(\overline{z}\right)}{1+k{\overline{p}}_k\left(\overline{z}\right)/{\overline{p}}_1\left(\overline{z}\right)}$$

(61)

we find the gelation time as functions of x. Let us plot it in our standard form

$$\Gamma \left(x\right)\mathsf{\equiv }\frac{1}{\beta t_{\mathrm{g}}}=(1+k^2{\overline{z}}^{k^{\prime }})ln\left(\frac{f^{\prime }{k}^{\prime }k{\overline{z}}^k}{f^{\prime }{k}^{\prime }k{\overline{z}}_x^k-1}\right)$$

(62)

where $\overline{z}=\overline{z}\left(x\right)$ is the solution of the conservation law (60).

Figure 4 (b) plots $1/\beta t_{\mathrm{g}}$, $1/\beta t_{\mathrm{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_{\mathrm{g}}=f^{\prime }{k}^{\prime }/k^{1/k^{\prime }}{(f^{\prime }{k}^{\prime }-1)}^{k/k^{\prime }}$ 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_{\mathrm{g}}$ in powers of the reverse rate constant $\beta$. We find

$$t_{\mathrm{g}}=t_{\mathrm{g}}^{\left(\mathrm{IR}\right)}\left\{1+\frac{R_{f,k}}{x^{k^{\prime }}}+O\left(\frac{1}{x^{2k^{\prime }}}\right)\right\}$$

(63)

with

$$R_{f,k}=\frac{{(1-p_{\mathrm{g}})}^{2k-1}+(2k-1)p_{\mathrm{g}}-1}{2k(2k-1)[1-{(1-p_{\mathrm{g}})}^{k^{\prime }}]{(1-p_{\mathrm{g}})}^k}$$

(64)

where $t_{\mathrm{g}}^{\left(\mathrm{IR}\right)}$ 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:

$${\mathrm{A}}_1+{\mathrm{A}}_{k-1}\leftrightarrows {\mathrm{A}}_k(k=2,3,\cdots s)$$

(65)

$$\frac{d{n}_k\left(t\right)}{dt}=J_k-J_{k+1}$$

(66a)

$$=-\{{\beta }_k+{\alpha }_k{n}_1\left(t\right)\}{n}_k\left(t\right)+{\beta }_{k+1}{n}_{k+1}\left(t\right)+{\alpha }_{k-1}{n}_1\left(t\right)n_{k-1}\left(t\right)$$

(66b)

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

$$\frac{d{n}_k\left(t\right)}{dt}=J_k-J_{k+1}$$

(66a)

$$=-\{{\beta }_k+{\alpha }_k{n}_1\left(t\right)\}{n}_k\left(t\right)+{\beta }_{k+1}{n}_{k+1}\left(t\right)+{\alpha }_{k-1}{n}_1\left(t\right)n_{k-1}\left(t\right)$$

(66b)

for $k\ge 2$, and

$$\frac{d{n}_1\left(t\right)}{dt}=-2J_2-\sum _{k=3}^s{J}_k$$

(67a)

$$=-2{\alpha }_1{n}_1{\left(t\right)}^2+2{\beta }_2{n}_2\left(t\right)+\sum _{k=3}^s\{{\beta }_k{n}_k\left(t\right)-{\alpha }_{k-1}{n}_1\left(t\right)n_{k-1}\left(t\right)\}$$

(67b)

for $k=1$. Here,

$$J_k\left(t\right)\equiv {\alpha }_{k-1}{n}_1\left(t\right)n_{k-1}\left(t\right)-{\beta }_k{n}_k\left(t\right)$$

(68)

for $k\ge 2$ is the flux between $k-1$ and k state ($J_1\equiv 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 ${\alpha }_k$ and ${\beta }_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\left(t\right)$ in order to find the equilibrium solution easily. Using equation (3), we find

$$\frac{d{p}_k\left(t\right)}{dt}=-\{{\beta }_k+{\alpha }_k\psi p_1\left(t\right)\}{p}_k\left(t\right)+\frac{k}{k+1}{\beta }_{k+1}{p}_{k+1}\left(t\right)+\frac{k}{k-1}{\alpha }_{k-1}\psi p_1\left(t\right)p_{k-1}\left(t\right)$$

(69)

for $k\ge 2$, and

$$\frac{d{p}_1\left(t\right)}{dt}=-2{\alpha }_1\psi p_1{\left(t\right)}^2+{\beta }_2{p}_2\left(t\right)+\sum _{k=3}^s\left\{\frac{{\beta }_k}{k}{p}_k\left(t\right)-\frac{{\alpha }_{k-1}}{k-1}\psi p_1\left(t\right)p_{k-1}\left(t\right)\right\}$$

(70)

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

$$({\beta }_k{K}_k-{\alpha }_{k-1}{K}_{k-1})-({\alpha }_k{K}_k-{\beta }_{k+1}{K}_{k+1})\overline{z}=0$$

(71)

must be fulfilled. We therefore propose the detailed balance condition such that the equlibrium constants satisfy the relation

$$K_k=\frac{{\alpha }_{k-1}}{{\beta }_k}{K}_{k-1}$$

(72)

By repeated use of this relation, we have a well known relation

$$K_k=\frac{{\alpha }_{k-1}{\alpha }_{k-2}}{{\beta }_k{\beta }_{k-1}}{K}_{k-2}=\cdots ={\lambda }_k{\lambda }_{k-1}\cdots {\lambda }_2$$

(73)

between the equilibrium constants $K_k$ and the stepwise association constants

$${\lambda }_k\equiv \frac{{\alpha }_{k-1}}{{\beta }_k}$$

(74)

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

$$\frac{d{p}_1}{dt}=-2{\alpha }_1{p_1}^2+({\beta }_2-\frac{1}{2}{\alpha }_2\psi p_1)p_2+\frac{{\beta }_3}{3}{p}_3$$

(75a)

$$\frac{d{p}_2}{dt}=2{\alpha }_1\psi {p_1}^2-({\beta }_2+{\alpha }_2\psi p_1)p_2+\frac{2{\beta }_3}{3}{p}_3$$

(75b)

$$\frac{d{p}_3}{dt}=\frac{3{\alpha }_2}{2}\psi p_1{p}_2-{\beta }_3{p}_3$$

(75c)

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 $\psi {\overline{p}}_1=\overline{z},\psi {\overline{p}}_2=2K_2{\overline{z}}^2,\psi {\overline{p}}_3=3K_3{\overline{z}}^3$. Detailed balance conditions lead to the relation $K_2={\alpha }_1/{\beta }_2\equiv {\lambda }_2$, and $K_3={\alpha }_1{\alpha }_2/{\beta }_2{\beta }_3\equiv {\lambda }_2{\lambda }_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 ${\xi }_1,{\xi }_2$ from the equilibrium state defined by $p_k\left(t\right)\equiv {\overline{p}}_k(1+{\xi }_k\left(t\right))$. After a lengthy calculation, we find

$$\frac{d}{dt}\left[\begin{array}{c}{\xi }_1\\ {\xi }_2\end{array}\right]=-\left[\begin{array}{cc}{a}_{1,1}& a_{1,2}\\ a_{2,1}& a_{2,2}\end{array}\right]\left[\begin{array}{c}{\xi }_1\\ {\xi }_2\end{array}\right]+\left[\begin{array}{c}{f}_1({\xi }_1,{\xi }_2)\\ f_2({\xi }_1,{\xi }_2)\end{array}\right]$$

(76)

where

$$a_{1,1}\equiv \frac{{\beta }_3}{3}+4{\alpha }_1\overline{z}+{\alpha }_2{K}_2{\overline{z}}^2$$

(77a)

$$a_{1,2}\equiv K_2\overline{z}\left(\frac{2{\beta }_3}{3}-2{\beta }_2+{\alpha }_2\overline{z}\right)$$

(77b)

$$a_{2,1}\equiv \frac{1}{K_2\overline{z}}\left(\frac{{\beta }_3}{3}-2{\alpha }_1\overline{z}+{\alpha }_2{K}_2{\overline{z}}^2\right)$$

(77c)

$$a_{2,2}\equiv {\beta }_2+{\beta }_3\left(\frac{2}{3}+{\lambda }_3\overline{z}\right)$$

(77d)

for the linear terms, and

$$f_1({\xi }_1,{\xi }_2)\equiv -(2{\alpha }_1{\xi }_1+{\alpha }_2{K}_2\overline{z}{\xi }_2)\overline{z}{\xi }_1$$

(78a)

$$f_2({\xi }_1,{\xi }_2)\equiv ({\alpha }_1{\xi }_1-{\alpha }_2{K}_2\overline{z}{\xi }_2){\xi }_1/K_2$$

(78b)

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 $\gamma$ of the matrix $\widehat{A}\equiv \left(a_{i,j}\right)$. The equation to find them is

$${\gamma }^2-A\left(\overline{z}\right)\gamma +B\left(\overline{z}\right)=0$$

(79)

where

$$A\left(\overline{z}\right)\equiv \mathrm{Tr}\left(\widehat{A}\right)={\beta }_2+{\beta }_3+(4{\alpha }_1+{\alpha }_2)\overline{z}+{\alpha }_2{K}_2{\overline{z}}^2$$

(80a)

$$B\left(\overline{z}\right)\equiv \mathrm{Det}\left(\widehat{A}\right)={\beta }_2{\beta }_3+4{\alpha }_1{\beta }_3\overline{z}+9{\alpha }_1{\alpha }_2{\overline{z}}^2$$

(80b)

Note that $B\left(\overline{z}\right)$ can be written as

$$B\left(\overline{z}\right)={\beta }_2{\beta }_3{\left(\overline{z}u\left(\overline{z}\right)\right)}^{\prime }={\beta }_2{\beta }_3\left(1+2\frac{{\overline{p}}_2}{{\overline{p}}_1}+3\frac{{\overline{p}}_3}{{\overline{p}}_1}\right)$$

(81)

so that it is regarded as an extension of KHDS decay rate (55) to the three-state model. (Prime ${}^{\prime }$ indicates differentiation.) The linear relaxation has two modes

$${\gamma }^{(\pm )}\left(\overline{z}\right)\equiv \frac{A\left(\overline{z}\right)}{2}\left\{1\pm \sqrt{1-\frac{4B\left(\overline{z}\right)}{A{\left(\overline{z}\right)}^2}}\right\}$$

(82)

The gel point is given by the condition

$$1-f^{\prime }\{p_2\left(t\right)+2p_3\left(t\right)\}=0$$

(83)

Using the deviation ${\xi }_k\left(t\right)$, this condition is transformed to

$${\xi }_1\left(t\right)+K_2\overline{z}{\xi }_2\left(t\right)=K_2+3K_3\overline{z}-1/2f^{\prime }\overline{z}$$

(84)

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 ${\beta }_2<<{\beta }_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

$$\frac{1}{t_{\mathrm{R}}}\simeq \frac{{\beta }_2{\beta }_3{\left(\overline{z}u\left(\overline{z}\right)\right)}^{\prime }}{{\beta }_2{b}_2\left(\overline{z}\right)+{\beta }_3{b}_3\left(\overline{z}\right)}=\overline{\beta }\left(\overline{z}\right)\left(1+2\frac{{\overline{p}}_2}{{\overline{p}}_1}+3\frac{{\overline{p}}_3}{{\overline{p}}_1}\right)$$

(85)

where $\overline{\beta }$ is the average of $\beta$ defined by

$$\frac{1}{\overline{\beta }\left(\overline{z}\right)}\equiv \frac{b_2\left(\overline{z}\right)}{{\beta }_2}+\frac{b_3\left(\overline{z}\right)}{{\beta }_3}$$

(86)

The relaxation time apparently takes an extended KHDS form. Here, we have eliminated ${\alpha }_k$ in favor of ${\beta }_k$ and ${\lambda }_k$ from $A\left(\overline{z}\right)$, and write

$$A\left(\overline{z}\right)={\beta }_2{b}_2\left(\overline{z}\right)+{\beta }_3{b}_3\left(\overline{z}\right)$$

(87)

by using

$$b_2\left(\overline{z}\right)\equiv 1+4{\lambda }_2\overline{z}$$

(88a)

$$b_3\left(\overline{z}\right)\equiv 1+{\lambda }_3\overline{z}+{\lambda }_2{\lambda }_3\overline{z}$$

(88b)

If the first association (pairing) is the rate-determining step, the relaxation time is approximately

$$\frac{1}{t_{\mathrm{R}}}\simeq \frac{{\beta }_2}{b_3\left(\overline{z}\right)}{\left(\overline{z}u\left(\overline{z}\right)\right)}^{\prime }$$

(89)

and goes back to the pairwise association in the limit of ${\lambda }_3\to 0$. If the second associaition (triple junction) is the rate-determining step, the relaxation time is approximately given by

$$\frac{1}{t_{\mathrm{R}}}\simeq \frac{{\beta }_3}{b_2\left(\overline{z}\right)}{\left(\overline{z}u\left(\overline{z}\right)\right)}^{\prime }$$

(90)

From these results, we may expect the extended KHDS form for the multiplicity more than three

$$\frac{1}{t_{\mathrm{R}}}\simeq \beta \left(1+\frac{1}{{\overline{p}}_1}\sum _{k\ge 2}k{\overline{p}}_k\right)=\beta \frac{{\langle k\rangle }_{\mathrm{w}}}{{\overline{p}}_1}$$

(91)

with some average $\beta$, 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{\xi }_2/dt\simeq 0$, the kinetic equation reduces to

$$\frac{d{\xi }_1}{dt}\simeq -\frac{B\left(\overline{z}\right)}{a_{2,2}}{\xi }_1$$

(92)

The rexation time is then given by

$$\gamma \left(\overline{z}\right)=\frac{{\beta }_2{\beta }_3}{a_{2,2}}{\left(\overline{z}u\left(\overline{z}\right)\right)}^{\prime }$$

(93)

The front factor can be regarded as an average backward rate constant $\overline{\beta }$ defined by

$$\frac{1}{\overline{\beta }\left(\overline{z}\right)}\equiv \frac{1}{{\beta }_3}+\frac{2/3+{\lambda }_3\overline{z}}{{\beta }_2}$$

(94)

Hence we find again the same form

$$\frac{1}{t_{\mathrm{R}}}=\overline{\beta }\left(\overline{z}\right)\left(1+2\frac{{\overline{p}}_2}{{\overline{p}}_1}+3\frac{{\overline{p}}_3}{{\overline{p}}_1}\right)$$

(95)

for the relaxation time, with slightly different average of $\beta$.

3.1.3. Slow-Mode Approximation

The general solution for the linear equation is given by

$$\left[\begin{array}{c}{\xi }_1\\ {\xi }_2\end{array}\right]=C_1{\mathbf{e}}^1{e}^{-{\gamma }^{(-)}t}+C_2{\mathbf{e}}^2{e}^{-{\gamma }^{(+)}t}$$

(96)

where ${\mathbf{e}}^1$ is the normalized eigen-vector for the eigen-value ${\gamma }^{(-)}$ (slow mode) and ${\mathbf{e}}^2$ for the fast mode. The constants $C_1,C_2$ can be fixed by the initial conditions $p_1\left(0\right)=1,p_2\left(0\right)=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_{\mathrm{R}}={\gamma }^{(-)}$, and the gelation time is $t_{\mathrm{g}}=t_{\mathrm{R}}Q\left(\overline{z}\right)$ with a thermodynamic factor

$$Q\left(\overline{z}\right)=ln\left\{\frac{2f^{\prime }\overline{z}{D}_1\left(\overline{z}\right)}{f^{\prime }({\langle k\rangle }_{\mathrm{w}}-1)-1}\right\}$$

(97)

from (84), where $D_1\left(\overline{z}\right)\equiv C_1(e_1^1+K_2\overline{z}{e}_2^1)$ is a function of $\overline{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 $\mu$ times larger than the first one:

$${\lambda }_2\equiv \lambda ,{\lambda }_3=\mu \lambda$$

(98)

Then the forward rate constants are given by

$${\alpha }_1={\beta }_2\lambda ,{\alpha }_2=\mu \lambda {\beta }_3$$

(99)

The conservation law takes the form

$$x\equiv \lambda \psi =\overline{z}(1+\overline{z}+\mu {\overline{z}}^2)$$

(100)

where $\lambda \overline{z}$ is written as $\overline{z}$ for simplicity. Now we have

$$A\left(\overline{z}\right)={\beta }_2(1+4\overline{z})+{\beta }_3(1+\mu \overline{z}+\mu {\overline{z}}^2)$$

(101a)

$$B\left(\overline{z}\right)={\beta }_2{\beta }_3(1+4\overline{z}+9\mu {\overline{z}}^2)$$

(101b)

In particular, for a special case of ${\beta }_2={\beta }_3\equiv \beta$ the relaxation time is given by

$$\frac{1}{\beta t_{\mathrm{R}}}=\frac{1}{2}A\left(\overline{z}\right)\left\{1-\sqrt{1-\frac{4B\left(\overline{z}\right)}{A{\left(\overline{z}\right)}^2}}\right\}$$

(102)

with

$$A\left(\overline{z}\right)=2+(\mu +4)\overline{z}+\mu {\overline{z}}^2$$

(103a)

$$B\left(\overline{z}\right)=1+4\overline{z}+9\mu {\overline{z}}^2$$

(103b)

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 $\mu$, the surface curve of the triple junction changes from A to D. For the limit of $\mu =0$ (curve D), the model reduces to the pairwise one. Figure 6 (b) plots The reciprocal gelation time $1/\beta t_{\mathrm{g}}$, the reciprocal relaxation time $1/\beta t_{\mathrm{R}}$, and the thermodynamic factor $Q\left(x\right)$ plotted against a scaled concentration $x\equiv \lambda \left(T\right)\psi$ for three ratios $\mu =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

$$\frac{d{\xi }_k}{dt}=-({\beta }_k+{\alpha }_k\overline{z})+{\beta }_k{\xi }_{k-1}+{\alpha }_k\overline{z}{\xi }_{k+1}+({\beta }_k-{\alpha }_k\overline{z}){\xi }_1(k\ge 2)$$

(104a)

$$\frac{d{\xi }_1}{dt}=-2{\alpha }_1\overline{z}(2{\xi }_1-{\xi }_2)-\sum _{k\ge 3}{\beta }_k{K}_k{\overline{z}}^k({\xi }_{k-1}-{\xi }_k+{\xi }_1)$$

(104b)

The gel-point condition (48) turns into

$$f^{\prime }\sum _{k\ge 2}(k-1){\overline{p}}_k{\xi }_k\left(t_{\mathrm{g}}\right)=1-f^{\prime }\sum _{k\ge 2}(k-1){\overline{p}}_k$$

(105)

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\le k\le k_2$ in the most time during the approach to equilibrium (assumption of quasi-stationary flow). We then have

$$J_k\simeq \left\{\begin{array}{cc}0\hfill & (1\le k\le k_1)\hfill \\ \mathrm{independent}\mathrm{of}k\equiv J\hfill & (k_1+1\le k\le k_2)\hfill \\ 0\hfill & (k_2+1\le k)\hfill \end{array}\right.$$

(106)

In Appendix C, we show that the flux J is proportional to the probability of free groups as

$$RJ\simeq \frac{m_0+M_1}{M_0}{\xi }_1$$

(107)

where $m_j$ and $M_j$ are the j-th moments for small aggregates and large micelles defined by

$$m_j\equiv \sum _{k=1}^{k_1}{k}^j{\overline{p}}_k$$

(108a)

$$M_j\equiv \sum _{k=k_2+1}^{\infty }{k}^j{\overline{p}}_k$$

(108b)

for $j=0,1,2\cdots$. Here the resistance R in the heat flow analogy is given by

$$R\equiv \sum _{k=k_1+1}^{k_2}\frac{1}{{\beta }_k{K}_k{\overline{z}}^k}=\frac{1}{\psi }\sum _{k=k_1+1}^{k_2}\frac{k}{{\beta }_k{\overline{p}}_k}$$

(109)

We also find that the probability of large micelles for $k\ge k_2+1$ is approximately given by

$${\xi }_k\simeq k{\xi }_1-RJ$$

(110)

The kinetic equation for $k=1$ then takes the form

$$\frac{d{\xi }_1}{dt}=-2J_2-\sum _{k=3}^{\infty }{J}_k\simeq -\frac{m_1+M_1}{R{M}_0}(k_2-k_1){\xi }_1$$

(111)

Hence, the relaxation time is

$$\frac{1}{t_{\mathrm{R}}}=(k_2-k_1)\frac{m_1}{R{M}_0}\left(1+\frac{M_1}{m_1}\right)$$

(112)

which is similar to AW formula. The solution of (111) is given by

$${\xi }_1\left(t\right)=C{e}^{-t/t_{\mathrm{R}}}$$

(113)

with $C\equiv (1-{\overline{p}}_1)/{\overline{p}}_1$.

Substituting (110) into the gel-point condition (105), we have

$$f^{\prime }\left\{\sum _{k=2}^{k_1}(k-1)k{\overline{p}}_k+\sum _{k\ge k_2+1}(k-1)k{\overline{p}}_k\right\}{\xi }_1\left(t_{\mathrm{g}}\right)-f^{\prime }RJ\sum _{k\ge k_2+1}{\overline{p}}_k=1-f^{\prime }\sum _{k\ge 2}(k-1){\overline{p}}_k$$

(114)

and hence we finally find the form

$$t_{\mathrm{g}}=t_{\mathrm{R}}Q\left(\overline{z}\right)$$

(115)

for the gelation time with the thermodynamic factor

$$Q\left(\overline{z}\right)=ln\left\{\frac{f^{\prime }CM}{f^{\prime }{\sum }_{k\ge 2}(k-1){\overline{p}}_k-1}\right\}$$

(116)

where

$$M\equiv \frac{M_1}{M_0}(m_1+M_1)-m_2-M_2$$

(117)

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,\cdots 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 ${\xi }_k\left(t\right)$ at time t. Since the moments are reduced to $m_1=m_2={\overline{p}}_1$ and $M_0={\overline{p}}_s,M_1=s{\overline{p}}_s$, we have

$$\frac{1}{t_{\mathrm{R}}}=(s-2)\frac{{\overline{p}}_1}{R{\overline{p}}_s}\left(1+s\frac{{\overline{p}}_s}{{\overline{p}}_1}\right)$$

(118)

This is again KHDS form with an effective rate constant

$${\beta }_{\mathrm{eff}}=(s-2)\frac{{\overline{p}}_1}{R{\overline{p}}_s}$$

(119)

written in terms of the resistance constant

$$R\equiv \sum _{k=2}^{s-1}\frac{k}{{\beta }_k{\overline{p}}_k}=\sum _{k=2}^{s-1}\frac{1}{{\beta }_k{K}_k{\overline{z}}^k}$$

(120)

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 $\beta$ is replaced by an effective one ${\beta }_{\mathrm{eff}}$. Relaxation is much slower, in particular, for high concentration region.

Figure 8 shows the gelation time $t_{\mathrm{g}}$ , the relaxation time $t_{\mathrm{R}}$, and the thermodynamic factor Q plotted against the scaled concentration $x\equiv \lambda \left(T\right)\varphi$ for the functionality $f=4$ with the varied upper bound $s=3,5,7,9,11$ of the junction multiplicity. For simplicity, uniform ${\beta }_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_{\mathrm{R}}\sim 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\le t_{\mathrm{g}}$ before the solution reaches the gel point. In this region, the conservation law (6) of the functional groups, equilibrium reactivities ${\overline{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_{\mathrm{g}}\le 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_{\mathrm{g}},1/t_{\mathrm{R}}$ and Q against the scaled concentration variable $x\equiv \lambda \left(T\right)\psi$.

From the theoretical modeling, the following conclusions can be drawn:

(1)

The gelation time $t_{\mathrm{g}}$, the relaxation time $t_{\mathrm{R}}$, and the thermodynamic factor Q are all functions of a single variable $x\equiv \lambda \left(T\right)\psi$ (scaled concentration), where $\lambda \left(T\right)$ 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_{\mathrm{g}}\left(x\right)=t_{\mathrm{R}}\left(x\right)Q\left(x\right)$. The thermodynamic factor $Q\left(x\right)$ is logarithmically singular at the equilibrium gel point $x=x_{\mathrm{g}}$, while the relaxation time $t_{\mathrm{R}}\left(x\right)$ 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_{\mathrm{g}}\left(x\right)$ of reversible cross-linking approaches that of irreversible one in the asymptotic region of large x (either high concentration $\varphi$ or high values of the association constant $\lambda \left(T\right)$). The power index of $1/t_{\mathrm{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.