Fumihiko Tanaka
Gelation Time of Network-Forming Polymer Solutions with Reversible Cross-Link Junctions of Variable Multiplicity
thermoreversible gelation,gelation time,relaxation time,multiple cross-links,stepwise association
The gelation time 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 on the basis of the kinetic equation for the stepwise cross-linking reaction as a function of the concentration, temperature, functionality f of the molecules, and multiplicity k of the cross-link junctions. It is shown that quite generally can be decomposed into the product of the relaxation time and a thermodynamic factor Q. They are functions of a single scaled concentration , where is the association constant and is the concentration. Therefore, the superposition principle holds with as a shift factor of the concentration. Additionally, they all depend on the rate constants of the cross-link reaction, and hence it is possible to estimate these microscopic parameters from macroscopic measurements of . The thermodynamic factor Q is shown to depend on the quench depth. It generates a singularity of logarithmic divergence as the temperature (concentration) approaches the equilibrium gel point, while the relaxation time changes continuously across it. Gelation time obeys a power law in the high concentration region, whose power index n is related to the multiplicity of the cross-links. 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 to find the rate-controlling steps in order for the minimization of the gelation time to be easier in the gel processing. For a micellar cross-linking covering a wide range of the multiplicity, as seen in hydrophobically-modified water-soluble polymers, is shown to obey a formula similar to the Aniansson–Wall law.
Gelation Time of Network-Forming Polymer Solutions with Reversible Cross-Link Junctions of Variable Multiplicity

The gelation time 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 on the basis of the kinetic equation for the stepwise cross-linking reaction as a function of the concentration, temperature, functionality f of the molecules, and multiplicity k of the cross-link junctions. It is shown that quite generally can be decomposed into the product of the relaxation time and a thermodynamic factor Q. They are functions of a single scaled concentration , where is the association constant and is the concentration. Therefore, the superposition principle holds with as a shift factor of the concentration. Additionally, they all depend on the rate constants of the cross-link reaction, and hence it is possible to estimate these microscopic parameters from macroscopic measurements of . The thermodynamic factor Q is shown to depend on the quench depth. It generates a singularity of logarithmic divergence as the temperature (concentration) approaches the equilibrium gel point, while the relaxation time changes continuously across it. Gelation time obeys a power law in the high concentration region, whose power index n is related to the multiplicity of the cross-links. 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 to find the rate-controlling steps in order for the minimization of the gelation time to be easier in the gel processing. For a micellar cross-linking covering a wide range of the multiplicity, as seen in hydrophobically-modified water-soluble polymers, is shown to obey a formula similar to the Aniansson–Wall law.

Thermoreversible gelation in solutions of polymers, as well as of low molecular weight molecules, has been attracting researchers’ 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. Some original studies, reviews and conceptual works have appeared with relation to responsive gels [8,9,10,11], hydrogels for biomedical applications [6,7,12], and hydrogen-bonding supramolecular gelators [13,14,15]. 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 is defined by the time necessary for the network-forming solutions of functional molecules to reach their gel points after a cross-linking reaction is started. It is related to the relaxation time 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, the food industry, the adhesion industry, etc. Gelation time is also very important from a 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 on the microscopic structure and formation process of the network junctions, in particular the rate-controlling steps in cross-linking. The gelation time critically depends on the thermodynamic conditions under which the 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 and , and find them as functions of the functionality of the molecules, multiplicity of the junctions, and 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 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. [16,17] reported the data on the gelation time of poly(vinyl alcohol) (PVA) in a mixed solvent of water and dimethylsulphoxide. They measured the time by a ball-drop method in the temperature quenching experiments for various polymer concentrations. The results were summarized in the form where is the concentration of PVA measured from the overlap concentration. The critical concentration at which becomes infinite was found to be independent of the temperature. (The log–log plot of this equation is referred to as Oakenfull plot in p.80 of the reference [2].) The exponent n was found to be 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 where B is a temperature-independent constant. Investigation was done by Mal et al. [18] on the solutions of crystalline polymer poly(vinylidene fluoride) (PVF) in two different solvents. The 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 0.45–0.60, much smaller than the value for PVA. Hong et al. [19] studied the same polymer PVF 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 to as the critical gel point is approached. Based on these data, a new model of gel formation in crystalline polymers was proposed. A similar study was done by Tobitani et al. [20,21] 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 took an activation type. They also found the exponent of the concentration to be 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 [22,23,24]. 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 [25], silica alkoxides [26,27], etc. Investigation on the gelation time was reported in the literature for the binary cross-linking in the mixture of functional molecules [28,29,30,31,32]. The main problem of binary gels is to find the optimal mixing ratio for minimizing the gelation time, so that plot of against the mixing ratio for different quenching temperatures is fundamental. From the 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 [33,34] 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. [35,36,37,38] extended Stockmayer’s model on the basis of a 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 [36]. The scaling and universality of gel formation was also discussed [38]. Dongen et al. [39] considered reversibility of the cross-linking reactions, and studied retardation effects due to the existence of reverse reaction. However, unfortunately, the effects of temperature and concentration were not explicitly considered in these studies. In spite of this diversity in experimental and theoretical studies, there has been no attempt to construct a unified picture of the gelation time of polymer solutions by starting from the fundamental kinetics of the cross-linking reaction. In this paper, we focus on the reversibility and variable cross-link multiplicity of a stepwise rate equation for the cross-linking reaction. Solving the rate equation for the probability of cross-linking with multiplicity k, we derive the gelation time as a function of the polymer concentration and temperature. In contrast to (1), it includes a temperature–concentration cross-term. The results are expressed in the common form in terms of the relaxation time 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. The gelation time for the reversible binary cross-linking will be reported in our forthcoming paper.

