# Microscale Fragmentation and Small-Angle Scattering from Mass Fractals

Authors: E. M. Anitas
Journal: Advances in Condensed Matter Physics
 (2015)

---

## Abstract

Using the small-angle scattering method, we calculate here the mono- and polydisperse structure factor from an idealized fragmentation model based on the concept of renormalization. The system
consists of a large number of fractal microobjects which are randomly oriented and whose positions are uncorrelated. It is shown that, in the fractal region, the monodisperse form factor is characterized
by a generalized power-law decay (i.e., a succession of maxima and minima superimposed on a simple power-law decay) and whose scattering exponent coincides with the fractal dimension of the scatterer.
The present analysis of the scattering structure factor allows us to obtain the number of fragments resulted at a given iteration. The results could be used to obtain additional structural information about systems obtained through microscale fragmentation processes.


---

## Body

## 1. Introduction

The number of fragments (such as those resulted from rocks weathering and explosions or produced by earth’s crust) as a function of their sizes over a wide range of scales usually can be described by a fractal distribution [1–3], which is responsible for various physical properties, such as hydraulic conductivity or moisture characteristics in soils [4]. Although a quantification of these processes using the renormalization group approach has been suggested in [5, 6], an important issue concerns the distribution of fragments at microscales obtained by various methods such us in ultrasonic fragmentation using optoacoustic lens [7].Small-angle scattering (SAS; neutrons, X-ray, light) [8, 9] is one of the most important techniques for investigating the microstructure of various types of systems which addresses the issue of size distribution, including the smallest and largest components. It yields the differential elastic cross section per unit solid angle as a function of the momentum transfer. The main advantage consists in its ability to differentiate between mass and surface fractals [10, 11], and it has been successfully used in studying the property of self-similarity across nano- and microscales [2], such as various types of membranes [12–15], cements [16], semiconductors [17], magnetic structures [18, 19], or biological structures [20–22]. Thus the concept of fractal geometry coupled with SAS technique can give new insights concerning the structural characteristics of such complex systems [10, 11, 23–30]. One of the main parameters which can be obtained is the fractal dimension D
m [1]. For a mass fractal it is given by the scattering exponent of the power-law SAS intensity I
(
q
)
∝
q
-
D
m where 0
<
D
m
<
3. For deterministic fractals, additional information can be obtained such as scaling factor (from the period of I
(
q
)
q
D
m in the logarithmic scale), the number of fractal iteration m (which equals the number of periods of function I
(
q
)
∝
q
-
D
m), and the total number N
m of structural units of which the fractal is composed.In this paper, we develop a theoretical model based on renormalization group approach which could describe various microscale fragmentation processes. We calculate analytically the mono- and polydisperse structure factor and show how to obtain the main structural parameter such as the fractal dimension and the number of fragments at a given iteration and show how to estimate the smallest and largest radius of the fractal from SAS data.


## 2. Fragmentation Model: Construction and Properties

The construction process of three-dimensional mass fractals is similar to that of mass generalized Cantor and Vicsek fractals (GCF and GVF) [26] in the sense that one follows a top-down approach in which an initial structure is repeatedly divided (by a single scaling factor) into a set of smaller structures of the same type according to a given rule which is changed randomly from one iteration to the next one. Thus, we first consider a cube of edge length l
0 (m
=
0; initiator) and in the first iteration (m
=
1) the number of cubes which are kept is related to the probability f in which the initiator will be fragmented into 27 cubes of edge length l
1
=
l
0
/
3. Figure 1 illustrates the construction of a generic model of fragmentation, where the sizes of remaining cubes at mth iteration are given by(1)
l
m
=
l
0
3
m
.At mth iteration the total number of particles will be given by(2)
N
m
=
1
-
f
1
+
27
f
+
27
f
2
+
⋯
+
27
f
m
,and in a good approximation (m
≫
1 and 27
f
>
1) it can be shown that the fractal dimension can be written as [5](3)
D
m
=
3
log
⁡
27
f
log
⁡
27
,with 1
/
27
<
f
<
1, and thus 0
<
D
m
<
3. Figure 2 shows the variation of the fractal dimension D
m with probability f in which the initiator will be fragmented into 27 cubes of edge length l
1
=
l
0
/
3.Figure 1
A two-dimensional projection of the fragmentation model. Note that at a given iteration all the fragments have the same size. The probability that the cube with edge lengthl
0 (m
=
0; initiator) will be fragmented into 27 cubes of edge length l
1 is f.Figure 2
Dependence of the fractal dimensionD
m of the probability f from (3).


