Next Article in Journal
Phase Separation and Pairing Fluctuations in Oxide Materials
Previous Article in Journal
Unconventional Transport Properties of Reduced Tungsten Oxide WO2.9
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Mechanics of Discrete Multicomponent Fragmentation

Department of Chemical Engineering, Pennsylvania State University, University Park, PA 16802, USA
Submission received: 7 September 2020 / Revised: 13 October 2020 / Accepted: 15 October 2020 / Published: 18 October 2020
(This article belongs to the Section Condensed Matter Theory)

Abstract

:
We formulate the statistics of the discrete multicomponent fragmentation event using a methodology borrowed from statistical mechanics. We generate the ensemble of all feasible distributions that can be formed when a single integer multicomponent mass is broken into fixed number of fragments and calculate the combinatorial multiplicity of all distributions in the set. We define random fragmentation by the condition that the probability of distribution be proportional to its multiplicity, and obtain the partition function and the mean distribution in closed form. We then introduce a functional that biases the probability of distribution to produce in a systematic manner fragment distributions that deviate to any arbitrary degree from the random case. We corroborate the results of the theory by Monte Carlo simulation, and demonstrate examples in which components in sieve cuts of the fragment distribution undergo preferential mixing or segregation relative to the parent particle.

1. Introduction

Objects disintegrate into fragments via impact, detonation, degradation, or cleavage of the bonds that hold the structure together. The object in question may range from the sub-atomic [1,2] and molecular [3,4,5] to living organisms [6], social structures [7,8], and celestial bodies [9], a diversity of scale and physics that is united by a common mathematical formalism. At its core, fragmentation is a branching process in which a parent object (“particle”) produces a set of offspring. The evolution of a population that undergoes splitting of this form is given by the fragmentation equation, an integro-differential equation that accounts for the generation and depletion of sizes due to fragmentation [10,11]. The primary input to this formulation is a breakup model that specifies the distribution of fragments produced by a given parent size and the relative rate at which different sizes break up. This population balance approach forms the basis for the mathematical treatment and numerical modeling of fragmentation in granular, colloid, and polymeric systems [11,12,13,14,15,16,17,18]. The mathematical literature of the fragmentation equation is rich and focuses on analytic solutions, existence criteria, and stability. Of particular interest is the emergence of “shattering”, a process akin to a phase transition that is demonstrated through the appearance of a finite population of particles with zero mass [19,20,21,22]. An alternative approach views fragmentation as the disintegration of bonds between the constitutive units of the particle and uses percolation theory to model and simulate the breakup of systems with topological structure. In contrast to the population balance method, which is a mean-field method, percolation treats fragmentation at the discrete probabilistic level [7,8,23]. Other treatments view fragmentation in a more abstract way as a partitioning of a discrete event space and use combinatorial and probabilistic methods to obtain the partition function and the mean distribution in this space [24,25,26,27].
A central question in fragmentation is the distribution of fragments per fragmentation event. The most common theoretical model is that of random binary fragmentation. In this model a parent cluster produces two fragments with uniform probability [22]. Empirical models have been proposed for the breakage of a single component into multiple pieces of unequal size and typically require a set of parameters that control the shape of the fragment distribution [28,29]. Systems of practical interest are almost always multicomponent. Pharmaceutical granulation is a case in point: granules contain an active pharmaceutical ingredient, an inert excipient, binder, and are characterized by additional attributes such as porosity or shape factor that behave as pseudo components [30]. Nonetheless, no generalized approaches exist to treat the problem of multicomponent fragmentation into arbitrary number of fragments. The mathematical treatment of multicomponent fragmentation into any number of pieces cannot be simply obtained as an extension of the one-component problem. In addition to the size distribution of fragments one must consider the distribution of components, provide rules for apportioning components to the fragments, offer a definition of what is meant by “random fragmentation” when both size and composition are distributed, and provide the means for constructing models that deviate from the random case to any extent desired.
The purpose of this paper is to address these questions by formulating the statistics of a single fragmentation event in the discrete domain for arbitrary number of fragments and components in a way that is general and not bound by the details of the particular application. Our interest is not in the physics behind the splitting of an object into smaller parts, but rather in the probabilistic treatment of the partitioning itself under the constraint of a conservation law, in this case conservation of mass. The main idea is this. We start with a multicomponent particle that is made of discrete units of any number of components, subject it to one fragmentation event with fixed number of fragments, and construct the set of all fragment distributions that can be obtained in this manner. We calculate the partition function of this ensemble of fragments, assign probabilities in proportion to the multiplicity of each distribution, and obtain the mean distribution in terms of the partition function. We then introduce a bias functional that biases the distribution away from that of random fragmentation. We present results from Monte Carlo simulations to corroborate the theory and show that components may preferentially mix or unmix in the fragments depending on the choice of the bias functional.

2. Random Fragmentation

2.1. One-Component Random Fragmentation

In discrete fragmentation, a particle composed of M integer units breaks up into N fragments, { m 1 , m 2 m N } that satisfy the mass balance condition
i = 1 N m i = M .
The distribution of fragments is given by vector n = ( n 1 , n 2 ) whose element n i is the number of fragments that contain i units of mass. We suppose that N is fixed but n is not; that is, if the fragmentation event is repeated with an identical parent particle the distribution of fragments may be different but the total number of fragments is always N. We refer to this process as N-nary fragmentation. All fragment distributions produced by this mechanism satisfy the following two conditions:
i = 1 n i = N ,
i = 1 i n i = M .
The first condition states that the number of fragments is N; the second that their mass is equal to the mass of the parent particle. Conversely, any distribution that satisfies the above two equations is a feasible distribution of fragments by N-nary fragmentation of mass M. Thus the set E M ; N of all distributions that satisfy Equation (2) forms the ensemble of fragment distributions produced from M. We assign a probability P ( n ) on the distributions in E M ; N , normalized over all distributions that satisfy Equations (2) and (3). Our goal is to obtain this distribution under various fragmentation models.
We will call the process random fragmentation if all ordered lists of N fragments produced by the same mass are equally probable. This views the ordered list of fragments, which we call configuration, as the elementary stochastic variable in this problem.

