Next Article in Journal
Influence of Interpenetrating Chains on Rigid Domain Dimensions in Siloxane-Based Block-Copolymers
Next Article in Special Issue
Knot Factories with Helical Geometry Enhance Knotting and Induce Handedness to Knots
Previous Article in Journal
Numerical Studies on Failure Mechanisms of All-Composite Sandwich Structure with Honeycomb Core under Compression and Impact Loading Conditions
Previous Article in Special Issue
Molecular Insights into the Wall Slip Behavior of Pseudoplastic Polymer Melt in Nanochannels during Micro Injection Molding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Silicic Acid Polymerization and SiO2 Nanoparticle Growth in Hydrothermal Solution

1
Research Geotechnological Center, Far Eastern Branch of Russian Academy of Sciences, 30, Severo-Vostochny Highway, 683002 Petropavlovsk-Kamchatsky, Russia
2
Chemical Department of Moscow State University, Leader Scientist, Russian Federation, Leninskie Gory, 1, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(19), 4044; https://doi.org/10.3390/polym14194044
Submission received: 15 July 2022 / Revised: 9 September 2022 / Accepted: 19 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Molecular Simulation and Modeling of Polymers)

Abstract

:
The approach of numerical simulation of orthosilicic acid OSA polymerization and SiO2 nanoparticle formation in hydrothermal solution have been developed based on the model of the homogeneous stage of nucleation and the subsequent growth of particles. The influence of surface tension on the interface of SiO2–water, the rate of molecular deposition, and Zeldovich factor Z were evaluated. Temperature dependence on time, pH, initial OSA concentration, and ionic strength are the main parameters that determine the kinetics of colloid phase formation, the final average size of SiO2 nanoparticles, and the particle size distribution and its polydispersity index. The results of the numerical simulation were verified with experimental data on OSA polymerization and measurement of nanoparticles sizes using the method of dynamic light scattering in a wide range of temperatures of 20–180 °C, pH = 3–9, SiO2 content Ct of 300–1400 mg/kg, and ionic strength Is of 0.0001–0.42 mol/kg. The results obtained can be used in the technology of hydrothermal synthesis of sols, gels, and nanopowders to regulate the kinetics of OSA polymerization and SiO2 nanoparticle growth, particle size distribution, morphology, and structure of products.

1. Introduction

Hydrothermal synthesis of SiO2 nanoparticles is the new technological method with a low production cost, high density of silanol groups, low toxicity, excluding the use of chemical reagents and other green chemistry methods of nanoparticle synthesis [1,2,3], and allows to obtain different forms of nanosilica—sols, gels, powders, products of the sol-gel process applied in concrete, agricultural plants, medicines, and other materials [4,5,6,7]. The technology of hydrothermal synthesis includes the stages of OSA polymerization, ultrafiltration membrane concentration of nanoparticles, the sol-gel process, vacuum sublimation, and others [4,5,6]. The temperature regime, the stage of polymerization, pH, initial OSA concentration Cs, and ionic strength influenced the rate of SiO2 nanoparticle formation, the average nanoparticle size, and the index of polydispersity of nanoparticles in hydrothermal sols. Therefore, the density, morphology, agglomerate dimensions and fractal structures, volume, surface, diameters, and structure of pores of nanopowders produced from hydrothermal sols are also affected [4,5]. There is a need to regulate the parameters at the stage of OSA polymerization of silicic acid by means of numerical simulation to achieve the necessary characteristics of sols, gels, and nanopowders.
The kinetics of OSA polymerization with different mechanisms has been studied by many researchers under various conditions: pH, initial concentration of dissolved silica Cs, solution temperature, ionic strength, and sizes and molecular weight of SiO2 particles formed.
At sufficiently high concentrations of silica in 2–6 wt.% sodium silicate salt solutions with the ratio SiO2:Na2O = 3.25 and with acidification at pH = 1.7, when the increase in molecular weight was accompanied mainly by the aggregation of small particles the degree of polymerization increased proportionally to the square root of the polymerization time [8].
At the monomer concentration of 0.24 % in the absence of salts [9,10], the average molecular weight at pH = 2.0 increased linearly depending on the square root of time; at pH = 3.2, 3.8, it increased proportionally with time; at pH = 4.36—linearly from the square of time. At pH = 2.2, temperature 90 °C, the average particle diameter was 1.3 nm. At a higher silica concentration of 8 g/L after about 150 h at pH = 3.0 or 75 h at pH = 2.0 the average molecular weight reached 300 SiO2 with a smaller resulting particle size of 0.9 nm [11]. The average molecular weight of the newly formed polymerized acid prepared by removing Na from a solution of sodium silicate at pH= 2–3 was about 660 SiO2 atoms [12]. At pH = 2.1, temperature 25 °C, in 0.24 wt.% SiO2 solution after 60 days after the start of polymerization, the proportion of monomeric silica was 14% and the particle diameter was 1.4 nm [13]. At pH = 8.8, growth of particles occurred during 1 min, that is 3∙104 times faster than at pH = 2.0. In studies performed at 25 °C, pH = 2.0, in the solution of monosilicic acid 2.4% SiO2 for 24 days the particle diameter was 1.65 nm; after 8 days it was about 3.4 nm [14,15,16]. In a wide range of SiO2 contents at pH =3.0 and pH = 6.1, the initial average diameter of the nuclei was 1.5 nm and the final particle size was about 3 nm [17]. In [18], after 30 days of polymerization at pH = 7.0 and 25 °C with a total concentration of SiO2 in a solution of 0.4 g/L the diameter of nanoparticles corresponding to the solubility was 3.7 nm.
In the pH = 7.0–10.0 region, when colloid particles are formed, but the aggregation process cannot proceed, the reaction order was 3 [13,19]. However, in [20], at different initial concentrations of OSA the apparent “order” of the polymerization reaction was determined as equal to 2 at a pH above 2 and was equal to 3 at a pH below 2 [20], whereas the second-order reaction rate constant decreased.
It was found that the rate of polymerization was the greatest at pH = 8.3 and close to the maximum in the pH = 7.0–9.0 region, and at pH < 7.0 and pH > 9.0 it noticeably decreased [21].
The detailed study of polymerization in the wide area of concentrations of 0.2–1.8 g/L of silica and pH = 4.0–10.0 at 25 °C confirmed that there is an induction period during which the polymerization of the monomer is only slight or nonexistent. A decrease in the monomer content as it polymerizes at pH = 7.0 and 25 °C at an initial content of 0.05 wt. % occurred according to the law 0.37∙(t − 0.35)−0.333 [22]. The induction period was 0.35 h. After 1 h after the decrease in the monomer content, silica particles began to grow rapidly without the formation of any intermediate product. The presence of a long induction period at pH = 4.5–5.5 and a temperature of 95 °C is shown in [23]. The induction period characterizes the metastability of the system and for aqueous silicic acid at pH = 7.75, t = 20 °C, Cs = 1450 mg/kg under mixing conditions it decreased from 3.25 to 2.75 h [20,21]. In the strongly acidic region, the rate of the process is proportional to [H+]1.2 at pH < 1.8, and in the weakly acidic region at pH > 3.4 –[OH-]0.9 [24,25]. The rate increases initially with an increase in temperature from 20 to 90 °C and decreases with a decrease in the initial OSA concentration [24,25]. The results obtained in [26] for OSA solutions with a SiO2 content of 300–1200 mg/kg and low salinity in the temperature range from 5 to 180 °C confirmed that the temperature changes at a constant silica concentration do not significantly affect the increase in induction time at temperatures from 90 to 180 °C.
In a dimer–monomer mixture at pH = 2.0 [27] the polymerization mechanism and particle sizes differed from the results of monosilicic acid polymerization obtained for polycyclic polymer formation in [14,15,16]: at the beginning the discrete particles of less than 4 nm in size formed, then such particles aggregated [27].
Thus, the experimental data obtained by the authors showed significant differences in the values of the polymerization reaction order, the rate constant, the induction period, the sizes of the formed nanoparticles, and the different effects of temperature, pH, and ionic strength on these parameters. Taking this into account, it is justified to use the method of numerical simulation of polymerization based on a certain mathematical model.
One of the mathematical models of orthosilicic acid polymerization is Fleming’s model [28]. Fleming suggested that the process under study is characterized by two main areas of flow. In the first region, the initial concentration of silicic acid exceeds the pseudo-equilibrium concentration of Cx, and the polycondensation reaction is a first-order reaction both with respect to the difference (Cs-Cx) and with respect to the surface concentration of CSiO-ionized groups of SiO, that is, the surface charge of colloid particles. The second area is defined by the condition Cs < Cx.
The decrease in the rate of transition of monomeric silica to colloid is explained using the Weres-Yi-Tsao model [29]. The authors studied two processes: the formation of colloid particles at the homogeneous stage of nucleation and the subsequent growth of particles due to OSA molecular deposition. From this model, developed on the basis of classical concepts of nucleation by Lose-Pound, it follows that the supersaturation of SN (T), equal to Cs/Ce, Ce (solubility of amorphous silica in water at temperature T), and pH are the main factors determining the rate of nucleation IN of silicic acid in an aqueous solution.
The authors of [30] used two methods to study the OSA polymerization, transfer in nanocolloids, and nanosilica precipitation in synthetic solutions based on the Na2SiO3 precursor: (1) the concentration model, when the rate of monomeric silica concentration Cs decrease is proportional to the 4-th order of k1∙Cs4 and the rate of nanosilica transfer to precipitated silica is proportional to the 1-st order of k2∙Cnano; and (2) the supersaturated model, when the rate of monomeric silica concentration Cs decrease is proportional to the 4-th order of the difference k1∙(Cs–Ce)4, Ce—solubility of amorphous silica in the water at 25 °C. In the wide range of initial Cs 240–1200 mg/kg, ionic strength 0.01–0.24 mol/kg, and pH = 3–7 it was obtained that constants k1 and k2 did not differ for concentration and supersaturation models. The models cannot predict the sizes and distribution of SiO2 particles.
The method of molecular dynamic simulation for OSA polymerization was developed in [31,32] at the temperature up 2700 K and can provide the structure of polymerized clusters which was in accordance with the 29Si NMR method.
The models with different orders of polymerization reaction [8], taking into account the concentration of OH groups [33] and ionic strength [28] have been used for analysis of experimental results of orthosilicic acid polymerization in hydrothermal solutions [34].
In [35], the model of Weres-Yi-Tsao [29] has been applied to experimental results in hydrothermal solutions because it can predict not only time dependence of the OSA concentration but also evaluate how the average size of SiO2 nanoparticles depends on time. The main purpose of this work was to investigate nanoparticle size distribution dependence on the temperature, ionic strength, and other parameters of the polymerization process needed to regulate the characteristics of nanosilica sol, gel, and powders produced using the technology of hydrothermal synthesis [4,5,6,7]. For this purpose, it is necessary: (1) to study the dependence of surface tension σsw, rate of molecular deposition Rmd, and Zeldovich factor Z on the temperature and pH to explain the mechanisms of their influence on the results of the polymerization process; (2) to examine the dependence of the average nanoparticle size and induction period on temperature, pH, and ionic strength; (3) to study the nanoparticle size distribution under conditions of time variable temperatures; (4) to verify the results of the numerical simulation with the experimental data in wide range of T, pH, and Is; and (5) to compare the results of the numerical simulation under time-variable temperature conditions at different technological stages of the hydrothermal sol production with the results of determining the average size of SiO2 nanoparticles and their size distribution using the DLS method.