## 3. Theory

Here, we restrict ourselves to two-phase systems, which are composed of homogeneous units of mass densityρ
m. The units are immersed into a solid matrix of pore density ρ
p. Then [8, 9] we can consider the system as if the units were frozen in a vacuum and had the density Δ
ρ
=
ρ
m
-
ρ
p. The density Δ
ρ is called scattering contrast, and thus the scattering intensity is given by(4)
I
q
=
n
Δ
ρ
2
V
2
F
q
2
,where n is the fractal concentration, V is the volume of each fractal, and F
(
q
) is the normalized form factor(5)
F
q
=
1
V
∫
V
e
-
i
q
·
r
d
r
,obeying F
(
0
)
=
1. The brackets ⋯ stand for the ensemble averaging over all orientations of q. Once a deterministic fractal is composed of the same objects, say, N cubes of edge length l, then(6)
F
q
=
ρ
q
F
0
q
l
N
,with ρ
q
=
∑
j
‍
e
-
i
q
·
r
j being the Fourier component of the density of the cubes, and r
j are the center-of-mass positions of cubes. Here, the cube form factor of unit edge length is given by [8](7)
F
0
t
=
sin
⁡
t
x
/
2
t
x
/
2
sin
⁡
t
y
/
2
t
y
/
2
sin
⁡
t
z
/
2
t
z
/
2
.Therefore, the scattering intensity becomes(8)
I
q
=
I
0
S
q
F
0
q
l
2
N
,with I
(
0
)
=
n
|
Δ
ρ
|
2
V
2 and(9)
S
q
≡
ρ
q
ρ
-
q
Nbeing the structure factor, and it is intimately connected to the pair distribution function. In the fractal region the form factor tends to one, and this implies that the scattering intensity is proportional to the structure factor in the fractal region.The monodisperse form factor atmth iteration can be written as [26](10)
F
m
q
=
F
0
q
l
m
∏
i
=
1
m
G
i
q
,with m
=
1,2
,
… and G
(
q
) being the generative function of the fractal and it specifies the positions of the scattering cubes inside the fractal. Thus, the scattering intensity becomes [26] (11)
I
m
q
I
m
0
=
F
m
q
2
.


## 4. Results and Discussion

