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.
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 [1]. For a mass fractal it is given by the scattering exponent of the power-law SAS intensity where . For deterministic fractals, additional information can be obtained such as scaling factor (from the period of in the logarithmic scale), the number of fractal iteration (which equals the number of periods of function ), and the total number 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 (; initiator) and in the first iteration () the number of cubes which are kept is related to the probability in which the initiator will be fragmented into 27 cubes of edge length . Figure 1 illustrates the construction of a generic model of fragmentation, where the sizes of remaining cubes at th iteration are given byAt th iteration the total number of particles will be given byand in a good approximation ( and ) it can be shown that the fractal dimension can be written as [5]with , and thus . Figure 2 shows the variation of the fractal dimension with probability in which the initiator will be fragmented into 27 cubes of edge length .
3. Theory
Here, we restrict ourselves to two-phase systems, which are composed of homogeneous units of mass density . The units are immersed into a solid matrix of pore density . Then [8, 9] we can consider the system as if the units were frozen in a vacuum and had the density . The density is called scattering contrast, and thus the scattering intensity is given bywhere is the fractal concentration, is the volume of each fractal, and is the normalized form factorobeying . The brackets stand for the ensemble averaging over all orientations of . Once a deterministic fractal is composed of the same objects, say, cubes of edge length , thenwith being the Fourier component of the density of the cubes, and are the center-of-mass positions of cubes. Here, the cube form factor of unit edge length is given by [8]Therefore, the scattering intensity becomeswith andbeing 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 at th iteration can be written as [26]with and 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]
4. Results and Discussion
For well-known systems such as Cantor sets or Vicsek fractals 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 of is shown in Figure 3. It is characterized by the presence of three main regions on a double logarithmic scale: at low () a plateau (Guinier region) which gives information about the overall size of the fractal, a fractal (generalized power law-decay) region at , and an asymptotic regionat high () which gives information about the number of fragments in the fractal (see also Figure 4).
The value of the fractal dimension can be obtained by considering the polydispersity in size of the fractals. Here, we use a distribution function of sizes in such a way that gives the probability of finding a fractal whose size falls within . We consider a log-normal distributionwhere . The mean length and its relative variance are and where . Thus, the polydisperse structure factor is given by [26]where is the total volume at the th iteration.
Figure 4 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 the structure factor is proportional to , 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 from the log-periodicity of the quantity versus in the fractal region. For the former model, the choice of the scattering units (here, cubes) according to a probability 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 . However, the asymptote of the fractal structure factor gives the number of particles inside the fractal and this can be used in (2), together with probability from (3) (with being known from the exponent of the power-law decay), to determine the iteration number of the fractal. Thus, by knowing , , and and using the relation [26] , 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 () 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 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.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The author acknowledges financial support from JINR-IFIN-HH projects awarded by the Romanian Plenipotentiary Representative at JINR.