Next Article in Journal
Microscopic, Spectroscopic, and Electrochemical Characterization of Novel Semicrystalline Poly(3-hexylthiophene)-Based Dendritic Star Copolymer
Next Article in Special Issue
Dispersion Homogeneity of Silicon Anode Slurries with Various Binders for Li-Ion Battery Anode Coating
Previous Article in Journal
Moisture Content Prediction in Polymer Composites Using Machine Learning Techniques
Previous Article in Special Issue
Rheological Properties and Melt Spinning Application of Controlled-Rheology Polypropylenes via Pilot-Scale Reactive Extrusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Peroxide-Based Crosslinking of Solid Silicone Rubber, Part I: Insights into the Influence of Dicumylperoxide Concentration on the Curing Kinetics and Thermodynamics Determined by a Rheological Approach

by
Maurício Azevedo
1,
Anna-Maria Monks
1,
Roman C. Kerschbaumer
1,*,
Sandra Schlögl
1 and
Clemens Holzer
2
1
Polymer Competence Center Leoben GmbH, Roseggerstrasse 12, 8700 Leoben, Austria
2
Department of Polymer Engineering and Science, Polymer Processing, Montanuniversitaet Leoben, Otto Gloeckel-Strasse 2, 8700 Leoben, Austria
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(20), 4404; https://doi.org/10.3390/polym14204404
Submission received: 12 September 2022 / Revised: 10 October 2022 / Accepted: 11 October 2022 / Published: 18 October 2022
(This article belongs to the Special Issue Rheology and Processing of Polymer Materials)

Abstract

:
Predicting the curing behaviour of industrially employed elastomeric compounds under typical processing conditions in a reliable and scientifically driven way is important for rubber processing simulation routines, such as injection moulding. Herein, a rubber process analyser was employed to study the crosslinking kinetics of solid silicone rubber based on the concentration of dicumylperoxide. A model was proposed to describe the optimal cure time variation with peroxide concentration and temperature, based on the analysis of processing parameters applying kinetic and thermodynamic judgments. Additionally, the conversion rate was described with the aid of a phenomenological model, and the effect of dicumylperoxide concentration on the final crosslink state was investigated using kinetic and thermodynamic explanations. Optimal curing time was affected both by temperature and dicumylperoxide concentration. However, the effects were less pronounced for high temperatures (>170 C) and high concentrations (>0.70 phr). A limit on the crosslink state was detected, meaning that the dicumylperoxide capacity to crosslink the silicone network is restricted by the curing mechanism. Curing restrictions were presumed to be primarily thermodynamic, based on the proton abstraction mechanism that drives the crosslinking reaction. In addition to providing more realistic crosslinking models for rubber injection moulding simulation routines, the results of this study may also explain the chemical behaviour of organic peroxides widely used for silicone crosslinking.

Graphical Abstract

1. Introduction

Poly(siloxane)s of high molecular weight, known as silicone rubbers, are materials based on macromolecules with successive silicon–oxygen bonds as a backbone, with the Si atom having two monovalent organic side groups. This alternating inorganic backbone of Si and O atoms provides strong and long chain bonds 106 kJ mol 1 and 1.64 Å, compared to 85 kJ mol 1 and 1.53 Å for C–C bond [1]), backbone flexibility (Si–O–Si bond angle of 143°, compared to 109.5° of the usual tetrahedron bond angles [2]), ease of side-group rotation, and low inter- and intramolecular forces [3]. As a result of the longer Si–O bond, compared to a C–C bond, there is more space for organic substituents, like methyl and phenyl, without an increase in steric hindrance or molecular congestion. Insertion of bulky substituents does not compromise the linearization of the molecules due to the long and flexible Si–O bond, allowing crystallization. The flexibility of the Si–O–Si bond is due to the oxygen atom delocalized lone electron pair, which is spread over the covalent bond region between Si and O atoms, specifically into the silicon vacant d orbital, allowing for a wider angle for the s p 3 hybridization. These molecular properties, linked to the inorganic nature of the backbone and the organic nature of the functional side groups, constitute an important bridge between inorganic and well-known organic polymers. Nowadays, silicone rubber has been used widely to produce medical devices, optical components, seals, and electronic parts, for example.
Silicone rubber, despite its very high molecular weight (>100,000 g mol 1 ), must be crosslinked to be employed in engineering applications [4]. The crosslinking reaction turns the soluble and highly viscous silicone rubber into an insoluble and viscoelastic network, following a specific crosslinking mechanism, kinetics, and thermodynamics. Industrially-applied elastomers, like the silicone rubber, are usually crosslinked by the use of specific chemicals (here denoted as crosslinkers), which along with temperature trigger the crosslinking reaction. Many are the chemical reactions that are able to join two adjacent macromolecules, being the most well-known ones the vulcanization (a mixed ionic-radical reaction at high temperature and applying sulphur, accelerators, and other additives) and the radical crosslinking (at high temperature usually employing organic peroxides) [5,6,7]. These two crosslinking methods are currently widely applied for producing rubber parts, employing different polymer processing techniques, such as compression and injection moulding. With the advent of numerical simulation, polymer processing has become a precise and highly technological industrial step. Before actual rubber injection moulding, for example, the process can be computer-aided simulated, in order to establish the most suitable processing conditions for a specific application [8,9]. During simulation, several models are applied to describe the injection moulding process, from the filling phase to the crosslinking phase of the product [10,11]. These models must account not only to describe the processing parameters, such as the time necessary to crosslink the silicone rubber product [12] or the pressure inside the cavity, but also to consider the kinetics and the thermodynamics of the related phenomena.
There are numerous scientific contributions in the field of silicone rubber crosslinking kinetics [13,14,15,16,17], like the one by Hernández-Ortiz and Osswald [18]; however, most of them do not approach processing-related parameters. Additionally, it is rare to find studies that cover the kinetics and the thermodynamics of silicone rubber crosslinking by analysing rheological experimental data. Thus, this study aims to cover these gaps in two different, but converging, fronts: first, to propose a processing-related point-of-view for the crosslinking kinetics by analysing important processing parameters; and second, to describe how the thermodynamics of dicumylperoxide-based silicone rubber crosslinking affects the overall crosslink state of the rubber product. The properties of the final silicone rubber network, which corroborate to the findings of the present study, are currently under investigation and will be thoroughly discussed in a separate publication.

2. Materials and Methods

This section briefly describes the elastomer under investigation and the methodology for the rheological measurements. For further details about the rubber process analyser device and its versatility, the following standards (ASTM, American Society for Testing and Materials, West Conshohocken, PA, USA) cover the main topics related to this kind of measurement: ASTM D5289 [19], ASTM D6204 (A, B, and C) [20], ASTM D6601 [21], and ASTM D8059 [22]. It is highly advised that the Appendix A to this publication is consulted for practical aspects of the described experiments.

2.1. Materials

High consistency poly(dimethylsiloxane), or solid PDMS (Xiameter™ RBB-2100-50, M w = 660 kg mol 1 , M w / M n = 1.8), containing approximately 26 wt% of an inorganic filler and no significant concentration of vinyl side groups (as confirmed by 1 1H-NMR studies, see Appendix A) was supplied by Dow Inc. (Midland, MI, USA). The PDMS characterization in terms of molecular weight, filler, and vinyl contents is described in the Appendix A. Dicumylperoxide (DCP) 99.9% (Peroxan DC) was supplied by Pergan GmbH (Bocholt, Germany), with active oxygen content of 5.91 wt%. Both components were thoroughly mixed in a 2-roll-mill varying the peroxide concentration: 0, 0.21, 0.35, 0.49, 0.70, 1.00, and 1.50 phr of DCP. These concentrations represent fractions of the suggested peroxide concentration (0.70 phr) by the silicone supplier, i.e., 0, 30, 50, 75, 100, 140, and 215%, respectively. All silicone compounds were stored at low temperatures (<5 C) prior to testing, aiming to avoid premature pre-curing.

2.2. Rotational Rheometry