2. Method of Numerical Simulation of OSA Polymerization

The initial concentration of molecules of orthosilicic acid formed due to dissolution of minerals of rocks [4,5,6,7] is determined using quartz solubility Ccr (mol/kg) of α-quartz in water at the temperature T(K) [36]:
lgC cr = 1.468 + 252.9 T 3.217 · 10 5 T 2 .  
After the hydrothermal solution reaches the surface pressure and the temperature decreases, the supersaturation of the solution depends on solubility Ce of amorphous silica [37]:
lgC e = 0.1185 1.126 · 10 3 T + 2.3305 · 10 5 T 2 3.6784 · 10 7 T 3 .  
Polymerization of orthosilicic molecules in oversaturated solution is provided by the formation of siloxane bonds and dehydration:
H 4 SiO 4 + H 4 SiO 4     Si 2 O 7 H 6 + H 2 O
The total silica concentration Ct in the solution is equal to the sum of the concentrations of colloid silica Ccol, soluble silicic acid Cs, and ions of silicic acids Cin:
C t = C col + C s + C in .  
The fraction of dimers does not exceed 1.0% from Ce, the fraction of trimers is no more than 0.1% [26], the fraction of tetramers and low-molecular cyclic polymers is less than 0.1%, and the fraction of Cin does not exceed 14.0%.
The physical and chemical characteristics of the hydrothermal solution from the Mutnovsky geothermal electric plant are presented in the Table 1.
We developed the mathematical model [29] which is based on the calculation of the rate of nucleation IN (nucl/(kg∙s)) dependent on oversaturation SN (T) = Cs/Ce:
I N = Q LP · Z · ( R md · A cr · N A · ρ · M Si 1 ) · e Δ Fcr k B · T .  
where ΔFcr is the change in free energy due to the formation of a nucleus of critical radius Rcr with a surface area of Acr = 4·π·Rcr2 and the number of SiO2 molecules ncr, σsw is the surface tension at the silica–water interface, Rmd is the rate of molecular deposition of SiO2 on a solid surface, kg/(m2∙s), kB = 1.38∙10−23 J/K, MSi = 0.060 kg/mol, NA =6.02∙1023 mol−1, QLP is the Lothe–Pound factor, QLP = 3.34∙1025 kg−1, and Z is the Zeldovich factor, calculated as follows:
R cr = 2 · σ SW · M Si ρ · N A · k B · T · lnS N .  
n cr = ( 4 π 3 ) · ( ρ · N A M Si ) · R cr 3 .  
Δ F cr = σ SW · A cr 3 = ( 16 · π 3 · σ sw 3 · M Si ρ · N A · k B · T · lnS N ) 2   .  
Z = [ ( 2 Δ F cr / n cr 2 ) / ( 2 · π · k B · T ) ] 0.5 = ( 2 3 ) · ( 3 · M Si 4 · π · ρ · N A · n cr 2 ) 1 3 · (   σ sw k B · T ) 0.5 .  
with ρ = 2200 kg/m3—the density of amorphous silica.
The function Rmd, which is the rate of molecular deposition of SiO2 as a function of temperature and pH of the solution, is expressed as follows in the model [29]:
R md = F ( pH ,   pH nom ) · k OH ( T ) · f f ( S a ) · ( 1 S N 1 ) .  
lg k oH = 3.11 4296.6 T .
f f = S t 5 ,   S a < S t .
f f = S t 5 + 5 · S t 4 · ( S a S t ) ,   S a > S t .  
lg S t = 0.0977 + 75.84 T .  
F ( pH , pH nom ) = h f · f ( pH ) + ( 1 h f ) · f ( pH nom ) .  
f ( pH nom ) = f ( pH nom ) f ( 7.0 ) .  
lg f = pH pK i + lg [ Na + ] .
The coefficient of surface tension σsw depends on temperature and pH via the function I(pH, pHnom):
σ sw = H σ T · S σ 2.302 · 10 3 · n 0 · k B · T · I ( pH , pH nom ) .
I ( pH ,   ph nom ) = 0.119 pH F ( pH , pH nom ( pH ) ) dp H .  
where Sa = (1 − αi)·SN, αi is a fraction of silicic acid in ionic form, pHnom = pH + lg([Na+]/0.069), [Na+] is an ion activity [Na+], mol/kg, pKi = 6.4, f(7.0) = 0.119, hf = 0.45, Hσ and Sσ are specific enthalpy and entropy of silica surface in water, Hσ = 63.68·10−3 J/m2, Sσ = 0.049·10−3 J/m2·K, and no = 6.84∙1018 m−2 [29].
Time dependence IN (t) was calculated considering the induction time τin that is necessary to grow and form a stable population of the particles with sizes close to critical [29]:
I N ( t ) = I N · ( 1 e Tp τ in )   .  
τ in = 1.08 · 10 6 · ( R md ) 1 · ( Q LP · Z · R cr 2 · e Δ Fcr k B · T ) 0.25 .  
Initial values were entered: the temperature T, pH, concentration of ions and ionic strength of the solution, total silica Ct content, time step DTp, initial radii R(I), and particle quantity MP(I) of every class “I”. The Runge–Cutt method was used for the numerical simulation.
The total content of colloid, monomeric silica, and the current value of supersaturation were estimated:
M col = I   = 1 I   =   N 4 · π · ρ · R 3 ( I ) · MP ( I ) 3 .  
C s = C t M col .  
Solution degree of supersaturation SN was:
S N = C s C e ( T ) .  
Then, the calculation of values of σsw, Rmd, Rcr, Z, τin, and IN(t) corresponding to the current values of SN, pH, and T using the equations was undertaken.
The quantity of new particles appearing during the time DTp in accordance with the current value of nucleation rate IN on a given program step (N+1) was calculated and it was equal to the quantity of particles in the new class N+1:
MP ( N + 1 ) = I N · DT p
The summary of the concentration of particles CONP was calculated in a program cycle:
CONP = I   = 1 I   =   N + 1 MP ( I )   .  
The increment in mass DPM of every class “I” during the time DTp was calculated:
DPM ( I ) = 4 · π · R 2 ( I ) R md · DT p
The calculation of particle radius R(I) corresponding to a new value of particle mass in every class I was:
R ( I ) = ( R 3 ( I ) + 3 · DPM ( I ) 4 · π · ρ ) 1 3 .  
The mean values of Ra, R2a, and R3a in all classes of particles were estimated using the program cycle:
R a = I + 1 I   =   N + 1 R ( I ) · MP ( I ) CONP   .  
R a 2 = I + 1 I   =   N + 1 R 2 ( I ) · MP ( I ) CONP   .  
R a 3 = I + 1 I   =   N + 1 R 3 ( I ) · MP ( I ) CONP   .  

