Next Article in Journal
Effect of Doxycycline Microencapsulation on Buccal Films: Stability, Mucoadhesion and In Vitro Drug Release
Previous Article in Journal
Theory of Charged Gels: Swelling, Elasticity, and Dynamics
Previous Article in Special Issue
Equilibrium Swelling of Biocompatible Thermo-Responsive Copolymer Gels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cellular Automata Modeling of Silica Aerogel Condensation Kinetics

Faculty of Chemical and Process Engineering, Warsaw University of Technology, 00-645 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Submission received: 29 January 2021 / Revised: 16 April 2021 / Accepted: 19 April 2021 / Published: 21 April 2021
(This article belongs to the Special Issue Modeling and Simulation of Gel-Based Materials)

Abstract

:
The formation of silica aerogels and the kinetics of condensation were investigated numerically. The influence of the reaction-limited to the diffusion-limited aggregation (RLA to DLA) transition on the reaction kinetics curves and the evolution of the aggregate size distribution during condensation were examined. The 2D cellular automaton was developed and applied to reflect the process of secondary particle aggregation. Several tendencies were observed due to the adjustment of the model parameters: the probability of condensation reaction and the particles’ concentration. The final wet-gel structures’ visualizations proves that the structure becomes more dense and compact due to entering the RLA regime. The simulation time (associated with the gelation time) decreased along with the increase in both model parameters. The lower the collision probability, the slower reaction becomes, and particles are more likely to penetrate the structure deeper until they finally join the aggregate. The developed model reflects the condensation process’s nature and its mechanisms properly and indicates a significant potential for further aerogel synthesis investigations and for the prediction of wet-gel properties according to condensation parameters.

1. Introduction

Silica aerogels have gained increasing attention since they were first synthesized in 1931 [1]. Utilizing alkoxysilanes as precursors, instead of silica oxides, provided control over the final material structure, which is the source of the aerogel’s unique properties. Nonetheless, the formation of aerogel depends on plenty of parameters: temperature, silane/solvent/water molar ratio, catalyst type and its concentration [2,3,4,5,6], to name just a few. The sol-gel procedure conditions lead to changes in the final product features due to differences in the kinetics parameters. The influence of a specific set of parameters remains unclear, and often structures are obtained through trial and error by each group of authors separately. The sol-gel method’s complexity, together with numerous process variables, stems the need for further investigation to better understand aerogel formation’s nature.
In parallel with experimental research, a theoretical and numerical study is being developed, allowing prediction of aerogels’ structure, mechanical and thermal properties and the dynamics of changes in this structure. To date, most of the attention has been directed to the numerical modeling of the aerogels’ structure (the exhaustive review can be found in [7]). Furthermore, much interest was paid to the mechanical properties of aerogels [8]. Among the used numerical methods, it is worth paying attention to molecular dynamics simulations [9,10], which were used to model the aerogels’ mechanical properties. The links between aerogel structures and their mechanical properties have also been analyzed utilizing hard-sphere aggregation models by Ma et al. [11,12]. Other methods used here are coarse-grained models [13,14,15] and multiscale models [7]. The thermal properties of aerogels have been modeled through, e.g., atomistic simulations [16]. The silica aerogels are mainly obtained by the sol–gel process that assembles small and nearly uniform particles into clusters or networks. This mechanism caused cellular automata (CA) [17] or other particle-based algorithms to model aerogel structure and properties. The latter includes, e.g., diffusion-limited aggregation (DLA), reaction-limited aggregation (RLA), diffusion-limited cluster–cluster aggregation (DLCA), reaction-limited cluster–cluster aggregation (RLCA) and a combination of the above. The highest popularity has been reached by the DLCA model [18,19]. DLCA modeling is used to describe, e.g., the structure of an aerogel network [20] and its mechanical deformation [21].
However, despite significant progress in aerogel modeling in recent years, there is still a lack of work describing the kinetics of an aerogel formation. One of a few investigations in this field is our last paper, where we have elaborated a simple kinetics model of the gelation process [22]. The results obtained by employing this model appeared to be in good agreement with experimental results. However, this model gave no insight into an aerogel structure. Other works coped with the kinetics of simple macromolecular structures formation, such as dimers, chains or rings [23,24]. The modelling results are in good agreement with the experimental results of NMR and UV Raman spectroscopy, but they do not allow one to find the complete dependence of gel mass on time.
In the current paper, we propose the two-dimensional particle-based RLA/RLCA and DLA/DLCA modeling of gelation dynamics. We also identify the relations between parameters of the model and parameters of the actual gelation process, and, thus, show the influence of these parameters on the dynamics of gelation kinetics, i.e., the shape of the kinetics curve and the time of gelation. Presented results also give some insight into the morphology and structure of the final product; however, the quantitative information about aerogel structure will provide a 3D model. The future model iteration has a promising potential for further evolution and alterations, taking into account data from experimental techniques, such as SAXS [25] or DLS.

