Evaluation of Changes in Structure of Modiﬁed Cement Composite Using Fractal Analysis

: The need for research using modern methods of physical and chemical analysis in combination with methods for processing the graphical results arises with an increase of developments aimed at changing the structure and properties of hydration-hardening materials. Currently, the method of scanning electron microscopy is widely used for the qualitative and quantitative analysis of materials. Additional information, including fractal analysis of the image, determines the information-entropy levels and the state of the information system in the form of the microstructure based on the generalized structural entropy that increases the volume and quality of analytical information extracted from the obtained images. Images with a magniﬁcation of 20,000 times were used to analyze the microstructure of a control cement composite and a carbon black modiﬁed composite. The fractal analysis of the images revealed a change in the relative organization of the system, which justiﬁedthe nature of the change in the physical and technical properties of the modiﬁed cement composite relative to the additive-free cement composite from the perspective of self-organization.


Introduction
Currently, a significant number of research and industrial developments have been accumulated on the introduction of nanotechnology in the production of building materials-in particular, the use of nano-structuring additives for composites based on hydrationhardening binders [1][2][3]. However, the number of studies attempting to determine the main factor that leads to significant changes in the composite structure is insufficient [4,5] compared with the number of applied projects aimed at obtaining the optimal material compositions [6][7][8][9]. Over the past decade, research and testing has been carried out by the authors on the modification of cement and gypsum composites with aqueous suspensions of multilayer carbon nanotubes [10][11][12][13]. The nature of the effect of these additives allows confirmation of the non-additive nature of the effect of the active components of these additives. A number of studies have been performed to determine the leading role of additives in changing the characteristics of composites as a result of the complex impact of the various components. At the same time, it was proved that the separate addition of the components of a complex additive does not have such an effect on the final properties of the material [4,14,15]. The synergistic effect occurs due to the complex impact of the surface-active substance (superplasticizer) with carbon nanotubes on the structure of hardening composite [10,11]. However, modern methods of materials research do not allow the phenomenon of a synergistic effect to be proved by direct methods; mainly, the phenomenon is confirmed by a combination of a complex of indirect methods [16,17].
Systems of diverse nature are considered in synergetics as a complex of subsystems that perform certain technological functions. Subsystems are interconnected by processes of intensive dynamic interaction and exchange of energy, matter, and information [18][19][20]. Self-organization of dynamic systems is a complex mechanism of self-management in natural and artificial systems in which the corresponding properties of the systems are formed without external influence. The effects of self-organization in natural systems arise as a result of changes in the "control parameters", while the effects of self-organization in artificial systems are formed as a result of the construction of feedback. Knowledge of the laws of the behavior of synergetic systems allows one to think deeply about the creation of artificial systems in which the processes of purposeful self-organization form the desired structures [21][22][23].
Materials based on hydration binders that are modified with amulticomponent additive containing carbon nanotubes can serve as an example of artificial systems capable of self-directed organization. To date, the optimal ratios of components in anadditive composition have been stated, leading to a significant enhancementof the physical and mechanical characteristics of matrices based on hydration-hardening binders [24][25][26]. At the same time, significant changes in structure of the material are noted by many researchers [27][28][29].
The complex of modern methods of physical and chemical studies indirectly allows stating the causes and mechanisms that justify a significant improvement of composite properties [30][31][32][33]. Studies of physical and mechanical properties confirm the non-additive nature of the effect of multicomponent additives [34][35][36]. Modern physical and chemical studies including scanning electron microscopy with X-ray microanalysis, X-ray phase analysis, IR spectral analysis, computed tomography, and differential thermal analysis are able to determine the changes occurring in the composite material at the level at which hydration products are formed [37][38][39]. However, the information extracted from the analysis of the material microstructure is often used by researchers only to determine the morphology and size of hydration products [40,41]. This does not allow full use of the information of the bitmap image [42][43][44].
The determinationof links among the elements of system andthe causality of movement of each element when interacting with other elements is the main property of information. The existence of a certain relationship between the production of entropy in the system and the exchange of entropy with the external environment is a prerequisite for structure formation in open systems. The effect of information, in contrast to entropy, is expressed in a tendency to increase orderliness and reduce uncertainty. Currently, the concept of entropy as a fundamental physical quantity is expanding and deepening in connection with the ability to evaluate information through its analysis. The analysis of the entropy characteristics of amicrostructure of concrete and mortar allows evaluation ofthe changes in itsintegral characteristics, such as the electrical resistance and the configuration entropy of the structural elements. However, physical entropy and information entropy may not coincide, and therefore it is necessary to introduce abstract concepts such as "entropy force", "information capacity", and"information potential" in order to remain within the framework of thermodynamic concepts when working with information.
Information transformations are the basis for thefractal analysis of structures of various natures. On this basis, the methods and computer programs for multifractal parametrization of flat images have been developed to solve problems for optimizing the fractal structures of a material in technological processes [45][46][47] as well as to study their stability and adaptive properties [48][49][50].
There is research showing that fractal analysis is used to study the structure of materials for construction purposes [51][52][53]. Thus, onpaper [54], the statistical and geometric characteristics of microstructure of the pore space of rock samples were studied on the basis of a3D micromodel. In the study [55], the variability (statistical stability) of the fractal characteristics (fractal dimension D and lacunarity Λ of the microstructure) from the action of various nanomodifiers was evaluated using the example of nano-structured concrete. It was stated that the value of the fractal dimension was more invariant when introducing the modifiers than the lacunarity. It was noted that the lacunarity significantly depends on the magnification at which the microstructure was obtained.
Analyzing the studies in this area, it can be concluded that fractal analysis is almost not used to assess the degree of self-organization of hardening cement systems. Thus, the aim of our research was to evaluate the changes in the structure of the modified cement composite using fractal analysis.