3. Results

3.1. Dependence of the Surface Tension σsw, Rate of Molecular Deposition Rmd, and Zeldovich Factor Z on the Temperature and pH

The main factors influencing the rate of the nucleation process IN and the final average particle size Ra are the degree of temperature T, the degree of supersaturation SN(t), and pH. In addition, the final result of the numerical simulation is influenced by the functions f(pH, pHnom) and i(pH, pHnom), specified in the Equations (10)–(19). The results of the calculation of the function f(pH) at different pH in the range from 4.0 to 8.7 are presented in Table 1.
With an increase in pH, the value of the f(pH) function increases (Table 2). Consequently, with an increase in pH, the rate of molecular deposition of Rmd increases, and the Rmd parameter will strengthen the tendency to decrease the final average particle size. With an increase in the initial concentration of Cs, the rate of molecular deposition also increases (Figure 1 and Figure 2). An increase in temperature has an inhomogeneous effect on the rate of molecular deposition (Figure 1 and Figure 2): first, due to the increase in the kinetics of the process with increasing temperature, the rate of molecular deposition increases, then Rmd decreases, since with increasing temperature, the solubility of Ce increases, the supersaturation of SN decreases, and the value of the rate of molecular deposition decreases. At a fixed pH, Rmd reaches a maximum at a certain temperature and the value of this maximum increases with an increase in the initial concentration (Figure 2).
An increase in the rate of molecular deposition with an increase in Cs and temperature will contribute to an increase in the rate of nucleation IN of OSA and a decrease in the final average particle size.
The values of the function i (pH) in the range from 6.0 to 8.7 are presented in Table 3. With an increase in pH, the value of the function i (pH) increases. Consequently, the surface tension coefficient and the critical radius Rcr reduced. The surface tension coefficient σsw decreases with increasing temperature and with increasing pH in accordance with Equations (18) and (19) (Figure 3a,b). Consequently, σsw will be a factor that, with an increase in temperature and pH, will contribute to a decrease in the critical radius Rcr, an increase in the nucleation rate IN, and a decrease in the final average particle size Ra, competing with the factor of reducing supersaturation SN.
The physical meaning of the Zeldovich factor Z, calculated using Equation (9), is that only a certain proportion of particles whose sizes have reached a critical radius continue to grow further. The remaining particles decrease in size. At the same time, the concentration of particles whose sizes are higher than the critical radius is less than the equilibrium concentration. With an increase in temperature at a constant pH, the value of the Zeldovich factor decreases, and at the same time its effect on the rate of nucleation IN competes with the opposite effect associated with a decrease in the supersaturation of SN. An increase in pH leads to an increase in the values of the Zeldovich factor (Table 4).

3.2. Results of the Numerical Simulation at Different Temperatures and pH: Average Particle Size, Induction Period, and the Particle Size Distribution

