Theoretical Investigation of the Hydrolytic Mechanism of α-Functionalized Alkoxysilanes as Effective Crosslinkers and the Difficulty of Deep Vulcanization in RTV Silicone Rubber.

The reactions between α-, γ-ethylenediaminemethyl trimethyl-ketoxime silane (α-, γ-EAMOS) and H₂O were investigated on the geometries of stationary points, the reaction pathway (IRC), thermodynamic and kinetic analysis by density functional theory (DFT) at the B3LYP/6-311G (d, p) level. Interestingly, the results showed that the hydrolysis activity of α-EAMOS is higher than that of γ-EAMOS, due to the influence of an amino substituent in position α-C on silicon. α-EAMOS can be used as a superior crosslinker for room temperature vulcanized (RTV) silicone rubber to achieve rapid crosslinking without a toxic catalyst. Besides, compared with the reaction between α-EAMOS and H₂O, the reactivity between α-EAMOS and hydroxy siloxane (HO⁻Si(CH₃)₂⁻OSiH₃) was discussed. Particularly, it revealed that the deep vulcanization of RTV silicone rubber is difficult.

1. Introduction

Carbofunctionalized trialkoxysilanes (X–(CH2)–Si(OR)3, where X is functional amines, chlorine, carboxyl etc., n = 1, 2, 3 is called α-, β-, γ-functionalized alkoxysilane, respectively) play a great role in organosilicon chemistry. They are widely used as coupling agents to connect inorganic and organic materials possessing various structures and properties [1,2]. And they are also used as crosslinkers for room temperature vulcanized (RTV) silicone rubber. The OR groups of trialkoxysilanes react with H2O to form silanols, which then condense with the polydimethylsiloxane (PDMS) to produce an elastomer with three-dimensional network structure [3,4]. The corresponding RTV silicone rubbers with excellent properties are widely used in sealants, adhesives, coatings, and spray foams [5,6,7].

Commercial crosslinkers used on a large scale have often been γ-functionalized alkoxysilanes (n = 3) because they are considered more stable than α- or β-functionalized alkoxysilanes (n = 1 or 2, β-effect makes β-functionalized alkoxysilanes unstable) [8,9]. However, the influence of the functional group on Si–OR of γ-functionalized alkoxysilanes is much smaller than that of α-functionalized alkoxysilanes because the γ-group lies far away from Si. When it comes to crosslinkers for silicone rubber, it is usually necessary to add a toxic catalyst, such as organotin, to complete the rapid vulcanization at room temperature.

Now, some α-silane-terminated prepolymers crosslink faster than the corresponding prepolymers produced from conventional γ-functionalized alkoxysilanes [10]. So, α-functionalized alkoxysilanes get more attention in silicon chemistry. For example, chloromethyltriethoxysilane was used as a self-catalytic crosslinker in water-based silicone rubber, which then had low modulus and high elongation [11]. The use of α-amine alkoxysilanes bearing piperazidine improved the stability of silicone rubber [12]. Meanwhile, various α-amine functionalized alkoxysilanes as crosslinkers accelerated the curing rate of the silicone rubber [13]. Just now, a series of α-amine ketoximesilanes [14] were prepared to be used as auto-catalyzed crosslinkers in RTV silicone rubber [15]. Thereby, it was simply proved by using density functional theory (DFT) calculations that α-amine ketoximesilanes have fast hydrolysis rates [16].

Although there are many studies on α-functionalized alkoxysilanes, the effect of the α-functional amino group is not fully understood yet. The geometries of α-, γ-functionalized silanes and their reaction processes with H2O need further investigation to find the hydrolytic mechanism, which is related to the vulcanization of RTV silicone rubber. Besides, as for the deep vulcanization of RTV silicone rubber, it is that the crosslinker first transforms to silanol (Si–OH). When there is rare water involved, the reactions between the crosslinker and H2O and between the crosslinker and HO–PDMS are competitive. Particularly, the deep vulcanizing process of RTV silicone rubber is difficult [17,18,19]. But, there are rare studies accounting for this.

In this regard, we studied the reactions between α-, γ-ethylenediaminemethyl trimethyl-ketoxime silane (EAMOS) and H2O, respectively. They were simulated by DFT calculations for the geometries of stationary points, and thermodynamic and kinetic analysis. To reveal the difficulty of deep vulcanization of RTV silicone rubber, the reaction between α-EAMOS and H2O and the reaction between α-EAMOS and HO–Si(CH3)2–OSiH3 were compared.