2. Results and Discussion

Figure 1 presents an examination of the influence of the DLA/DLCA to RLA/RLCA transition on the numerically designated condensation reaction kinetics curves. The model parameters were (i) concentration of aerogel secondary particles, c, and (ii) probability of condensation reaction, P. Simulations were conducted for three chosen concentration values: c ≈ 1%, 5% and 12.5% and four probability values: P = 0.001, 0.01, 0.1 (reaction-limited (cluster) aggregation, RLA/RLCA) and 1 (diffusion-limited (cluster) aggregation, DLA/DLCA). Simulations were conducted three times for every case to verify the results’ repeatability (blue, red and yellow curves, as shown in Figure 1).
The time of a simulation (associated with the moment when 100% of secondary particles are connected within one aggregate, i.e., the structure is completely cross-linked) decreases as the system approaches the DLA/DLCA regime due to the higher probability of reaction between two adjacent particles/clusters. The same tendency occurs with an increase in secondary particle concentrations. A higher amount of particles within the computing domain is also responsible for a higher probability of collision between particles/clusters for the same reaction probability value.
The final number of aggregates above the threshold size (0.5% of the total number of cells) is associated with a mass of condensing aerogel. The two alternative methodologies for the experimental designation of condensation kinetics was thoroughly described in our previous paper [22]—via the filtration of gelling solution and determination of the resulting product mass or via the UV–Vis spectrophotometry. Both methods provided the sigmoidal shape of kinetics curves for both samples with or without macroscopic phase separation and “particle aggregates” type of microstructure [26]. We identified two stages and the prevailing mechanisms of aerogel condensation kinetics (Figure 2a). During the first stage of condensation kinetics, the dominant mechanism is the reaction between free molecules of the precursor or aggregates containing few molecules. In the second stage, a surface reaction appears between single precursor molecules and already formed aerogel spheres. The merging of bigger (already “detectable”) clusters was not considered, as it do not affect the mass of the condensing product.
Comparing the previously obtained results (experimental and numerical) with the ones obtained by cellular automata modeling, one can notice that the first condensation stage is not especially evident. It may result from too small a computing domain (more information in Section 4: Materials and Methods). Otherwise, to capture the initial stage of condensation kinetics a different method must be proposed, reflecting primary particles’ formation due to the polycondensation reaction between precursor molecules (e.g., the coarse-grained molecular dynamic). The superposition of two numerical methods shows the potential to thoroughly reflect the reaction kinetics and provide an insight (visualization) into the process.
Figure 2b presents preliminary comparison of the numerically obtained kinetics curve (model parameters: P = 0.01, c = 12.5%) and experimental measurements. The chemical composition of the sample was 1:2:0.6:0.6 (vol. MTMS: MeOH: oxalic acid 0.01 M: ammonia 1 M) [22]. The results were normalized to the maximum value of condensing aerogel mass and total measurement time. The comparison indicates a similarity between the shape of experimental and numerical curves.
Figure 3 presents the visualizations of the final structures obtained by the cellular automata modeling for the chosen values of reaction probability and concentration. For each concentration value (c ≈ 1%, 5% and 12.5%), one can notice that the structure becomes more dense and compact along with the decrease in the condensation reaction probability value. This is a consequence of reaction-limited aggregation—the probability of merging particles/clusters is significantly lower. Thus, the particles are more likely to penetrate the structure deeper until they join the aggregate at some point.
Figure 4 shows exemplary SEM (scanning electron microscopy) structures obtained experimentally during our previous research [21]. The particle aggregates type of structure [26] resembles the character of the numerically obtained wet-gel structures (Figure 3).
The developed model also provides the aggregates’ size distribution (dependent on the number of particles within the aggregates) during the condensation time. The exemplary results are presented in Figure 5. At the system’s initial state (t = 0), the vast majority of particles are single (not connected). However, some particles already merge at the stage of random generation of their initial location. During the condensation, the total number of aggregates decreases, and the distribution moves towards bigger aggregates until the system is entirely cross-linked (t = 1427 dt).
As the presented studies’ main goal is to develop a model reflecting the condensation kinetics, the gelation time value was compared for both experimental and numerical methods.
Figure 6a shows the condensation kinetics obtained by the UV–Vis spectrophotometry for samples with variable ammonia concentration (0.11, 0.17 and 0.22 [mol/dm3]) and constant MTMS/MeOH and acid/base molar ratios (more information in Section 4: Materials and Methods). The cross points marked on the kinetics curves stand for the gelation time values, estimated by the tilting test tube method. It was noticed that gelation occurs approximately right after the inflexion of a kinetics curve.
A similar tendency was observed for numerical results (Figure 6b) when the gelation time was assumed at the point when the structure is cross-linked in 99% (which means that 99% of secondary particles are a part of a detectable aggregate).
Figure 7 presents the analysis of the reaction probability’s influence—assumed to be dependent on base catalyst concentration. Experimentally designated values of gelation time for samples with a variable base/precursor ratio are presented in Figure 7a. The dependence can be approximated by the power function:
t g e l = 36.2 · n N H 4 O H / n M T M S 1.5
With coefficient R2 = 0.85.
Power function can also approximate the dependence of gelation time on reaction probability (Figure 7b):
t g e l = A · P n
The values of parameters in the above formula are given in Table 1.
By comparing Equations (1) and (2), we may observe that reaction probability P as a model parameter is a growing function of molar ammonia to MTMS ratio. More precisely, it reads:
P n N H 4 O H / n M T M S α
where the exponent α = 1.5 n equals between 2 and 3, dependent on the concentration of the secondary particles. By comparing (3) with (8), we may deduce how the ammonia to MTMS molar ratio lowers the activation energy:
E a l n n N H 4 O H / n M T M S