The model solution we consider is a polymer solution in which the number N of reactive (associative) molecules (denoted as R{A}) with degree of polymerization n are dissolved in the number of solvent molecules (S). Molecules can be any type, such as high molecular weight linear polymers, star polymers, 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,22,23]. We are based on the lattice-theoretical picture of polymer solutions [40], 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 in the unit of the cell volume. We assume incompressibility of the solution, so that we have for the total volume. The volume fraction of each component is then given by for the reactive molecule, for the solvent. In terms of the functional groups, the number concentration of A-groups on the reactive molecules is . The fundamental picture of our research problem is summarized in Figure 1. We study the time development of the solution after its temperature is changed suddenly from the initial one to the final one (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 1a). The temperature quench depth is defined by , while the concentration quench depth is , where and 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 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 for an arbitrarily chosen functional group to belong to the junction of multiplicity k is related to the number of junctions by the equation After a long time, the solution reaches its equilibrium state with equilibrium reactivity . Figure 1b schematically shows the time development of , its linearized approximation , irreversible counterpart , reactivity of the sol part 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 [24], we studied in detail thermoreversible gelation and phase separation in solutions capable of unary (self) cross-linking. We started from the equilibrium condition where is the equilibrium constant, is the probability for a functional group to remain free from reaction, and is its equilibrium value. This equilibrium condition leads to for the reactivity in terms of the number of free groups . From the normalization condition of , we find the conservation law where