2. Computational Methods

Here, the computational methods referred to the default convergence limit of B3LYP/6-311G (d, p) of DFT in Gaussian 03 (C. 02, Jinan, China) [20]. DFT calculations were carried out using Becke’s three-parameter hybrid function (B3LYP) level using the 6-311G (d, p) basis set. The geometry optimization was applied for the ground states and these ground states were assumed to be a singlet state. The theoretical calculations were performed by this quantum set. All the above-mentioned calculations for geometry optimization were carried out with the Gaussian 03 quantum chemistry program-package.

The reactions between α-, γ-EAMOS and H2O were studied by using DFT at the B3LYP/6-311G (d, p) level. The molecules were optimized. The geometries of various stationary points and harmonic vibrational frequencies were calculated. Only one imaginary frequency for each transition state could be found. The reaction pathways were investigated and confirmed by intrinsic reaction coordinate (IRC) calculations (Figure 1). The changes of thermodynamic functions including entropy, enthalpy, free energies, and equilibrium constant were calculated through reactants, transition states, and products in Eyring transition state theory. Zero-point energy correction was also carried out at the B3LYP/6-311G (d, p) level.

Figure 1

Geometries of stationary points of the reaction between α-ethylenediaminemethyl trimethyl-ketoxime silane (EAMOS) and H2O at B3LYP/6-311G** level.

Having tried different reaction modes, the reaction path here needs the lowest energy. In order to guarantee the reality of the reaction path, the structures of transition states (the imaginary frequency is 1) were confirmed by vibration analysis. IRC of all the structures of transition states were calculated (Figure S1). The optimization was further performed along both sides of the IRC curve. One is reactants, and the other is products, which proved that the reaction path was reliable.

3. Results and Discussion

3.1. Stationary Points and Pathways with α-EAMOS and H2O as an Example

IRC calculation (Figure 1, Tables S1–S3) showed that in the reaction between α-EAMOS and H2O, the first step is the strong electronegative 32O (NBO charge: −0.918) of H2O, coordinated with two hydrogens (30H, 33H), moving to the electropositive 1Si (distance of 1Si–32O decreased from 2.038 Å to 1.803 Å, NBO charge: 2.156). Then, 32O and 1Si form a five-coordination intermediate by using the 3d orbital of Si. Afterwards, 3O (NBO charge: −0.732) leaves from 1Si (the 3O–1Si bond length increases from 1.710 Å to 2.065 Å), and 30H leaves from 32O (the 30H–32O bond length increases from 0.985 Å to 1.099 Å). This process crosses the barrier potential (78.8 kJ·mol−1), resulting in the transition state (TS-α-C). After that, 32O keeps close to 1Si, (the 32O–1Si bond length increases from 1.803 Å to 1.661 Å) resulting in an Si–O bond. 30H (NBO charge: 0.510) leaves from 32O and keeps close to 3O, (the 30H–32O bond length increases from 1.099 Å to1.814 Å and the 30H–3O bond length decreases from 1.368 Å to 0.975 Å). Finally, the product P-α-C is obtained. Here, TS-α-C was confirmed as first order saddle point by IRC analysis and the only imaginary frequency of transition is −841.9 HZ.

Besides, the reaction between γ-EAMOS and H2O follows a similar pathway shown in Figure S2 and Tables S4–S6.

3.2. The Highest Occupancy Orbit (HOMO) of Each Stationary Point of Reactions between α-, γ-EAMOS and H2O

The HOMO of the stationary points of the reactions between α-, γ-EAMOS and H2O in Figure 2 and Figure 3 were compared. In Re-α-C (H2N(CH2)2NHCH2Si(O–N=CEtMe)3), when H in the α position is replaced by an electronegative ethylenediamine group, ‘p’ electrons in N enter the 3d orbital of Si, resulting in p–π conjugation. Thus, the electronegativity of Si increases; it is easier to be attacked by nucleophiles. While electronegative OH in H2O comes close to Si to form a complex, the Si–O bond will break easily and separate.

Figure 2