3. Conclusions

The developed model, despite its simplicity, reflects the nature of the condensation process properly. The final structures become denser while the system is approaching the RLA/RLCA regime. The shape of the kinetics curve obtained from a numerical model was similar to that obtained from experimental investigations. The time of gelation decreased with the increase in reaction probability value and the increase in secondary particle concentrations, which itself increased the likelihood of a collision.
The gelation time dependence on the reaction probability (numerical results) and the concentration of the base catalyst (experimental) were approximated by power functions. The gelation time happened right after the inflexion of the experimental kinetics curve. This observation stays in good correspondence with numerical results with gelation time assumption as the 99% extent of cross-linking. Furthermore, it was deduced that the ammonia to MTMS molar ratio lowers activation energy, according to the Equation (4).
The model indicates a significant potential for further aerogel synthesis investigations and predictions of wet-gel properties according to condensation parameters. Further iteration requires the transition into the 3D system, while the present 2D model does not describe the complexity of the 3D aerogel structure.
The next step should be correlating the size of an aggregate with its diffusivity (utilizing, e.g., Einstein–Stokes equation), thus the frequency of motion in the CA system. Another factor that should be considered is the condensation temperature, affecting the reaction mixture’s viscosity and particles’ diffusivity.

4. Materials and Methods

4.1. Cellular Automata Studies

