
MBW:CoagulationFragmentation ModelsFrom MathBioAuthors: Jonathan A.D. Wattis Reviewed by: Dustin Keck All work shown below can be assumed to originate from Wattis unless otherwise cited. ContentsExecutive SummaryThis paper focuses on fundamental discrete deterministic meanfield models for both coagulation and fragmentation processes. On the other hand, the paper does not describe in any detail stochastic modeling for these processes nor does it cover a host of other ways of modeling coagulation and/or fragmentation. Throughout the paper, the authors make several major assumptions. First, the systems described are composed of a large number of fundamental and identical units ('atoms') that cannot be broke down any farther. Second, these fundamental units can join together to form clusters. Third, the authors can then use law of mass action to describe the interactions of these clusters. Fourth, in each of the models discussed, the authors assume a discretization of the size of the 'atoms', while they assume time is a continuous variable. The paper then takes on three distinct sections. It first reviews Smoluchowski's coagulation equation (done in 1917) and describes some updated results. Next, this paper explains the BeckerDoring equations that account for both coagulation and fragmentation. The paper then concludes by describing a few, more complex models than either Smoluchowski's equation or the BeckerDoring model.
Context/Biological Phenomenon under ConsiderationCoagulationfragmentation processes occur across of breadth of areas. This paper mentions the formation of aerosols, colloidal aggregates (such as shaving cream or mayonnaise)^{[1]}, polymers, which are large molecules composed of repeating structural units^{[2]}, and celestial bodies on astronomical scales. Of particular interest to us, coagulationfragmentation processes also occur in bacterial pathogens such as Klebsiella pneumoniae, see APPM4390:Bacterial Aggregates in the Bloodstream. The coagulation process consists of the dynamics occurring when 'atoms' of a particular substance join together to form larger clusters of those 'atoms'. Conversely, fragmentation entails the dynamics involved when pieces of these clusters break away from their original cluster to form smaller clusters of their own. Again, as discussed in APPM4390:Bacterial Aggregates in the Bloodstream, understanding these dynamics can lead to better treatments for conditions such as bloodstream infections. HistoryMarian Smoluchowski was an ethnic Polish scientist in the AustroHungarian Empire. He was a pioneer of statistical physics and an avid mountaineer ^{[3]}. He actually described Brownian motion independently of Albert Einstein. Smoluchowski's work, inculding that on coagulation, had dramatic consequences in the early 1900's. Fulinski writes, "Let us remember that at the beginning of our century the 'molecular hypothesis' was just a hypothesis, which needed proofs, and that the first experimental proofs — honoured in 1926 by Nobel prizes in physics (for Perrin) and in chemistry (for Svedberg) — were based on Smoluchowski’s theory of Brownian motion, criticalopalescence, sedimentation and coagulation." ^{[4]} Becker and Doring were both physists. Their equations used in this paper added complexity to the model used by Smoluchowski by making coagulation a reversible process which introduces fragmentation. On the other hand, their model greatly simplified Smoluchowski by assuming that reaction of clusters only included interactions with monomers. In other words, a cluster of size 7 only includes a coagulations between clusters of size 6 and monomers (clusters of size 1) to form clusters of size 7, and clusters of size 7 joining with monomers to form clusters of size 8. The reverse process, fragmentation, of each of these are also included. Of course, there have been modifications made to these models as time has evolved. Fournier and Laurencot have considered models with homogeneous kernels of massconserving type. Also, Menon and Pego have shown a number of interesting results for the continuous Smoluchowski coagulation equations with a constant kernel. Regarding the BeckerDoring equations, Dreyer and Duderstadt, as well as Hermann et al. take a different approach with respect to the available free energy of the system. The last section of the paper discuss a few more even more general approaches taken to these types of problems. Full model with results from the paperThe first set of models discussed in this paper are derived from Smoluchowski's coagulation equation, which models the process of binary aggregation. Aggregation and coagulation are used interchangeably throughout the paper. If we want to know how cluster size distributions change over time, we can describe the process as:
where, represents a cluster of size fundamental units. Denoting the reaction rate constant as and the concentration at time of , the law of mass action gives
Specifically, this paper uses this model to answer the following questions. Is the total mass in the system conserved? Do selfsimilar solutions exist and what is their form? The authors emphasize that for a specific application, the form of is not necessarily easy to determine. While many different forms of the aggregation kernel have been used, this paper describes three specific cases in detail. In all three cases, this paper uses a couple of important concepts and their resultant equations. First they use a generating function approach. With this approach, they define
where is the total number of cluster in the system, is the total mass of the system, and the th moment of the cluster size distribution function is
All three cases also use the following initial conditions:
The first specific form of models Brownian coagulation in which case . Here the rate constant does not depend on the size of the cluster. Specifically for this example, the authors use . In combination with the generating function this results in the single partial differential equation (PDE)
which after substituting the initial conditions leads to a selfsimilar solution of the form . A selfsimilar solution means that the system eventually "forgets" its initial conditions, so that given the initial mass and compactly supported initial conditions, the longterm solution is not too complex. The second specific form of models a limited case of gravitational coagulation in which case . Specifically, the authors use . The PDE in this case is . which after substituting the initial conditions leads to a selfsimilar solution of the form . The third specific form of models branched chain polymerisation in which case . Specifically, the authors use . The resulting PDE is . which after substituting the initial conditions leads to a selfsimilar solution of the form . In this situation, the solution introduces the concept of gelation, which will be discussed further in the Analysis section of this review. The paper then goes on to show that systems with more general assumptions behave somewhat like those discussed above. Considering only what happens at large times and large cluster sizes, kernels of the following form are considered: , where with . Four types of behavior are observed under these assumption, all of which are discussed in the Analysis section of this review. The next set of models described in this paper are derived from the BeckerDoring equations. These equations allow for both aggregation and fragmentation of a single particle at a time. In this case the reaction looks like
where the forward rate is denoted and the backward rate is denoted . The paper provides a helpful example of what this means. For instance, the only reaction that involve are and where the rate of change is . (1) They also define the general rate of change of concentrations as
where is defined as . Finally, while there are a number of choices for , this paper considers two formulations of the BeckerDoring equations. In the first case, the authors assume is a constant, so . For the second case, total mass of the system is fixed, so the monomer concentration, varies. As with the Smoluchowski equations, ideally we would like explicit solutions to the BeckerDoring equations, but here again, the authors emphasize in actual application this can be extremely difficult, so studying long term behavior of systems takes precedence. They also make an important distinction between equilibrium and steadystate. Equilibrium implies the the forward and backward rates of reaction are equal vs. steadystate where the time derivatives are zero. The following highlights some of the more important models/parts of the models for the two different assumptions on the monomer concentration with the results explained in the next section. First, we must address some notation. At equilibrium, all of the fluxes, are zero, so the equilibrium concentrations are . (2) On the other hand, at steadystate, for some constant . We then have . Now in the case of constant monomer concentration we also define which shows which process, aggregation or fragmentation, is dominating the reaction. When , fragmentation dominates aggregation. After some derivation, the authors show that for , the steadystate solution is approached at large times, whereas for , the equilibrium solution is approached at large times. The paper then examines these solutions for three cases. When . (3) When we have . (4) Finally, when we have . (5) With the other approach, the authors assume constant total mass in the system, so (6) At this point, the paper summarizes the results of Ball, Carr, and Penrose as they sought solutions to (6). There are four cases, based on specific combinations of aggregation and fragmentation rate coefficients, in which metastability occurs in the last two. In this context, metastability is when the system is resistant to many perturbations away from that state, but still evolves on some trajectory. For case 1, we have , where is the radius of convergence of . where is the total mass of the equilibrium solution and . Then choosing rates such that fragmentation dominates, . For case 2, we have . Here for any we have for any intial mass . For case 3, we have . Here has a maximum value defined as . Now, if we get similar results as case 2, but if we get weak convergence. Then with very specific aggregation and fragmentation rates, we have . Then for case 4, we have . In this case, the only equilibrium solution is which means for all as , so that for every we have . The models in the final section of the paper reflect more generalized processes of coagulation and fragmentation. They mention three processes in particular. First, the authors consider binary fragmentation which generalizes the stepwise fragmentation of the BeckerDoring equations. In this case, we have (7) where under the assumption that , we get . (8) Second, the authors briefly describe Smoluchowski coagulationfragmentation, in which case . Here, they are still only considering binary fragmentation where both the aggregation rate, and the fragmentation rate, are symmetric. Then applying the law of mass action, (9) . (10) Third, the authors examine multicomponent coagulationfragmentation, where multiple species form heterogeneous clusters. A couple of examples are offered. In one such example, if each cluster, has units of type A and units of type B, we have , (11) where (12) defined on the region Analysis/interpretationThe Smoluchowski equations only allow for coagulation, and from the three specific cases of aggregation kernels discussed, we learn a little about the form of the distribution of large clusters at large times. In the first two cases where and we had a conservation of mass and the shape of the distribution decays exponentially in . In the third case where we had gelation. When gelation occurs, there is a finite time at which points all of the clusters form a single mass (or gel) and then continues to grow. For the more general assumptions, where , we had four types of behaviors. First, there was no gelation. This occurs when , and , which was like the first two specific cases above. The second behavior, which occurs when , and ,was delayed gelation just as we saw in the third specific case above. The third behavior was instantaneous gelation, which occurs when and max. The fourth behavior was nonexistence of solutions which occurs when . These four regions of behavior are illustrated in Fig. 1. The BeckerDoring equations add complexity to the modelling in that now fragmentation is allowed, but only clustermonomer interactions are allowed, which actually simplifies the modelling as shown in equation 1. As mentioned in the modeling section, the BeckerDoring equations are divided under two different assumptions. First, the authors assume that the monomer concentration is constant. They then use as a measure of the relative strength of aggregation to fragmentation. In each of the three cases for values of we get different distributions of the concentrations of cluster sizes (equations 3, 4, and 5). When , the evolution of the evolution of the distribution function to its equilibrium state can be seen in Fig. 2. Three regions can be seen. First, when is small, is already at its equilibrium. Then there is a transition region where goes to its equilibrium, equation 2. Finally, in the region where r is large, stays close to its initial conditions. The authors highlight that if one wants to study the evolution of for some large , the concentration of clusters size remains low for a long time, before it goes to and remains at its equilibrium value. When , evolution of the distribution function is a stationary diffusive wave which can be seen in Fig. 3. When the equilibrium solution diverges, but the system converges to a steadystate solution as shown in Fig. 4.
In the first two cases, given any there is a unique solution . The difference is that in case 1, as , but in case 2, as . On the other hand, analysis gets more complicated in cases 3 and 4. In case 3, has a maximum and at that maximum the total mass is finite. When the initial mass is less than the mass at maximum , results are similar to case 2. Conversely, when the initial mass is greater than the mass at maximum , we have weak convergence to a metastable process. The authors offer this comparison. In the Smoluchowski model, when mass was not conserved we had gelation (a kind of superparticle). With the BeckerDoring model, the kinetics are slower and gelation can't occur at some finite time. This is where they make an important distinction. Metastable states last a long time and evolve quite slowly, but they are not equilibrium states. In the final section of the paper, the authors show how the previous models can be generalized. In the first generalized model, binary fragmentation is considered described by equations 7 and 8. They point out that solutions to those equations have been shown to exist, but they do not conserve mass. In the next generalized model, the authors study Smoluchowski coagulationfragmentation. They use the law of mass action, symmetric rate constants, and derive a system of equations (9 and 10) which behave much like the BeckerDoring equations. They point out that under a variety of different aggregation and fragmentation kernels, many authors have proven existence and uniqueness of solutions. Finally, the authors examine multicomponent coagulationfragmentation where multiple species form heterogeneous clusters. In the example defined by equations 11 and 12, they make several new modifications to the systems and expect to conserve quantities of both components A and B. With all of these models, the authors remind us of several challenges. One is determining appropriate rate coefficients. Another is the understanding of multicomponent systems. They conclude with the motivation that this area of study has many unanswered questions with a myriad of applications. References/External Links