2.1.1. Probability of Random Fragment Distribution

Proposition 1.
The probability of distribution n produced by random N-nary fragmentation of mass M is
P ( n ) = n ! M 1 N 1 ,
where n ! = ( n 1 , n 2 ) ! is the multinomial coefficient of distribution n ,
n ! = ( i n i ) ! i n i ! = N ! n 1 ! n 2 ! .
Proof. 
First we note that the number of ordered lists that can be formed by breaking integer M into N fragments is
Ω M ; N ( 1 ) = M 1 N 1 .
This is the number of ways to partition integer M into N parts and can be shown easily as follows [31]: thread M balls into a string and partition them into N pieces by cutting the string at N 1 points (Figure 1). There are M 1 points where we can cut and must choose N 1 of them. The number of ways to do this is the binomial factor on the RHS of Equation (6).
If all ordered lists of fragments are equally probable, the probability of ordered list m = ( m 1 , m 2 m N ) is
Prob ( m ) = 1 Ω M ; N ( 1 ) .
There are n ! ordered lists with the same distribution of fragments n . Accordingly, the probability of n is
P ( n ) = n !   Prob ( m ) = n ! Ω M ; N ( 1 ) .
This proves the proposition. □
The multinomial factor n ! is the multiplicity of distribution n , namely, the number of configurations (ordered lists of fragments) represented by n . Using ω ( n ) = n ! to notate this multiplicity, the probability of distribution is expressed as
P ( n ) = ω ( n ) Ω M ; N ( 1 ) ,
and Ω ( 1 ) satisfies
n ω ( n ) = Ω M ; N ( 1 ) .
The summation is over all distributions n E M ; N , namely, over all distributions that satisfy Equations (2) and (3). Accordingly, Ω M ; N ( 1 ) is the total multiplicity in the ensemble, equal to the number of ordered configurations of fragments that can be produced from integer mass M breaking into N fragments. We refer to Ω M ; N ( 1 ) as the partition function of the one-component ensemble of fragments.

2.1.2. Mean Fragment Distribution

Each distribution n appears in the ensemble of fragment distributions with probability P ( n ) ; the mean distribution of fragments is their ensemble average:
n = n n   P ( n ) ,
with P ( n ) from Equation (4) and with the summation going over all distributions that are produced by N-nary fragmentation of integer mass M.
Proposition 2.
The mean distribution in N-nary random fragmentation is
n k N = Ω M k ; N 1 ( 1 ) Ω M ; N ( 1 ) = M k 1 N 2 M 1 N 1 ,
with M N 2 and k = 1 , M N + 1 .
Proof. 
First we write the probability of distribution in the form
P ( n ) = N ! Ω M ; N ( 1 ) i = 1 α i n i n i !
with α i > 0 and note that this reverts to Equation (4) when α i = 1 . We will retain the factors α i and will set them equal to 1 at the end. The normalization condition on the probability P ( n ) reads
Ω M ; N ( 1 ) = N ! n i = 1 α i n i n i ! .
The derivative of log Ω M ; N ( 1 ) with respect to α k is
log Ω M ; N ( 1 ) α k = N ! α k Ω M ; N ( 1 ) n n k i α i n i n i ! = n k α k ,
where n k is the mean value of n k in the ensemble of fragments. We also have
Ω M ; N ( 1 ) α k = N ( N 1 ) ! n n k 0 α i n k 1 ( n k 1 ) ! = N   Ω M k ; N 1 ( 1 ) .
The summand in the expression in the middle amounts to removing one fragment of mass k from all distributions of the ensemble; accordingly, the quantity in braces is the partition function Ω M k ; N 1 ( 1 ) . Combining Equations (15) and (16) and setting α k = 1 we obtain Equation (12). A similar proof was given by Durrett et al. [25] for a closely related system. □
Equation (12) was previously obtained by Montroll and Simha [12] via a combinatorial derivation. Notably, it is the same distribution as in discrete binary aggregation (the reverse process of binary fragmentation) with constant kernel, derived by Hendriks et al. [24] who also credit older unpublished work by White. It also appears outside the context of fragmentation when the probability distribution is of the form of Equation (13) (see for example [32]).
For large M the fragment distribution becomes
n k N ( N 1 ) M 1 k M N 2 .
This is the continuous limit of random fragmentation of a straight line into N segments, an elementary result that has been derived multiple times in the literature. The earliest report known to us is by Feller [33] who corrected an earlier approximation by Ruark [34].

2.2. Two-Component Random Fragmentation

2.2.1. Representations of Bicomponent Populations

We now consider a particle that is made of two components. The particle contains M A units of component A, M B units of component B, and its mass is M = M A + M B . The distribution of fragments is given by the two-dimensional vector n = { n a , b } where n a , b is the number of fragments that contain a units of A and b units of B. This distribution satisfies the conditions
  a = 0 b = 0 n a , b = N ,
  a = 0 b = 0 a   n a , b = M A ,
  a = 0 b = 0 b   n a , b = M B .
The set E M A , M B ; N of all distributions that satisfy the above conditions constitutes the set of feasible distributions in bicomponent fragmentation. Strictly speaking, the upper limit in these summations is constrained by a M A , b M B . Under the convention that n a , b = 0 outside the meaningful range of a and b, we may set the upper limit to .
The color-blind size distribution or simply “size distribution” n A + B = { n k } is the distribution of the mass of the fragments k = a + b regardless of composition:
n k = a = 0 k n a , k a ,   k = 1 , 2
and satisfies the conditions
  k = 1 n k = N ,
  k = 1 k n k = M A + M B = M .