Samples Preparation
Sample-beams with dimensions of 4 × 4 × 16 cm were made to study the features of composite structure formation. Ordinary Portland cement of CEM I 42.5 N strength class was used as a binder. The chemical and mineral composition of cement wasrepresented by the following minerals as a percentage: tricalcium silicate C 3 S-64.6, bicalcium silicate C 2 S-10.7, tricalcium aluminate C 3 A-7.0, andtricalciumaluminoferrite C 4 AF-14.7, MgO-1.4.
Cement-sand mortars (mortar modified with an additive and the control mortar) were made to determine theirmechanical characteristics. Soot was used as a modifying additive. It acted as a control parameter for changing the functional quality of the information system of the mortar microstructure. Soot with aparticle size less than 80nm and with aspecific surface area of 65-100 m 2 /g was used. The amount of carbon black used to modify the samples was dosed as a percentage of the binder weight. A change in the configuration entropy of the system is expected when using technical soot as a modifier.
The water-to-binder ratio was 0.52 for all compositions. The samples were stored in wet conditions for 28 days prior to the test.

Testing of Samples
Tests for determiningthe tensile strength in bending and compressive strength of samples were carried out using aPGM-100 MG4-A hydraulic press with amaximum load of 100 kN and a loading speed of 0.5 MPa/s. Samples of control and modified compositions were selected for microstructure study. Preliminary deposition of a gold-based conductive film was applied to samples. The scanning electron micrographs were obtained using aMIRA 3 TUSCAN microscope at the AdMaS Research Center located at the Technical University of Brno. A statistically significant number of images were taken in the study, and the most characteristic results were analyzed and are presented. Differential thermogravimetry (DTG) and differential scanning calorimetry (DSC) methods were used. Laboratory studies were conducted on a TGA/DSC1 thermal analyzer produced by Mettler-Toledo Vostok. Shooting conditions: measurement interval 50-1100 • C, heating rate 10 deg/min, platinum crucibles, working medium-air. Obtained results were discussed in previously published data [56].

Calculation of Multifractal Parameters
Multifractal analysis of concrete surface images was carried out in accordance with the methodology described in papers [9,10] and includedthe following steps:

-
Obtaining the digital image of the concrete surface under study using an electron microscope according to the scheme in Figure 1; -Selecting the sections of scanner images of the microstructure of concrete samples; -Calculating the multifractal characteristics for the scanning area; -Studying the obtained multifractal characteristics. The calculation of the multifractal parameters was performed on the basis of parametrization of electron microscopic images of microstructures using a computer program [57]. For this purpose, an equallycellular square partition of the studied image structure wasused. A measure in the form of a positive number was mapped to each cell. It characterized the area of the structure in thecell.
The processing principle wasas follows. Images of fractal objects weredivided into square cells with aside δ that variedfrom 3 to 50 pixels, digitizing the image size of r × c, The calculation of the multifractal parameters was performed on the basis of parametrization of electron microscopic images of microstructures using a computer program [57]. For this purpose, an equallycellular square partition of the studied image structure wasused. A measure in the form of a positive number was mapped to each cell. It characterized the area of the structure in thecell.
The processing principle wasas follows. Images of fractal objects weredivided into square cells with aside δ that variedfrom 3 to 50 pixels, digitizing the image size of r × c, where r is the number of rows and c is the number of columns. The minimum cell is the pixel of the image x ij , i = 1, r and j = 1, c, which characterizes the grayscale from 0 to 255, where 0 means black and 255 means white. The technique of multifractal parametrization of the image is known and it uses the calculation of moments p q k -the probability of distribution of points of the set in the cells [46,54]; the developed program provides the choice of the gray scale interval when digitizing the image. The number of single pixels is calculated in each "non-empty" cell by formula: where r(k) and c(k) are, respectively, the number of the row and column of the pixel with which the k-th cell begins.
The number of single pixels of the image is calculated by the formula: The "population" of the k-th cell can be determined: The normalization condition must be met: The moment p q k -th order determines the value of D q , which corresponds to the degree of sparsity of the multifractal multitude: where q ∈ (−∞, ∞). In the literature, the minimum and maximum sparsity are limited to the consideration of D +40 and D −40 .
For integer values q(−40, −39, . . . , 39, 40) the discrete values D q,δ depending on the cell size of δ are: The proposed method differs from the binary image because the image is analyzed for the grayscale interval (the optical density region can be set arbitrarily).
The multifractal parametrization of the material structure image describes its statistical properties by the spectrum of generalized dimensions of Renee of D q . For q = 1, we get the information dimension of D 1 (the Shannon index). The value of K s = D 1 − D 40 is used as an order measure. The quantitative measure of chaos in the system under study is taken as The states of static stability of the structure, the so-called maxima of the information dimension, are stated from the graph of dependence of the information dimension of D 1 on the scale of d of structure elements. States with the maximum value of the information dimension of D 1i = const can be considered as a memory matrix written as the static stability of the active elements of structure, where D 1i is the information dimension of the i-th level. The generalized information entropy is equal to the sum of the information dimensions of the stationary levels of H = ∑D 1i , and it characterizes the static state of the structural entropy.
Information transformations of the image are accompanied by configuration changes of the active elements of structure, causing the increase of the measure of ∆D randomness. The states of the dynamic stability of structure-namely, the maxima of ∆D* and the information dimension of D1*-are stated from the graph of dependence of ∆D on the scale of d of the structure elements. The information entropy of H* = ∑D* 1i characterizes the dynamic state of the structural entropy, where D* 1i is the information dimension of the i-th level of randomness.
An increase of the system complexity occurs due to an increase of the number of heterogeneous elements and connections between them, creating the N stable stationary states with the set of probabilities of p 1 , p 2 , . . . , p N . The amount of information corresponding to this state of the system is expressed by the K.Shannon's formula H = −∑ N j P j log 2 P j . The number of N determines the information capacity-namely, the maximum amount of information of H max that the system contains-and it is determined under the equally probable condition of P = 1/N. The maximum structural entropy of the information system of H max = −log 2 P = log 2 N corresponds to the equilibrium state of physical chaos in the system, and it is measured in bits.
The difference of H max − H * and R = 1 − H * /H max express the absolute and relative organization of the system. Relative randomness of the system is S = H * /H max . The structural entropy of the dynamic state of the information system as H* = H max corresponds to the state of physical chaos, and the loss of the functional quality of structure, herewith the relative organization, is R = 0.