The two-dimensional computing domain representing the condensing mixture was divided into a regular grid of cells with periodic boundary conditions (Figure 8). Each cell can be described by one of a finite number of states—in this particular case, we only had two possible states: (i) a cell was occupied by a secondary particle or (ii) the cell was empty. According to the set percentage concentration value (assumed to be a ratio of occupied to total number of cells), the number of secondary particles was given as one of the model parameters. An initial state of the system (for t = 0) was pseudo-randomly generated. According to the set reaction probability value, if two particles meet in the adjacent cells (the Moore type of neighborhood was applied), the reaction of condensation can occur. The reaction probability is a significant model parameter deciding whether aggregation is limited by diffusion (P = 1, DLA/DLCA) or reaction (P < 1, RLA/RLCA). If the reaction between adjacent cells does indeed occur, the two particles start to act as one aggregate. Particles belonging to different aggregates may also join each other if they appear in neighboring cells. The probability of reaction is assumed to be the same as for free particles. The reaction between particles that already belong to aggregates results in the aggregate growth according to the DLCA or RLCA mechanism, depending on P value.
The probability P is dependent on the kinetics constant of the condensation reaction. Indeed, assuming that two clusters may join in every time step Δt with probability P, we conclude that in every time step there appears to be a number of newly formed aggregates Nnew given as:
N n e w = P · N c o l l
where N c o l l is a number of colliding pair of aggregates and is proportional to the square of aggregates (or secondary particles) concentration N c o l l c 2 . Thus,
N n e w = α · P · c 2
where α is a proportional coefficient, and the condensation rate is given as:
v N n e w t = α P t · c 2
On the other hand, the kinetics theory of chemical reaction reads as follows:
v = k × c 2
Thus, we have k P t .
Following the above, probability P depends on the activation energy, temperature and collision effectiveness. In the same manner as k, P is the increasing function of collision effectiveness and temperature and a decreasing function of activation energy in the form of
P e x p E a k B T
where kB is a Boltzmann constant.
In our parametric study, P was the first parameter of our model.
The second parameter of the model was a relative concentration of occupied cells (which does not change during the simulation) and should be interpreted as a volumetric concentration of precursor particles.
During the simulation, the particles started to diffuse in random directions, with the total probability of the movement equal to one in every time step. The direction of the motion was chosen randomly with uniform probability 0.125, which means that the sedimentation was not considered. Indeed, for objects as small as secondary particles or small aggregates, diffusion was the dominant mechanism of motion. If the particle was a part of an aggregate, the rest of the particles forming that aggregate moved in the same direction as the first one. This means the clusters’ diffusivity was assumed to be not related to their size, which is significant but (at this point of research) an intended simplification. Another simplification arising from these rules of motion is that the aggregates did not change their shape (e.g., as a result of deflection). For relatively small aggregates, this assumption seems to be valid. The computation procedure stops at the point of complete cross-linking, i.e., when all secondary particles were connected, forming one giant aggregate as presented on the video provided in the Supplementary Materials (Video S1). The animation shows an exemplary simulation for DLA mechanism (concentration = 1%, probability of reaction = 100%). The particles represent aerogel secondary particles. The animation shows their motion and aggregation, until the point when the structure is 100% cross-linked (all particles are connected).
For the experimental and numerical kinetics curves comparison, it was assumed that an aggregate could be detected when it reached some threshold size. This threshold corresponds to the minimum weight that is visible during the measurements. Three sizes (aggregates consisting of 0.1, 0.5 or 2.5% of the total number of cells) were tested (Figure 9). The curve for 2.5% was not smooth enough to obtain results corresponding to the experimental data. This work shows the results obtained for the value of 0.5%.
Three sizes of the computing domain were tested: 100 × 100, 200 × 200, 1000 × 1000 cells. The condensation kinetics curves and final structure visualizations are presented in Figure 10. As far as the final structures are concerned, there is no significant impact of the domain size on the gel morphology. However, the kinetics curves stronger resemble the sigmoidal characters for bigger computing domain sizes. To reflect the kinetics curves’ character, observed with UV–Vis spectroscopy, and to limit the computing power, the 200 × 200 grid was chosen.
For the numerical description of the condensation step of silica aerogels’ synthesis, the C++ programming language was applied, and the MATLAB software was used for the graphs’ preparation. The code was included in the Supplementary Materials (Code S1).