These are the same as Equations (2) and (3) in the one-component case for a particle with mass M A + M B . Accordingly, the feasible set of the color-blind distribution is E M ; N with M = M A + M B .
The sieve-cut distribution n A | k = { n a | k } is the number of fragments with size k that contains a units of component A:
n a | k = n a , k a ,   ( a = 1 k , k = 1 M ) .
and satisfies the normalizations
  k = 1 a = 0 k n a | k = N ,   k = 1 a = 0 k a n a | k = M A ,   k = 1 a = 0 k k n a | k = M A + M B .
We divide the sieve-cut distribution by the number of fragments of size k to obtain the compositional distribution of component A within fragments of fixed size k,
c a | k = n a | k n k .
The compositional distribution is normalized to unity and may be interpreted as the conditional probability to obtain a fragment with a units of A, given that the fragment has mass k. The bicomponent distribution may now be expressed in terms of the color-blind distribution n A + B and the compositional distribution c a | k in the form
n a , k a = n k   c a | k .
If we divide both sides by the total number of fragments, the result reads as a joint probability: the probability n a , k a / N to obtain a fragment with mass k that contains a units of component A is equal to the probability n k / N to obtain a fragment of mass k times the probability c a | k to obtain a fragment with a units of component A given that the mass of the fragment is k.

2.2.2. The Ensemble of Random Fragment Distributions

Random fragmentation is implemented by analogy to the one-component case: we line up the unit masses in the particle into a string and cut at N 1 places. Every cut is equally probable and so is every permutation in the order of the beads.
Proposition 3.
The probability of fragment distribution n in random bicomponent fragmentation is
P ( n ) = n ! Ω M A , M B ; N ( 2 ) a = 0 b = 0 a + b a n a , b ,
where n ! is the multinomial coefficient of the bicomponent distribution,
n ! = N ! a = 0 b = 0 n a , b !
and Ω M A , M B ; N ( 2 ) is the two-component partition function, given by
Ω M A , M B ; N ( 2 ) = M A + M B M A Ω M A + M B ; N ( 1 ) .
Proof. 
First we count the number of ordered sequences of fragments (configurations). Configurations are distinguished by the order the fragments and by the order of components within fragments (Figure 1). We color the components and place them in a line in some order. There are M A As and M B Bs; the number of permutations is M A + M B M A . Each permutation produces Ω M A + M B ; N ( 1 ) configurations with Ω ( 1 ) given in Equation (6). The total number of configurations, therefore, is their product and proves Equation (29).
Since all configurations are equally probable, the probability of fragment distribution n is proportional to the number of configurations with that distribution. This is equal to the number of permutations in the order of the fragments times the number of permutations in the order of components within the fragments. The number of permutations in the order of fragments is given by the multinomial factor of bicomponent distribution in Equation (28). The number of permutations of components within a fragment that contains a units of A and b units of B is a + b a and since there are n a , b such fragments, the total number of internal permutations in distribution n is
a = 0 b = 0 a + b a n a , b .
The probability of distribution n is equal to the product of Equations (28) and (30) divided by the total number of configurations, given by Equation (29). The result is Equation (27) and proves the proposition. □
As a corollary, we obtain the multiplicity of the bicomponent distribution,
ω ( n ) = n ! a = 0 b = 0 a + b a n a , b .
Thus, we write
P ( n ) = ω ( n ) Ω M A , M B ; N ( 2 )
with Ω ( 2 ) = n ω ( n ) .
An alternative result for P ( n ) is obtained by expressing the bicomponent distribution n in terms of the color-blind distribution n A + B and all sieve-cut distributions n A | k . The result is
P ( n ) = n A + B ! Ω M A , M B ; N ( 2 ) k = 1 n A | k ! a = 0 k k a n a | k
and is based on the identity
n ! a = 0 b = 0 k a n a , b = n A + B ! k = 0 n A | k ! a = 0 k k a n a | k .
Here, n A + B ! is the multinomial coefficient of the color-blind distribution,
n A + B ! = N ! n 1 !   n 2 !
and n A | k ! is the multinomial coefficient of the sieve-cut distribution,
n A | k ! = n k ! n 0 | k !   n 1 | k !   n k | k ! .

2.2.3. Mean Fragment Distribution