Highest occupied molecular orbital (HOMO) of stationary points of the reaction between α-EAMOS and H2O at the B3LYP/6-311G** level.

Figure 3

HOMO of stationary points of the reaction between γ-EAMOS and H2O at the B3LYP/6-311G** level.

However, in Re-γ-C (H2N(CH2)2NH–CH2CH2CH2Si(O–N=CEtMe)3), the electron cloud of N will not conjugate with Si because the group lies far away from the Si atom. Thus, there is rare electron interaction on Si. As a result, the reaction between α-EAMOS and H2O is easier than the reaction between γ-EAMOS and H2O.

3.3. The Reaction between α-EAMOS and HO–Si(CH3)2–OSiH3

IRC calculation (Figure 4, Tables S7–S9) showed that 32O in HO–Si(CH3)2–OSiH3 comes close to 1Si (distance of 1Si–32O from 4.127 Å to 1.827 Å). Then, 32O and 1Si share a five-coordination bond by using the 3d orbital of Si. Afterwards, 2O leaves from 1Si, (the 1Si–2O bond length increases from 1.684 Å to 2.035 Å), and 30H leaves from 32O, (the 30H–32O bond length increases from 0.969 Å to 1.076 Å). This process crosses the barrier potential, 118.9 kJ·mol–1, resulting the transition state (TS-Si–OH). After that, 32O keeps close to 1Si (distance of 32O–1Si decreases from 1.827 Å to 1.611 Å), resulting in a Si–O bond. 30H leaves from 32O and keeps close to 2O (the 30H–32O bond length increases from 1.076 Å to 3.004 Å and the 30H–2O bond decreases from 1.387 Å to 0.995 Å). Finally, the product P–Si–OH is obtained. Here, TS-Si–OH was confirmed as the first order saddle point by IRC analysis and the only imaginary frequency of transition is −468.9 HZ.

Figure 4

Geometries of stationary points of the reaction between α-EAMOS and HO–Si(CH3)2–OSiH3 at the B3LYP/6-311G** level.

3.4. Thermodynamic and Kinetic Parameters

The entropy changes from reactant to transition state were calculated using vibration frequency. They were based on the equilibrium geometries of the transition states (TS-α, TS-γ-C, TS-Si–OH) and the products (P-α, P-γ-C, P–Si–OH). Then, the reaction rates were obtained by transition state theory, while entropy changes, enthalpy changes, Gibbs’ free-energy changes, and the equilibrium constants were calculated by using statistical thermodynamic method. Kinetic analysis further proved the possibility and practicality of these reactions. For calculating rate constants, some insertion points were selected in both sides of the transition states along the steepest descent energy curve calculated by IRC. Then, the Hessian matrices of these points, reactants, products, and the transition states were calculated.