Let us first consider the simplest case of pairwise cross-linking reaction for which k takes only two values: (free), and (bound). Cross-linking by a covalent bond is included as the irreversible limit of such a 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 . Probability p is the conventional reactivity of the cross-linking reaction. The time development kinetic equation for the number of cross-link junctions can be written as and hence we have where and are the rate constant of the forward reaction and backward reaction. (Throughout this paper, we avoid conventional symbols and to prevent confusion with multiplicity k.) By using the scaled time , we have a simple second order equation where Here, a new constant is the equilibrium constant written in terms of the rate constants. The solution of this equation with the initial condition is given by where and The equilibrium reactivity is given by . Therefore, has the meaning of relaxation time of the reversible reaction, i.e., the time necessary for the system to reach its chemical equilibrium. A detailed derivation of (14) is given in Appendix A. The result is drawn in Figure 2a for the functionality . In Figure 2b, the weight- and number-average molecular weight of the three-dimensional cross-linked polymers ( and in the unit of the molecular weight of the primary molecules) are plotted against time. They are explicitly given by [33] 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. However, there are other possibilities [35,36,37,38]. In this paper, we focus on the process of approaching the gel point, and prevent discussion on the post-gel regime. In Figure 2b, 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) . The number-average remains finite at the gel point. Because the relaxation time is defined by the speed to reach the equilibrium state, it is given by . In terms of the bare time, we have In the limit of irreversible reaction where with finite , is small, and , so that we go back to the Stockmayer’s result [33,34] (see in Figure 1b) 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 [33,34,40,41,42,43] give the gel point for the reactivity. Substituting this value for in (14), we find for the gelation time where is a function of a combined variable of the temperature and concentration. Since this factor drastically depends upon the quench depth , 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 , so that goes logarithmically to infinity as x approaches the critical value at the gel point. (In what follows, we use abbreviated notations for simplicity.) Near the gel point, we have . A lengthy calculation leads to the form 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 can be expressed in terms of the equilibrium reactivity. The solution of the conservation law (6) for a pairwise cross-linking is given by We therefore have , and hence 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 and take limit while is kept finite. Because x takes a large value, we expand the right hand side of this equation in powers of . After a quite a lengthy calculation, we find where the front factor is the gelation time of an irreversible reaction (covalent bonding). It can directly be obtained from Stockmayer’s form (19) by fixing . The first correction due to reversible reaction is . Its coefficient is found to be Let us refer to it as retardation coefficient. In the next section, we shall derive the retardation coefficient for the cross-linking with arbitrary multiplicity k. To compare with experimental data, we plot in Figure 3 the reciprocal of the gelation time as a function of the combined variable x for varied functionality f. Near the sol–gel transition point , it goes logarithmically to zero as 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 [16,17] 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 a much lower concentration than the overlap concentration. For sufficiently high x, it approaches the irreversible limit In this limit, the gelation time is separable into a product of the temperature factor and the concentration factor . 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. The only assumption is that the gel point is assumed to be given by the conventional tree approximation. We now examine the validity of the 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 and consider only linear terms of the deviation in (11). We have where is the derivative of . Substituting the equilibrium value of p leads to Hence, the relaxation time in the linear approximation agrees with the exact one. If we assume , the reactivity in the linear approximation takes the form (see Figure 1b) The gel-point condition then leads to (21), but with of slightly different form 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.

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 before we move onto a more complex case of stepwise association. A functional group is either free () or reacted (). Let us refer to it as a fixed-multiplicity model (see Figure 4). We have a rate equation for the fixed-multiplicity reversible reaction Or, equivalently (). Hence, we have where . Time is scaled as and The coefficient of the reverse reaction is changed to where is the equilibrium constant. In the study of heat-setting -Lactoglobulin, Tobitani and Ross-Murphy [20,21] solved this equation in the irreversible limit of . The rate equation can easily be integrated into this limit. We find The relaxation is not exponential, but obeys a power law. Let us assume that all functional groups are free in the initial state (). In our previous study of gelation with cross-links of variable multiplicity [23,24], we showed that the gel point is in general given by the condition where is the average multiplicity of the junctions. The average was referred to as a branching number [23,24]. In the case of the fixed-multiplicity model, the gel point is found by the condition Substituting this value into (47), we find Therefore, it turns out that temperature and concentration are separable for irreversible gelation. The temperature factor is It goes back to (28) for the irreversible pairwise cross-linking of . 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 where and hence the relaxation time is by using equilibrium values of the reactivities. They can be written as functions of the solution of the conservation Equation (6) with The relaxation time similar form to (55) was first proposed by Kresheck et al. [44] (referred to as KHDS), and later by Muller [45] 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 The equilibrium gel point is given by the condition and hence 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 [23] where is a binding constant per one functional group. The conservation law becomes where is the scaled concentration of the free groups. Solving this equation for , and substituting the solution to we find the gelation time as functions of x. Let us plot it in our standard form where is the solution of the conservation law (60). Figure 4b plots , , and Q as functions of x for varied multiplicity k for tetra-functional molecules as a typical example. Multiplicity goes back to the pairwise association studied in the previous section. All curves diverge logarithmically as x approaches the critical value from above. For , we know that the gel point is not a monotonically decreasing function of k, but takes a minimum at (see Figure 7.5 (a) in the reference [4]). The retardation coefficient for f-functional primary molecules with cross-link multiplicity k can be found in a similar method of expanding in powers of the reverse rate constant . We find with where is the irreversible gelation time (51) by Tobitani–Ross–Murphy [20]. Detailed calculation is given in Appendix B.