Results and Discussion
Studies of physical and mechanical properties of the control (additive-free) and modified compositions show that the compressive strength of the control sample was 20.8 MPa, and the tensile strength in bending was 4.26 MPa. The compressive strength of the modified sample with 0.02% of soot was 26.9 MPa; the tensile strength in bending was 5.2 MPa.
The results of X-ray phase analysis and differential thermal analysis show that the number of changes in the composition of the modified sample wasinsignificant, despite a significant increase of its strength characteristics. It can be assumed that the effect of the additive on the conditions of structure formation does not change the composition of hydration products. The phenomenon of the synergistic effect is an example of nonlinear phenomena; in the analysis of the synergistic effect, forecasting is primary, since the variance of the results becomes comparable with the mathematical expectation, which affects the requirement for the accuracy of the results.
The results of scanning electron microscopy show changes in the shape and size of the structural elements ( Figure 2). additive on the conditions of structure formation does not change the composition of hydration products. The phenomenon of the synergistic effect is an example of nonlinear phenomena; in the analysis of the synergistic effect, forecasting is primary, since the variance of the results becomes comparable with the mathematical expectation, which affects the requirement for the accuracy of the results.
The results of scanning electron microscopy show changes in the shape and size of the structural elements ( Figure 2). Morphological differences in structure, differences in phase and chemical composition, and local differences in physical and mechanical properties by volume were observed in the cement matrix. Modification of matrices of hydraulic and non-hydraulic binders by dispersed additives also ledto significant changes in the morphology of the hydration products [10][11][12].The transformations of substances and the transformation of energy occurredin separate elementary volumes, which were interconnected by the diffusion transfer of atoms and energy. This supports the non-equilibrium state of the open system and the processes of self-organization accompanied by changes in the thermodynamic variables of the system.
The concrete structure was heterogeneous due to the technological features and features of the structure formation. The concrete structure was characterized by a chaotic distribution of defects. Hardening was characterized by simultaneous ordering processes-the crystal grewover time and degradation processes-the formation and change of the pore space that ledto a non-equilibrium state of the system.
Thus, the influence of additives on the early hydration of the hardening cement matrix that ledto the increase of thermodynamic fluctuations of hereditary disorder in a chaotic system wasone of the possible mechanisms of the influence of additives on the processes of structure formation. Hence, it was possible to manifest the self-organization by acquiring its own nonlinear dynamics, which reducedthe structural entropy and increasedthe adaptation. The order processed and formed local regions with short-range order developing as a structuring information principle long before the appearance of the center of formation of a new phase, as a process of production and materialization of in- Morphological differences in structure, differences in phase and chemical composition, and local differences in physical and mechanical properties by volume were observed in the cement matrix. Modification of matrices of hydraulic and non-hydraulic binders by dispersed additives also ledto significant changes in the morphology of the hydration products [10][11][12].The transformations of substances and the transformation of energy occurredin separate elementary volumes, which were interconnected by the diffusion transfer of atoms and energy. This supports the non-equilibrium state of the open system and the processes of self-organization accompanied by changes in the thermodynamic variables of the system.
The concrete structure was heterogeneous due to the technological features and features of the structure formation. The concrete structure was characterized by a chaotic distribution of defects. Hardening was characterized by simultaneous ordering processes-the crystal grewover time and degradation processes-the formation and change of the pore space that ledto a non-equilibrium state of the system.
Thus, the influence of additives on the early hydration of the hardening cement matrix that ledto the increase of thermodynamic fluctuations of hereditary disorder in a chaotic system wasone of the possible mechanisms of the influence of additives on the processes of structure formation. Hence, it was possible to manifest the self-organization by acquiring its own nonlinear dynamics, which reducedthe structural entropy and increasedthe adaptation. The order processed and formed local regions with short-range order developing as a structuring information principle long before the appearance of the center of formation of a new phase, as a process of production and materialization of information in the form of a specific topological scheme (nature, shape, size of additive particles, and their aggregation into clusters with strong associative connections).
The results of the multifractal analysis of the electron microscopic image of microstructure of the control sample are shown in Figure 3. The entropy parameters of the static and dynamic stabilities of the microstructure are given in Table 1.
formation in the form of a specific topological scheme (nature, shape, size of additive particles, and their aggregation into clusters with strong associative connections).
The results of the multifractal analysis of the electron microscopic image of microstructure of the control sample are shown in Figure 3. The entropy parameters of the static and dynamic stabilities of the microstructure are given in Table 1.    The deterministic function of the information system of the image of the control sample microstructures has the relative randomness S = 0.55 − 0.69 and the relative organization R = 0.31 − 0.45. One can note the value of 70% of the renewal of active elements of the microstructure in the spectrum of levels of the dynamic state of microstructures of the AB and BB cells. The entropy of the equilibrium levels of the AB and BB cells is D1 eq = 0.87534 and D1 eq = 0.8805, respectively. For the BA cell (Table 1), the update is 50%, and the structural entropies of the static and dynamic states are close in magnitude. The invariant variety is observed in the change ∆D (Figure 3) in the form of increasing the resonant maxima in a wide range of changes in the scale of d. The relative organization of the functional state of the information system of the BA microstructure is R = 0.31. The structural entropy of the static state exceeds the entropy of the dynamic state in the information system of the AA microstructure, so the information capacity of K = exp(H) = 4.7 × 10 5 is higher than the capacity of K* = exp(H*) = 3.4 × 10 4 .