Rate constant: k (T) = g(kbT/h) exp (ΔS#/R − ΔH#/RT)

g: Wigner correction factor, g = 1 + (hν#/kbT)2/24;

kb: Boltzmann constant (1.3806505(24) × 10−23 J·K−1);

T: thermodynamic temperature (K);

h: Planck constant (6.6260693(11) × 10−34 J·s);

R: molar gas constant (8.31441 ± 0.00026 J·mol−1·K−1);

V#: imaginary frequency of transition (cm3·mol−1);

ΔS#: standard molar activation entropy (J·mol−1·K−1);

ΔH#: standard molar activation enthalpy (kJ·mol−1).

Equilibrium constant: K (T) = exp (−ΔG × 1000/ (8.314T))

ΔG: Gibbs free energy (kJ·mol−1).

As shown in the kinetic calculations in Table 1 and Table 2, the rate constants will increase and the equilibrium constants will decrease gradually as the temperature increases. Interestingly, the reaction between α-EAMOS and H2O has high rate constants and equilibrium constants. In detail, the reaction between α-EAMOS and H2O has the highest rate constant which is three times that of γ-EAMOS and H2O, and ten orders of magnitude higher than that of α-EAMOS and HO–Si(CH3)2–OSiH3 at 298 K. This difference will decrease gradually as the temperature increases.

Table 1

Rate Constants k (T) calculated at the B3LYP/6-311G (d, p) level.

Rate Constants k(T) (s−1)
T/K	α-C	γ-C	Si–OH
298	3.23 × 10−1	1.09 × 10−1	4.80 × 10−11
400	1.48 × 103	3.66 × 102	8.81 × 10−6
500	2.10 × 105	4.12 × 104	1.04 × 10−2
600	5.79 × 106	9.45 × 105	1.16 × 100
700	6.29 × 107	8.76 × 106	3.35 × 101
800	3.79 × 108	4.61 × 107	4.19 × 102
900	1.54 × 109	1.66 × 108	2.99 × 103
1000	4.78 × 109	4.63 × 108	1.45 × 104
1100	1.21 × 1010	1.06 × 109	5.26 × 104
1200	2.63 × 1010	2.12 × 109	1.55 × 105
1300	5.09 × 1010	3.78 × 109	3.86 × 105
1400	8.98 × 1010	6.19 × 109	8.45 × 105
1500	1.47 × 1011	9.46 × 109	1.67 × 106

Table 2

Equilibrium Constants K (T) calculated at the B3LYP/6-311G (d, p) level.

Equilibrium Constants K(T)
T/K	α-C	γ-C	Si–OH
298	1.04 × 1014	5.79 × 1013	9.28 × 1011
400	1.71 × 1011	1.06 × 1011	4.56 × 108
500	4.12 × 109	2.73 × 109	5.10 × 106
600	3.48 × 108	2.41 × 108	2.49 × 105
700	6.01 × 107	4.26 × 107	2.84 × 104
800	1.61 × 107	1.17 × 107	5.50 × 103
900	5.80 × 106	4.26 × 106	1.53 × 103
1000	2.56 × 106	1.91 × 106	5.45 × 102
1100	1.32 × 106	9.87 × 105	2.34 × 102
1200	7.55 × 105	5.71 × 105	1.15 × 102
1300	4.71 × 105	3.59× 105	6.34 × 101
1400	3.15 × 105	2.41 × 105	3.79 × 101
1500	2.22 × 105	1.71 × 105	2.43 × 101

As the thermodynamic calculations in Table S10 show, the reaction between α-EAMOS and H2O is entropy-increasing, exothermic, and spontaneous. The reaction between α-EAMOS and HO–Si(CH3)2–OSiH3 is exothermic, spontaneous, but entropy-decreasing. This is due to the polyreaction of two reactants from disorder to order. Besides, the reactions between α-EAMOS and H2O exhibit similar aviation energy (78.8 kJ·mol−1, Tables S2 and S5), while that of α-EAMOS and HO–Si(CH3)2–OSiH3 (118.9 kJ·mol−1, Table S8) is higher.

In conclusion, α-functionalized silanes have a fast hydrolysis rate. Using α-functionalized ketoxime silanes as effective autocatalytic crosslinkers in RTV silicone rubber by decreasing the tack-free time and improving thermal stability (Figure S3) [15] or mechanical properties has been reported [16,21,22]. Besides, α-functionalized silanes are easier to react with H2O on the surface of silicone rubber but difficult to react with HO-PDMS (silicone rubber gum) in the vulcanizing system. It proved that the corresponding silicone rubber is difficult to deep vulcanize at room temperature.

4. Conclusions

According to DFT calculations, the amine group on α-EAMOS is only one C atom away from Si, and the lone pair of electrons on the amine can be conjugated into the 3d orbital of Si, thus increasing the electronegativity of Si. Then, the hydrolyzable group on Si was attacked easily by H2O. But, there is no conjugation effect in γ-EAMOS, whose hydrolysis rate was lower than that of α-EAMOS. Interestingly, it proved that α-amino substituted alkoxy silanes could be used as efficient crosslinkers in RTV silicone rubber.

Besides, the reaction between α-EAMOS and H2O had a faster reaction rate, lower activation energy, and higher reaction trend than those of the reaction between α-EAMOS and HO–Si(CH3)2–OSiH3. So, α-functionalized ketoxime silanes could react with H2O quickly on the surface of silicone rubber but they are difficult to react with HO-PDMS. Then, there is rare cross-linking effect inside. As a result, the deep vulcanization of the corresponding silicone rubber is difficult.

Here, DFT calculations were explanative and markedly consistent with our experimental results. Not only is this study important for basic research, it may also facilitate an increase in the use of α-functionalized ketoxime silanes to improve the properties of RTV silicone rubber in the silicon industry.