In most physical gels, we expect that cross-linking proceeds via step-by-step association of the free functional groups: 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 for , and for . Here, for is the flux between and k state (), i.e., the number of junctions changing to k junctions per unit time as a result of forward and backward reaction with the rate constant and (see Figure 5). Similar rate equations were proposed by Aniansson and Wall [46,47,48] 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 the micellization problem they move freely. However, 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 in order to find the equilibrium solution easily. Using Equation (3), we find for , and for . Let us next confirm that the equilibrium distribution (5) satisfies these kinetic equations. On substitution of (5) into (69), we find that the relation must be fulfilled. We therefore propose the detailed balance condition such that the equlibrium constants satisfy the relation By repeated use of this relation, we have a well-known relation between the equilibrium constants and the stepwise association constants for 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).

First, we study the three-state model of cross-linking. In this model, we have free (), double (), and triple () cross-link junctions. Some biopolymers form cross-links of either double helix or triple helix [2,49,50]. Depending on the environmental condition, there is a competition between them. We can study such competition in forming cross-links by using a three-state rate equation. The kinetic equations are The conservation law of the total number of functional groups is described by the normalization condition . The equilibrium solutions are given by . Detailed balance conditions lead to the relation , and . Let us first eliminate triple association by using the normalization condition. Then, we transform the kinetic equation in terms of the deviation from the equilibrium state defined by . After a lengthy calculation, we find where for the linear terms, and for the nonlinear (quadratic) terms. We can numerically solve these equations for a given set of rate constants. However, 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 . The equation to find them is where Note that can be written as so that it is regarded as an extension of the KHDS decay rate (55) to the three-state model. (Prime indicates differentiation.) The linear relaxation has two modes The gel point is given by the condition Using the deviation , this condition is transformed to

We first look into the cases where the two rate constants of the backward reaction are largely different. If , the cross-linking is controlled by the first step association (pair formation), and vice versa. In both cases, we have the situation , so that the relaxation time is approximately given by where is the average of defined by The relaxation time apparently takes an extended KHDS form. Here, we have eliminated in favor of and from , and write by using If the first association (pairing) is the rate-determining step, the relaxation time is approximately and goes back to the pairwise association in the limit of . If the second associaition (triple junction) is the rate-determining step, the relaxation time is approximately given by From these results, we may expect the extended KHDS form for the multiplicity more than three with some average , because the average multiplicity is related to the average branching number of the cross-link junctions.

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 , the kinetic equation reduces to The rexation time is then given by The front factor can be regarded as an average backward rate constant defined by Hence, we find again the same form for the relaxation time, with slightly different average of .

The general solution for the linear equation is given by where is the normalized eigen-vector for the eigen-value (slow mode) and for the fast mode. The constants can be fixed by the initial conditions . 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 , and the gelation time is with a thermodynamic factor from (84), where is a function of , which is non-singular across the gel point. Therefore, we have again logarithmic divergence of the gel 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: Then, the forward rate constants are given by The conservation law takes the form where is written as for simplicity. Now we have In particular, for a special case of the relaxation time is given by with Figure 6a shows a schematic free energy surface of association drawn against the reaction coordinates for the three-state model. With a decrease in the ratio , the surface curve of the triple junction changes from A to D. For the limit of (curve D), the model reduces to the pairwise one. Figure 6b plots the reciprocal gelation time , the reciprocal relaxation time , and the thermodynamic factor plotted against a scaled concentration for three ratios of the stepwise constants for the tetra-functional molecules (). Due to the existence of state, relaxation is slower than the fixed-multiplicity model of .