The results of numerical simulation showed (Table 5) that when the temperature increases at a constant pH the final average size of SiO2 nanoparticles increases because the degree of oversaturation and IN become higher. When the pH decreases at a constant temperature the final average size increases (Table 5) due to the increase in the surface tension and the increase in the critical size of the nucleus Rcr and the decrease in IN. If the initial concentration of orthosilicic acid increases the final average size decreases because the degree of oversaturation increases, IN increases, and Rcr decreases (Table 6).
After nucleation and polymerization completion the fraction of dimers is not more than 1.0 % of solubility Ce(T), the concentration of trimers is about 0.1 %, and the concentration of tetramers and polymers with circular structures of less than 6 units of SiO2 is lower than 0.1 % [26].
The concentration of ions H3SiO4 and H2SiO42− is not more than 14.0 % calculated from the dependences of constants of orthosilicic acid ionization of the first K1 and second K2 stages:
lg K 1 = 2549 T 15.36 · 10 6 · T 2  
lg K 2 = 5.37 3320 T 20 · 10 3 · T  
The function of oversaturation SN vs. time Tp showed the long period of induction τin and homogeneous nucleation which becomes much longer at low values of pH when SN changes slightly and the period of rapid growth of SiO2 nanoparticles when SN decreases sharply (Figure 4 and Figure 5, Table 7). The number of new particles PNS depending on the duration of nucleation Tp reached a constant value when the rate of nucleation Rnuc rapidly decreased (Figure 4 and Figure 5). The period of homogeneous nucleation increases when pH decreases and the initial degree of oversaturation SN decreases at higher temperatures and low initial Cs concentrations (Table 7). The relation (thom/thetg) increases when the temperature increases and pH decreases (Table 7).
In Figure 6, the distributions of the number of particles vs. radius formed at constant T = 50–60 °C, pH = 4 and 8, and initial OSA concentrations Cs = 700 and 800 mg/kg are shown. The form of the function of SiO2 particle number distribution by radius has a maximum corresponding to the average radius Ra. Figure 6 shows an asymmetrical distribution, that is the right branch of the silica particle size distribution formed at constant temperature is convex and the left branch is concave.
The average radius of SiO2 nanoparticles and the distribution by radius were determined using the method of dynamic light scattering (DLS). Values Ra were in the range from 5.0 up to 20.0 nm and values of particles radii were in the range from 1.0 up to 50.0 nm [4,5,6,7].

3.3. Simulation of OSA Polymerization at Different Ionic Strengths

The dependence of Cs (Tp) is also affected by the ionic strength of the Is solution. The simulation results were obtained at t = 20 °C, pH = 8.0, and Cs = 700 mg/kg and various values of Is: 0.014, 0.07, 0.14, 0.28, 0.4, 0.8, 1.0, 1.2, 1.4, 7.1, 10.6, and 14.2 mol/kg. With an increase in Is from 0.014 to 0.8 mol/kg, the Cs(Tp) dependence shifts to the left and acquires a concave shape (Figure 7a,b). In the range of ionic strength values from 0.8 to 1.4 mol/kg, with an increase in Is, the rate of decrease of Cs(Tp) was lower and the curves shift to the right. The shape of the curves becomes convex. At values of Is > 1.4 mol/kg, the curves Cs (Tp) are shifted to the left again.
The dependence of the final average particle size on the ionic strength of the Is solution is non-monotonic also (Table 8, Figure 8). With an increase in Is from 0.014 to 0.8 mol/kg, Ra increases. Starting from the value Is = 0.8 mol/kg, the final average particle size decreases. This non-monotonic dependence is explained by the fact that with an increase in Is, the rate of molecular deposition of Rmd increases due to an increase in the degree of ionization of orthosilicic acid αSIL. Therefore, the rate of nucleation of OSA IN increases (Equation (5)). At the same time, the growth rate of the formed particles increases (Equation (10)). Thus, the final average size of silica particles is influenced by two competing factors—growth of IN and Rmd which determines the non-monotonic behavior of the Ra (Is) dependence.

3.4. Simulation of OSA Polymerization under Time-Variable Temperatures

A significant influence of the variable temperature regime on the rate of the OSA nucleation process and the final average size of silica particles was revealed.
To simulate the process of nucleation of OSA under variable temperature conditions, the following law of temperature t reduction with polymerization time Tp was chosen:
t = t 2 + ( t 1 t 2 ) · e T p taun 0
Dependence (1.34) corresponds to the rate of heat q, W/(m2∙s) exchange with an environment with a constant temperature to according to Newton’s law q = α∙F∙(t − t0), t1, t2 are the initial and final temperatures of the solution; α—coefficient of heat exchange W/(m2∙s∙°C); F—area of heat exchange, m2; taun0—the characteristic cooling time of the solution. Calculations were performed at t1 = 100 °C and t2 = 20 °C for solutions of the Mutnovsky hydrothermal field at Cs = 700 mg/kg and pH = 8.0; taun0 was taken to be 120, 240, 360, and 480 min.
Figure 9 shows that the smaller the tun0, the higher the level of supersaturation of SN. Consequently, the rate of OSA nucleation increases. This leads to a decrease in the final average particle size.
The effect of the initial concentration on the rate of the nucleation process and the final average size of silica particles in the variable temperature regime turned out to be ambiguous. The temperature was lowered from 100 to 20 °C, pH = 8.0, and taun0 = 360 min. The initial concentration assumed various constant values—500, 600, 700, 800, and 1000 mg/kg. Figure 10 shows a comparison of the dependencies obtained.
At variable temperatures, the supersaturation functions of time have a non-uniform behavior. The supersaturation of the solution SN = Cs/Ce depends on the solubility of Ce, which decreases, and faster than the value of the concentration of Cs.
Figure 9c and Figure 10c show a more symmetrical distribution, that is both branches of the silica particle size distribution are convex in contrast to the case of constant temperatures when the right was convex and the left was concave. Thus, the temperature regime at the stage of polymerization strongly influenced the index of polydispersity of hydrothermal sols and the density, morphology, agglomerate dimensions, volume, surface, diameters, and structures of pores of nanopowders produced from hydrothermal sols [5].

3.5. Comparison of the Experimental Data with Results of the Numerical Simulation

Experimental data Cs(Tp) obtained in sodium silicate solutions with Na+ cations removed by ion exchange and low ionic strength [26] were in agreement with the results of numerical simulation in the wide range of temperatures and initial OSA concentration Cs (Figure 11) [35]: t = 30, 90, 120,150, 180 °C, Cs = 383, 620, 110, 1200 mg/kg, pH = 7.0, 7.4, 7.85; Is = 0.000132–0.00134 mol/kg.
The comparison of the experimental data on OSA polymerization and the results of the numerical simulation for the solution of the Pauzhetsky hydrothermal field showed that the discrepancy of about 12–14.2% was within the uncertainty of the yellow molybdate method for determining the concentration of dissolved silicic acid for the following conditions [35]: t = 20 °C; Cs = 315.6–338 mg/kg; pH = 8.2–8.35; and Is = 0.04284 mol/kg.
The results of the simulation of OSA polymerization at the temperatures of 84–95 °C and high ionic strength for solutions with different chemical compositions showed an agreement with the experimental data obtained at Wairakei and Broadlands hydrothermal fields [35]: t = 84 °C, pH = 8.1, and Cs = 570 mg/kg (Wairakei); t = 95 °C, pH = 8.0, and Cs = 620 mg/kg (well 11 at Broadlands); t = 95 °C, pH = 7.7, and Cs = 900 mg/kg (well 22 at Broadlands); Is = 0.05–0.06 mol/kg.
The results of the simulation of the distribution of the number of SiO2 nanoparticles vs. radius for solutions in which OSA polymerized at the temperatures 33–43 °C, pH = 6.7–8.1, Cs = 366–546 mg/kg, and Is = 0.043 mol/kg corresponded to the results of particles radii measured using the DLS method [35,38] in the range of 50–225 nm.
The experimental results of [24,25] on the maximum nucleation rate vmax determined as (Cs1 − Cs2)/(Tp2 − Tp1) were in agreement with the results of the numerical simulation in the region of pH = 4–6, in which the model used [35] was applicable (pH = 4–8).
The comparison of the average diameter measured using helium chromatography in a synthetic solution with an initial Cs = 880 mg/kg at 20 °C and pH = 7.0 [39] were in good agreement with results of the numerical simulation [35].
The values of the k1 constant for concentration and supersaturation models obtained in [30] for synthetic solutions in the wide range of initial Cs 240–1200 mg/kg, ionic strength 0.01–0.24 mol/kg, and pH = 3–7 corresponded to the times Tp of Cs decrease calculated in our numerical simulations for the same conditions.
The experimental results [35] with Cs(Tp) curves characterized by kp constants of the OSA reaction of the polymerization at different temperatures 20–50 °C, different pH = 4–8.5, and ionic strengths 0.01–1.1 mol/kg corresponded to the models [8,28,33].