Rubber process analyser (D-RPA 3000 Montech Werkstoffprüfmaschinen GmbH, Buchen, Germany) equipment, or moving die rheometer (MDR), was employed as a rotational rheometer to characterize the curing behaviour of the PDMS/DCP compounds at different temperatures: 140, 150, 160, 170, and 180 C for 1 h. The rotational deformation was set to 0.5° at a constant frequency of 1.667 Hz. Processing-related parameters were calculated from the torque vs. time curves: minimum torque ( M L ), maximum torque ( M H , as the torque value at t = 1 h, since no reversion or marching modulus was observed), scorch time ( t s 1 ), and optimum cure time ( t 90 ). The scorch time was defined as the time correspondent to increase M L by 1 dN m (torque unit), and the optimum cure time as the time correspondent to 90% of M H . All experiments were performed in triplicate for each sample (7 dicumylperoxide concentrations and 5 temperatures).
To describe the effect of temperature and dicumylperoxide concentration on the optimum cure time, fitting of the data was performed using the Levenberg–Marquardt algorithm [23] or damped least-squares method. A general equation to define t 90 was proposed, into which the experimental data was fitted:
t 90 = A 0 e x p E a R T ( D C P ) β
In Equation (1), the optimum cure time is described as a function of two factors. To the temperature was assigned an Arrhenius relation, since time and temperature were under study. The influence of dicumylperoxide concentration, on the other hand, was described by a power law with order β .
For the crosslinking kinetics study, only the data related to the curing temperature of 160 C was used, since this is the common and the manufacturer-advised curing temperature for the silicone rubber under investigation. However, the usual m o d u s   o p e r a n d i is to fit the data related to several thermal programs, those being either different temperatures or different heating rates. Three ASTM standards currently describe how to estimate kinetic parameters utilizing thermal analysis: ASTM E2041-13 (Standard test method for estimating kinetic parameters by differential scanning calorimeter using the Borchardt and Daniels method) [24], ASTM E698-18 (Standard Test Method for Kinetic Parameters for Thermally Unstable Materials Using Differential Scanning Calorimetry and the Flynn/Wall/Ozawa Method) [25], and ASTM E2781-11 (Standard Practice for Evaluation of Methods for Determination of Kinetic Parameters by Thermal Analysis) [26]. Along with these standards, the publications by Vyazovkin [27,28,29,30,31] give important information related to determining the kinetics of several phenomena, such as crystallization, glass transition, and crosslinking.

3. Theoretical Background

The description of curing kinetics in terms of experimental values is detailed below by applying a phenomenological model. In contrast to the mechanistic models, which are based on solving numerous differential equations that describe mass action kinetics, these are formulated based on the curing degree, or conversion [32]. For engineering and processing simulation purposes, phenomenological models produce parameters that are easy to implement, making them more convenient than mechanistic models. As a result, the phenomenological model considers the curing reaction as a single chemical step, which is obviously a simplification of the whole crosslinking phenomenon that includes several sequential reactions, which will be discussed in the following section.

3.1. Crosslinking Kinetic Modelling

Concerning rheological experiments, the crosslink conversion rate d α d t can be calculated considering the difference between the torque before and after the curing reaction and the actual torque at a certain time t. Specifically for RPA, the torque before curing is defined as the minimum torque ( M L ) and the torque after curing as the maximum torque ( M H ), being their difference expressed as Δ M = M H M L . Using Equation (2), the conversion rate at a given point in time can be calculated based on the following mathematical relation:
d α d t = 1 Δ M d ( M M L ) d t
Considering that crosslink formation is a temperature-induced process, the same conversion rate can be represented as a function of temperature k ( T ) , and as a function of the appropriate kinetic model f ( α ) , as reviewed by Vyazovkin [28]:
d α d t = k ( T ) f ( α )
It is well-known and reported that k ( T ) is an Arrhenius-like rate constant, i.e., it is a correlation between temperature T (K), a pre-exponential factor A (time 1 ), an activation energy E (J mol 1 ), and the molar gas constant R (8.3145 J mol 1 K 1 ) as follows:
k ( T ) = A e x p E R T
Substituting Equation (4) into Equation (3), and taking a logarithm followed by differentiation against 1 T gives rise to the equation that defines the isoconversional principle: considering two samples, if the curing rates d α d t are the same and the conversions α are equal, then the curing rate is only a function of temperature. The isoconversional principal equation is then written as:
ln ( d α d t ) T 1 α = ln ( A ) T 1 α + ln ( f ( α ) ) T 1 α E R α
Analysing Equation (5) enables one to make the following conclusions: the term l n ( A ) T 1 α is zero, since differentiating a constant results zero; and l n ( f ( α ) ) T 1 α also equals zero, since the kinetic model f ( α ) is only a function of the conversion. Finally, one can conclude that the activation energy for a certain conversion value α (now referred as E α ) does not depend on the chosen kinetic model, i.e., it is model-free:
E α = R ln ( d α d t ) T 1 α
The assumption that the kinetic model does not change with temperature is reasonable, since the temperature range related to a given conversion is so narrow that no change would be detected [28]. Obtaining the values for activation energy in terms of conversion without taking the kinetic model into consideration is very convenient and already gives light to important details about the crosslinking mechanism. Thus, this principle was employed to study the crosslinking kinetics during the present study.
There are many different approaches that apply the isoconversional principle to determine the activation energy of thermally induced processes. Vyazovkin et al. [27] and Zhang [33] give a broad overview about these methods, which include the integral isoconversional one used in this research: the Friedman-like [34], for isothermal experiments, represented by Equation (7).
ln ( t α , i ) = ln [ f ( α ) A α ] E α R T α , i
The subscripts α , i indicate a given conversion α and a given thermal program i (temperature for isothermal experiments or heating rate for dynamic experiments) [29]. It is possible to notice after careful analysis that Equation (7) can be written in the linear form y = a x + b , assuming 1 T as x. Thus, the activation energy E α can be determined for each conversion value (here represented by t α , i ) by plotting the left-hand side of Equation (7) over 1 T . This strategy was used in the present research to determine an initial guess for the activation energy, to be implemented in the subsequent fitting procedure.
Regarding the kinetic model function f ( α ) , it is widely accepted [13,14,16,28] that the crosslinking reaction of silicone rubbers follows the autocatalytic model, first proposed by Šesták-Berggren [35] in the form of Equation (8), and further modified by Kamal [36] to include the temperature dependence, as seen in Equation (9):
f ( α ) = α m ( 1 α ) n
d α d t = ( k 1 + k 2 α m ) ( 1 α ) n
where m and n are the reaction orders and k 1 and k 2 are the Arrhenius rate constants, as previously defined in Equation (4). For the present study, a third version of the autocatalytic model was assumed, where k 1 and k 2 are taken as constants with activation energies E α , 1 and E α , 2 , in a way that the final kinetic model equation can be written, after proper substitution using Equations (4) and (9), as:
d α d t = A 1 e x p E α , 1 R T + A 2 e x p E α , 2 R T α m ( 1 α ) n
The subscript 1 denotes the n t h order contribution to the crosslinking model, while the subscript 2 denotes the auto-catalytic contribution. Similarly, n is the reaction order for the n t h model, and m is the reaction order for the auto-catalytic part.
The RPA data (calculated conversion rate d α d t values based on experimentally-determined torque data, according to Equation (2)) for T = 160 C was fitted to Equation (10) for conversion α values between 0.1 and 0.9, as normally recommended [13,27]. Initially, it was assumed [13] that E α , 1 = E α , 2 = E α as calculated applying the isoconversional approach. This E α value was then implemented as the activation energy parameters initial guess at the preliminary fitting routine and the parameters A 1 , A 2 , m, n, and E α were determined employing the Levenberg–Marquardt algorithm [23] or damped least-squares method. After the first preliminary fitting, the previously calculated parameters were employed as initial guesses for the final calculations of the kinetic parameters A 1 , A 2 , m, n, E α , 1 , and E α , 2 . All fitting procedures were performed utilizing Python coding language (Python 3.8, PyCharm Community Edition 2021.2.2, Jet Brains, Prague, Czechia).

3.2. Peroxide-Based Crosslinking Mechanism of PDMS