Proposition 4.
The mean distribution of fragments in random bicomponent fragmentation is
n a , b N = a + b a Ω M A a , M B b ; N 1 ( 2 ) Ω M A , M B ; N ( 2 )
Proof. 
The proof follows in the steps of the one-component problem. We express the multiplicity and the partition function in the form
ω ( n ) = N ! a = 0 a = b α a , b n a , b n a , b ! ,
Ω M A , M B ; N ( 2 ) = N ! n a = 0 a = b α a , b n a , b n a , b !
With α a , b = a + b a we recover the result for random fragmentation but for the derivation we treat α a , b as a variable. Following the same procedure that led to Equation (12) we now obtain
n a , b N = α a , b Ω M A a , M B b ; N 1 ( 2 ) Ω M A , M B ; N ( 2 ) .
To arrive at this result, we note that differentiation of the partition function with respect to α a , b by analogy to Equation (16) amounts to removing one cluster that contains a units of A and b units of B, thus producing the partition function Ω M A a , M B b ; N 1 ( 2 ) in the numerator of Equation (40). Setting α a , b = a + b a we obtain Equation (37). □
Alternative Proof.
An alternative proof will be obtained by a mean-field argument. First we write the mean distribution in the form
n a | k = n k c a | k ¯ ,
where n a | k and n k are the ensemble averages of n a | k and n k , respectively, and c a | k ¯ = n a | k / n k is the compositional distribution within the mean distribution. We begin with the observation that the mean color-blind distribution is the same as in the one-component case. This follows from the fact that the choice of the points at which the string of beads is cut is independent of the compositional makeup of the particle (Figure 1). Thus, n k is given by Equation (12) with M = M A + M B :
n k N = Ω M A + M B k ; N 1 ( 1 ) Ω M A + M B ; N ( 1 )
We obtain the compositional distribution by the following construction. Imagine that all possible distributions are stacked vertically to form a table so that column 1 contains the first fragment in all distributions, column 2 contains all second fragments, and so on. All columns are permutations of each other (this follows from the construction of the fragments illustrated in Figure 1) and since all permutations are equally likely (this follows from the condition of random fragmentation), all columns have the same fragment and compositional distribution; therefore, we only need to consider one of them. The equivalent problem now is this: count the number of ways to select a beads from a pool of M A As and k a beads from a pool of M B Bs and take its ratio over the total number of ways to pick k beads:
c a | k ¯ = M A a M B k a M A + M B k .
The mean distribution then is the product of the size and compositional distributions:
n a | k = n k   c a | k ¯ ,
or
n a , b N = M A a M B b M A + M B a + b Ω M A + M B a b ; N 1 ( 1 ) Ω M A + M B ; N ( 1 )
It is straightforward algebra to show that this is equivalent to Equation (37).
For M A a , M B b , the compositional distribution goes over to the binomial:
M A a M B b M A + M B a + b a + b a ϕ A a ϕ B b ,
with ϕ A = M A / ( M A + M B ) , ϕ B = 1 ϕ A . It is a more compact expression than Equation (43), but is valid only in the asymptotic limit.
Figure 2 shows compositional distributions for a bicomponent particle with M A = 4 units of A and M B = 3 units of B. As a means of a demonstration we show the results of a Monte Carlo simulation, which are seen to be in excellent agreement with theory. The binomial distribution, also shown for comparison, is only in qualitative agreement because the fragment masses are small and the conditions for asymptotic behavior are not met in this case.

2.3. Any Number of Components

Extension to any number of components follows in a straightforward manner from the bicomponent case but the notation becomes less transparent. Suppose the parent particle consists of K components A , B and contains M A units of A, M B units of B, and so on. The distribution of fragments is now expressed by the K-dimensional vector n = { n a , b } that gives the number of fragments that contain a units of component A, b units of B, etc. This distribution satisfies
  a , b n a , b = N
  a , b z n a , b = M Z ;   z = a , b
where M Z is the mass of component z = a , b in the parent particle. The set of all distributions that satisfy the above conditions constitutes the ensemble of all distributions that are produced by the fragmentation of the parent particle into N fragments.
Random fragmentation is once again implemented as shown in Figure 1: Given a string of colored beads, we cut it at N 1 random points to produce N fragments. All permutations of the beads are equally probable. Accordingly, configurations are equally probable. The number of configurations is
Ω M ; N ( K ) = M !   Ω M ; N ( 1 ) = M ! M A ! M B ! M 1 N 1 ,
where M = M A + M B + is the total mass of the particle and M ! = ( M A , M B ) ! . The multiplicity ω ( n ) of distribution n is the number of configurations with that distribution and is given by the number of permutations in the order of fragments and in the order of components within each fragment:
ω ( n ) = n ! a , b ( a + b + ) ! a ! b ! n a , b = n ! c c ! n c
with c = ( a , b ) . The probability of distribution n is
P ( n ) = ω ( n ) Ω M ; N ( K )
and the partition function is the sum of multiplicities of in the ensemble:
Ω M ; N ( K ) = n ω ( n ) .
The mean distribution of fragments is
n c | k N = c ! Ω M c ; N 1 ( K ) Ω M ; N ( K )
and is the generalization of (37). Alternatively, the mean distribution can be expressed by analogy to Equation (44) as the product of the color blind distribution with a mean compositional distribution:
n c | k = n k   c c | k ¯ .
The mean color-blind size distribution n k / N is the same as in one-component fragmentation,
n k N = Ω M k ; N 1 ( 1 ) Ω M ; N ( 1 )
with M = M A + M B + , k = a + b + , and c c ¯ is the conditional probability that the compositional vector of fragment size k in the mean distribution is c = ( a , b ) :
c c | k ¯ = M A a M B b M a + b + .
This is the generalization of Equation (43).

3. Nonrandom Bicomponent Fragmentation

In random fragmentation, all permutations are equally probable. We now bias the probability of the permutation by a functional W ( n ) of the fragment distribution such that the probability of fragment distribution n is
P ( n ) = ω ( n ) W ( n ) Ω ˜ M A , M B ; N ( K )
with
Ω ˜ M A , M B ; N ( K ) = n ω ( n ) W ( n ) .
and ω ( n ) from Equation (50). Here, ω ( n ) is the intrinsic multiplicity of n in the ensemble of fragments, while the product ω ( n ) W ( n ) ω ˜ ( n ) is its apparent (biased) multiplicity as weighted by the bias functional and distinguished by the tilde. Similarly, the partition function Ω ˜ is the summation of the apparent (biased) multiplicities of all distributions in E M A , M B ; N . All permutations in the same configuration of fragments are equally probable under this formulation, as they all have the same distribution n . With W = 1 , we recover the random case (all permutations in all configurations are equally probable). Accordingly, “random” and “unbiased” both refer to uniform bias W = 1 .

3.1. Linear Ensemble