AABA
The high potential of choice (information capacity) of the static state forms a gradient force of influence on the dynamic state, probably providing the increase of the relative organization of R = 0.45 of the microstructure of the AA cell. A single information maximum is observed in the near-order change of ∆D in the scale change interval of d. It can be assumed that the development of the crystal intergrowth is provided due to the high potential of choice (information capacity) of the static state of the information system. This confirms the nonequilibrium of states due to the presence of gradient-phase components in the microstructure of the control sample.
The results for determining the entropy parameters of structure of the control additivefree sample are presented in Table 1, where H is the general structural information entropy, H max is the maximum information, H* is the configuration entropy, R is the value of relative organization of the system, D1equil is the value of the information dimension, K is the information capacity, and ∆D is the measure of randomness.
The results of the multifractal analysis of the electron microscopic image of the microstructure of the modified sample are shown in Figure 4. The entropy parameters of the static and dynamic stability of microstructure are given in Table 2. The deterministic function of the information system of the image of the co sample microstructures has the relative randomness S = 0.55 − 0.69 and the relativ ganization R = 0.31 − 0.45. One can note the value of 70% of the renewal of active elem of the microstructure in the spectrum of levels of the dynamic state of microstructur the AB and BB cells. The entropy of the equilibrium levels of the AB and BB cells is =0.87534 and D1eq =0.8805, respectively. For the BA cell (Table.1), the update is 50%, the structural entropies of the static and dynamic states are close in magnitude. Th variant variety is observed in the change ∆D (Figure 3) in the form of increasing the onant maxima in a wide range of changes in the scale of d. The relative organizatio the functional state of the information system of the BA microstructure is R = 0.31. structural entropy of the static state exceeds the entropy of the dynamic state in th formation system of the AA microstructure, so the information capacity of K = exp( 4.7· 10 5 is higher than the capacity of K* = exp(H*) = 3.4·10 4 .
The high potential of choice (information capacity) of the static state forms a g ent force of influence on the dynamic state, probably providing the increase of the tive organization of R = 0.45 of the microstructure of the AA cell. A single informa maximum is observed in the near-order change of ∆D in the scale change interval of can be assumed that the development of the crystal intergrowth is provided due to high potential of choice (information capacity) of the static state of the information tem. This confirms the nonequilibrium of states due to the presence of gradient-p components in the microstructure of the control sample.
The results for determining the entropy parameters of structure of the control a tive-free sample are presented in Table 1, where H is the general structural informa entropy, Hmax is the maximum information, H* is the configuration entropy, R is the v of relative organization of the system, D1equil is the value of the information dimen K is the information capacity, and∆D is the measure of randomness.
The results of the multifractal analysis of the electron microscopic image of the crostructure of the modified sample are shown in Figure 4. The entropy parameters o static and dynamic stability of microstructure are given in Table 2.  The functional state of the microstructure of the modified sample is determined by the deterministic and stochastic functions of the information system.
The spectrum of information-entropy levels of dynamic deterministic states of the information systems of the control sample for the BB section and the modified sample for  The functional state of the microstructure of the modified sample is determined by the deterministic and stochastic functions of the information system.
The spectrum of information-entropy levels of dynamic deterministic states of the information systems of the control sample for the BB section and the modified sample for the AA section areshown in Figure 4. The influence of the scale on the changes in the information entropy of D 1 and the ordering of K S × 10 −3 in the spectrum of levels of dynamic state of the information system of the modified and control microstructures are analyzed in Figures 5-7.        Based on the fractal analysis of the images, it can be concluded that soot particles act as centers with a high energy potential during the modification of cement composite. Polydispersity clusters are formed around these centers, consisting of various more stable and energetically advantageous hydration products from the point of view of crystallization under these conditions. It is possible to form a cluster network under the condition of a uniform distribution of such soot particles in the hardening cement matrix. Stochastic fractals act as a self-similar unit in this cluster network.

Conclusions
Studies of physical and mechanical properties of the control (additive-free) and modified compositions have shownthat the compressive strength of the control composition is 20.8 MPa, and the tensile strength in bending is 4.26 MPa. The compressive strength of the modified sample with 0.02% of soot is 26.9 MPa, the tensile strength in bending is 5.2 MPa. The results of X-ray phase analysis and differential thermal analysis showthat the number of changes in the composition of the modified sample wasinsignificant, despite a significant increase of its strength characteristics. It can be assumed that the effect of the additive on the conditions of structure formation did not change the composition of hydration products. The results of scanning electron microscopy show changes in the shape and size of structural elements. The analysis of the spectrum of information-entropy levels of the dynamic state of the information system of the microstructures of control and modified samples ( Figures 5 and 6) shows an increase of the information entropy of D 1i and a decrease of the orderliness of K S levels with a scale change (the increase of the cell size of d).
It should be noted that the modification of cement composite increases the thermodynamic fluctuations of chemical and phase compositions, affectsthe parameters of K s orderliness and the entropy of D 1 in the information system of microstructure, and leads to disorder both of the scale invariance and the monotony of dependencies. The break points of the dependencies correspond to non-equilibrium phase transitions. The modification is accompanied by the increase of the orderliness of K s × 10 −3 and the decrease of the information entropy in spectrum of dynamic deterministic state of the composition, with the additive due to the processes of dynamic self-organization in the course of information transformations. The processes of thermodynamic self-organization are observed as nonequilibrium phase transitions in the form of violations of scale invariance in Figures 5 and 6 at the d scale up to 45 pixels.
Based on the fractal analysis of the images, it can be concluded that soot particles act as centers with a high energy potential during the modification of cement composite. Polydispersity clusters are formed around these centers, consisting of various more stable and energetically advantageous hydration products from the point of view of crystallization under these conditions. It is possible to form a cluster network under the condition of a uniform distribution of such soot particles in the hardening cement matrix. Stochastic fractals act as a self-similar unit in this cluster network.

Conclusions
Studies of physical and mechanical properties of the control (additive-free) and modified compositions have shownthat the compressive strength of the control composition is 20.8 MPa, and the tensile strength in bending is 4.26 MPa. The compressive strength of the modified sample with 0.02% of soot is 26.9 MPa, the tensile strength in bending is 5.2 MPa. The results of X-ray phase analysis and differential thermal analysis showthat the number of changes in the composition of the modified sample wasinsignificant, despite a significant increase of its strength characteristics. It can be assumed that the effect of the additive on the conditions of structure formation did not change the composition of hydration products. The results of scanning electron microscopy show changes in the shape and size of structural elements.
Electronic micrographs with amagnification of 20,000 times were used for fractal analysis of microstructures of the control and soot-modified cement composites. Soot with aparticle size less than 80 nm and with aspecific surface area of 65-100 m 2 /g was used as a modifier. The calculation of the multifractal parameters was performed on the basis of parametrization of electron microscopic images of microstructures using acomputer program. The fractal analysis of the images revealed a change in the relative organization of the system, which justified the nature of the change of mechanical properties of the modified material relative to the additive-free systemfrom the perspectiveof self-organization.
Based on the fractal analysis of the images it can be concluded that soot particles acted as centers with a high energy potential during the modification of the cement composite. Modification of composition with dispersed soot led to a change in the configuration entropy of the mixture ∆S* due to the relationship between entropy and the measure of ordering of the microstructure by scale invariance. It can be argued that there is a direct and inverse relationship when controlling the structural state through modification.
Polydispersity clusters are formed around these centers, consisting of various more stable and energetically advantageous hydration products from the point of view of crystallization under these conditions. It is possible to form a cluster network under the condition of a uniform distribution of such soot particles in the hardening cement matrix. Stochastic fractals act as a self-similar unit in this cluster network.