Crosslinks formation between two adjacent PDMS macromolecules, when aided by an organic peroxide, is based on the generation of radicals by the peroxide decomposition reaction. The interaction between peroxide radicals and PDMS macromolecules can occur via hydrogen abstraction from a methyl and/or vinyl groups or via addition to the double bond on the vinyl group, if this is present on the chain. If the energies of the resulting radicals are compared and considering that hydrogen abstraction is more favourable if the radical energy is lowered, hydrogen abstraction from the vinyl group is unfavourable, due to the similar energy of alkoxy and vinylic radicals [37].
When dicumylperoxide is concerned, the main chemical reactions involved during PDMS crosslinking via hydrogen abstraction from methyl side groups are:
i
C 18 H 22 O 2 Δ 2 C 9 H 11 O ·
ii
C 9 H 11 O · C 8 H 8 O + C H 3 ·
iii
C 9 H 11 O · + [ S i O S i ( C H 3 ) ] n C 9 H 11 O H + [ S i O S i ( C H 2 · ) ] n
iv
C H 3 · + [ S i O S i ( C H 3 ) ] n C H 4 + [ S i O S i ( C H 2 · ) ] n
v
2 [ S i O S i ( C H 2 · ) ] n [ S i O S i ( C H 2 ) ( C H 2 ) S i O S i ] n
vi
C 9 H 11 O H C 9 H 10 + H 2 O
Radical generation is accomplished at reactions i and ii, via thermal decomposition of dicumylperoxide and beta scission of the cumyloxy radical, generated in the reaction i, respectively. There are two radicals responsible for silicone crosslinking: the cumyloxy and the methyl radicals, which are able to abstract a proton from the methyl side group of PDMS (reactions iii and iv), creating two macromolecular radicals. It was shown by Baquey et al. [38] that the main responsible for PDMS proton abstraction is the radical directly formed after peroxide decomposition, which in the present case is the cumyloxy. However, the methyl radical is also considered in this study as able to abstract protons from the silicone polymer chain. Finally, the macromolecular radicals, when adjacent to each other, can react with one another, creating a covalent bond between the carbon atoms from the methyl side groups (reaction v).
For the case where double bond addition occurs, the cumyloxy (denoted as R and the methyl (denoted as R’) radicals can add to the less-substituted secondary carbon (due to its higher stability as intermediate radical) of the vinyl group, as shown in reactions vii and viii, leading to the formation of macromolecular radicals. Similar to what happens in the hydrogen abstraction case, adjacent macromolecular radicals are able to react with each other by recombination, creating a crosslink point (reaction ix), or the macromolecular radical can continue adding to another adjacent vinyl group. It is worth mentioning that additional reactions result in products with fragments from the peroxide initiator, i.e., the fragment R or R’ is attached to the polymer chain.
vii
R · + [ S i O S i ( H 2 C = C H ) ] n C 9 H 12 O + [ S i O S i ( H 2 C R C H · ) ] n
viii
R · + [ S i O S i ( H 2 C = C H ) ] n C 9 H 12 O + [ S i O S i ( H 2 C R C H · ) ] n
viiii
2 [ S i O S i ( H 2 C R C H · ) ] n [ S i O S i ( H 2 C R C H ) ( C H R C H 2 ) S i O S i ] n
When vinyl moieties are present in the PDMS backbone, hydrogen abstraction and vinyl addition happen simultaneously, even though, considering alkoxy (from dicumylperoxide) and methyl radicals (both 105 kcal mol 1 ), the formation of secondary R 2 C H · radicals is more favourable than primary R C H 2 · radicals (97 against 100 kcal mol 1 , respectively) [37], due to the alkyl substitution stability effect. However, vinyl addition extension is limited to the vinyl concentration on the PDMS chain, while hydrogen abstraction has virtually no concentration limitation in terms of methyl presence; it is only limited due to steric hindrance and/or thermodynamical reasons.
Concurrently to dicumylperoxide-mediated crosslink formation in PDMS, by-products are formed due to hydrogen abstraction or decomposition of radicals. In reaction iii, cumyl alcohol C 9 H 11 O H is formed after hydrogen abstraction by the cumyloxy radical, and further dehydrated and converted into α -methylstyrene (reaction vi). In reaction ii, acetophenone C 8 H 8 O is the by-product of cumyloxy beta scission. All these by-products do not contribute to the network formation, but are treated and understood as low molecular weight and volatile compounds that may migrate, causing mainly the very characteristic odour of peroxide-crosslinked silicone products. Cumene may also be formed after reduction of α -methylstyrene [39].
The aforementioned crosslinking kinetic model accounts for the overall crosslink formation. This model is not able to differentiate between abstraction or addition reactions due to its phenomenological approach [40], but clarifies the rate of crosslink formation. A complete description of the complex set of reactions that underline peroxide crosslinking of PDMS can be performed considering mechanistic approaches, like population balance equations, which is not the main focus of the present study, where a phenomenological kinetic model was proven to be sufficient.

4. Results

This section is divided into three major subsections. The first, comprehending Section 4.1, discusses the curing kinetics and its association with processing conditions, proposing a mathematical relation for the optimum cure time based on the experimental data. The second, including Section 4.2, reports the crosslinking kinetics and presents the kinetic parameters, along with the reasoning behind their trend as the dicumylperoxide concentration increases. Finally, the third, Section 4.3, presents thermodynamic considerations regarding the dicumylperoxide-based crosslinking of silicone rubber.

4.1. Curing Characteristics

Figure 1 displays the curing curves for the PDMS/DCP systems for different dicumylperoxide concentrations and at different curing temperatures. Curing curves are usually analysed in rubber industry and research to identify the crosslinking profile in terms of minimum detected torque, induction time, velocity, curing time, and maximum detected torque. A torque increase is present for all samples with dicumylperoxide, regardless of the concentration, due to the formation of a denser polymer network, promoted by crosslinking. Connecting adjacent PDMS macromolecules forms a stiffer and more elastic network, which imposes a higher torque to the device, in order to keep the set deformation.
For all dicumylperoxide concentrations, the temperature increase shortens the induction or scorch time and lowers the minimum detected torque, as also depicted in Table 1. Since dicumylperoxide decomposition is initiated and promoted by heat (reaction i in Section 3.2), a higher curing temperature delivers more heat to the rubber compound, leading to a shorter crosslinking induction time, represented by the t s 1 values. For a constant curing temperature, the higher the dicumylperoxide concentration, the lower the induction time. The effect of dicumylperoxide concentration on the induction time is explained by the higher concentration of radicals that are formed after decomposition and by the autocatalytic nature of this process, as also reported by Kruželák et al. [41,42]. Industrially, this value is related to process safety and connected to the mouldability of the rubber compound, which is limited after t s 1 is reached. As a matter of example, short scorch times are preferred for extrusion of hollow rubber profiles, to avoid extruded collapse after passing through the matrix. On the other hand, for rubber injection moulding, scorch time must be long enough to guarantee the complete filling of the mould before curing.
The minimum detected torque M L was also slightly lowered by temperature increase, which is mainly an effect of temperature over the rubber compound’s viscosity. The effect of temperature over rubber viscosity is well known, and the processability is related to this value when RPA is employed in the rubber industry. Along with the values for Mooney viscosity, the minimum detected torque can clarify aspects about rubber processing, including mixing and moulding. No effects of dicumylperoxide concentration were detected at the minimum torque due to the low concentration, the low molecular weight, and the dispersion of the peroxide into the rubber matrix.
The curing rate is enhanced by temperature, as the cure curves, for all dicumylperoxide concentrations, become steeper for higher curing temperatures. Dicumylperoxide decomposition kinetics is favoured by temperature [43,44], leading to a faster crosslinking reaction for the PDMS/DCP systems as well. Crosslinking velocity is also promoted by higher concentrations of dicumylperoxide, due to the higher concentration of radicals and the autocatalytic nature of this process. Table 2 describes the optimal cure time t 90 , which decreases as temperature and dicumylperoxide concentration increase. This decrease is an effect of temperature over the dicumylperoxide decomposition and rubber crosslinking kinetics, and of the quantity of radicals that are formed. Industrially, the selection of the most appropriate curing time is a compromise between processing time, processing temperature, and part quality in terms of mouldability.
As an indicative of the final crosslink density, the maximum detected torque represents the maximum shear force resistance of the rubber compound at a given time and temperature, for the set deformation [4]. The maximum torque values, as shown in Table 2, tend to increase with temperature and with the dicumylperoxide concentration, suggesting that the effective number of crosslink points is also enhanced by temperature and dicumylperoxide amount. Even though the maximum torque values increase with dicumylperoxide, a threshold of approximately 11–12 dN m is reached for the highest temperatures/dicumylperoxide concentrations.
Regarding the samples without dicumylperoxide, an increase in the detected torque is observed for all samples, with a steeper increase for higher temperatures. Even though a decrease in the torque within the first minutes due to viscosity reduction is detected, torque increases with time, which may indicate some crosslinking phenomenon not associated to the presence of dicumylperoxide, but also triggered by temperature. This phenomenon may be associated to peroxide or metal contamination, that ultimately would lead to crosslinking when the temperature is increased.
When cycle time for moulding is evaluated for industrial applications, the optimum cure time t 90 is a factor usually taken into account. An estimation for the optimum cure time within the temperature and dicumylperoxide ranges that were investigated is shown in Figure 2, taking Equation (1) as the mathematical model for the parameter. The proposed equation can fit almost 99% of the experimental data, as represented by the R 2 factor.
Temperature has a major influence over the optimum cure time for lower peroxide concentrations, while for concentrations above 0.70 phr the effect is not prominent. Similarly, the peroxide concentration has a stronger influence on the optimum cure time for lower temperatures, becoming weaker as the temperature increases. The less distinguished effect on the optimum cure time for higher temperatures and higher DCP concentrations is related to two distinguished reasons. For the temperature, the higher range (>150 C) is well above the dicumylperoxide decomposition temperature, which is around 120 C [44]. This means that the decomposition and the preceding crosslinking reaction will take place almost instantly and rapidly, and the different peroxide concentrations will not impact the optimum cure time in a notable way. Regarding the DCP concentration, the higher concentration range (>0.70 phr) shows a weaker effect of temperature on the optimum cure time, possibly due to the autocatalytic nature of peroxide-based crosslinking reactions. In these cases, the number of generated radicals is enough to compensate the temperature effect and overlap it. Hou et al. [45] recently reported the same effect of temperature on the optimum cure time for dicumylperoxide-crosslinked silicone rubber.
Here in this paper, an equation is proposed to model the optimum cure time for peroxide-cured high consistency silicone elastomers. The parameters for this equation are also described in Figure 2. The activation energy is associated to the amount of energy necessary to reach the optimum cure time, which is also associated to the amount of energy necessary to decompose the dicumylperoxide molecules as stated in the chemical equation i in Section 3.2. Thus, the value of 115 kJ mol 1 is in accordance with the values reported in the literature [43,44] for DCP decomposition, and with the activation energy values that will be presented further in the crosslinking kinetics discussion. The exponent for the dicumylperoxide concentration has a negative value due to the decrease of the optimum cure time with the peroxide concentration and is less than 1 possibly due to its limitation in crosslinking the polymer, which will also be further explained. The pre-exponential factor, similarly to the activation energy, presents a value that will also similarly reported when the crosslinking kinetics is discussed further in this paper. It is important to state here that these parameters and this specific equation fit to the polymer and to the crosslinking system that were investigated in this paper, at the same time. For other polymeric systems, with different crosslinking mechanisms and kinetics, another set of equations and parameters must be proposed and evaluated.

4.2. Crosslinking Kinetics

Prediction of how a polymer crosslinks when exposed to heat during typical processing conditions, like compression or injection moulding, is a key step to control the overall process and the part quality. The crosslinking kinetics, i.e., the evolution of conversion α over time, was evaluated in terms of the torque increase detected at the rubber process analyser at 160 C, which was the temperature set for compression moulding.
Based on the curves shown in Figure 1, the conversion rate d α d t was calculated employing Equation (2) and is plotted versus time in Figure 3a and versus conversion in Figure 3b. The effect of increasing the dicumylperoxide concentration is seen as an acceleration of the reaction, i.e., higher cure velocity, when analysing the conversion rate evolution on time. Additionally, the maximum cure rate is shifted to lower times. These findings are related to the concentration of peroxide radicals that are formed at 160 C and its catalytic effect on the formation of further radical species. Regarding the evolution of cure velocity over conversion, or reaction path, the maximum cure rate remains constant at a fixed conversion value between 0.2 and 0.4 for the dicumylperoxide concentrations up to 0.70 phr. This means that the highest rate of crosslink formation occurs at a constant conversion range, regardless of the dicumylperoxide concentration below 0.7 phr. It is possible to speculate that the formation of the characteristic three-dimensional crosslinked network occurs mostly within the 0.2 and 0.4 conversion range for DCP concentration equal or lower than 0.70 phr, as a result of the decomposition of the majority of dicumylperoxide molecules. Within this conversion range, the concentration of radicals likely decreases, either by transfer to the polymer chain or by termination, resulting in the end of the crosslinking reaction. For the samples with higher dicumyperoxide concentrations (1.00 phr and 1.50 phr), the abruptly fast crosslinking reaction, as seen by the conversion rate values, shifts the maximum reaction speed to lower conversion values, probably due to the high amount of radicals that are formed.
As extensively described by Vyazovkin et al. [27,28,31], the determination of the kinetic parameters starts with calculating the activation energy without defining the most appropriate kinetic model, a strategy that is called model-free. The activation energy calculation following a Friedman-like approach as described by Equation (7) is exemplified for the 0.21 phr sample in Figure 4a. The slopes of the linear-fitted lines (correlation factor for the linear fittings were higher than 0.960) were used to determine the activation energy variation with conversion, as shown in the diagram of Figure 4b. A decrease of the activation energy is realized as the dicumylperoxide content increases, representing thermodynamically easier crosslinking reactions. This behaviour was already observed before, but in a kinetic way: the higher peroxide concentration accelerated the crosslinking reaction. Both thermodynamic and kinetic phenomena are behind the auto-catalytic nature of dicumylperoxide decomposition, which precedes the polymer crosslinking. Bianchi et al. [46] observed the same behaviour when studying the crosslinking of poly(ethylene-vinyl acetate) with dicumylperoxide, similarly reporting activation energy values between 85 and 105 kJ mol 1 .
Duh et al. [43] reported an activation energy for dicumylperoxide thermal decomposition of 110–150 kJ mol 1 , while Lv et al. [44] reported values between 150 and 200 kJ mol 1 , which agree with the activation energy determined in the present study, i.e., varying between 85 and 140 kJ mol 1 .
The activation energy is fairly constant in the conversion range 0.1–0.9, which is reasonable for a system that does not undergo drastic molecular diffusion changes during crosslinking, which would impose mass transfer limitations to the process [28]. Systems that experience a greater change in viscosity, like epoxy thermoset resins, are more sensitive to the diffusion hindrance caused by the crosslinking reaction (due to increase in viscosity), representing an increase of the activation energy [47]. Another consideration that deserves commenting is the heat transfer change due to the formation of a crosslinked layer at the top and bottom of the specimen inside the rubber process analyser during the first instances of curing. As studied by Cheheb et al. [48], the thermal conductivity of a crosslinked rubber sample is 10% higher than a non-crosslinked one, but still much lower than the metallic RPA cavity. This means that during curing, the heat transfer changes. However, this limitation was not further studied in the present work due to the low thickness of the RPA specimen (which imposes a fast heating of the whole sample) and the lack of scope of this investigation. Since the activation energy exhibited a relatively constant trend with conversion, the experimental data were fitted to the Kamal model using the mean values to calculate the kinetic parameters. Activation energy deviations that are observed at higher conversion values ( α > 0.9 ), which may indicate a decrease of the activation energy, arise from the small variation of torque at the final stages of crosslinking, which leads to high standard deviation of data.
Conversion rate data were fitted applying Equation (10), resulting in the kinetic parameters that are shown in Table 3. The calculated kinetic energies are in accordance with the values previously determined via the Friedman approach, and also showing a decrease with the dicumylperoxide concentration. Reaction orders m and n do not follow a specific trend and are considered here as being rather constant. Considerations regarding the reaction orders are associated with the crosslinking reaction mechanism, which is understood as the same regardless of the peroxide’s quantity. The pre-exponential values A 1 and A 2 have a tendency of decrease as the dicumylperoxide concentration increases.
However, analysing the pre-exponential factor alone may lead to misunderstandings [31], so the Arrhenius rate constants k 1 and k 2 are also reported in Table 3. The constants k 1 and k 2 are measurements of the system’s reactivity and are enhanced as the dicumylperoxide concentration increases. This is an additional evidence that corroborates the catalytic effect of the dicumylperoxide concentration.
Analysing the pre-exponential factors solely may lead to misinterpretations since it is often related to the activation entropy, meaning that the higher the pre-exponential factor, the more accelerated the chemical reaction is, as stated in Equation (11) [49]. This mathematical relation is a different way of expressing the Arrhenius rate constant previously denoted as Equation (4).
k ( T ) = k B T h e x p Δ S R e x p Δ H R T
The factor that contains the activation entropy Δ S , the Boltzmann constant k B , and the Planck constant h was previously assigned as the pre-exponential factor A; while the inverse dependence on temperature is represented by the activation enthalpy Δ H factor. In the case of the present PDMS/DCP systems, it is not correct to assume that due to the decrease of the pre-exponential factor, the crosslinking reaction is decelerated as the dicumylperoxide concentration increases. Actually, the effect on the system’s reactivity, or the overall free energy of the activated reaction intermediate Δ G , is dominated by the enthalpy factor Δ H via the decrease of the activation energy, which is a pure effect of the peroxide catalyst’s nature, as explained by the usual thermodynamic relationship, derived from the second law of thermodynamics:
Δ G = Δ H T Δ S
The free energy Δ G is thus lowered strongly by the decrease of the activation enthalpy Δ H in comparison with the not so evident entropic effect. It is beyond the scope of this paper to explain why the activation entropy decreases as the dicumylperoxide increases, thus justifying further investigations. However, one rationalization may be that the active sites for radical reaction, which ultimately leads to a crosslink point, become limited as the dicumylperoxide concentration increases, possibly justifying the crosslinking threshold that is observed on the final crosslinked elastomer specimen. A second reasoning is that as the dicumylperoxide concentration increases, the radical concentration is also enhanced, increasing the probability of radical recombination before the reaction with the polymer chain, likewise limiting the crosslinking process. The last reasoning was already studied by Parks and Lorenz [50] when investigating the efficiency of dicumylperoxide to react with dimethyloctadiene.
After time integration, the conversion values calculated according to the model parameters were compared to the experimental ones, which are plotted in Figure 5. Within the conversion range 0.1–0.9, the calculated conversion values accordingly predict the experimental values with correlation factors higher than 0.99. For the specific temperature (160 C) and the dicumylperoxide concentrations that were evaluated, it is possible to predict how the crosslinking conversion behaves with time, meaning that this model can be possibly integrated into processing simulation routines. Along with the proposed mathematical relationship depicted at Equation 1 for the optimum cure time, these suggestions may strongly enhance the quality of processing simulations that deal with reactive systems, such as elastomers.
Incorporating the variable dicumylperoxide concentration into the crosslinking kinetics model was not included in the scope of this paper, since it would demand an effort beyond the selected scientific targets. However, the ideal crosslinking kinetics model should embody not only the temperature of curing, but also a factor that correlates with the elastomer recipe in terms of the crosslinking systems. This becomes a suggestion for further studies.

4.3. Thermodynamic Reasoning behind the Dicumylperoxide-Based Crosslinking Reaction

Since no significant vinyl concentration was detected at the PDMS matrix via high-field solution 1 H-NMR (as described in the Appendix A), considerations regarding the crosslinking mechanism relies on analysing the radical energies associated to the involved species. As reported by Dluzneski [37], the radical ability to abstract protons obeys the following ranking, where the bond dissociation energy of the first proton abstraction is expressed inside parenthesis and in kcal mol 1 ( P h stands for C 6 H 5 ):
P h O · ( 87 ) < P h C · ( 88 ) C = C C · ( 88 ) < R 3 C · ( 91 ) < R 2 H C · ( 97 ) < R H 2 C · ( 100 ) < C = C · ( 104 ) < H 3 C · ( 105 ) R O · ( 105 )
R O · and H 3 C · radicals are generated with the higher bond dissociation energies, meaning they are the ones more prone to abstract a proton from the organic side group at the polymer chain. Considering reactions iii and iv as previously described in Section 3.2, these are characterized by the radical stabilization energy (RSE) [51] of the newly formed radical, which is the macromolecular radical [ S i O S i ( C H 2 · ) ] n , relative to the methyl and cumyloxy radicals. A macromolecular radical of this type behaves like a secondary carbon radical R H 2 C · , where R is the polymer backbone, meaning it has energy of about 100 kcal mol 1 . For the proton abstraction reaction to happen, the RSE has to be negative, i.e., a decrease in the radical energy has to occur. For the related radical species, a decrease of only 5 kcal mol 1 takes place, which is probably enough for the reaction to thermodynamically proceed. This net decrease of around 5 kcal mol 1 is confirmed by an analogy by Bordwell et al. [52], where they report that the RSE for M e 3 S i C H 2 · when compared to the methyl radical H 3 C · is −6 kcal mol 1 .
Furthermore, this thermodynamic aspect related to the decrease in the radical energy can be understood as an enthalpic impediment, different from an entropic boundary that rises if steric hindrance or diffusion limitations takes place. By analysing once again Equation (12), one can realize that the factor affecting the reaction free energy is the enthalpic contribution, by decreasing the activation energy. This thermodynamic “boost” is associated to the crosslinking reaction and accomplished by the catalytic effect, i.e., decreasing the activation energy of the crosslinking process and allowing the reaction to take place in a thermodynamically favoured scenario, given the thermodynamic circumstances of proton abstraction. Besides, increasing the temperature also enhances the entropic factor T Δ S , which also leads to a decrease of the free energy.
Thus, one can understand that the crosslinking reaction is limited by the proton abstraction capacity of the radical species R O · and H 3 C · , which, even though it exists, is considered as having a low thermodynamic impulse (RSE∼5 kcal mol 1 ). This is perhaps the factor that impedes the dicumylperoxide concentration to improve the crosslink density, as detected by the limited maximum detected torque. One evidence that corroborates this reasoning is the fact that Verheyn et al. [53] reported a substantial increase of silicone rebound resilience as the concentration of either di(2,4-dichloro-benzoyl)peroxide or di(4-methylbenzoyl)peroxide increased. Both peroxides decompose, producing aryl carboxylic radicals, which are even more stable than cumyloxy and methyl radicals, leading to more negative values for RSE.

5. Conclusions

This work reported a processing-oriented crosslinking kinetics study of a solid silicone rubber, or high consistency silicone, which is a widely employed material for manufacturing of technical rubber products. The investigation was carried out applying a rheological approach, by employing an industrially well-known device, the rubber process analyser (RPA). Different dicumylperoxide concentrations were studied, in order to determine their effect on processing parameters and on kinetic parameters.
Rheological experimental data showed that dicumylperoxide concentration enhancement accelerates the crosslinking reaction by shortening the induction time and increasing the reaction velocity. Besides, it results in higher but limited detected maximum torque, with distinguished extents for lower and higher temperatures. A correlation between the optimum cure time, the curing temperature, and the dicumylperoxide concentration was given, aiming to predict the cycle time during compression or injection moulding, for example.
Determination of the kinetics parameters showed that even though there was a decrease of the pre-exponential factor, the overall reactivity (represented by the Arrhenius rate constant) increased as the dicumylperoxide concentration rose. The overall reactivity was boosted due to the activation energy decrease caused by the increase on the dicumylperoxide concentration. This effect was associated to the catalytic effect of the radicals responsible for crosslinking.
For the first time, a thermodynamic and kinetic-based explanation was given for the limited effect of dicumylperoxide crosslinking efficiency on solid poly(dimethylsiloxane) rubber, considering the hydrogen abstraction mechanism. The findings of this work are useful for simulating processing cycles, enhancing the simulation similarity with real moulding. The effect of this specific crosslinking mechanism on the final silicone network and resulting properties is currently under investigation and will be disclosed in a separate study.

Author Contributions

Conceptualization, M.A. and R.C.K.; methodology, M.A. and R.C.K.; formal analysis, M.A. and A.-M.M.; investigation, M.A. and A.-M.M.; data curation, M.A. and A.-M.M.; writing—original draft preparation, M.A.; writing—review and editing, M.A., A.-M.M., R.C.K., S.S. and C.H.; supervision, R.C.K. and C.H.; project administration, R.C.K. and C.H.; funding acquisition, R.C.K. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

The research work presented at this publication was part of the COMET project “Innovative material characterization methods for reactive systems” number 1072920 (ReaSys project) and accomplished within the framework of the COMET-program of the Austrian Federal Ministry for Climate Action, Environment, Energy, Mobility, Innovation and Technology and the Federal Ministry for Digital and Economic Affairs.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data presented in this publication are only available upon request to the corresponding author, assuming a formal approval of the involved companies and co-authors.

Acknowledgments

The research work of this paper was performed at the Polymer Competence Center Leoben GmbH (PCCL, Austria) within the framework of the COMET-program of the Federal Ministry for Climate Action, Environment, Energy, Mobility, Innovation and Technology and the Federal Ministry for Digital and Economic Affairs, with contributions from the Montanuniversitaet Leoben (Department of Polymer Engineering and Science, Polymer Processing and Materials Science and Testing of Polymers). The authors kindly acknowledge Bernard Lechner for experimental support, and Johannes Macher for the help with fitting the experimental data using Python. PCCL is funded by the Austrian Government and the State Governments of Styria, Lower Austria and Upper Austria.

Conflicts of Interest

M.A., A.-M.M., R.C.K., and S.S. are employed by Polymer Competence Center Leoben GmbH. The company had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript and in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
A i Preexponential factor
DCPDicumylperoxide
DSCDifferential scanning calorimetry
hPlanck constant
Δ H Enthalpy change
E a Activation energy
1 H NMRProton-1 nuclear magnetic resonance
k i Arrhenius rate constant
k B Boltzmann constant
mAutocatalytic reaction order
MHMaximum detected torque
MLMinimum detected torque
M i Molecular weight of a chain
M n Number average molecular weight
M w Weight average molecular weight
M z Z-average molecular weight
M z + 1 Z+1-average molecular weight
n n t h -order reaction order
N i Number of chains of a specific molecular weight
PDMSPoly(dimethylsixolane)
phrParts per hundred rubber
ppmParts per million
RMolar gas constant
RPARubber process analyser
RSERadical stabilization energy
Δ S Entropy change
TTemperature
tTime
t s 1 Scorch time
t 90 Optimum cure time
TGAThermogravimetric analysis
α Conversion
δ Chemical shift

Appendix A

Appendix A.1. Determination of the Molecular Weight

The molecular weight distribution was determined via gel permeation chromatography (GPC) by employing a GPC analyser (Shimadzu Corp., Kyoto, Japan) coupled with a refractive index detector RID-20 and a MZ-Gel SDplus Linear 5 μ m column. Sample preparation was accomplished by producing three 10 mg mL 1 silicone rubber solution samples in triplicate (samples A, B, and C) in tetrahydrofuran by solubilization. Only the sample 0 phr dicumylperoxide was analysed. These solutions were subsequently filtered through a 0.45 μ m PTFE membrane filter, aiming to get rid of the filler incorporated to the silicone rubber gum. Calibration was made using poly(styrene) references and a calibration curve was built.
Several averaged molecular weights are reported in Table A1, along with the polydispersity indexes M w M n and M z M w . These values were determined after molecular weight calculation using the calibration curve, originating the molecular weight distribution curves shown in Figure A1.
Table A1. Averaged molecular weights (g mol 1 ) and polydispersity indexes for three silicone rubber samples (0 phr) as determined by GPC.
Table A1. Averaged molecular weights (g mol 1 ) and polydispersity indexes for three silicone rubber samples (0 phr) as determined by GPC.
Sample M n M w M z M z + 1 M w M n M z M w
A373,887677,0371,043,5061,397,8821.810811.54128
B365,736650,0121,007,8521,356,4951.777271.55051
C358,198647,9471,015,9851,377,6671.808911.56801
Mean365,940658,3301,022,4501,377,3501.7991.553
± std. dev±6406±13,253±15,256±16,896±0.015±0.011
Figure A1. Molecular weight distribution curves for three (ac) silicone rubber samples (0 phr). The notations along the curves are, from left to right, the number average, the weight average, and the higher average molecular weights.
Figure A1. Molecular weight distribution curves for three (ac) silicone rubber samples (0 phr). The notations along the curves are, from left to right, the number average, the weight average, and the higher average molecular weights.
Polymers 14 04404 g0a1
Values for averaged molecular weights and polydispersity indexes are in accordance with typical values for solid poly(siloxane) elastomers [54]. The polydispersity M w M n of around 2 is also typical for polycondensation polymerization reactions, one of the normally employed synthesis route for poly(siloxane)s [55].

Appendix A.2. Determination of the Filler Content

The silicone rubber grade under investigation is filled with an inorganic filler (as described in Section 2.1) and its content was determined by thermogravimetric analysis (TGA). TGA was performed in a TGA/DSC 1 device (Mettler Toledo, Greifensee, Switzerland) under 20 mL min 1 of N 2 at 10 K min 1 from 25 C to 900 C. Filler content was determined as the mass residue at 900 C, since pure silicone rubber gum experiences no residue at this temperature according to the literature [56]. The thermogram for the sample without dicumylperoxide is shown in Figure A2, where it is possible to observe the high thermal stability of silicone rubbers [57], and the residue at 900 C.
Figure A2. Thermogram for the 0 phr dicumylperoxide silicone rubber under inert atmosphere and 10 K min 1 heating rate.
Figure A2. Thermogram for the 0 phr dicumylperoxide silicone rubber under inert atmosphere and 10 K min 1 heating rate.
Polymers 14 04404 g0a2

Appendix A.3. Determination of the Vinyl Content

Determination of the vinyl presence or absence in the silicone rubber under study was vital to correctly describe the crosslinking mechanism taking the thermodynamic and kinetic parameters that were calculated into consideration. The vinyl content was then determined via high-field solution proton nuclear magnetic resonance spectroscopy ( 1 H NMR), by analysing the vinyl signals at about 6 ppm. The silicone rubber gum without dicumylperoxide (0 phr sample) was solubilized in deuterated chloroform ( C D C l 3 ) and analysed in a 300 MHz Bruker Advance III spectrometer (Bruker, Billerica, MA, USA) with 64 scans and reference as the residual solvent peak at 7.26 ppm.
Figure A3 shows the 1 H NMR full spectrum, as well as a detailed 1 H NMR spectrum with focus on the 5.5–6.5 ppm chemical shift ( δ ) region (usually where vinyl peaks show up [58]), for the 0 phr dicumylperoxide solid silicone sample. The full spectrum shows characteristic poly(dimethylsiloxane) signals, mainly at −0.20–0.20 ppm. Peaks at the 1.5 ppm region can be associated to water. No significant vinyl signals were present at the 5.5–6.5 ppm region. Normally, even at low concentrations (down to 2%, for example), vinyl groups appear as a set of multiple peaks at about 6 ppm. Thus, it is possible to state that this silicone grade has no significant amount of vinyl moieties, being mainly composed of poly(dimethylsiloxane).
Figure A3. Full (a) 1 H NMR spectrum in C D C l 3 for the 0 dicumylperoxide solid silicone, and detailed spectrum (b) for the vinyl region.
Figure A3. Full (a) 1 H NMR spectrum in C D C l 3 for the 0 dicumylperoxide solid silicone, and detailed spectrum (b) for the vinyl region.
Polymers 14 04404 g0a3aPolymers 14 04404 g0a3b

References

  1. Shit, S.; Shah, P. A Review on Silicone Rubber. Natl. Acad. Sci. Lett. 2013, 36, 355–365. [Google Scholar] [CrossRef]
  2. Calderón, J.; López, D.; Pérez, E.; Montesinos, J. Polysiloxanes as polymer matrices in biomedical engineering: Their interesting properties as the reason for the use in medical sciences. Polym. Bull. 2020, 77, 2749–2817. [Google Scholar] [CrossRef]
  3. Fink, J. Liquid Silicone Rubber: Chemistry, Materials, and Processing; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  4. Kim, D.; Park, J.; Lee, D.; Seo, K. Correlation between the Crosslink Characteristics and Mechanical Properties of Natural Rubber Compound via Accelerators and Reinforcement. Polymers 2020, 12, 2020. [Google Scholar] [CrossRef]
  5. Nardelli, F.; Calucci, L.; Carignani, E.; Borsacchi, S.; Cettolin, M.; Arimondi, M.; Giannini, L.; Geppi, M.; Martini, F. Influence of Sulfur-Curing Conditions on the Dynamics and Crosslinking of Rubber NEtworks: A Time-Domain NMR Study. Polymers 2022, 14, 767. [Google Scholar] [CrossRef] [PubMed]
  6. Kruželák, J.; Hložeková, K.; Kvasnicáková, A.; Tomanová, K.; Hudec, I. Application of Sulfur and Peroxide Curing Systems for Cross-Linking of Rubber Composites Filled with Calcium Lignosulfonate. Polymers 2022, 14, 1921. [Google Scholar] [CrossRef]
  7. Heideman, G.; Datta, R.; Noodermeer, J.; van Baarle, B. Activators in Accelerated Sulfur Vulcanization. Rubber Chem. Technol. 2004, 77, 512–541. [Google Scholar] [CrossRef]
  8. Traintinger, M.; Kerschbaumer, R.; Lechner, B.; Friesenbichler, W.; Lucyshyn, T. Temperature Profile in Rubber Injection Molding: Application of a Recently Developed Testing Method to Improve the Process Simulation and Calculation of Curing Kinetics. Polymers 2021, 13, 380. [Google Scholar] [CrossRef]
  9. Stieger, S.; Mitsoulis, E.; Walluch, M.; Ebner, C.; Kerschbaumer, R.; Haselmann, M.; Mostafaiyan, M.; Kämpfe, M.; Kühnert, I.; Wießner, S.; et al. On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber. Polymers 2021, 13, 2323. [Google Scholar] [CrossRef]
  10. Bont, M.; Barry, C.; Johnston, S. A review of liquid silicone rubber injection molding: Process variables and process modeling. Polym. Eng. Sci. 2021, 61, 331–347. [Google Scholar] [CrossRef]
  11. Ou, H.; Sahli, M.; Barrière, T.; Gelin, J.C. Multiphysics modelling and experimenal investigations of the filling and curing phases of bi-injection moulding of thermoplastics polymer/liquid silicone rubbers. Polym. Eng. Sci. 2017, 92, 3871–3882. [Google Scholar]
  12. Lucyshyn, T.; D’Avernas, L.V.; Holzer, C. Influence of the Mold Property on the Injection Molding Cycle Time and Warpage Depending on the Polymer Processed. Polymers 2021, 13, 3196. [Google Scholar] [CrossRef]
  13. Bardelli, T.; Marano, C.; Vangosa, F. Polydimethylsiloxane crosslinking kinetics: A systematic study on Sylgard184 comparing rheological and thermal approaches. J. Appl. Polym. Sci. 2021, 138, 51013. [Google Scholar] [CrossRef]
  14. Hong, I.K.; Lee, S. Cure kinetics and modeling the reaction of silicone rubber. J. Ind. Eng. Chem. 2013, 19, 42–47. [Google Scholar] [CrossRef]
  15. Jalkanen, T. Isoconversional kinetic analysis for determining the rate of cross-linking for Pt and peroxide cure silicone rubbers. Thermochim. Acta 2021, 703, 178982. [Google Scholar] [CrossRef]
  16. Karami, Z.; Jazani, O.; Navarchian, A.; Saeb, M. Cure Kinetics of Silicone/Halloysite Nanotube Composites. J. Vinyl Addit. Technol. 2020, 26, 548–565. [Google Scholar] [CrossRef]
  17. Ke, Q.; Wu, C.; Chen, X. Model-free cure kinetics of additional liquid silicone rubber. Thermochim. Acta 2020, 688, 178584. [Google Scholar] [CrossRef]
  18. Hernández-Ortiz, J.; Osswald, T. Modeling processing of silicone rubber: Liquid versus hard silicone rubbers. J. Appl. Polym. Sci. 2011, 119, 1864–1871. [Google Scholar] [CrossRef]
  19. ASTM D5289-19a; Standard Test Method for Rubber Property—Vulcanization Using Rotorless Cure Meters. American Society for Testing and Materials: West Conshohocke, PA, USA, 2019.
  20. ASTM D6204-19a; Standard Test Method for Rubber—Measurement of Unvulcanized Rheological Properties Using Rotorless Shear Rheometers. American Society for Testing and Materials: West Conshohocke, PA, USA, 2019.
  21. ASTM D6601-19; Standard Test Method for Rubber Properties—Measurement of Cure and After—Cure Dynamic Properties Using a Rotorless Shear Rheometer. American Society for Testing and Materials: West Conshohocke, PA, USA, 2019.
  22. ASTM D8059-19; Standard Test Method for Rubber Compounds—Measurement of Unvulcanized Dynamic Strain Softening (Payne Effect) Using Sealed Cavity Rotorless Shear Rheometers. American Society for Testing and Materials: West Conshohocke, PA, USA, 2019.
  23. Marquardt, D. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Indust. App. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  24. ASTM E2041-13(2018); Standard Test Method for Estimating Kinetic Parameters by Differential Scanning Calorimeter Using the Borchardt and Daniels Method. American Society for Testing and Materials: West Conshohocke, PA, USA, 2018.
  25. ASTM E698-18; Standard Test Method for Kinetic Parameters for Thermally Unstable Materials Using Differential Scanning Calorimetry and the Flynn/Wall/Ozawa Method. American Society for Testing and Materials: West Conshohocke, PA, USA, 2018.
  26. ASTM E2781-11; Standard Practice for Evaluation of Methods for Determination of Kinetic Parameters by Thermal Analysis. American Society for Testing and Materials: West Conshohocke, PA, USA, 2016.
  27. Vyazovkin, S.; Burnham, A.; Criado, J.; Pérez-Maqueda, L.; Popescu, C.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta 2011, 520, 1–19. [Google Scholar] [CrossRef]
  28. Vyazovkin, S. A time to search: Finding the meaning of variable activation energy. Phys. Chem. Chem. Phys. 2016, 18, 18643–18656. [Google Scholar] [CrossRef]
  29. Vyazovkin, S. Kissinger Method in Kinetics of Materials: Things to Beware and Be Aware of. Molecules 2020, 25, 2813. [Google Scholar] [CrossRef] [PubMed]
  30. Vyazovkin, S.; Burnham, A.; Favergeon, L.; Koga, N.; Moukhina, E.; Pérez-Maqueda, L.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for analysis of multi-step kinetics. Thermochim. Acta 2020, 689, 178597. [Google Scholar] [CrossRef]
  31. Vyazovkin, S. Determining Preexponential Factor in Model-Free Kinetic Methods: How and Why? Molecules 2021, 26, 3077. [Google Scholar] [CrossRef] [PubMed]
  32. Lang, M.; Hirner, S.; Wiesbrock, F.; Fuchs, P. A Review on Modeling Cure Kinetics and Mechanisms of Photopolymerization. Polymers 2022, 14, 2074. [Google Scholar] [CrossRef]
  33. Zhang, X. Application of Kinetic Methods in Thermal Analysis: A Review. Eng. Sci. 2021, 14, 1–13. [Google Scholar] [CrossRef]
  34. Friedman, H. Kinetics of thermal degradation of char-forming plastics from thermogravimetry. Application to a phenolic plastic. J. Polym. Sci. C Polym. Symp. 1964, 6, 183–195. [Google Scholar] [CrossRef]
  35. Šesták, J. Šesták–Berggren equation: Now questioned but formerly celebrated—what is right. J. Therm. Anal. Calorim. 2017, 127, 1117–1121. [Google Scholar] [CrossRef]
  36. Kamal, M. Thermoset characterization for moldability analysis. Polym. Eng. Sci. 1974, 14, 231–239. [Google Scholar] [CrossRef]
  37. Dluzneski, P. Peroxide vulcanization of elastomers. Rubber Chem. Technol. 2000, 74, 451–492. [Google Scholar] [CrossRef]
  38. Baquey, G.; Moine, L.; Babot, O.; Degueil, M.; Maillard, B. Model study of the crosslinking of polydimethylsiloxanes by peroxides. Polymer 2005, 46, 6283–6292. [Google Scholar]
  39. Wiśniewska, P.; Zedler, Ł.; Formela, K. Processing, Performance Properties, and Storage Stability of Ground Tire Rubber Modified by Dicumyl Peroxide and Ethylene-Vinyl Acetate Copolymers. Polymer 2021, 13, 4014. [Google Scholar] [CrossRef] [PubMed]
  40. Likozar, B.; Krajnc, M. Simulation of chemical kinetics of elastomer crosslinking by organic peroxides. Polym. Eng. Sci. 2009, 49, 60–72. [Google Scholar] [CrossRef]
  41. Kruželák, J.; Hakošová, S.; Kvasničáková, A.; Hudec, I. Dicumyl Peroxide used as curing Agent for different Typer of Rubber Matrices Part I: Effect of Temperature. KGK Kautsch. Gummi Kunststoffe 2020, 73, 36–42. [Google Scholar]
  42. Kruželák, J.; Hakošová, S.; Kvasničáková, A.; Hudec, I. Dicumyl Peroxide used as curing Agent for different Typer of Rubber Matrices Part II: Effect of peroxide Content. KGK Kautsch. Gummi Kunststoffe 2020, 73, 36–40. [Google Scholar]
  43. Duh, Y.S.; Kao, C.S.; Lee, W.L.W. Chemical kinetics on thermal decompositions of dicumyl peroxide studies by calorimetry. J. Therm. Anal. Calorim. 2017, 127, 1089–1098. [Google Scholar] [CrossRef]
  44. Lv, J.; Chen, L.; Chen, W.; Gao, H.; Peng, M. Kinetic analysis and self-accelerating decomposition temperature (SADT) of dicumyl peroxide. Thermochim. Acta 2013, 571, 60–63. [Google Scholar] [CrossRef]
  45. Hou, Q.; Yin, L.; Xu, L.; Tan, J. Effects of composite reinforcing filler, vulcanizing temperature, and pressure on mechanical properties of gasket material for proton exchange membrane fuel cells. J. Appl. Polym. Sci. 2022, 139, 52298. [Google Scholar] [CrossRef]
  46. Bianchi, O.; Martins, J.; Fiorio, R.; Oliveira, R.; Canto, L. Changes in activation energy and kinetic mechanism during EVA crosslinking. Polym. Test. 2011, 30, 616–624. [Google Scholar] [CrossRef]
  47. Granado, L.; Tavernier, R.; Foyer, G.; David, G.; Caillol, S. Comparative curing kinetics study of high char yield formaldehyde- and terephthalaldehyde-phenolic thermosets. Thermochim. Acta 2018, 667, 42–49. [Google Scholar] [CrossRef]
  48. Cheheb, Z.; Mousseau, P.; Sarda, A.; Deterre, R. Thermal Conductivity of Rubber Compounds Versus the State of Cure. Macromol. Mater. Eng. 2012, 297, 228–236. [Google Scholar] [CrossRef]
  49. Glasstone, S.; Laidler, K.; Eyring, H. The Theory of the Rate Processes; McGraw-Hill Book Co.: New York, NY, USA, 1941. [Google Scholar]
  50. Parks, C.; Lorenz, O. Crosslinking efficiency in the reaction of dicumyl peroxide with dimmethyloctadiene. J. Polym. Sci. 1961, 50, 287–298. [Google Scholar] [CrossRef]
  51. Hioe, J.; Zipse, H. Radical stability and its role in synthesis and catalysis. Org. Biomol. Chem. 2010, 8, 3609–3617. [Google Scholar] [CrossRef]
  52. Bordwell, F.; Zhang, X.M.; Alnakkar, S. Effects of Adjacent Acceptors and Donors on the Stabilities of Carbon-Centerd Radicals. J. Am. Chem. Soc. 1992, 447, 7623–7629. [Google Scholar] [CrossRef]
  53. Verheyen, F.; Giesen, R.U.; Heim, H.P. Mechanical Properties of Polydimethylsiloxane as a Function of the Amount and Type of Crosslinking Agent. Int. Polym. Process. 2017, 32, 337–342. [Google Scholar] [CrossRef]
  54. Soroceanu, A.; Stiubianu, G. Siloxane Matrix Molecular Weight Influences the Properties of Nanocomposites Based on Metal Complexes and Dielectric Elastomer. Materials 2021, 14, 3352. [Google Scholar] [CrossRef] [PubMed]
  55. Yokoyaza, T.; Yokoyama, A. Chain-Growth Polycondensation: Living Polymerization Nature in Polycondensation and Approach to Condensation Polymer Architecture. Polym. J. 2004, 36, 65–83. [Google Scholar] [CrossRef] [Green Version]
  56. Gao, X.; Hiu, H.; Wei, H.; Zheng, J.; Huang, G. Effect of incompletely condensed tri-silanol-phenyl-POSS on the thermal stability of silicone rubber. Polym. Bull. 2019, 76, 2835–2850. [Google Scholar] [CrossRef]
  57. Shi, S.; Lei, B.; Li, M.; Cui, X.; Wang, X.; Fan, X.; Tang, S.; Shen, J. Thermal decomposition behaviour of a thermal protection coating composite with silicone rubber: Experiment and modelling. Prog. Org. Coat. 2020, 143, 105609. [Google Scholar] [CrossRef]
  58. Saalwächter, K.; Reichert, D. Polymer Applications of NMR in Encyclopedia of Spectroscopy and Spectrometry, 3rd ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2017. [Google Scholar]
Figure 1. Rotational rheometry (RPA, 0.5° deformation at 1.667 Hz) curves for silicone rubber compounds with varied dicumylperoxide concentrations (bg) at different crosslinking temperatures. The plot (a) (torque vs. time curve) is related to the pure solid silicone rubber without dicumylperoxide.
Figure 1. Rotational rheometry (RPA, 0.5° deformation at 1.667 Hz) curves for silicone rubber compounds with varied dicumylperoxide concentrations (bg) at different crosslinking temperatures. The plot (a) (torque vs. time curve) is related to the pure solid silicone rubber without dicumylperoxide.
Polymers 14 04404 g001aPolymers 14 04404 g001b
Figure 2. Optimum cure time t 90 as a function of the dicumylperoxide concentration and the temperature. The dots represent the experimental data and the surface indicates the model described at Equation (1), with fitting parameters as follows: A 0 (pre-exponential factor) = 2.02 × 10 14 min 1 , E a (activation energy) = 1.182 × 10 5 J mol 1 , and β = —0.603. The correlation factor R 2 was 0.988.
Figure 2. Optimum cure time t 90 as a function of the dicumylperoxide concentration and the temperature. The dots represent the experimental data and the surface indicates the model described at Equation (1), with fitting parameters as follows: A 0 (pre-exponential factor) = 2.02 × 10 14 min 1 , E a (activation energy) = 1.182 × 10 5 J mol 1 , and β = —0.603. The correlation factor R 2 was 0.988.
Polymers 14 04404 g002
Figure 3. Conversion rate d α d t at 160 C for silicone rubber compounds with different dicumylperoxide concentrations as a function of time (a) and conversion (b).
Figure 3. Conversion rate d α d t at 160 C for silicone rubber compounds with different dicumylperoxide concentrations as a function of time (a) and conversion (b).
Polymers 14 04404 g003
Figure 4. Friedman-like isoconversional approach for calculating the activation energy, regarding the sample 0.21 phr (a) and the activation energy for different dicumyperoxide concentration as function of conversion (b).
Figure 4. Friedman-like isoconversional approach for calculating the activation energy, regarding the sample 0.21 phr (a) and the activation energy for different dicumyperoxide concentration as function of conversion (b).
Polymers 14 04404 g004
Figure 5. Comparison between the experimental and the calculated conversion values for the crosslinking at 160 C for silicone rubber compounds with dicumylperoxide concentrations. The calculated conversion values are plotted from 0.1 to 0.9 of conversion.
Figure 5. Comparison between the experimental and the calculated conversion values for the crosslinking at 160 C for silicone rubber compounds with dicumylperoxide concentrations. The calculated conversion values are plotted from 0.1 to 0.9 of conversion.
Polymers 14 04404 g005
Table 1. Curing parameters for silicone rubber compounds concerning minimum torque ( M L ) and induction or scorch time ( t s 1 ) for different crosslinking temperatures and varied dicumylperoxide concentrations.
Table 1. Curing parameters for silicone rubber compounds concerning minimum torque ( M L ) and induction or scorch time ( t s 1 ) for different crosslinking temperatures and varied dicumylperoxide concentrations.
Temperature140  C150  C160  C170  C180  C
ML t s 1 ML t s 1 ML t s 1 ML t s 1 ML t s 1
(dN m)(min)(dN m)(min)(dN m)(min)(dN m)(min)(dN m)(min)
0 phr0.58n/a0.57n/a0.56n/a0.54n/a0.41n/a
0.21 phr0.6119.980.606.160.582.120.560.960.530.53
0.35 phr0.6211.130.623.630.581.290.550.640.550.39
0.49 phr0.609.400.593.020.571.160.550.580.540.36
0.70 phr0.618.060.602.630.581.010.560.530.540.35
1.00 phr0.613.760.581.310.560.630.550.570.560.28
1.50 phr0.582.620.561.050.530.570.570.370.530.28
Table 2. Curing parameters for silicone rubber compounds concerning maximum torque ( M H ) and optimum cure time ( t 90 ) for different crosslinking temperatures and varied dicumylperoxide concentrations.
Table 2. Curing parameters for silicone rubber compounds concerning maximum torque ( M H ) and optimum cure time ( t 90 ) for different crosslinking temperatures and varied dicumylperoxide concentrations.
Temperature140  C150  C160  C170  C180  C
MH t 90 MH t 90 MH t 90 MH t 90 MH t 90
(dN m)(min)(dN m)(min)(dN m)(min)(dN m)(min)(dN m)(min)
0 phr0.72n/a0.77n/a0.91n/a1.07n/a1.28n/a
0.21 phr7.7645.459.8920.9210.488.9111.254.5411.542.92
0.35 phr8.7231.1510.2612.9410.865.4311.572.9911.761.64
0.49 phr8.9327.079.9710.1810.874.7211.572.4411.641.25
0.70 phr9.424.379.958.9310.984.2111.572.1311.821.19
1.00 phr11.4016.7411.928.7512.394.512.562.9312.631.22
1.50 phr11.0512.5811.556.9711.784.1412.074.5312.061.04
Table 3. Kinetic parameters determined after fitting of the experimental conversion rate d α d t to the Kamal model (Equation (10)) for silicone rubber compounds crosslinked at 160 C with different dicumylperoxide concentrations.
Table 3. Kinetic parameters determined after fitting of the experimental conversion rate d α d t to the Kamal model (Equation (10)) for silicone rubber compounds crosslinked at 160 C with different dicumylperoxide concentrations.
A 1 A 2 E α , 1 E α , 2 mn k 1 k 2
Sample(min 1 )(min 1 )(kJ mol 1 )(kJ mol 1 )(min 1 )(min 1 )
0.21 phr 1.1 × 10 15 6.4 × 10 15 131.35131.351.1671.6170.15910.9261
0.25 phr 6.4 × 10 14 3.8 × 10 15 127.47127.471.2981.6210.27201.6150
0.49 phr 4.4 × 10 14 2.4 × 10 15 125.75125.751.2101.5520.30151.6444
0.70 phr 3.4 × 10 14 1.9 × 10 15 123.92123.921.3501.6670.38722.1638
1.00 phr 3.2 × 10 10 4.2 × 10 11 87.6787.602.0482.7600.871511.438
1.50 phr 1.0 × 10 10 1.1 × 10 11 83.6183.611.8202.6100.82699.0960
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Azevedo, M.; Monks, A.-M.; Kerschbaumer, R.C.; Schlögl, S.; Holzer, C. Peroxide-Based Crosslinking of Solid Silicone Rubber, Part I: Insights into the Influence of Dicumylperoxide Concentration on the Curing Kinetics and Thermodynamics Determined by a Rheological Approach. Polymers 2022, 14, 4404. https://doi.org/10.3390/polym14204404

AMA Style

Azevedo M, Monks A-M, Kerschbaumer RC, Schlögl S, Holzer C. Peroxide-Based Crosslinking of Solid Silicone Rubber, Part I: Insights into the Influence of Dicumylperoxide Concentration on the Curing Kinetics and Thermodynamics Determined by a Rheological Approach. Polymers. 2022; 14(20):4404. https://doi.org/10.3390/polym14204404

Chicago/Turabian Style

Azevedo, Maurício, Anna-Maria Monks, Roman C. Kerschbaumer, Sandra Schlögl, and Clemens Holzer. 2022. "Peroxide-Based Crosslinking of Solid Silicone Rubber, Part I: Insights into the Influence of Dicumylperoxide Concentration on the Curing Kinetics and Thermodynamics Determined by a Rheological Approach" Polymers 14, no. 20: 4404. https://doi.org/10.3390/polym14204404

APA Style

Azevedo, M., Monks, A. -M., Kerschbaumer, R. C., Schlögl, S., & Holzer, C. (2022). Peroxide-Based Crosslinking of Solid Silicone Rubber, Part I: Insights into the Influence of Dicumylperoxide Concentration on the Curing Kinetics and Thermodynamics Determined by a Rheological Approach. Polymers, 14(20), 4404. https://doi.org/10.3390/polym14204404

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