The bias functional W will remain unspecified. This allows us to choose the bias so as to produce any desired distribution of fragments. A special but important case is when W is of the product form
W ( n ) = a b ( w a , b ) n a , b ,
where w a , b are factors that depend on a and b but not on the fragment distribution itself. The log of the bias is then a linear functional of n :
log W ( n ) = a b n a , b log w a , b .
The result states that the log of the bias is homogeneous functional of n with degree 1, i.e., log ( λ n ) = λ log W ( n ) for any λ > 0 . We refer to this case as linear bias with the understanding that linearity actually refers to the log of W.
Proposition 5.
The mean distribution of fragments under the bias in Equation (59) is
n a , b N = w a , b a + b a Ω ˜ M A a , M B b ; N 1 ( 2 ) Ω ˜ M A , M B ; N ( 2 )
with
Ω ˜ M A , M B ; N ( 2 ) = N ! n a = 0 b = 0 w a , b n a , b n a , b ! a + b a n a , b
and the summation over all n that satisfy Equations (18)–(20).
Proof. 
We write the apparent multiplicity ω ˜ ( n ) of distribution n as
ω ˜ ( n ) = N ! a = 0 a = b ( α a , b ) n a , b n a , b ! ,
and the probability of distribution as
P ( n ) = N ! Ω ˜ M A , M B ; N a = 0 a = b ( α a , b ) n a , b n a , b ! ,
with
α a , b = w a , b a + b a ,
The claim of Proposition 5 then follows directly from Proposition 4. □

3.2. Composition-Independent Bias

If the bias factors are of the form w a , b = g a + b , where g k is a function of a single variable, the acceptance probability of a configuration of fragments depends on the mass k = a + b of the fragment but not on its composition. This leads to a simple expression for the mean distribution by the following argument. With reference to Figure 1, fix the points where the string is cut; this amounts to fixing the color blind distribution of the fragments. All permutations of components are equally probable because they have the same distribution. Accordingly, the compositional distribution is the same as in the random case and is given by Equation (43). The size distribution on the other hand is biased and is the same as when the same bias is applied to one-component distribution. The final result is
n a , b N = M A a M B b M A + M B a + b n a + b N ,
where n a + b = n k is the one-component size distribution under bias w a , b = g a + b ,
n k N = g k Ω ˜ M A + M B k ; N 1 ( 1 ) Ω ˜ M A + M B ; N ( 1 )
with
Ω ˜ M A + M B ; N ( 1 ) = N ! n k = 1 g k n k n k ! .
Except for special forms of g k the partition function will not be generally available in closed form. Table 1 summarizes three cases for which exact results are possible. All three cases are associated with distributions encountered in binary aggregation [35]. The partition functions in cases 1 and 2 refer to the constant and sum kernels, respectively, and are exact; case 3 is associated with the product kernel and gives the asymptotic limit of the partition function for M , N 1 , M / N < 2 , conditions that refer to the pre-gel state [36].
In the general case, w a , b depends on both a and b explicitly and affects both the size and compositional distributions. This case will be demonstrated by simulation in the next section.

4. Simulation of Biased Fragmentation

Except for certain special forms of the bias the mean fragment distribution cannot be calculated analytically and the only recourse is stochastic simulation. Here, we describe a Monte Carlo (MC) algorithm for sampling the ensemble of distributions. We will then use this method to demonstrate results for two cases of biased bicomponent fragmentation.

4.1. Monte Carlo Sampling by Exchange Reaction

Suppose m = ( ( a 1 , b 1 ) , ( a N , b N ) ) is a configuration of N bicomponent fragments, such that fragment i contains a i units of component A and b i units of component B. The probability of configuration is equal to the probability of its distribution, P ( n ) , divided by its multiplicity ω ( n ) ; from Equation (57) this probability is
Prob ( m ) = W ( m ) Ω M A , M B ; N ( 2 ) ,
where W ( m ) = W ( n ) is the bias of configuration m , equal to the bias of its distribution n . If the bias functional is of form in Equation (59), its value on configuration m is
W ( m ) = i = 1 N w a i , b i .
Here, the product is over the N fragments in the configuration, whereas in Equation (59) it is over all all units a and b in the distribution. Suppose that two fragments i and j exchange mass according to the reaction
( a i , b i ) + ( a j , b j ) ( a i , b i ) + ( a j , b j )
under the conservation conditions a i + a j = a i + a j and b i + b j = b i + b j . This amounts to a transition m m between configurations with equilibrium constant
K m m = Prob ( m ) Prob ( m ) = W ( m ) W ( m ) = w a i , b i w a j , b j w a i , b i w a j , b j .
The stationary distribution of this exchange reaction is the same as that in Equation (57) [35]. Accordingly, the ensemble of fragment distributions may be sampled via exchange reactions by tuning the equilibrium constant to the selection functional according to Equation (72).
To implement this sampling computationally, we represent fragments as a list of 1s (representing component A) and 0s (component B). We pick two clusters i and j at random, merge them into a single list, randomize the order of components, and break them into two new fragments by picking a break point at random. We accept the resulting configuration by the Metropolis criterion based on the equilibrium constant in Equation (72): we accept the result of the exchange if rnd K m m , where rnd is a random number uniformly distributed in ( 0 , 1 ) ; otherwise we reject. With W = 1 , every exchange reaction is accepted, which amounts to random fragmentation. The randomization of the order of components in the merged list ensures that all permutations are considered with equal probability. If W = 1 , the resulting configuration is always accepted and the distribution conforms to random fragmentation. This is how the MC results in Figure 2 were obtained and the agreement with the theoretical distribution serves as a validation of the numerical algorithm.

4.2. Two Examples

In random fragmentation ( w a , b = 1 ), the compositional distribution is given by Equation (43). We may choose the bias functional so as to produce deviations in either direction relative to the random case. It is possible to produce positive deviations (preferential segregation of components in the fragments relative to random mixing) or negative deviations (more intimate mixing than in random mixing). We demonstrate both behaviors using the two examples below:
  • Case I (positive deviations)
    w a , b = ( a + 1 ) α + ( b + 1 ) α
  • Case II (negative deviations)
    w a , b = ( a + 1 ) α ( b + 1 ) α