4. Simulation the Particle Size Distribution in Hydrothermal Sol Production Technology

The technological flow sheet for the solution of the Mutnovsky geothermal electric power plant GeoPP includes the well, separator, pipe silencer, tank for the solution cooling and aging to develop OSA polymerization, ultrafiltration membrane module, and volumes for concentrated sol, Figure 12a. The temperature decreases with the linear dependence in every element from the well to silencer from 300 to 96 °C at pH = 8.0 and with the dependence (1.33) in the aging tank in which the temperature decreases from 96 to 70 °C. When the duration of aging was 20 h and taun = 11 h, the average SiO2 nanoparticle radius calculated using the numerical simulation was 12.5 nm (Figure 12b) and it was 12.9 nm measured using the DLS method (Figure 13). When the duration of aging was 100 h and taun = 15 h, the final average diameters were 68.6 and 66 nm using the DLS method. If pH = 5.0, the temperature decrease during aging of 24 h was from 96 to 23 °C, taun = 12 h, and the calculated diameter was 160 nm and the diameter measured using the DLS method was 154 nm.
The temperature and aging duration can be varied to regulate the average particle diameter from 5 to 160 nm.
For the solution from wells at Cerro Prieto GeoPP with the values Cs = 950 mg/kg, pH = 7.3, Is = 0.4227 mol/kg, aging duration 65 min, taun = 60 min, and temperature decrease from 100 to 30 °C [23,24], the final average radius calculated in the present work was RA = 1.32 nm and it was 1.35 nm as a result of the calculation in [29].
For the solutions of the Wairakei GeoPP, New Zealand, at aging temperatures of 70 and 20 °C during 50 h, Cs = 490 mg/kg, and Is = 0.05 mol/kg, the calculated average radii were 60.7 and 9.5 nm, and the radii measured using DLS were 60 and 10 nm.

5. Conclusions

The mathematical model we adopted allowed us to investigate the influence of various physical and chemical factors in a wide range of values on the rate of polymerization of OSA in a hydrothermal solution and on the SiO2 nanoparticle size distribution. The results of the numerical simulation were verified using a comparison with experimental data. The temperature, pH, initial OSA concentration, and ionic strength are the main parameters that determine the kinetics of colloid phase formation, the final average size of SiO2 nanoparticles, and the form of the particle size distribution.
The final average radius decreases when the temperature decreases, and the degree of oversaturation becomes higher taking into account the surface tension, rate of molecular deposition, and Z-factor. The average radius increases when the pH decreases due to the increase in the surface tension and the critical size of the nucleus. The average radius decreases when the initial OSA concentration and degree of oversaturation increases.
The function of oversaturation SN depending on the time of polymerization Tp showed the long period of induction τin and homogeneous nucleation which becomes much longer at low values of pH. The period of homogeneous nucleation increases when pH decreases and the initial degree of oversaturation SN decreases at higher temperatures and low initial Cs concentrations. The relation of homogeneous and heterogenous durations (thom/thetg) increases when the temperature increases and pH decreases.
The dependence of the curve Cs vs. Tp and final average particle size Ra on the ionic strength Is is non-monotonic. With an increase in Is from 0.014 to 0.8 mol/kg, the kinetics of OSA polymerization accelerate, but Ra increases. Starting from the value of Is = 0.8 mol/kg, the final average particle size decreases. This non-monotonic dependence is explained by the fact that with an increase in Is, the rate of molecular deposition of Rmd increases due to an increase in the degree of ionization of orthosilicic acid αSIL. Therefore, the rate IN of OSA nucleation increases. At the same time, the growth rate of the formed particles increases. Thus, the final average size of silica particles is influenced by two competing factors—the increases in IN and Rmd which determine the non-monotonic behavior.
At variable temperatures, the supersaturation functions of time have a non-uniform behavior. The supersaturation of the solution SN = Cs/Ce depends on the solubility of Ce, which decreases, and faster than the value of the concentration of Cs. Both branches of the silica particle size distribution are concave in contrast to the case of constant temperatures. The variation in temperatures during the stage of OSA polymerization can influence the index of polydispersity of particles, its distribution regarding sizes, and the homogeneity and symmetrical form of the distribution. The temperature regime at the stage of polymerization highly influenced the index of polydispersity of hydrothermal sols and, thus, the density, morphology, agglomerate dimensions and fractal structures, volume, surface, diameters, and structures of pores of nanopowders produced from hydrothermal sols.
The results obtained can be used for regulating the kinetics of SiO2 nanoparticle formations, sizes, the polydispersity of the distribution of sizes, and the concentrations of nanoparticles in the technology of the hydrothermal synthesis of sol, gel, and powder.
The results can be applied to the technology of the production of synthetic water sols based on Na2SiO3 precursors and the technology of precipitated silica.

Author Contributions