4.2. Experimental

The investigated samples were prepared by a two-step (acid-base) sol–gel method. Trimethoxymethylsilane (MTMS) was used as a precursor for the developed synthesis (purchased from Sigma-Aldrich) and methanol (purchased from Stanlab) as a solvent. As acid and base catalysts, the ammonia water (Eurochem BGD) and oxalic acid (Sigma-Aldrich, 0.01 M aqueous solution) were applied.
The syntheses were conducted according to molar ratios provided in Table 2. The oxalic acid solution (0.01 M) was added to start the hydrolysis of MTMS dissolved in methanol. After 1 h, a certain amount of ammonia solution was added. The sample was stirred vigorously for another 30 s, and 2.5 mL of the reaction mixture was transferred into a polystyrene cuvette and placed into the spectrophotometer calibrated on pure methanol. The rest of the sample was used to observe and estimate the time of gelation value by a simple tilting test tube method.
For experimental designation of condensation kinetics (for samples No. 6, 7 and 8), a UV–Vis spectrophotometer (Genesys 10 S UV–Vis spectrophotometer, Thermo Scientific, Waltham, MA, USA) was used at wavelength λ = 633 nm with time intervals (the 20 s for a sample with 0.22 M ammonia solution concentration, 60 s for samples with 0.17 and 0.11 M), until the absorbance remained approximately constant. The absorbance increase in time is proportional to the increase in obtained gel mass during the condensation step of aerogel synthesis [22]. The time of gelation and kinetics for each sample were measured 3 times for each sample and averaged.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/gels7020050/s1, Video S1: animation DLA, conc = 1%, Code S1: program code.

Author Contributions

Conceptualization: J.M.G., N.H.B. and B.N.; methodology: J.M.G. and N.H.B.; investigation: R.P. (Rafał Pakuła), R.P. (Robert Przewodzki), J.M.G. and N.H.B.; writing: N.H.B., J.M.G. and B.N. All authors have read and agreed to the published version of the manuscript.

Funding