In case I, the fragment bias w a , b is an additive function of the amounts of the two components. Considering that a + b is constrained by mass balance, the fragment bias is large for fragments that are rich in either component but small for fragments that are relatively mixed. This ought to favor the formation of fragments in which the components are relative segregated. The fragment bias in case II is a multiplicative function of the amounts of the two components. It is large in fragments that contain both components but quite small if one component is present in excess of the other. This form ought to produce fragments that are better mixed than fragments produced by random fragmentation.
We test these behaviors in Figure 3 which shows results for α = 4 . In this example the particle contains an equal number of units of each component, M A = M B = 20 , and breaks into N = 4 pieces. In both cases the size distribution deviates from that in random fragmentation. Compositional distributions are shown for sieve-cut masses k = 2 , 4, and 8. The additive bias (case I) produces distributions that are more spread out relative to the random case. For k = 2 , in particular, the compositional distribution is inverted relative to the random case. This indicates strong segregation, as the majority of fragments contains either A or B, while only few fragments in this size contain both components. As the fragment size increases the segregation of components is less pronounced, though always present, as indicated by the fact that the random distribution is always narrower. The opposite behavior is observed in case II: distributions are narrower than those in random fragmentation, especially in the smaller fragment sizes. In this case, the fragments are better mixed relative to the parent particle. As a general trend in both cases, deviations from random mixing are most pronounced in small fragment sizes. Large fragments on the other hand are close to randomly mixed. There is simply not enough material to produce a large fragment that is highly enriched in one component; thus mixing prevails. This limitation is not present in small clusters.

5. Concluding Remarks