For well-known systems such as Cantor sets or Vicsek fractalsG
i
(
q
) has known analytical expressions [26]. Here, we use a random distribution of cubes at each iteration for writing the generative function.The monodisperse structure factor for various values off is shown in Figure 3. It is characterized by the presence of three main regions on a double logarithmic scale: at low q (q
≲
1
/
l
0) a plateau (Guinier region) which gives information about the overall size of the fractal, a fractal (generalized power law-decay) region at 1
/
l
0
≲
q
≲
1
/
l
m, and an asymptotic region(12)
S
q
≃
1
N
m
,at high q (1
/
l
m
≲
q) which gives information about the number of fragments in the fractal (see also Figure 4).Figure 3
The monodisperse structure factor for the fragmentation model at various values of probabilityf shows a superposition of maxima and minima on a simple power-law decay in the fractal region. The curves are shifted down vertically for clarity with a factor of 1
0
-
2 (red; middle curve) and, respectively, with 1
0
-
4 (green; lower curve). Iteration number is m
=
6.Figure 4
The polydisperse structure factor for the fragmentation model at various values of probabilityf shows a simple power-law decay in the fractal region with the scattering exponent equal to the fractal dimension of the fractal. The horizontal lines indicate the asymptotic values of the structure factor (see (12)). Iteration number is m
=
6 and relative variance σ
r
=
0.6.The value of the fractal dimension can be obtained by considering the polydispersity in size of the fractals. Here, we use a distribution functionD
N
(
l
) of sizes in such a way that D
N
(
l
)
d
l gives the probability of finding a fractal whose size falls within (
l
,
l
+
d
l
). We consider a log-normal distribution(13)
D
N
l
=
1
σ
l
2
π
1
/
2
exp
⁡
-
ln
⁡
l
/
l
0
+
σ
2
/
2
2
2
σ
2
,where σ
=
[
l
n
(
1
+
σ
r
2
)
]
1
/
2. The mean length and its relative variance are μ
≡
l
D and σ
r
≡
(
l
D
2
-
l
0
2
)
1
/
2
/
μ where ⋯
D
=
∫
0
∞
‍
⋯
D
N
(
l
)
d
l. Thus, the polydisperse structure factor is given by [26](14)
I
m
q
=
n
Δ
ρ
2
∫
0
∞
F
m
q
2
V
m
2
D
N
l
d
l
,where V
m is the total volume at the mth iteration.Figure4 shows the polydisperse structure obtained from Figure 3 together with (14). As expected, the influence of polydispersity is to smooth the scattering curves, and to preserve the scattering exponent in the fractal region. At high q the structure factor is proportional to 1
/
N
m, according to (12). Note that the value of fractal dimension is in agreement with (3).The results show that in an ideal SAS experiment (which shows the presence of Guinier, fractal, and Porod regions) the overall size and the smallest size in the fractal can be obtained from both the mono- and polydisperse intensity, through the position of the crossover between Guinier and fractal region, and, respectively, between fractal and Porod region. However, a reliable estimation of the fractal dimension can be obtained only in the polydisperse case.An important difference between the model developed here and those based on exact self-similar fractals is that for the latter models we have a periodic distribution for the minima of monodisperse intensity [26, 29, 31], which allows us to obtain the number of fractal iteration, as well as the fractal scaling factor β
s from the log-periodicity of the quantity I
(
q
)
q
D
m versus q in the fractal region. For the former model, the choice of the scattering units (here, cubes) according to a probability f leads to an asymmetric distribution with respect to the origin of the initiator. This difference is reflected in the polydisperse intensity by oscillations of various amplitudes (see Figure 4) around the simple power-law decay I
(
q
)
∝
q
-
D
m. However, the asymptote of the fractal structure factor gives the number of particles N
m inside the fractal and this can be used in (2), together with probability f from (3) (with D
m being known from the exponent of the power-law decay), to determine the iteration number m of the fractal. Thus, by knowing m, N
m, and D
m and using the relation [26] N
m
=
(
1
/
β
s
)
m
D
m, we can then still recover the scaling factor.Depending on the type and availability of the sample and on the experimental setup, various approaches could be followed in order to extract structural information. Thus, if the structure factor given by (12) is multiplied by a form factor corresponding to the basic scattering unit, then one obtains the SAS intensity in which the asymptotic region of the structure factor is replaced by a Porod decay (I
(
q
)
∝
q
-
4) and thus additional information can be obtained. For example, in [32], from the calculated small-angle X-ray scattering intensity of a fragmentation, fractal silica nanoagglomerates, the fractal structure, and aggregation number have been obtained by using the well-known Beaucage model [33], and the specific surface area has been obtained with the help of the integral parameter Q and the Porod constant [9].


## 5. Conclusions

We have developed an idealized statistical self-similar fractal model for describing fragmentation processes based on 3D fractal and have used the SAS technique to characterize the microstructural properties of such fragments. The model is based on the concept of renormalization, and the suggested approach has important advantages since SAS is a nondestructive technique, the obtained information is averaged over a macroscopic volume, and the model features allow us to obtain additional information about the fractal fragments, such as the number of structural units of which the fractal is composed, the scaling factor or the fractal iteration number.We have calculated the monodisperse and polydisperse structure factor and have shown how to obtain the smallest and largest sizes inside the fractal, from the position of the crossover between Guinier and fractal region, and, respectively, between fractal and Porod region. We have shown that the slope of SAS intensity is in agreement with the fractal dimension of the fractal, as given by (3).The model can be used in understanding various physical properties such as hydraulic conductivity or moisture characteristics for systems obtained through fragmentation and it can be easily extended to describe more complex structures, such as by keeping fragments of different sizes obtained at various iterations (surface fractals), by taking into account multiple scaling factors at a given iteration (multifractals), changing the scaling factors with the iteration number according to a certain rule (fat fractals) or a combination of them.



---

*Source: 501281-2015-08-05.xml*