Research was funded by POB Technologie Materiałowe of Warsaw University of Technology within the Excellence Initiative: Research University (IDUB) programme (IDUB-TM-2-3).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kistler, S.S. Coherent expanded aerogels and jellies. Nature 1931, 127, 741. [Google Scholar] [CrossRef]
  2. Borovin, E.; Callone, E.; Ribot, F.; Diré, S. Mechanism and kinetics of oligosilsesquioxane growth in the in situ water production sol-gel route: Dependence on water availability. Eur. J. Inorg. Chem. 2016, 2016, 2166–2174. [Google Scholar] [CrossRef]
  3. Lee, A.S.; Choi, S.S.; Baek, K.Y.; Hwang, S.S. Hydrolysis kinetics of a sol-gel equilibrium yielding ladder-like polysilsesquioxanes. Inorg. Chem. Commun. 2016, 73, 7–11. [Google Scholar] [CrossRef]
  4. Endo, H.; Takeda, N.; Unno, M. Synthesis and properties of phenylsilsesquioxanes with ladder and double-decker structures. Organometallics 2014, 33, 4148–4151. [Google Scholar] [CrossRef]
  5. Kessler, D.; Löwe, H.; Theato, P. Synthesis of defined poly(silsesquioxane)s: Fast polycondensation of trialkoxysilanes in a continuous-flow microreactor. Macromol. Chem. Phys. 2009, 210, 7. [Google Scholar] [CrossRef]
  6. Luo, Y.; Li, Z.; Zhang, W.; Yan, H.; Wang, Y.; Li, M.; Liu, Q. Rapid synthesis and characterization of ambient pressure dried monolithic silica aerogels in ethanol/water co-solvent system. J. Non-Cryst. Solids 2019, 503, 214–223. [Google Scholar] [CrossRef]
  7. Maximiano, P.; Duares, L.; Simoes, P. Overview of Multiscale Molecular Modeling and Simulation of Silica Aerogels. Ind. Eng. Chem. Res. 2019, 58, 18905. [Google Scholar] [CrossRef]
  8. Patil, S.P.; Parale, V.G.; Park, H.-H.; Markert, B. Mechanical modeling and simulation of aerogels: A review. Ceram. Int. 2020; in press. [Google Scholar] [CrossRef]
  9. Rivas Murillo, J.S.; Bachlechner, M.E.; Campo, F.A.; Barbero, E.J. Structure and mechanical properties of silica aerogels and xerogels modeled by molecular dynamics simulation. J. Non-Cryst. Solids 2010, 356, 1325–1331. [Google Scholar] [CrossRef]
  10. Patil, S.P.; Rege, A.; Itskov, M.; Markert, B. Mechanics of Nanostructured Porous Silica Aerogel Resulting from Molecular Dynamics Simulations. J. Phys. Chem. B 2017, 121, 5660–5668. [Google Scholar] [CrossRef] [PubMed]
  11. Ma, H.-S.; Roberts, A.P.; Prévost, J.-H.; Jullien, R.; Scherer, G.W. Mechanical structure-property relationship of aerogels. J. Non-Cryst. Solids 2000, 277, 127–141. [Google Scholar] [CrossRef]
  12. Ma, H.-S.; Prévost, J.-H.; Jullien, R.; Scherer, G.W. Computer simulation of mechanical structure-property. J. Non-Cryst. Solids 2001, 285, 216–221. [Google Scholar] [CrossRef]
  13. Gelb, L.D. Simulating Silica Aerogels with a Coarse-Grained Flexible Model and Langevin Dynamics. J. Phys. Chem. C 2007, 111, 15792–15802. [Google Scholar] [CrossRef]
  14. Ferreiro-Rangel, C.A.; Gelb, L.D. Investigation of the bulk modulus of silica aerogel using molecular dynamics simulations of a coarse-grained model. J. Phys. Chem. B 2013, 117, 7095–7105. [Google Scholar] [CrossRef]
  15. Ferreiro-Rangel, C.A.; Gelb, L.D. Computational study of uniaxial deformations in silica aerogel using a coarse-grained model. J. Phys. Chem. B 2015, 119, 8640–8650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Morthomas, J.; Gonçalves, W.; Perez, M.; Foray, G.; Martin, C.L.; Chantrenne, P. A novel method to predict the thermal conductivity of nanoporous materials from atomistic simulations. J. Non-Cryst. Solids 2019, 516, 89–98. [Google Scholar] [CrossRef]
  17. Menshutina, N.; Lebedev, I.; Lebedev, E.; Paraskevopoulou, P.; Chriti, D.; Mitrofanov, I. A cellular automata approach for the modeling of a polyamide and carbon aerogel structure and its properties. Gels 2020, 6, 35. [Google Scholar] [CrossRef]
  18. Meakin, P. Formation of fractal clusters and networks by irreversible diffusion-limited aggregation. Phys. Rev. Lett. 1983, 51, 1119. [Google Scholar] [CrossRef]
  19. Kolb, M.; Botet, R.; Jullien, R. Scaling of kinetically growing clusters. Phys. Rev. Lett. 1983, 51, 1123. [Google Scholar] [CrossRef]
  20. Hasmy, A.; Anglaret, E.; Foret, M.; Pelous, J.; Jullien, R. Small-angle neutron scattering investigation of long-range correlations in silica aerogels: Simulations and experiments. Phys. Rev. B 1994, 50, 6006–6016. [Google Scholar] [CrossRef] [PubMed]
  21. Ma, H.-S.; Prévost, J.-H.; Scherer, G.W. Elasticity of DLCA model gels with loops. Int. J. Solids Struct. 2002, 39, 4605–4614. [Google Scholar] [CrossRef]
  22. Borzęcka, N.H.; Nowak, B.; Gac, J.M.; Głaz, T.; Bojarska, M. Kinetics of MTMS-based aerogel formation by the sol-gel method—Experimental results and theoretical description. J. Non-Cryst. Solids 2020, 547, 120310. [Google Scholar] [CrossRef]
  23. Zhang, X.-Q.; Trinh, T.T.; van Santen, R.A.; Jansen, A.P.J. Mechanism of the Initial Stage of Silicate Oligomerization. J. Am. Chem. Soc. 2011, 133, 6613. [Google Scholar] [CrossRef] [PubMed]
  24. Depla, A.; Lesthaeghe, D.; van Erp, T.S.; Aerts, A.; Houthoofd, K.; Fan, F.; Li, C.; Van Speybroeck, V.; Waroquier, M.; Kirschock, C.E.A.; et al. 29SiNMR and UV-Raman Investigation of Initial Oligomerization Reaction Pathways in Acid-Catalyzed Silica Sol-Gel Chemistry. J. Phys. Chem. C 2011, 115, 3562. [Google Scholar] [CrossRef]
  25. Abdusalamabov, R.; Itskov, M.; Milow, B.; Reichenauer, G.; Rege, A. Investigation of the fractal properties of aerogels by diffusion-limited aggregation models. In Proceedings of the 8th GACM Colloquium on Computational Mechanics: For Young Scientists From Academia and Industry, Kassel, Germany, 28–30 August 2019. [Google Scholar] [CrossRef]
  26. Nakanishi, K.; Kanamori, K. Organic-inorganic hybrid poly(silsesquioxane) monoliths with controlled macro- and mesopores. J. Mater. Chem. 2005, 15, 3776–3786. [Google Scholar] [CrossRef]