We have presented a treatment of multicomponent fragmentation on the basis of random fragmentation in combination with a functional that biases the ensemble of feasible distributions. The two key notions in this treatment are the set of feasible distributions and the multiplicity of distribution within this set as established by the rules that define “random” fragmentation. In the random-fragmentation ensemble, distributions are proportional to their multiplicity. This problem is analytically tractable and we have presented its solution for any number of components and fragments. A third key notion is that of the bias functional that modulates the probability of feasible distributions and allows us to obtain fragment distributions other than that of random fragmentation with deviations in either direction.
Random fragmentation is not endowed with any special universality. In certain problems, such as the linear chain in Figure 1, selecting the bonds to break at random might be a reasonable physical model and random fragmentation applies; in general though, this will not be the case. The importance of random fragmentation is mathematical. Similar to the “fair coin” or the “ideal solution,” it provides an analytically solvable reference case from which to measure deviations in systems that do not conform to this model. The tool that quantifies these deviations is the bias functional. This functional, analogous to the activity coefficient in solution thermodynamics [37], permits the systematic construction of distributions that exhibit any degree of deviation from the random case. This is a key result of this formulation.
In single-component fragmentation, the quantity of interest is the mean size distribution of the fragments. In multicomponent systems we are additionally concerned with the compositional distribution of the fragments. This introduces a new dimension to the problem and raises questions of mixing and unmixing of components. Do fragments inherit the compositional characteristics of the parent particle? Do they become progressively more well mixed or less? Both behaviors are possible and are quantified via the bias functional W. This functional is where the mathematical theory of fragmentation presented here makes contact with the physical mechanisms that lead to the disintegration of material particles. To make this connection quantitatively, one must begin with the a physical model of fragmentation that assigns probabilities to all possible distributions of fragments that can be generated. This is a major undertaking and is specific to the particular problem that is being considered. The point we wish to make is that once such results are available, their reduction into a compositional distribution passes through the bias functional, which represents the contact point between physics and the mathematical formulation of fragmentation.
Lastly, the connection to statistical mechanics should not be lost. We have constructed an ensemble whose fundamental element (“microstate”) is a the ordered configuration of fragments; its total number in the ensemble is the partition function. The higher-level stochastic variable (the observable) is the distribution of fragments and its probability is determined by its multiplicity in the ensemble. The form of the probability in Equation (13), also known as Gibbs distribution [27], is encountered in time reversible processes as well as in population balances of aggregation and breakup [24,25,26,27,38]. The derivation of the mean distribution in the random case follows in the steps of the Darwin–Fowler method [39]. Additionally, the compositional distribution in random breakup is given asymptotically by the binomial distribution in Equation (46). This establishes a reference for compositional interactions analogous to that of the ideal solution in thermodynamics. In fact, the Shannon entropy of the binomial distribution is the ideal entropy of mixing when two pure components coalesce into a single particle that contains mass fraction ϕ A of component A. These connections are not coincidental. Biased sampling from a distribution generates a probability space of distributions and when the base distribution is exponential, this ensemble obeys thermodynamics [40]. In fragmentation the base distribution is a multicomponent exponential: the size distribution in Equation (12) goes over to the exponential distribution when M , N 1 . In this limit, the ensemble of fragments becomes mathematically equivalent to a thermodynamic ensemble of two components with interactions that produce positive or negative deviations relative to ideal solution.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References and Note

  1. Gross, D.H.E. Microcanonical thermodynamics and statistical fragmentation of dissipative systems. The topological structure of the N-body phase space. Phys. Rep. 1997, 279, 119–201. [Google Scholar] [CrossRef]
  2. Das, C.; Dasgupta, S.; Lynch, W.; Mekjian, A.; Tsang, M. The thermodynamic model for nuclear multifragmentation. Phys. Rep. 2005, 406, 1–47. [Google Scholar] [CrossRef] [Green Version]
  3. Désesquelles, P.; Pasquiers, S.; Blin-Simiand, N.; Magne, L.; Van-Oanh, N.T.; Thomas, S.; Domin, D. The statistical molecular fragmentation model compared to experimental plasma induced hydrocarbon decays. Phys. Chem. Chem. Phys. 2020, 22, 7586–7596. [Google Scholar] [CrossRef] [PubMed]
  4. Martinet, G.; Díaz-Tendero, S.; Chabot, M.; Wohrer, K.; Negra, S.D.; Mezdari, F.; Hamrita, H.; Désesquelles, P.; Padellec, A.L.; Gardés, D.; et al. Fragmentation of Highly Excited Small Neutral Carbon Clusters. Phys. Rev. Lett. 2004, 93, 063401. [Google Scholar] [CrossRef] [PubMed]
  5. Hervieux, P.A.; Zarour, B.; Hanssen, J.; Politis, M.F.; Martín, F. Fragmentation in collisions of Na9+clusters with Cs atoms. J. Phys. B At. Mol. Opt. Phys. 2001, 34, 3331–3355. [Google Scholar] [CrossRef]
  6. Gueron, S.; Levin, S.A. The dynamics of group formation. Math. Biosci. 1995, 128, 243–264. [Google Scholar] [CrossRef]
  7. Chen, Y.; Paul, G.; Cohen, R.; Havlin, S.; Borgatti, S.P.; Liljeros, F.; Stanley, H.E. Percolation theory applied to measures of fragmentation in social networks. Phys. Rev. E 2007, 75, 046107. [Google Scholar] [CrossRef] [Green Version]
  8. Allard, A.; Hébert-Dufresne, L.; Noël, P.A.; Marceau, V.; Dubé, L.J. Exact solution of bond percolation on small arbitrary graphs. EPL Europhys. Lett. 2012, 98, 16001. [Google Scholar] [CrossRef] [Green Version]
  9. Brown, W.K.; Karpp, R.R.; Grady, D.E. Fragmentation of the universe. Astrophys. Space Sci. 1983, 94, 401–412. [Google Scholar] [CrossRef]
  10. Melzak, Z.A. A scalar transport equation. Trans. Am. Math. Soc. 1957, 85, 552–560. [Google Scholar] [CrossRef]
  11. Ziff, R.M.; McGrady, E.D. The kinetics of cluster fragmentation and depolymerisation. J. Phys. A Math. Gen. 1985, 18, 3027. [Google Scholar] [CrossRef]
  12. Montroll, E.W.; Simha, R. Theory of Depolymerization of Long Chain Molecules. J. Chem. Phys. 1940, 8, 721–726. [Google Scholar] [CrossRef]
  13. Simha, R. Kinetics of Degradation and Size Distribution of Long Chain Polymers. J. Appl. Phys. 1941, 12, 569–578. [Google Scholar] [CrossRef]
  14. Simha, R. On the Degradation of Branched Chain Molecules. J. Chem. Phys. 1956, 24, 796–802. [Google Scholar] [CrossRef]
  15. Ziff, R.M.; McGrady, E.D. Kinetics of polymer degradation. Macromolecules 1986, 19, 2513–2519. [Google Scholar] [CrossRef]
  16. McGrady, E.D.; Ziff, R.M. “Shattering” Transition in Fragmentation. Phys. Rev. Lett. 1987, 58, 892–895. [Google Scholar] [CrossRef]
  17. Ziff, R.M. New solutions to the fragmentation equation. J. Phys. A Math. Gen. 1991, 24, 2821. [Google Scholar] [CrossRef] [Green Version]
  18. Krapivsky, P.L.; Ben-Naim, E. Scaling and multiscaling in models of fragmentation. Phys. Rev. E 1994, 50, 3502–3507. [Google Scholar] [CrossRef] [Green Version]
  19. Connaughton, C.; Dutta, A.; Rajesh, R.; Siddharth, N.; Zaboronski, O. Stationary mass distribution and nonlocality in models of coalescence and shattering. Phys. Rev. E 2018, 97, 022137. [Google Scholar] [CrossRef] [Green Version]
  20. Krapivsky, P.L.; Ben-Naim, E. Shattering transitions in collision-induced fragmentation. Phys. Rev. E 2003, 68, 021102. [Google Scholar] [CrossRef] [Green Version]
  21. Ernst, M.H.; Szamel, G. Fragmentation kinetics. J. Phys. A Math. Gen. 1993, 26, 6085. [Google Scholar] [CrossRef]
  22. Ziff, R.M. An explicit solution to a discrete fragmentation model. J. Phys. A Math. Gen. 1992, 25, 2569. [Google Scholar] [CrossRef] [Green Version]
  23. Bauer, W.; Dean, D.; Mosel, U.; Post, U. New approach to fragmentation reactions: The nuclear lattice model. Phys. Lett. B 1985, 150, 53–56. [Google Scholar] [CrossRef]
  24. Hendriks, E.M.; Spouge, J.L.; Eibl, M.; Schreckenberg, M. Exact solutions for random coagulation processes. Z. Für Phys. B Condens. Matter 1985, 58, 219–227. [Google Scholar] [CrossRef]
  25. Durrett, R.; Granovsky, B.; Gueron, S. The Equilibrium Behavior of Reversible Coagulation-Fragmentation Processes. J. Theor. Probab. 1999, 12, 447–474. [Google Scholar] [CrossRef]
  26. Berestycki, J. Exchangeable Fragmentation-Coalescence Processes and their Equilibrium Measures. Electron. J. Probab. 2004, 9, 770–824. [Google Scholar] [CrossRef]
  27. Berestycki, N.; Pitman, J. Gibbs Distributions for Random Partitions Generated by a Fragmentation Process. J. Stat. Phys. 2007, 127, 381–418. [Google Scholar] [CrossRef] [Green Version]
  28. Austin, L.; Shoji, K.; Bhatia, V.; Jindal, V.; Savage, K.; Klimpel, R. Some Results on the Description of Size Reduction as a Rate Process in Various Mills. Ind. Eng. Chem. Process Des. Dev. 1976, 15, 187–196. [Google Scholar] [CrossRef]
  29. Rozenblat, Y.; Grant, E.; Levy, A.; Kalman, H.; Tomas, J. Selection and breakage functions of particles under impact loads. Chem. Eng. Sci. 2012, 71, 56–66. [Google Scholar] [CrossRef]
  30. Abberger, T. Chapter 24 Population balance modelling of granulation. In Handbook of Powder Technology; Granulation, A.D., Salman, M.H., Seville, J., Eds.; Elsevier Science B.V.: Amsterdam, The Netherlands, 2007; Volume 11, pp. 1109–1186. [Google Scholar] [CrossRef]
  31. Bóna, M. A Walk Through Combinatorics—An Introduction to Enumeration and Graph Theory, 2nd ed.; World Scientific Publishing: Singapore, 2006. [Google Scholar]
  32. Evans, M.R.; Hanney, T. Nonequilibrium statistical mechanics of the zero-range process and related models. J. Phys. A Math. Gen. 2005, 38, R195–R240. [Google Scholar] [CrossRef] [Green Version]
  33. Feller, W. On the Time Distribution of So-Called Random Events. Phys. Rev. 1940, 57, 906–908. [Google Scholar] [CrossRef]
  34. Ruark, A. The Time Distribution of So-Called Random Events. Phys. Rev. 1939, 56, 1165–1167. [Google Scholar] [CrossRef]
  35. Matsoukas, T. Generalized Statistical Thermodynamics: Thermodynamics of Probability Distributions and Stochastic Processes; Springer International Publishing: Cham, Switzerland, 2019. [Google Scholar] [CrossRef]
  36. Matsoukas, T. Statistical Thermodynamics of Irreversible Aggregation: The Sol-Gel Transition. Sci. Rep. 2015, 5, 8855. [Google Scholar] [CrossRef] [PubMed]
  37. Matsoukas, T. Fundamentals of Chemical Engineering Thermodynamics; Pearson Education, Inc.: London, UK, 2013. [Google Scholar]
  38. Kelly, F.P. Reversibility and Stochastic Networks; Wiley: Cambridge, UK, 2011. [Google Scholar]
  39. Schrödinger, E. Statistical Thermodynamics; Dover Publications: Mineola, NY, USA, 1989; Reprint of the 2nd edition originally published by Cambridge University Press, 1952, under the subtitle A Course Seminar Lectures Delivered in January-March 1944 at the School of Theoretical Physics, Dublin Institute for Advanced Studies. [Google Scholar]
  40. Matsoukas, T. Thermodynamics Beyond Molecules: Statistical Thermodynamics of Probability Distributions. Entropy 2019, 21, 890. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Random fragmentation of integer mass M into N pieces is equivalent to breaking a string with M beads at N 1 random points. With M = 10 , N = 3 the number of possible partitions is 9 2 = 36 . If the mass is made up of two colors, every permutation of the beads is equally probable; with M A = 6 , M B = 4 the number of partitions increases by the factor 6 + 4 4 = 210 and the total number of permutations is 7560.