Conceptualization, V.V.P. and A.A.C.; methodology, V.V.P. and A.A.C.; investigation, V.V.P., A.A.C. and D.S.G.; data curation, V.V.P. and A.A.C.; writing—original draft preparation, V.V.P., A.A.C. and D.S.G.; writing—review and editing, V.V.P. and D.S.G.; supervision, V.V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meng, L.; Duwal, S.; Lane, J.M.D.; Ao, T.; Stoltzfus, B.; Knudson, M.; Park, C.; Chow, P.; Xiao, Y.; Fan, H.; et al. Pressure Induced Assembly and Coalescence of Lead Chalcogenide Nanocrystals. J. Am. Chem. Soc. 2021, 143, 2688–2693. [Google Scholar] [CrossRef] [PubMed]
  2. Dubey, R.; Rajesh, Y.; More, M. Synthesis and Characterization of SiO2 Nanoparticles via Sol-gel Method for Industrial Applications. Mater. Today Proc. 2015, 2, 3575–3579. [Google Scholar] [CrossRef]
  3. Duan, H.; Wang, D.; Li, Y. Green chemistry for nanoparticle synthesis. Chem. Soc. Rev. 2015, 44, 5778–5792. [Google Scholar] [CrossRef] [PubMed]
  4. Potapov, V.; Fediuk, R.; Gorev, D. Obtaining sols, gels and mesoporous nanopowders of hydrothermal nanosilica. J. Sol-Gel Sci. Technol. 2020, 94, 681–694. [Google Scholar] [CrossRef]
  5. Potapov, V.; Fediuk, R.; Gorev, D. Hydrothermal SiO2 Nanopowders: Obtaining Them and Their Characteristics. Nanomaterials 2020, 10, 624. [Google Scholar] [CrossRef]
  6. Potapov, V.; Fediuk, R.; Gorev, D. Membrane concentration of hydrothermal SiO2 nanoparticles. Sep. Purif. Technol. 2020, 251, 117290. [Google Scholar] [CrossRef]
  7. Potapov, V.; Efimenko, Y.; Fediuk, R.; Gorev, D. Effect of hydrothermal nanosilica on the performances of cement concrete. Constr. Build. Mater. 2020, 269, 121307. [Google Scholar] [CrossRef]
  8. Iler, R.K. Polymerization of polycilicic acid derived from 3.3 ratio sodium silicate. J. Phys. Chem. 1953, 57, 604–607. [Google Scholar] [CrossRef]
  9. Alexander, G.B. The Reaction of Low Molecular Weight Silicic Acids with Molybdic Acid. J. Am. Chem. Soc. 1953, 75, 5655–5657. [Google Scholar] [CrossRef]
  10. Alexander, G.B. The Polymerization of Monosilicic Acid. J. Am. Chem. Soc. 1954, 76, 2094–2096. [Google Scholar] [CrossRef]
  11. Schwarz, R.; Knauff, K.G. Zur Kenntnis der Kieselsäuren. X. Über Alkoxysilane und Oligokieselsäuren. Z. Für Anorg. Und Allg. Chem. 1954, 275, 176–192. [Google Scholar] [CrossRef]
  12. Bechtold, M.F. Polymerization and Properties of Dilute Aqueous Silicic Acid from Cation Exchange. J. Phys. Chem. 1955, 59, 532–541. [Google Scholar] [CrossRef]
  13. Goto, K.J. States of silica in aqueous solutions. II. Chem. Soc. Jap. Pure Chem. Sect. 1955, 76, 1364–1366. [Google Scholar]
  14. Hoebbel, D.; Wieker, W. Über Kondensationsreaktionen der Monokieselsäure. Z. Für Anorg. Und Allg. Chem. 1973, 400, 148–160. [Google Scholar] [CrossRef]
  15. Hoebbel, D.; Wieker, W. Über die dünnschichtchromatographische Trennung trimethylsilylierter Kieselsäuren. Z. Anorg. Allgem. Chem. 1974, 405, 163–166. [Google Scholar] [CrossRef]
  16. Hoebbel, D.; Wieker, W.; Hoebbel, D.; Wieker, W. 29Si-NMR-Spektroskopie an Silicatlösungen. IV. Untersuchungen zur Kondensation der Monokieselsäure. Z. Für Anorg. Und Allg. Chem. 1977, 428, 43–52. [Google Scholar]
  17. Baumann, H. Polymerization and depolymerization de kieselsaure unter verschiedenen begingungen. Kolloid-Z. 1959, 162, 28–35. [Google Scholar] [CrossRef]
  18. Ginzburg, F.L.; Sheidina, L.D. Radiochemistry. 1973, 15, 410. Available online: http://sciencejournals.ru/journal/radkhim/ (accessed on 18 September 2022).
  19. Goto, K. Effect of pH on Polymerization of Silicic Acid. J. Phys. Chem. 1956, 60, 1007–1008. [Google Scholar] [CrossRef]
  20. Okkerse, C. Submicroporous and Macroporous Silica. Ph.D. Thesis, Delftsche Uitgevers Maatschappif N.V., Delft, The Netherlands, 1961. [Google Scholar]
  21. Richardson, E.; Waddams, J.A. Use of the silico-molybdate reaction to investigate the polymerization of low molecular weight silicic acids in dilute solutions. Res. Corresp. 1954, 7, 43. [Google Scholar]
  22. Tarutani, T. Chromatographic behavior of silicic acid on sephadex columns. J. Chromatogr. 1970, 50, 523–526. [Google Scholar] [CrossRef]
  23. Makrides, A.K.; Turner, M.; Slaughter, J. Condensation of silica from super-saturated silicic acid solutions. J. Colloid Interface Sci. 1980, 73, 345–367. [Google Scholar] [CrossRef]
  24. Frolov, Y.G.; Shabanova, N.A.; Popov, V.V. Influence of temperature and pH on polycondensation of silicic acid in an aqueous medium. Colloid. J. 1983, 45, 179–182. [Google Scholar]
  25. Frolov, Y.G.; Shabanova, N.A.; Popov, V.V. Polycondensation silicic acid in an aqueous medium. Influence of the concentration of silicic acid. Colloid. J. 1983, 45, 382–386. [Google Scholar]
  26. Rothbaum, H.; Rohde, A. Kinetics of silica polymerization and deposition from dilute solutions between 5 and 180 °C. J. Colloid Interface Sci. 1979, 71, 533–559. [Google Scholar] [CrossRef]
  27. Coudurier, M.; Baudru, R.; Donnet, J.B. Etude de la polycondensation de l’acide disilique. Bull. Soc. Chem. Fr. 1971, 9, 3154–3160. [Google Scholar]
  28. Fleming, B.A. Kinetics of reaction between silicic acid and amorphous silica surfaces in NaCl solutions. J. Colloid Interface Sci. 1986, 110, 40–64. [Google Scholar] [CrossRef]
  29. Weres, O.; Yee, A.; Tsao, L. Kinetics of silica polymerization. J. Colloid Interface Sci. 1981, 84, 379–402. [Google Scholar] [CrossRef]
  30. Conrad, C.F.; Icopini, G.A.; Yasuhara, H.; Bandstra, J.Z.; Brantley, S.L.; Heaney, P.J. Modeling the kinetics of silica nanocolloid formation and precipitation in geologically relevant aqueous solutions. Geochim. Et Cosmochim. Acta 2007, 71, 531–542. [Google Scholar] [CrossRef]
  31. Garofalini, S.H.; Martin, G. Molecular Simulations of the Polymerization of Silicic Acid Molecules and Network Formation. J. Phys. Chem. 1994, 98, 1311–1316. [Google Scholar] [CrossRef]
  32. Rao, N.Z.; Gelb, L.D. Molecular Dynamics Simulations of the Polymerization of Aqueous Silicic Acid and Analysis of the Effects of Concentration on Silica Polymorph Distributions, Growth Mechanisms, and Reaction Kinetics. J. Phys. Chem. B 2004, 108, 12418–12428. [Google Scholar] [CrossRef]
  33. Greenberg, S.A. Polymerization of silicic acid in alkaline solutions a kinetics study. J. Polym. Sci. 1958, 27, 523–527. [Google Scholar] [CrossRef]
  34. Potapov, V.V.; Serdan, A.A.; Kashpura, V.N.; Gorbach, V.A.; Tyurina, N.A.; Zubakha, S.V. Polycondensation of orthosilicic acid in hydrothermal solutions at different temperatures, pH values, and ionic strengths. Glas. Phys. Chem. 2007, 33, 44–49. [Google Scholar] [CrossRef]
  35. Potapov, V.V.; Cerdan, A.A.; Kashutina, I.A. Numerical Simulation of the Polycondensation of Orthosilicic Acid and of the Formation of Silica Particles in Hydrothermal Solutions. J. Volcanol. Seism. 2019, 13, 216–225. [Google Scholar] [CrossRef]
  36. Crerar, D.A.; Anderson, G.M. Solubility and solvation reactions of quartz in dilute hydrothermal solutions. Chem. Geol. 1971, 8, 107–122. [Google Scholar] [CrossRef]
  37. Marshall, W.L. Amorphous Silica Solubilities: I. Behaviour in Aqueous Sodium Nitrate Solutions: 25–300 °C, 0–6 Molal. Geochim. Cosmochim. Acta 1980, 44, 907–913. [Google Scholar] [CrossRef]
  38. Ohsawa, S.; Kawamura, T.; Nakamatsu, N.; Yusa, Y. Geothermal blue water colored by colloidal silica. In Proceedings of the World Geothermal Congress, Kyushu, Tohoku, Japan, 28 May–10 June 2000; pp. 663–668. [Google Scholar]
  39. Crerar, D.A.; Axtmann, E.V.; Axtmann, R.C. Growth and ripening of silica polymers in aqueous solutions. Geochim. Cosmochim. Acta 1981, 45, 1259–1266. [Google Scholar] [CrossRef]