Let us next consider cross-link junctions whose multiplicity is stabilized around a large number 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 both chain ends (telechelic polymers), such as hydrophobic ethoxylated urethane (called HEUR) [51,52,53], and hydrophobic poly(N-isopropylacrylamide) [54,55,56,57]. 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 [58,59,60], are another important example whose cross-links are micelles of various shapes; spherical, cylindrical, planar, etc.; consisting of the block segments. The size distribution of the cross-links changes as the reaction proceeds step by step, from the 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 [44,45,46,47]. The functional groups treated here are, however, attached to the polymer main chains. Therefore, their association takes place under the constraint of chain movement. However, here we assume the equal reactivity, and start from the stepwise kinetic Equations (69) and (70). Their linearized forms are The gel-point condition (48) turns into Following the idea by Aniansson and Wall [46] (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 judiciously choose a multiplicity , below which the population distribution is dominant in the initial stage, and , above which the population is dominant in the final stage. A stationary state is approximately retained between them for in the most time during the approach to equilibrium (assumption of quasi-stationary flow). We then have In Appendix C, we show that the flux J is proportional to the probability of free groups as where and are the j-th moments for small aggregates and large micelles defined by for . Here, the resistance R in the heat flow analogy is given by We also find that the probability of large micelles for is approximately given by The kinetic equation for then takes the form Hence, the relaxation time is which is similar to AW formula. The solution of (111) is given by with . Substituting (110) into the gel-point condition (105), we have and hence we finally find the form for the gelation time with the thermodynamic factor where 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 are allowed. Then, the last step from to s is expected to be slower than other steps, and hence we can choose and . Figure 7 shows the quasi-stationary approximation schematically. Circles connected by dotted lines indicate the probability deviation at time t. Since the moments are reduced to and , we have This is again KHDS form with an effective rate constant written in terms of the resistance constant 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 . Relaxation is much slower, in particular, for the high concentration region. Figure 8 shows the gelation time , the relaxation time , and the thermodynamic factor Q plotted against the scaled concentration for the functionality with the varied upper bound of the junction multiplicity. For simplicity, uniform is assumed. The case of goes back to the three-state model studied above. With an 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 as . We can study these results in more detail by comparing the experimental data of associating polymer solutions measured by using a temperature-jump technique.

On the basis of the kinetic equation of cross-linking reaction, we have found the gelation time in the temperature/concentration jump experiments. In the high concentration region, we find the power law with the power index n: for pairwise cross-linking, for simultaneous association of k functional groups, and for stepwise association with the maximum multiplicity s. These results can be compared with the experimental data for (1) to study the formation process of the cross-links. For the amplitude, we find where is the enthalpy of association. The factor is associated because there are bonds in a cross-link of multiplicity k. In the vicinity of the gel point , on the other hand, we have found a logarithmic singularity . The index corresponds to the critical exponent of a phase transition. It is generally a fractional number. However, since our kinetic equation is based on a mean-field picture (equal reactivity), we have reached a mean-field exponent . As for the temperature quench, we can do a similar analysis in the critical region by the plot We have again for the index. For the deep quench temperature jump experiment (), however, we exponential law (121). For simplicity, we assumed cold-setting gelation, i.e., gelation at low temperature. In nature, however, we often observe heat-setting gelation phenomena. Some examples are biopolymer hydrogels and temperature-responsive polymers in water. Polymers are dissolved in water due to hydration. When the temperature is raised, polymers are dehydrated, and exposed to water, followed by the intermolecular association of the hydrophobic segments interfering with LCST phase separation. For heat-setting gelation, therefore, we have to consider the dehydration process (dissociation of bound water molecules from the polymer chain segments) before the cross-linking reaction. We will treat them in a separate paper. Throughout this paper, we have confined our study to the pre-gel region before the solution reaches the gel point. In this region, the conservation law (6) of the functional groups, equilibrium reactivities (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 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 a 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 [33], and the latter by Flory [41,42,43]. Later, Ziff and Stell [35,36,37,38,39] 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.

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—simultaneous plot of and Q against the scaled concentration variable in a single graph. From the theoretical modeling, the following conclusions can be drawn: (1) The gelation time , the relaxation time , and the thermodynamic factor Q are all functions of a single variable (scaled concentration), where 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. Data for different concentrations can be superimposed onto a single curve by using an appropriate temperature shift factor. (2) These three factors obey a fundamental relation . The thermodynamic factor is logarithmically singular at the equilibrium gel point , while the relaxation time 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 of reversible cross-linking approaches the power law of the irreversible one in the asymptotic region of large x (either high concentration or high values of the association constant ). The power index of lies at somewhere between (simultaneous cross-linking) and (stepwise cross-linking). Hence, the reaction kinetics, simultaneous or stepwise, can be inferred by measuring the power. (4) For 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, and enables a comparison with experimental data. 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.