Figure 1. Random fragmentation of integer mass M into N pieces is equivalent to breaking a string with M beads at N 1 random points. With M = 10 , N = 3 the number of possible partitions is 9 2 = 36 . If the mass is made up of two colors, every permutation of the beads is equally probable; with M A = 6 , M B = 4 the number of partitions increases by the factor 6 + 4 4 = 210 and the total number of permutations is 7560.
Condensedmatter 05 00064 g001
Figure 2. The compositional distribution c a | k ¯ in particles of mass k = 2 , 3 , and 4. The parent particle contains M A = 4 units of A, M B = 3 units of B, and breaks into N = 4 pieces. Lines are from Equation (37) and points are from MC simulation after 20,000 fragmentation events. The excellent agreement between MC and theory demonstrates the exact nature of Equation (37) and validates the MC method. The binomial distribution is only asymptotically valid and in this case is a poor approximation because the size of the fragments is small.
Figure 2. The compositional distribution c a | k ¯ in particles of mass k = 2 , 3 , and 4. The parent particle contains M A = 4 units of A, M B = 3 units of B, and breaks into N = 4 pieces. Lines are from Equation (37) and points are from MC simulation after 20,000 fragmentation events. The excellent agreement between MC and theory demonstrates the exact nature of Equation (37) and validates the MC method. The binomial distribution is only asymptotically valid and in this case is a poor approximation because the size of the fragments is small.
Condensedmatter 05 00064 g002
Figure 3. Size and compositional distributions at fragment sizes k = 2 , 4 , and 8 for two bias functionals: (ad): w a , b = ( 1 + a ) 4 ( 1 + b ) 4 ; (eh): w a , b = ( 1 + a ) 4 + ( 1 + b ) 4 . In both cases, the particle contains M A = 20 units of A, M B = 20 units of B, and breaks into N = 4 fragments.
Figure 3. Size and compositional distributions at fragment sizes k = 2 , 4 , and 8 for two bias functionals: (ad): w a , b = ( 1 + a ) 4 ( 1 + b ) 4 ; (eh): w a , b = ( 1 + a ) 4 + ( 1 + b ) 4 . In both cases, the particle contains M A = 20 units of A, M B = 20 units of B, and breaks into N = 4 fragments.
Condensedmatter 05 00064 g003
Table 1. Closed form results for three composition-independent bias functionals.
Table 1. Closed form results for three composition-independent bias functionals.
w a , k a Ω M ; N ( 1 )
Case 11 M 1 N 1
Case 2 k k 1 k ! m M N N ! M ! M 1 N 1
Case 3 2 ( 2 k ) k 2 k ! m M N N ! M ! 2 M 1 N 1
asymptotically for M , N 1 , M / N < 2 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Matsoukas, T. Statistical Mechanics of Discrete Multicomponent Fragmentation. Condens. Matter 2020, 5, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat5040064

AMA Style

Matsoukas T. Statistical Mechanics of Discrete Multicomponent Fragmentation. Condensed Matter. 2020; 5(4):64. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat5040064

Chicago/Turabian Style

Matsoukas, Themis. 2020. "Statistical Mechanics of Discrete Multicomponent Fragmentation" Condensed Matter 5, no. 4: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat5040064

Article Metrics

Back to TopTop