Figure 1. Dependence of the rate of molecular deposition Rmd (g/(cm2∙min)) on the temperature t (°C) and concentration Cs (mg/kg). (a) pH = 2.0; (b) pH = 4.0; (c) pH = 6.0; (d) pH = 8.0. Cs = 300, 400, 500, 600, 700, 800, 900, 100, 110, and 1200 mg/kg.
Figure 1. Dependence of the rate of molecular deposition Rmd (g/(cm2∙min)) on the temperature t (°C) and concentration Cs (mg/kg). (a) pH = 2.0; (b) pH = 4.0; (c) pH = 6.0; (d) pH = 8.0. Cs = 300, 400, 500, 600, 700, 800, 900, 100, 110, and 1200 mg/kg.
Polymers 14 04044 g001aPolymers 14 04044 g001bPolymers 14 04044 g001c
Figure 2. The temperature t (°C) corresponded to the maximum rate of molecular deposition Rmd at different concentrations Cs, pH = 8.0.
Figure 2. The temperature t (°C) corresponded to the maximum rate of molecular deposition Rmd at different concentrations Cs, pH = 8.0.
Polymers 14 04044 g002
Figure 3. Dependence of the surface tension σsw (mJ/m2) on (a) temperature t (°C) and (b) on pH. Initial concentration Cs = 700 mg/kg.
Figure 3. Dependence of the surface tension σsw (mJ/m2) on (a) temperature t (°C) and (b) on pH. Initial concentration Cs = 700 mg/kg.
Polymers 14 04044 g003
Figure 4. The number of new particles PNS (nucl/kg) versus time Tp; T = 50 °C, the initial concentration is Cs = 700 mg/kg, pH = 8.
Figure 4. The number of new particles PNS (nucl/kg) versus time Tp; T = 50 °C, the initial concentration is Cs = 700 mg/kg, pH = 8.
Polymers 14 04044 g004
Figure 5. The rate of nucleation Rnuc (nucl/(kg∙s)) versus time Tp; T = 50 °C, the initial concentration is Cs = 700 mg/kg, pH = 8.
Figure 5. The rate of nucleation Rnuc (nucl/(kg∙s)) versus time Tp; T = 50 °C, the initial concentration is Cs = 700 mg/kg, pH = 8.
Polymers 14 04044 g005
Figure 6. The particle size distribution at different pH and initial OSA concentrations Cs. (a) T = 50 °C, Cs = 700 mg/kg, pH = 4; (b) T = 50 °C, Cs = 700 mg/kg, pH = 8; (c) T = 60 °C, Cs = 800 mg/kg, pH = 8. PNS—total quantity of the formed SiO2 particles.
Figure 6. The particle size distribution at different pH and initial OSA concentrations Cs. (a) T = 50 °C, Cs = 700 mg/kg, pH = 4; (b) T = 50 °C, Cs = 700 mg/kg, pH = 8; (c) T = 60 °C, Cs = 800 mg/kg, pH = 8. PNS—total quantity of the formed SiO2 particles.
Polymers 14 04044 g006
Figure 7. The dependence of the concentration of dissolved silicic acid (a) and supersaturation (b) on time. t = 20 °C, pH = 8.0, initial concentration Cs = 700 mg/kg, and various values of ionic strength (mol/kg).
Figure 7. The dependence of the concentration of dissolved silicic acid (a) and supersaturation (b) on time. t = 20 °C, pH = 8.0, initial concentration Cs = 700 mg/kg, and various values of ionic strength (mol/kg).
Polymers 14 04044 g007
Figure 8. Particle size distribution with different ionic strength values. t = 20 °C, pH = 8.0, and initial concentration Cs = 700 mg/kg, Is (mol/kg).
Figure 8. Particle size distribution with different ionic strength values. t = 20 °C, pH = 8.0, and initial concentration Cs = 700 mg/kg, Is (mol/kg).
Polymers 14 04044 g008
Figure 9. The time dependences of the concentration of dissolved silicic acid Cs (a), supersaturation (b), and particle size distribution (c). The temperature decreases from 100 to 20 °C, pH = 8.0, and Cs = 700 mg/kg.
Figure 9. The time dependences of the concentration of dissolved silicic acid Cs (a), supersaturation (b), and particle size distribution (c). The temperature decreases from 100 to 20 °C, pH = 8.0, and Cs = 700 mg/kg.
Polymers 14 04044 g009aPolymers 14 04044 g009b
Figure 10. Dependence of Cs on Tp (a), SN on Tp (b), and particle size distribution (c). The temperature decreases from 100 to 20 °C and taun0 = 360 min.
Figure 10. Dependence of Cs on Tp (a), SN on Tp (b), and particle size distribution (c). The temperature decreases from 100 to 20 °C and taun0 = 360 min.
Polymers 14 04044 g010aPolymers 14 04044 g010b
Figure 11. Comparing the experimental data [19] with the results of the numerical simulation. (a) t = 30 °C, Cs = 383 mg/kg, pH 7.0; (b) t = 150 °C, Cs = 1008 mg/kg, pH 7.85. Is = 0.000132–0.00134 mol/kg. ▬—results of the simulation; ▲—experimental data.
Figure 11. Comparing the experimental data [19] with the results of the numerical simulation. (a) t = 30 °C, Cs = 383 mg/kg, pH 7.0; (b) t = 150 °C, Cs = 1008 mg/kg, pH 7.85. Is = 0.000132–0.00134 mol/kg. ▬—results of the simulation; ▲—experimental data.
Polymers 14 04044 g011
Figure 12. Parameters involved in the technological flow sheet of silica extraction. (a) Temperature versus time; (b) average silica particle radius as a function of time.
Figure 12. Parameters involved in the technological flow sheet of silica extraction. (a) Temperature versus time; (b) average silica particle radius as a function of time.
Polymers 14 04044 g012
Figure 13. The particle size distribution based on DLS results. d is the particle diameter, nm. Well 054 at the Mutnovsky field.
Figure 13. The particle size distribution based on DLS results. d is the particle diameter, nm. Well 054 at the Mutnovsky field.
Polymers 14 04044 g013
Table 1. Characteristics of the hydrothermal solution from the Mutnovsky field, ionic strength Is = 14.218 mmol/kg, Is = ∑cj∙zj2/2, cj, zj (concentration (mol/kg) and charge of the j-th ion in solution), (zNa = 1, zK = 1, zCa = 2, zMg = 2, zAl = 3), pH = 9.35, elec trical conductivity σel = 1.1–1.3 mS/cm. nc—not calculated.
Table 1. Characteristics of the hydrothermal solution from the Mutnovsky field, ionic strength Is = 14.218 mmol/kg, Is = ∑cj∙zj2/2, cj, zj (concentration (mol/kg) and charge of the j-th ion in solution), (zNa = 1, zK = 1, zCa = 2, zMg = 2, zAl = 3), pH = 9.35, elec trical conductivity σel = 1.1–1.3 mS/cm. nc—not calculated.
ComponentConcentration
mg/Lmg equ/L
Na+239.410.413
K+42.01.074
Ca2+1.60.0798
Mg2+0.720.0592
Fe2, 3+<0.1<0.0053
Al3+0.270.033
NH4+1.10.0609
Li+0.710.102
Total cation concentration285.911.827
Cl198.55.591
HCO381.01.327
CO32−19.90.663
SO42−192.13.9995
HS4.950.15
H2S05.92-
Fn/dn/d
Total anion concentration496.511.73
H3BO3106.9nc
(H4SiO4)t1190nc
(H4SiO4)s222nc
Mineralization Mh1638.9nc
Table 2. Values of function f at different pH.
Table 2. Values of function f at different pH.
pHf(pH)pHf(pH)pHf(pH)
4.00.000255.90.0177.50.21
4.50.000796.00.0228.00.32
5.00.00246.50.0558.50.44
5.50.00757.00.1198.70.49
Table 3. Values of function i at different pH.
Table 3. Values of function i at different pH.
pHi (pH)
6.00.01
6.50.03
7.00.07
7.50.15
8.00.29
8.50.48
8.70.56
Table 4. Dependence of Zeldovich factor Z on the temperature and pH (-).
Table 4. Dependence of Zeldovich factor Z on the temperature and pH (-).
t, °CpH
4.04.55.05.56.06.5
203.8∙10−23.8∙10−23.8∙10−23.8∙10−23.8∙10−23.9∙10−2
502.1∙10−22.1∙10−22.1∙10−22.1∙10−22.2∙10−22.2∙10−2
1005.8∙10−3-5.8∙10−35.8∙10−36.0∙10−36.1∙10−3
1301.4∙10−31.3∙10−31.3∙10−31.3∙10−31.4∙10−31.4∙10−3
t, °CpH
7.07.58.08.58.78.99
204.2∙10−24.8∙10−26.4∙10−21.1∙10−11.5∙10−13.3∙10−1
502.4∙10−22.8∙10−24∙10−27.6∙10−21.2∙10−13.5∙10−1
1006.7∙10−38.2∙10−31.2∙10−23.1∙10−26.2∙10−28.6∙10−3
1301.6∙10−32∙10−33.3∙10−31∙10−22.6∙10−4-
Table 5. The final average radius of silica particles at a constant temperature and varying pH (nm). Cs(t = 0) = Ct = 700 mg/kg.
Table 5. The final average radius of silica particles at a constant temperature and varying pH (nm). Cs(t = 0) = Ct = 700 mg/kg.
pHt, °C
2040506080100
45.2819.2448.06156.556342.1-
55.2919.0547.30153.155858.7-
62.9816.8240.69128.014684.7-
72.977.6615.2838.17701.4-
80.981.211.501.773.4818.7
8.50.640.730.780.841.011.26
Table 6. Values of the final mean size of silica particles (nm) under different temperatures t and initial concentrations of solution Cs, pH = 8; (-) means that it was not estimated.
Table 6. Values of the final mean size of silica particles (nm) under different temperatures t and initial concentrations of solution Cs, pH = 8; (-) means that it was not estimated.
CS,
mg/kg
t, °C
20405060
30093.5---
35012.41019.9--
4004.3549.1586.3-
5001.794.088.9632.7
6001.2461.84-3.80
7000.981.21-1.77
8000.871.05-1.31
10000.730.83-0.95
CS,
mg/kg
t, °C
70808590
300----
350----
400----
500----
6007.5624.1456.14169.8
700-3.48-6.49
800-1.81-2.30
1000-1.10-1.21
Table 7. Durations of homogeneous thom and heterogeneous thetg nucleation at different temperatures, pH, Cs.
Table 7. Durations of homogeneous thom and heterogeneous thetg nucleation at different temperatures, pH, Cs.
Temperature, °CpHInitial OSA Concentration Cs, mg/kgthom. minthetg, minthom/thetg
20570080,385.34,369,1950.018398
2067008255.45454,617.50.018159
207700749.6381,119.240.009241
20870050.205133.770.009779
505700118,281.13,027,5790.039068
50670011,123.75250,812.40.044351
507700715.8720,328.970.035215
50870013.47776.680.01735
8057004,352,9387.15·1070.0608
806700373,867.65,574,7100.067065
80770010,519.41143,413.70.07335
80870019.635328.120.05984
10057007.07·1095.32·10100.133
10067004.89·1083.53·1090.139
100770032,828.602.19·1070.15
1008700116.61638.400.182
12057001.08·10194.94·10190.219
12067003.93·10171.78·10180.22
12077004.71·10131.87·10140.252
120870026,324.9273,435.340.3584
Table 8. The final average size of silica particles in solutions with different ionic strength values. t = 20 °C, pH = 8.0, and Cs = 700 mg/kg.
Table 8. The final average size of silica particles in solutions with different ionic strength values. t = 20 °C, pH = 8.0, and Cs = 700 mg/kg.
IS, mol/kgRa, nmIS, mol/kgRa, nmIS, mol/kgRa, nm
0.0140.990.41.0481.41.042
0.071.0290.81.0497.10.9
0.141.0361.01.04710.60.79
0.281.0441021.04514.20.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Potapov, V.V.; Cerdan, A.A.; Gorev, D.S. Silicic Acid Polymerization and SiO2 Nanoparticle Growth in Hydrothermal Solution. Polymers 2022, 14, 4044. https://doi.org/10.3390/polym14194044

AMA Style

Potapov VV, Cerdan AA, Gorev DS. Silicic Acid Polymerization and SiO2 Nanoparticle Growth in Hydrothermal Solution. Polymers. 2022; 14(19):4044. https://doi.org/10.3390/polym14194044

Chicago/Turabian Style

Potapov, Vadim V., Angel A. Cerdan, and Denis S. Gorev. 2022. "Silicic Acid Polymerization and SiO2 Nanoparticle Growth in Hydrothermal Solution" Polymers 14, no. 19: 4044. https://doi.org/10.3390/polym14194044

APA Style

Potapov, V. V., Cerdan, A. A., & Gorev, D. S. (2022). Silicic Acid Polymerization and SiO2 Nanoparticle Growth in Hydrothermal Solution. Polymers, 14(19), 4044. https://doi.org/10.3390/polym14194044

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