Figure 1. Condensation kinetics curves obtained by cellular automata modeling for different values of reaction probability (P) and concentrations of secondary particles (c).
Figure 1. Condensation kinetics curves obtained by cellular automata modeling for different values of reaction probability (P) and concentrations of secondary particles (c).
Gels 07 00050 g001
Figure 2. (a) The assumed division of condensation kinetics curves into two stages and (b) qualitative comparison of experimental [22] and numerical condensation kinetics curve.
Figure 2. (a) The assumed division of condensation kinetics curves into two stages and (b) qualitative comparison of experimental [22] and numerical condensation kinetics curve.
Gels 07 00050 g002
Figure 3. Final wet–gel structure visualizations obtained by cellular automata modeling for different values of reaction probability and concentrations of secondary particles.
Figure 3. Final wet–gel structure visualizations obtained by cellular automata modeling for different values of reaction probability and concentrations of secondary particles.
Gels 07 00050 g003
Figure 4. SEM images of exemplary aerogel structures. Chemical compositions (vol. MTMS: MeOH: oxalic acid 0.01 M: ammonia 1 M) (a) 1:2:0.6:0.6 and (b) 1:2:0.8:0.8. Magnification 2000× [22].
Figure 4. SEM images of exemplary aerogel structures. Chemical compositions (vol. MTMS: MeOH: oxalic acid 0.01 M: ammonia 1 M) (a) 1:2:0.6:0.6 and (b) 1:2:0.8:0.8. Magnification 2000× [22].
Gels 07 00050 g004
Figure 5. The evolution of the aggregate size distribution during condensation (c = 5%, P = 1). Colors of aggregates were randomly generated for easy distinction.
Figure 5. The evolution of the aggregate size distribution during condensation (c = 5%, P = 1). Colors of aggregates were randomly generated for easy distinction.
Gels 07 00050 g005
Figure 6. The condensation kinetics curves obtained (a) experimentally (UV–Vis spectrophotometry) and (b) numerically. The gelation time values (cross points) estimated by (a) tilting test tube method and (b) assuming the gelation time point occurs for structure cross-linked in 99%.
Figure 6. The condensation kinetics curves obtained (a) experimentally (UV–Vis spectrophotometry) and (b) numerically. The gelation time values (cross points) estimated by (a) tilting test tube method and (b) assuming the gelation time point occurs for structure cross-linked in 99%.
Gels 07 00050 g006
Figure 7. The dependence of (a) gelation time on base/precursor molar ratio (for different base/acid ratios, marked with different point colors) and (b) logarithm of gelation time on the reaction probability logarithm.
Figure 7. The dependence of (a) gelation time on base/precursor molar ratio (for different base/acid ratios, marked with different point colors) and (b) logarithm of gelation time on the reaction probability logarithm.
Gels 07 00050 g007
Figure 8. Scheme of periodic boundary conditions.
Figure 8. Scheme of periodic boundary conditions.
Gels 07 00050 g008
Figure 9. The influence of the threshold aggregate size on the shape of kinetics curves. Secondary particles concentration: 10%, reaction probability P = 10%, grid size: 200 × 200 cells.
Figure 9. The influence of the threshold aggregate size on the shape of kinetics curves. Secondary particles concentration: 10%, reaction probability P = 10%, grid size: 200 × 200 cells.
Gels 07 00050 g009
Figure 10. The influence of the threshold aggregate size on the shape of kinetics curves. Secondary particles concentration: 1%, reaction probability P = 10%, grid size: 200 × 200 cells.
Figure 10. The influence of the threshold aggregate size on the shape of kinetics curves. Secondary particles concentration: 1%, reaction probability P = 10%, grid size: 200 × 200 cells.
Gels 07 00050 g010
Table 1. The parameters describing the power dependence of gelation time on the reaction probability.
Table 1. The parameters describing the power dependence of gelation time on the reaction probability.
Secondary Particles Concentration (%)A Parametern ParameterR2
13901.8−0.560.94
5761.8−0.610.95
12.590.1−0.771
Table 2. The chemical composition and the values of gelation time for the investigated samples.
Table 2. The chemical composition and the values of gelation time for the investigated samples.
Sample NrMolar RatioConc. NH4OH [mol/dm3]Gel. Time [s]S.dev [s]
MTMSMeOHH2OOxalic AcidNH4OH
117.3614.871.340.130.19913.67129.62
217.3614.931.340.070.092445.00141.07
317.3616.521.490.150.20737.67329.82
417.3618.171.640.160.21540.67147.78
517.3619.781.790.220.28634.6749.00
617.3619.821.790.180.22366.00214.00
717.3619.861.790.130.17966.6715.00
817.3619.911.790.090.111411.67128.00
917.3621.471.940.190.23328.00129.90
1017.3623.082.090.260.30324.007.94
1117.3623.122.090.210.24301.67158.66
1217.3623.172.090.160.18640.0017.32
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Borzęcka, N.H.; Nowak, B.; Pakuła, R.; Przewodzki, R.; Gac, J.M. Cellular Automata Modeling of Silica Aerogel Condensation Kinetics. Gels 2021, 7, 50. https://doi.org/10.3390/gels7020050

AMA Style

Borzęcka NH, Nowak B, Pakuła R, Przewodzki R, Gac JM. Cellular Automata Modeling of Silica Aerogel Condensation Kinetics. Gels. 2021; 7(2):50. https://doi.org/10.3390/gels7020050

Chicago/Turabian Style

Borzęcka, Nina H., Bartosz Nowak, Rafał Pakuła, Robert Przewodzki, and Jakub M. Gac. 2021. "Cellular Automata Modeling of Silica Aerogel Condensation Kinetics" Gels 7, no. 2: 50. https://doi.org/10.3390/gels7020050

APA Style

Borzęcka, N. H., Nowak, B., Pakuła, R., Przewodzki, R., & Gac, J. M. (2021). Cellular Automata Modeling of Silica Aerogel Condensation Kinetics. Gels, 7(2), 50. https://doi.org/10.3390/gels7020050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop