# Characterization of Dielectric Relaxation Process by Impedance Spectroscopy for Polymers: Nitrile Butadiene Rubber and Ethylene Propylene Diene Monomer

Authors: Jae Kap Jung; Young Il Moon; Gyung Hyun Kim; Nae Hyung Tak
Journal: Journal of Spectroscopy
 (2020)

---

## Abstract

We invented a dispersion analysis program that analyzes the relaxation processes from dielectric permittivity based on a combination of the Havriliak–Negami and conductivity contribution functions. By applying the created program to polymers such as nitrile butadiene rubber (NBR) and ethylene propylene diene monomer (EPDM), several relaxation processes were characterized: anα process due to segmental motions of the C-C bond, an α′ process attributed to fluctuations in the end-to-end dipole vector of the polymer chain, the conduction contribution by the filler observed above room temperature, and secondary relaxation processes β and γ of motion for the side group in NBR. In the EPDM specimen, the β process associated with the rotational motion of the side groups, the α process associated with the relaxation of local segmental motion, and the αβ process associated with the origin of the β process at high temperatures above 305 K were observed. The Maxwell–Wagner–Sillars effect and conduction contribution were also presented. The molecular chains responsible for the relaxation processes were assigned by building molecular models of the two polymers. The temperature dependence of the relaxation strength and the shape parameters that characterize the process were investigated. From the temperature-dependent relaxation analysis, the merged αβ process, activation energy, and glass transition temperature were determined and compared.


---

## Body

## 1. Introduction

Many kinds of rubbery polymers are used in O-rings for hydrogen gas sealing under high-pressure hydrogen infrastructure conditions; nitrile butadiene rubber (NBR), ethylene propylene diene monomer (EPDM), and fluoroelastomer (FKM) rubbers are used as sealants and gaskets in valves and pipelines [1–6]. Previous studies have investigated physical properties such as the glass transition, volume swelling, and thermal stability, mechanical properties such as tensile strength, compression, elongation, and shear behavior at high hydrostatic pressure, and transport properties such as permeability, diffusivity, and solubility of hydrogen.Dielectric spectroscopy provides a means to characterize the dielectric and electrical behaviors of polymers and has been widely employed to investigate the nature of molecular motion and related relaxation behaviors [7–13]. The dielectric response is usually controlled by α relaxation, whose relaxation time steeply increases when the temperature decreases toward the glass transition temperature, Tg. The α process corresponds to the segmental motions of the main chains. In addition to α relaxation, additional relaxations called secondary relaxations or β relaxations can be active on faster time scales [14–17].Only a few dielectric studies have been performed on NBR. Dielectric relaxation spectroscopy and the dynamic mechanical analysis of NBR composites have identifiedα and α′ relaxation processes. The observed phenomena and their temperature dependence have been investigated in view of the molecular dynamics of the chains, filler effects, and glass transition [18]. The observed relaxations follow the Havriliak–Negami (HN) approach for the distribution of relaxation times. However, relaxations related to conductivity and secondary relaxations have not been detected in NBR. The α relaxation peaks of EPDM have been detected in the graph of dielectric tan δ data versus temperature, where the height and position of the peaks changed with the thermal aging time [19]. The α process for neat EPDM nanocomposites was observed only in the measurement of the imaginary dielectric permittivity versus the frequency, in which the dielectric spectra were well fitted by the HN permittivity function. However, the β relaxation peak and other relaxations were not detected [20].Research on dielectric spectroscopy has been performed with limited observations of relaxation processes and has not revealed all the relaxations embedded in polymers. The interfacial or Maxwell–Wagner–Sillars (MWS) relaxation process occurs in heterogeneous systems, including the conduction contribution in two polymers, which has not been investigated until now. Moreover, relatively less attention has been given to the investigation of the dielectric and electrical behaviors of polymers, especially at temperatures higher than the glass transition temperatureTg. The dielectric spectrum is mainly governed by the relaxation time for relaxation processes in polymers. This spectrum occurs over a wide frequency range due to line broadening from distributions in the relaxation times and the superposition of corresponding relaxation processes [21, 22]. Therefore, it is difficult to fully discern and characterize the overall relaxation processes embedded in the studied polymers.This work decomposes the superimposed dielectric spectrum and assigns the peaks of the relaxation processes in the dielectric spectrum, resulting in the characterization of several relaxation processes with a developed dispersion analysis program [23]. The upgraded algorithm is used to analyze the relaxation processes by selectively satisfying either real or imaginary dielectric data based on three governing functions—the Cole and Cole [24], Davidson and Cole [25], and HN [26] functions—together with the conduction contribution. The upgraded program is applied to characterize the entire relaxation processes that exist in the two polymers—NBR and EPDM—under study. As the extension of dielectric spectroscopy research on NBR was investigated recently [23], research on EFDM has been completed. The present work was performed for NBR specimens produced from both different batches and compositions in [23]. Furthermore, this research is extended to lower frequencies reaching 0.01 Hz. These investigated results for both polymers are discussed and reviewed. The molecular chains of the two polymers responsible for the relaxation processes are assigned by molecular modeling. The temperature dependence of the relaxation strength and shape parameters for the relaxation spectra are investigated. Upon finding the α and β relaxations, conduction levels and MWS contributions in the two polymers, the activation energy, and glass transition temperature are obtained and discussed. The origin of the merged αβ process in EPDM was first revealed from the temperature dependence of the relaxation spectra and relaxation strength.


## 2. Materials and Methods

### 2.1. Sample Preparation

NBR is a butadiene-based synthetic rubber made by copolymerizing combinations of butadiene (CH2CH=CHCH2) and acrylonitrile (CH2CHCN). The NBR specimen used in this work was synthesized by a Korean domestic company and contained 50% carbon black as a filler. The acrylonitrile content amounts to 34% of the total amount of NBR. The glass transition temperature varies depending on the acrylonitrile content. NBR is widely used as a sealing material due to its excellent gas resistance [27]. Thus, NBR polymer materials have been used as O-ring seals for flange connections, threaded connectors, and various valves in high-pressure hydrogen infrastructures.EPDM is an ethylene-based synthetic rubber made by copolymerizing various combinations of ethylene (CH2CH2), propylene (CH2CHCH3), and ENB monomer. EPDM elastomers have excellent heat, ozone, weathering, and aging resistance [28]. EPDMs can be used in a wide range of applications, and typical applications include radiators, heater hoses, O-rings, and gaskets. The EPDM specimen used herein was also synthesized by a domestic Korean company, and 34% carbon black was included as a filler during the fabrication of the EPDM specimen. The chemical compositions of both the NBR and EPDM specimens are summarized in Table 1 [29].Table 1
Chemical compositions of NBR and EPDM specimens.
FunctionNBREPDMPolymerNBR (40)∗EPDM (58)Filler-reinforcingCarbon black (50)Carbon black (34)Processing aid1, 2-benzenedicarboxylic acid (6)—Antioxidant2-benzimidazolethiol (2)—Curing agentSulfur (2)Zinc oxide (3)Dicumyl peroxide (5)∗Numbers in () are weight ratios in %.


### 2.2. Impedance Spectroscopy System

The impedance spectroscopy system used to measure the complex dielectric permittivity (Figure1) is composed of a temperature-controlled chamber, a digital multimeter (Fluke 8842 A) for temperature monitoring, and an impedance analyzer (VSP-300) with a GPIB interface to a PC. The temperature is controlled and maintained by a heater and a circulating refrigerating medium varying at an interval of 5 K with temperature stability better than 0.2 K in a temperature range of 233 K to 404 K. The disc-shaped specimens (50 mm in diameter and 2.4 mm in thickness) are placed between two copper electrodes for this measurement.Figure 1
Configuration of the automatic impedance spectroscopy system for measuring the complex dielectric permittivity.Both the real and imaginary permittivities from the measurements using the impedance analyzer are calculated with the following relation:(1)ε′=dωε0AZ″Z′2+Z″2,ε″=dωε0AZ′Z′2+Z″2,where impedance Z=Z′−jZ″ and Z′ and Z″ are the real and imaginary impedance values, respectively. A is the area of the specimen, and d is the thickness of the cylindrically shaped specimen. ε0 is the vacuum permittivity of 8.854 × 10−12 F/m.


### 2.3. Dispersion Analysis Program

A dispersion analysis program for analyzing the complex dielectric permittivity in various polymers was developed. The program was written in C# in the Microsoft Visual Studio Net 2018 graphical user interface (GUI). It can be executed on a Microsoft Windows system. Figure2 shows the algorithm of the dispersion analysis program, which employs parallel Nelder-Mead optimization. The original data were composed of a set of dielectric spectra and the frequency measured at each temperature. Basically, one isothermal spectrum can be fitted with a single model function.Figure 2
Algorithm of the dispersion analysis program, which employs parallel Nelder-Mead optimization.


### 2.4. Analysis Function and Relaxation Process

#### 2.4.1. Analysis Function

To understand the characteristics of the relaxation processes in the polymers, the experimental complex dielectric permittivity,ε∗ω, can be fitted with an HN function. The HN function can be phenomenologically described as a combination of the conductivity term with the HN functional form as follows [26]:(2)ε∗ω=ε′−iε″=−iσdcε0ωN+ε∞+∑k=1,2,…Δεk1+iωτHNk1−akbk,where the first term σdc/ε0ωN corresponds to the conductivity, σdc is the DC conductivity, and N is an exponent that characterizes the conduction process. Its contribution has the strongest influence at the highest temperature and the lowest frequency. The last term is the HN function. The summation symbol indicates the existence of more than one relaxation process. τHNk is the average characteristic relaxation time of the corresponding process k. The exponents ak and bk (0≤ak≤1,bk≤1) are the HN fitting shape parameters describing the distributions and shape of the relaxation times. The values of ak and bk are related to the linewidth and asymmetry, respectively, of the corresponding relaxation loss peak. As the values of a and b are close to 1, the linewidth and shape become broader and more symmetric, respectively. The dielectric relaxation strength, εs−ε∞, is represented by Δεk.


#### 2.4.2. Relaxation Process

Theα relaxation responsible for the segmental motions of the main chains in the polymer is related to the glass transition [30]. The α relaxation reveals the following features: the asymmetric loss peak is often well described by the empirical HN function in the frequency domain, and the relaxation strength, Δεα, increases with decreasing temperature. α relaxation is located at lower frequencies or higher temperatures than β relaxation, indicating the stronger temperature dependency of the former. For cooperative processes, the temperature dependence of the relaxation rate obeys the Vogel–Fulcher–Tamman (VFT) law as follows [31–35]:(3)logfα=logf∞,α−AT−T0,where f∞,α is the relaxation rate at an infinite frequency, A is a material-specific parameter, T is the temperature, and T0 is the Vogel temperature, which is related to the dynamic glass transition process of the polymer and is usually less than Tg.β relaxation is the secondary relaxation process in amorphous polymers [36]. β relaxation originates from localized fluctuations of parts of the main chain or hindered rotations of polymer side groups. For most polymers, β relaxation has been found below the glass transition temperature (Tg). The β relaxation mechanism is a thermally activated process characterized by an Arrhenius temperature dependence of the relaxation rates [37, 38].(4)fβ=f∞,βexp−EakT.Ea represents the potential barrier (activation energy) between two possible equilibrium states. The activation energy Ea, related to the slope of logfβ versus 1/T, depends on both the internal rotation barriers and the environment of a moving unit. The dielectric strength (Δεβ) of β relaxation normally increases with temperature. In the frequency domain, dielectric β relaxation displays a broad and symmetrically shaped loss peak with a half-width that generally narrows with increasing temperature.The interfacial polarization effects are often called MWS effects [39, 40]. The effects in heterogeneous systems were first highlighted by Maxwell and later modified successively by Wagner and Sillars. Interfacial polarization processes occur in heterogeneous dielectrics, especially at low frequencies, as a result of the build-up of space charges at interfaces between two media having different permittivity and conductivity values. The relaxation is also characterized by an Arrhenius temperature dependence of the relaxation rates, as shown in equation (4).Moreover, the DC conductivity can be considered the same transport phenomena in polymers by considering the movement of different moving units (chain segment, dipole, and charge carrier). Similar to theα relaxation time, the temperature dependence of DC conductivity can be described by the VFT equation [41, 42]:(5)σdc=σ0exp−B′T−T0′,where B′ is a constant, σ0 is the conductivity at a very high temperature, and T0′ is the temperature at which the conductivity ideally reduces to zero.The temperature dependence of DC conductivity can also be described by the Arrhenius equation [43]:(6)σdc=σ0exp−EakT.


## 2.1. Sample Preparation

NBR is a butadiene-based synthetic rubber made by copolymerizing combinations of butadiene (CH2CH=CHCH2) and acrylonitrile (CH2CHCN). The NBR specimen used in this work was synthesized by a Korean domestic company and contained 50% carbon black as a filler. The acrylonitrile content amounts to 34% of the total amount of NBR. The glass transition temperature varies depending on the acrylonitrile content. NBR is widely used as a sealing material due to its excellent gas resistance [27]. Thus, NBR polymer materials have been used as O-ring seals for flange connections, threaded connectors, and various valves in high-pressure hydrogen infrastructures.EPDM is an ethylene-based synthetic rubber made by copolymerizing various combinations of ethylene (CH2CH2), propylene (CH2CHCH3), and ENB monomer. EPDM elastomers have excellent heat, ozone, weathering, and aging resistance [28]. EPDMs can be used in a wide range of applications, and typical applications include radiators, heater hoses, O-rings, and gaskets. The EPDM specimen used herein was also synthesized by a domestic Korean company, and 34% carbon black was included as a filler during the fabrication of the EPDM specimen. The chemical compositions of both the NBR and EPDM specimens are summarized in Table 1 [29].Table 1
Chemical compositions of NBR and EPDM specimens.
FunctionNBREPDMPolymerNBR (40)∗EPDM (58)Filler-reinforcingCarbon black (50)Carbon black (34)Processing aid1, 2-benzenedicarboxylic acid (6)—Antioxidant2-benzimidazolethiol (2)—Curing agentSulfur (2)Zinc oxide (3)Dicumyl peroxide (5)∗Numbers in () are weight ratios in %.


## 2.2. Impedance Spectroscopy System

The impedance spectroscopy system used to measure the complex dielectric permittivity (Figure1) is composed of a temperature-controlled chamber, a digital multimeter (Fluke 8842 A) for temperature monitoring, and an impedance analyzer (VSP-300) with a GPIB interface to a PC. The temperature is controlled and maintained by a heater and a circulating refrigerating medium varying at an interval of 5 K with temperature stability better than 0.2 K in a temperature range of 233 K to 404 K. The disc-shaped specimens (50 mm in diameter and 2.4 mm in thickness) are placed between two copper electrodes for this measurement.Figure 1
Configuration of the automatic impedance spectroscopy system for measuring the complex dielectric permittivity.Both the real and imaginary permittivities from the measurements using the impedance analyzer are calculated with the following relation:(1)ε′=dωε0AZ″Z′2+Z″2,ε″=dωε0AZ′Z′2+Z″2,where impedance Z=Z′−jZ″ and Z′ and Z″ are the real and imaginary impedance values, respectively. A is the area of the specimen, and d is the thickness of the cylindrically shaped specimen. ε0 is the vacuum permittivity of 8.854 × 10−12 F/m.


## 2.3. Dispersion Analysis Program

A dispersion analysis program for analyzing the complex dielectric permittivity in various polymers was developed. The program was written in C# in the Microsoft Visual Studio Net 2018 graphical user interface (GUI). It can be executed on a Microsoft Windows system. Figure2 shows the algorithm of the dispersion analysis program, which employs parallel Nelder-Mead optimization. The original data were composed of a set of dielectric spectra and the frequency measured at each temperature. Basically, one isothermal spectrum can be fitted with a single model function.Figure 2
Algorithm of the dispersion analysis program, which employs parallel Nelder-Mead optimization.


## 2.4. Analysis Function and Relaxation Process

### 2.4.1. Analysis Function

To understand the characteristics of the relaxation processes in the polymers, the experimental complex dielectric permittivity,ε∗ω, can be fitted with an HN function. The HN function can be phenomenologically described as a combination of the conductivity term with the HN functional form as follows [26]:(2)ε∗ω=ε′−iε″=−iσdcε0ωN+ε∞+∑k=1,2,…Δεk1+iωτHNk1−akbk,where the first term σdc/ε0ωN corresponds to the conductivity, σdc is the DC conductivity, and N is an exponent that characterizes the conduction process. Its contribution has the strongest influence at the highest temperature and the lowest frequency. The last term is the HN function. The summation symbol indicates the existence of more than one relaxation process. τHNk is the average characteristic relaxation time of the corresponding process k. The exponents ak and bk (0≤ak≤1,bk≤1) are the HN fitting shape parameters describing the distributions and shape of the relaxation times. The values of ak and bk are related to the linewidth and asymmetry, respectively, of the corresponding relaxation loss peak. As the values of a and b are close to 1, the linewidth and shape become broader and more symmetric, respectively. The dielectric relaxation strength, εs−ε∞, is represented by Δεk.


### 2.4.2. Relaxation Process

Theα relaxation responsible for the segmental motions of the main chains in the polymer is related to the glass transition [30]. The α relaxation reveals the following features: the asymmetric loss peak is often well described by the empirical HN function in the frequency domain, and the relaxation strength, Δεα, increases with decreasing temperature. α relaxation is located at lower frequencies or higher temperatures than β relaxation, indicating the stronger temperature dependency of the former. For cooperative processes, the temperature dependence of the relaxation rate obeys the Vogel–Fulcher–Tamman (VFT) law as follows [31–35]:(3)logfα=logf∞,α−AT−T0,where f∞,α is the relaxation rate at an infinite frequency, A is a material-specific parameter, T is the temperature, and T0 is the Vogel temperature, which is related to the dynamic glass transition process of the polymer and is usually less than Tg.β relaxation is the secondary relaxation process in amorphous polymers [36]. β relaxation originates from localized fluctuations of parts of the main chain or hindered rotations of polymer side groups. For most polymers, β relaxation has been found below the glass transition temperature (Tg). The β relaxation mechanism is a thermally activated process characterized by an Arrhenius temperature dependence of the relaxation rates [37, 38].(4)fβ=f∞,βexp−EakT.Ea represents the potential barrier (activation energy) between two possible equilibrium states. The activation energy Ea, related to the slope of logfβ versus 1/T, depends on both the internal rotation barriers and the environment of a moving unit. The dielectric strength (Δεβ) of β relaxation normally increases with temperature. In the frequency domain, dielectric β relaxation displays a broad and symmetrically shaped loss peak with a half-width that generally narrows with increasing temperature.The interfacial polarization effects are often called MWS effects [39, 40]. The effects in heterogeneous systems were first highlighted by Maxwell and later modified successively by Wagner and Sillars. Interfacial polarization processes occur in heterogeneous dielectrics, especially at low frequencies, as a result of the build-up of space charges at interfaces between two media having different permittivity and conductivity values. The relaxation is also characterized by an Arrhenius temperature dependence of the relaxation rates, as shown in equation (4).Moreover, the DC conductivity can be considered the same transport phenomena in polymers by considering the movement of different moving units (chain segment, dipole, and charge carrier). Similar to theα relaxation time, the temperature dependence of DC conductivity can be described by the VFT equation [41, 42]:(5)σdc=σ0exp−B′T−T0′,where B′ is a constant, σ0 is the conductivity at a very high temperature, and T0′ is the temperature at which the conductivity ideally reduces to zero.The temperature dependence of DC conductivity can also be described by the Arrhenius equation [43]:(6)σdc=σ0exp−EakT.


## 2.4.1. Analysis Function

To understand the characteristics of the relaxation processes in the polymers, the experimental complex dielectric permittivity,ε∗ω, can be fitted with an HN function. The HN function can be phenomenologically described as a combination of the conductivity term with the HN functional form as follows [26]:(2)ε∗ω=ε′−iε″=−iσdcε0ωN+ε∞+∑k=1,2,…Δεk1+iωτHNk1−akbk,where the first term σdc/ε0ωN corresponds to the conductivity, σdc is the DC conductivity, and N is an exponent that characterizes the conduction process. Its contribution has the strongest influence at the highest temperature and the lowest frequency. The last term is the HN function. The summation symbol indicates the existence of more than one relaxation process. τHNk is the average characteristic relaxation time of the corresponding process k. The exponents ak and bk (0≤ak≤1,bk≤1) are the HN fitting shape parameters describing the distributions and shape of the relaxation times. The values of ak and bk are related to the linewidth and asymmetry, respectively, of the corresponding relaxation loss peak. As the values of a and b are close to 1, the linewidth and shape become broader and more symmetric, respectively. The dielectric relaxation strength, εs−ε∞, is represented by Δεk.


## 2.4.2. Relaxation Process

Theα relaxation responsible for the segmental motions of the main chains in the polymer is related to the glass transition [30]. The α relaxation reveals the following features: the asymmetric loss peak is often well described by the empirical HN function in the frequency domain, and the relaxation strength, Δεα, increases with decreasing temperature. α relaxation is located at lower frequencies or higher temperatures than β relaxation, indicating the stronger temperature dependency of the former. For cooperative processes, the temperature dependence of the relaxation rate obeys the Vogel–Fulcher–Tamman (VFT) law as follows [31–35]:(3)logfα=logf∞,α−AT−T0,where f∞,α is the relaxation rate at an infinite frequency, A is a material-specific parameter, T is the temperature, and T0 is the Vogel temperature, which is related to the dynamic glass transition process of the polymer and is usually less than Tg.β relaxation is the secondary relaxation process in amorphous polymers [36]. β relaxation originates from localized fluctuations of parts of the main chain or hindered rotations of polymer side groups. For most polymers, β relaxation has been found below the glass transition temperature (Tg). The β relaxation mechanism is a thermally activated process characterized by an Arrhenius temperature dependence of the relaxation rates [37, 38].(4)fβ=f∞,βexp−EakT.Ea represents the potential barrier (activation energy) between two possible equilibrium states. The activation energy Ea, related to the slope of logfβ versus 1/T, depends on both the internal rotation barriers and the environment of a moving unit. The dielectric strength (Δεβ) of β relaxation normally increases with temperature. In the frequency domain, dielectric β relaxation displays a broad and symmetrically shaped loss peak with a half-width that generally narrows with increasing temperature.The interfacial polarization effects are often called MWS effects [39, 40]. The effects in heterogeneous systems were first highlighted by Maxwell and later modified successively by Wagner and Sillars. Interfacial polarization processes occur in heterogeneous dielectrics, especially at low frequencies, as a result of the build-up of space charges at interfaces between two media having different permittivity and conductivity values. The relaxation is also characterized by an Arrhenius temperature dependence of the relaxation rates, as shown in equation (4).Moreover, the DC conductivity can be considered the same transport phenomena in polymers by considering the movement of different moving units (chain segment, dipole, and charge carrier). Similar to theα relaxation time, the temperature dependence of DC conductivity can be described by the VFT equation [41, 42]:(5)σdc=σ0exp−B′T−T0′,where B′ is a constant, σ0 is the conductivity at a very high temperature, and T0′ is the temperature at which the conductivity ideally reduces to zero.The temperature dependence of DC conductivity can also be described by the Arrhenius equation [43]:(6)σdc=σ0exp−EakT.


## 3. Results and Discussion

### 3.1. NBR

Five relaxation processes,α, α′, β, γ and conductivity relaxations embedded in NBR, are identified for the first time from the imaginary dielectric spectra versus frequency from 0.01 Hz to 1 MHz for all the temperatures investigated. Figure 3 shows an example of the decomposition of the imaginary relaxation spectra, composed of two out of five relaxation processes at 275 K, by employing a dispersion analysis program. In Figure 3(a), a screenshot for the complex permittivity spectrum versus the frequency for the temperature range of 233 K∼404 K at 5 K intervals is shown. The spectrum can be decomposed by introducing two HN models with eight specific parameters, as shown in the bottom table of Figure 3(b). Several parameters involved in the optimization strategy, such as the weighting factor λ, can be adjusted in the dialog window.Figure 3
An example of the decomposition of two relaxation processes in NBR by the dispersion analysis program at 275 K. (a) Real (red line) and imaginary (blue line) fitted permittivity values versus the frequency for all temperatures and (b) the decomposition of the two relaxation processes by the dispersion analysis program. In (b), the measured real (red ×) and imaginary (blue ×) permittivity and the fitted real (red line) and imaginary (blue line) permittivity graphs at 275 K are shown, and their fitted parameters are listed at the bottom table.Theα′, α, β, γ and conductivity relaxations identified from the imaginary dielectric spectra for all temperatures investigated appear in the following temperature ranges.(i)
α′ process at temperatures ranging from 254 K to 404 K(ii)
α process at temperatures ranging from 233 K to 341 K(iii)
β process at temperatures ranging from 233 K to 269 K(iv)
γ process at temperatures ranging from 233 K to 269 K(v)
Conduction contribution at temperatures ranging from 323 K to 404 KTypical dielectric loss spectra of NBR corresponding to five relaxation processes at seven temperatures (404 K, 359 K, 326 K, 293 K, 272 K, 251 K, and 233 K) are presented in Figure4, which shows the evolution of the five relaxation peaks as the temperature increases. All the relaxation peaks shift toward higher frequencies with increasing temperature. The black solid lines through the measured data (open squares) are the sums of the relaxation processes for the α′ (black dashed line), α (red dashed line), β (blue dashed line), and γ (gray dashed line) processes fitted by the HN permittivity functions and conductivity (navy dashed line) process by the fitted conductivity function of the first term in equation (2). As an example, the loss spectra at 404 K are composed of conductivity and α′ processes. The loss spectra at other temperatures are also composed of two or three relaxation processes.Figure 4
Decomposed relaxation process of NBR for seven representative temperatures (a) 404 K (b) 359 K (c) 326 K (d) 293 K (e) 272 K (f) 251 K, and (g) 233 K. Five relaxation processes,α′ (black dashed line), α (red dashed line), β (blue dashed line), γ (gray dashed line), and conduction (navy dashed line), contribute to the imaginary permittivity data.
(a)(b)(c)(d)(e)(f)(g)Figure4 depicts the evolution of typical α′ dielectric spectra at five representative temperatures, excluding 251 K and 233 K, and it displays primary α dielectric spectra detected in a temperature range from 233 K to 341 K at five representative temperatures (326 K, 293 K, 272 K, 251 K, and 233 K). The α and α′ processes at these temperatures are well fitted by the HN equation. The α relaxation process induced by the small dipole moment is responsible for the segmental motions of the main chains from the butadiene and acrylonitrile structures caused by the rotation of the C-C bonds under the applied electric field. The α′ process is usually called normal-mode relaxation, which was observed in a Stockmayer type A polymer with a dipole moment vector parallel to the principal axis. The α′ process originates from the entire fluctuation of the end-to-end vector. The proposed α and α′ relaxation processes in the NBR molecular chain structure are shown in Figure 5.Figure 5
Proposedα and α′ relaxation processes in the NBR molecular structure.Additional relaxations called secondary relaxations, that is,β and γ processes with faster time scales, are observed at low temperatures from 233 K to 269 K. The β relaxation versus frequency displays a loss peak with a symmetrical shape, which is ascribed to the motion of the cyanide side group with a larger dipole moment. Moreover, conduction spectra are detected. Figure 4 depicts the imaginary dielectric spectra containing the α and α′ processes at three representative temperatures (404 K, 359 K, and 326 K). At these temperatures, the navy dashed lines are fitted by the conductivity contribution of the first term of equation (2). At high temperatures above 323 K and low frequencies less than a few hundred kHz, a sharp increase in the imaginary dielectric spectra with decreasing frequency is correlated with the conduction of charges, which mainly occurs in polymers due to highly conductive filler.According to the fitting of equation (2), the relaxation strength Δεk and shape parameters (a and b) for the four relaxation processes except conductivity are plotted as a function of temperature, as shown in Figure 6. In the α′ process (Figure 6(a)), Δεα′ decreases with increasing temperature. This temperature dependence is predicted by the Debye theory of dielectric relaxation, as improved by Bottcher et al. [44, 45]. The value of a lies in the range of 0.2∼0.5, and the value of b is ∼1.0 (Figure 6(b)). This finding implies that the linewidth is relatively sharp with a linewidth of 105 Hz∼106 Hz, and its shapes are symmetrical.Figure 6
Temperature dependence of the relaxation strength and shape parameter for the four relaxation processes in NBR.
(a)(b)(c)(d)(e)(f)(g)(h)For theα process (Figures 6(c) and 6(d)), Δεα decreases with increasing temperature. Similar to α′, this temperature dependence is also predicted by the Debye theory of dielectric relaxation. a and b are determined to be 0.0∼0.9 and 0.5∼1.0, respectively. This finding implies again that the linewidth is more or less broad (104 Hz∼1015 Hz) and the shape is symmetrical. Δεβ and Δεγ for secondary relaxation processes (β and γ), as shown in Figures 6(e) and 6(g), are smaller than those for the primary relaxation process, as expected. The a values (Figures 6(f) and 6(h)) for the β and γ processes are 0.25 and ∼0.5, respectively, corresponding to lines of 106 Hz and ∼1015 Hz. The b value for the β process is found to be 1, as shown in Figure 6(f), indicating that the process complies with the Cole-Cole function rather than the HN function. The b value for the γ process (Figure 6(h)) is found to be less than 0.5, and the relaxation loss spectra are asymmetrical.In the case of conductivity relaxation, the DC conductivityσdc and exponent N are determined as a function of temperature as a result of fitting the first term of equation (2). The determined N value for the temperature variations is nearly 1, implying that σdc is a constant independent of the frequency.The central frequency (f0) of the imaginary loss peak for the corresponding four relaxation processes, α′, α, β, and γ, and the conduction contribution are plotted as a function of the reciprocal temperature, as shown in Figure 7, which is the relaxation map. The f0 value for the four relaxation processes decreases with decreasing temperature, as also shown in Figure 4. Simultaneously, the fitting results (solid lines) obtained by assuming the Arrhenius temperature dependence in equation (4) for the α′ relaxations give activation energies Ea,α′ = 69.6 kJ/mol. The activation energy for the α′ process is comparable to the Ea = 83.1 kJ/mol value obtained from the thermogravimetric analysis [46]. Similar to the α′ process, the activation energies for the β and γ processes determined by assuming the Arrhenius temperature dependence are Ea,β = 32.0 kJ/mol and Ea,γ = 24.3 kJ/mol, respectively, which is the typical value of the activation energy.Figure 7
Reciprocal temperature dependence of the central frequency (f0) of the imaginary loss peak for the α′ (black circle), α (red circle), β (blue circle), γ (gray circle), and conductivity (navy circle) processes in NBR. The corresponding solid lines are the fitted results for the relaxation processes.On the other hand, the relaxation rate of theα relaxation process is normally fitted by the VFT law, as shown in equation (3), giving rise to the glass temperature, Tg = 242.4 K, as shown in Figure 7. The obtained glass transition temperature Tg is very similar to the 242.2 K measurement of the differential scanning calorimetry.As a result of fitting the first term of equation (2) for the conductivity contribution, the DC conductivity σdc is determined as a function of temperature; as shown in Figure 7σdc is found to increase with increasing temperature, which confirms the Arrhenius dependence as expected. The reciprocal temperature dependence of conductivity satisfies the Arrhenius temperature dependence of equation (6), which gives an activation energy of Ea,c = 38.7 kJ/mol.


### 3.2. EPDM

The complex dielectric permittivity for the EPDM specimen is measured in the 0.01 Hz to 1 MHz and 233 K to 404 K ranges. Similar to NBR, we also decompose the superimposed dielectric spectra by employing the dispersion analysis program. In Figure8(a), a screenshot for the complex permittivity spectrum for all temperatures with a 5 K interval is shown. Figure 8(b) shows an example of the decomposition of relaxation spectra at 272 K, which comprises three relaxation processes. The decomposition procedure is performed by the three HN functions, as shown in Figure 8(b).Figure 8
An example of the decomposition of relaxation spectra by the dispersion analysis program at 272 K in EPDM. (a) Real (red line) and imaginary (blue line) fitted permittivity line versus frequency for all temperatures and (b) the decomposition of three relaxation processes by the dispersion analysis program. In (b), the measured data of the real (red ×) and imaginary (blue ×) permittivity and the fitted real (red line) and imaginary (blue line) permittivity graphs at 272 K are shown, and their fitted parameters are listed in the bottom table.The five relaxation processes,β, α, αβ, MWS, and conductivity, are determined from imaginary dielectric spectra for all the temperatures investigated. The relaxation processes are observed following the temperature ranges.(i)
α relaxation process at temperatures ranging from 233 K to 302 K(ii)
β relaxation process at temperatures ranging from 233 K to 305 K(iii)
Theαβ merged relaxation at temperatures ranging from 305 K to 404 K originating from the β process(iv)
MWS relaxation process at temperatures ranging from 233 K to 404 K(v)
Conduction contribution at temperatures ranging from 308 K to 404 K.Typical dielectric loss spectra of EPDM at seven temperatures (404 K, 374 K, 344 K, 313 K, 281 K, 257 K, and 233 K) are displayed in Figure9, which shows the evolution of four relaxation peaks as the temperature increases. The black solid lines through the data (open squares) are the fitted sums of the MWS (orange dashed line), β (blue dashed line) processes by two HN fits, and the conductivity process (navy dashed line) by the fitted conductivity function (404 K, 374 K, 344 K, and 313 K) or the fitted sums of the MWS (orange dashed line), α (red dashed line), and β (blue dashed line) processes by three HN fits (281 K, 257 K, and 233 K). The αβ merged relaxation process from 305 K to 404 K originates from the β process. All the relaxation peaks shift toward a higher frequency with increasing temperature.Figure 9
Relaxation processes of the EPDM specimen at seven temperatures ranging from 233 K to 404 K. Four relaxation processes,α (red dashed line), β (blue dashed line), MWS (orange dashed line), and conductivity (navy dashed line), contribute to the imaginary permittivity data. (a) 404 K (b) 374 K (c) 344 K (d) 313 K (e) 281 K (f) 257 K, and (g) 233 K.
(a)(b)(c)(d)(e)(f)(g)Theα dielectric spectra of the EPDM specimen are observed below room temperature. Each dielectric loss spectrum versus frequency displays a wide loss peak with a symmetrical shape, which is attributed to the relaxation of local segmental motion of the main chain. The β relaxation spectra detected in the entire temperature range investigated display a narrower loss peak than that of the α process with a symmetrical shape. The β process is attributed to the relaxation of the rotational motion of the side groups. The α process is responsible for the rotation of the main chain (C-C bond) in the ethylene, propylene, and ENB units, as shown in Figure 10. The β relaxation process is based on the reorientation of the CH3 bond (methyl) in the side chain of the propylene structure.Figure 10
Proposedα and β relaxation processes in the EPDM molecular structure.Theα relaxation spectra are well fitted by the Cole-Cole function with b = 1 in equation (2). Figure 11(a) shows the temperature dependence of a for the α process obtained from the dispersion analysis program. The a values for the α process are independent of temperature. The value for a is found to be 0.85 with b = 1. This finding implies that relaxation spectra have a symmetrical shape with a linewidth above 1015 Hz.Figure 11
Temperature dependence of the (a), (b) shape parameter and (c) relaxation strength for theα and β processes.
(a)(b)(c)Theβ relaxation spectra detected in the entire temperature range investigated display a narrower loss peak compared to the α process with a symmetrical shape, fitted by the Cole-Cole equation in the case of b = 1 in equation (2). Figure 11(b) shows the temperature dependence of the shape parameter obtained from the dispersion simulation with equation (2). The a value for the β process remains constant for temperature variations within the measurement uncertainty. The values for a are found to be 0.50∼0.65, corresponding to a linewidth of 1012 Hz∼1014 Hz with a symmetrical shape.As a result of a simulation from the dispersion analysis program, theΔε values for both the α and β processes as a function of temperature are also determined, as shown in Figure 11(c). The relaxation strength shows Δεα>Δεβ for all the temperatures investigated. This is generally observed for most amorphous polymers having a dipole moment attached to the main chain. For the α process, Δεα decreases with increasing temperature. This temperature dependence is predicted by the Debye theory of dielectric relaxation. However, the Δεβ value for the β relaxation process increases with temperature, which is caused by an increase in the fluctuation angle of the dipole vector or by the number of moving dipoles increasing the temperature. These are also general behaviors for the β process. The peaks of the α and β processes, which are separated below 305 K (indicated by arrows), merge above this temperature, as shown in Figure 11(c). The merged α and β processes are called the αβ process. However, Figure 11(c) shows that the temperature behavior for the merged peak may follow that for the β process.MWS spectra are detected in the entire temperature range investigated. Each dielectric loss spectrum versus frequency displays a clear loss peak with a symmetrical shape above 290 K (Figure9), which may be attributed to the interfacial interaction between the filler and polymer matrix. This MWS relaxation may be due to the accumulation of charged ions at the interfaces of inhomogeneous media. This relaxation arises from the free ingredients added to EPDM, which are present at the stage of processing and are now immobilized in the specimen, which leads to conductivity, and consequently, the charges can migrate under the applied electric field. MWS relaxation is usually confounded with ionic conduction, as these two phenomena originate from charged carrier mobility.The MWS relaxation spectra are fitted with the HN function in equation (2). Figure 12(a) shows the temperature dependence of a and b obtained from the dispersion analysis program. The values for a for the MWS relaxation spectra range from 0∼0.4. The relaxation spectra have narrow linewidths in the range from 104 Hz∼106 Hz. The linewidth increases with increasing temperature. The shape parameter b is found to be 1 at all the temperatures investigated, implying that the imaginary permittivity for MWS relaxation is symmetrical with respect to the frequency. This finding again indicates that this behavior of the relaxation spectra also satisfies the Cole-Cole equation. Moreover, Figure 12(b) shows that ΔεMWS for the MWS relaxation process decreases with temperature, similar to that of the α process.Figure 12
Temperature dependence of the (a) shape parameter and (b) relaxation strength for the MWS relaxation spectra.
(a)(b)Conductivity spectra are detected above room temperature, which is shown in Figure9 at the four representative temperatures (404 K, 374 K, 344 K, and 313 K). The conductivity spectra are attributed to the conductive carbon black filler contained in EPDM. It is normally observed that the imaginary permittivity increases with decreasing frequency and with increasing temperature. Two relaxation processes corresponding to both the MWS polarization and conductivity are normally observed in the composite containing conducting filler. The conductivity (navy solid lines) through the data (open square) is the fitting result by the first term of equation (2), and the DC conductivity σdc and exponent N are determined as a function of temperature. The determined N value is nearly 0.85∼0.95, implying that σdc is almost constant and independent of the frequency.Thef0 values of the imaginary loss peak for the corresponding α, β and MWS relaxation processes are plotted as a function of the reciprocal temperature, as shown in Figure 13. The fitting results (solid line) of these relaxations are also displayed. According to equation (4), the Arrhenius temperature dependence of f0 in the relaxation rates for the β and MWS relaxation processes results in Ea,β = 16.4 kJ/mol and Ea, MWS = 68.0 kJ/mol, respectively. By assuming a non-Arrhenius temperature dependence of the α relaxation rate, the relaxation rate can be fitted by the VFT law, as shown in equation (3); this approach provides glass transition temperature values of Tg = 211.4 K (Figure 13).Figure 13
Reciprocal temperature dependence of the central frequency (f0) of the imaginary loss peak for α (red circles), β (blue circles), MSW (orange circles), and conductivity (navy circles) processes in EPDM. The corresponding solid lines are the fitted results for the relaxation processes.The navy solid lines through the data (navy open squares) are the fitting results by the first term of equation (2), and the DC conductivity σdc is determined as a function of the temperature, as shown in Figure 13σdc is found to increase with increasing temperature, which follows the Arrhenius dependence as expected. The reciprocal temperature dependence of the conductivity data fit by equation (6) gives activation energy of Ea = 10.9 kJ/mol.Theαβ peaks that merge in the temperature range from 305 K to 404 K are separated below approximately 305 K (Figure 11(c)). The intersection at 305 K in the relaxation strength for the α and β processes, as shown in Figure 11(c), is the same phenomenon as the relaxation peak indicated by the arrow in Figure 13. The merging and splitting of both the α and β relaxations are observed for the first time in our work. There are three ways to interpret the splitting of the α and β processes [8]. Among them, the merged α and β relaxations originate from the β relaxation because the merged relaxation is described by the Arrhenius dependence of the β relaxation (Figure 13), which is also found in the temperature dependence of the relaxation strength for the α and β relaxations, as shown in Figure 11(c).


## 3.1. NBR

Five relaxation processes,α, α′, β, γ and conductivity relaxations embedded in NBR, are identified for the first time from the imaginary dielectric spectra versus frequency from 0.01 Hz to 1 MHz for all the temperatures investigated. Figure 3 shows an example of the decomposition of the imaginary relaxation spectra, composed of two out of five relaxation processes at 275 K, by employing a dispersion analysis program. In Figure 3(a), a screenshot for the complex permittivity spectrum versus the frequency for the temperature range of 233 K∼404 K at 5 K intervals is shown. The spectrum can be decomposed by introducing two HN models with eight specific parameters, as shown in the bottom table of Figure 3(b). Several parameters involved in the optimization strategy, such as the weighting factor λ, can be adjusted in the dialog window.Figure 3
An example of the decomposition of two relaxation processes in NBR by the dispersion analysis program at 275 K. (a) Real (red line) and imaginary (blue line) fitted permittivity values versus the frequency for all temperatures and (b) the decomposition of the two relaxation processes by the dispersion analysis program. In (b), the measured real (red ×) and imaginary (blue ×) permittivity and the fitted real (red line) and imaginary (blue line) permittivity graphs at 275 K are shown, and their fitted parameters are listed at the bottom table.Theα′, α, β, γ and conductivity relaxations identified from the imaginary dielectric spectra for all temperatures investigated appear in the following temperature ranges.(i)
α′ process at temperatures ranging from 254 K to 404 K(ii)
α process at temperatures ranging from 233 K to 341 K(iii)
β process at temperatures ranging from 233 K to 269 K(iv)
γ process at temperatures ranging from 233 K to 269 K(v)
Conduction contribution at temperatures ranging from 323 K to 404 KTypical dielectric loss spectra of NBR corresponding to five relaxation processes at seven temperatures (404 K, 359 K, 326 K, 293 K, 272 K, 251 K, and 233 K) are presented in Figure4, which shows the evolution of the five relaxation peaks as the temperature increases. All the relaxation peaks shift toward higher frequencies with increasing temperature. The black solid lines through the measured data (open squares) are the sums of the relaxation processes for the α′ (black dashed line), α (red dashed line), β (blue dashed line), and γ (gray dashed line) processes fitted by the HN permittivity functions and conductivity (navy dashed line) process by the fitted conductivity function of the first term in equation (2). As an example, the loss spectra at 404 K are composed of conductivity and α′ processes. The loss spectra at other temperatures are also composed of two or three relaxation processes.Figure 4
Decomposed relaxation process of NBR for seven representative temperatures (a) 404 K (b) 359 K (c) 326 K (d) 293 K (e) 272 K (f) 251 K, and (g) 233 K. Five relaxation processes,α′ (black dashed line), α (red dashed line), β (blue dashed line), γ (gray dashed line), and conduction (navy dashed line), contribute to the imaginary permittivity data.
(a)(b)(c)(d)(e)(f)(g)Figure4 depicts the evolution of typical α′ dielectric spectra at five representative temperatures, excluding 251 K and 233 K, and it displays primary α dielectric spectra detected in a temperature range from 233 K to 341 K at five representative temperatures (326 K, 293 K, 272 K, 251 K, and 233 K). The α and α′ processes at these temperatures are well fitted by the HN equation. The α relaxation process induced by the small dipole moment is responsible for the segmental motions of the main chains from the butadiene and acrylonitrile structures caused by the rotation of the C-C bonds under the applied electric field. The α′ process is usually called normal-mode relaxation, which was observed in a Stockmayer type A polymer with a dipole moment vector parallel to the principal axis. The α′ process originates from the entire fluctuation of the end-to-end vector. The proposed α and α′ relaxation processes in the NBR molecular chain structure are shown in Figure 5.Figure 5
Proposedα and α′ relaxation processes in the NBR molecular structure.Additional relaxations called secondary relaxations, that is,β and γ processes with faster time scales, are observed at low temperatures from 233 K to 269 K. The β relaxation versus frequency displays a loss peak with a symmetrical shape, which is ascribed to the motion of the cyanide side group with a larger dipole moment. Moreover, conduction spectra are detected. Figure 4 depicts the imaginary dielectric spectra containing the α and α′ processes at three representative temperatures (404 K, 359 K, and 326 K). At these temperatures, the navy dashed lines are fitted by the conductivity contribution of the first term of equation (2). At high temperatures above 323 K and low frequencies less than a few hundred kHz, a sharp increase in the imaginary dielectric spectra with decreasing frequency is correlated with the conduction of charges, which mainly occurs in polymers due to highly conductive filler.According to the fitting of equation (2), the relaxation strength Δεk and shape parameters (a and b) for the four relaxation processes except conductivity are plotted as a function of temperature, as shown in Figure 6. In the α′ process (Figure 6(a)), Δεα′ decreases with increasing temperature. This temperature dependence is predicted by the Debye theory of dielectric relaxation, as improved by Bottcher et al. [44, 45]. The value of a lies in the range of 0.2∼0.5, and the value of b is ∼1.0 (Figure 6(b)). This finding implies that the linewidth is relatively sharp with a linewidth of 105 Hz∼106 Hz, and its shapes are symmetrical.Figure 6
Temperature dependence of the relaxation strength and shape parameter for the four relaxation processes in NBR.
(a)(b)(c)(d)(e)(f)(g)(h)For theα process (Figures 6(c) and 6(d)), Δεα decreases with increasing temperature. Similar to α′, this temperature dependence is also predicted by the Debye theory of dielectric relaxation. a and b are determined to be 0.0∼0.9 and 0.5∼1.0, respectively. This finding implies again that the linewidth is more or less broad (104 Hz∼1015 Hz) and the shape is symmetrical. Δεβ and Δεγ for secondary relaxation processes (β and γ), as shown in Figures 6(e) and 6(g), are smaller than those for the primary relaxation process, as expected. The a values (Figures 6(f) and 6(h)) for the β and γ processes are 0.25 and ∼0.5, respectively, corresponding to lines of 106 Hz and ∼1015 Hz. The b value for the β process is found to be 1, as shown in Figure 6(f), indicating that the process complies with the Cole-Cole function rather than the HN function. The b value for the γ process (Figure 6(h)) is found to be less than 0.5, and the relaxation loss spectra are asymmetrical.In the case of conductivity relaxation, the DC conductivityσdc and exponent N are determined as a function of temperature as a result of fitting the first term of equation (2). The determined N value for the temperature variations is nearly 1, implying that σdc is a constant independent of the frequency.The central frequency (f0) of the imaginary loss peak for the corresponding four relaxation processes, α′, α, β, and γ, and the conduction contribution are plotted as a function of the reciprocal temperature, as shown in Figure 7, which is the relaxation map. The f0 value for the four relaxation processes decreases with decreasing temperature, as also shown in Figure 4. Simultaneously, the fitting results (solid lines) obtained by assuming the Arrhenius temperature dependence in equation (4) for the α′ relaxations give activation energies Ea,α′ = 69.6 kJ/mol. The activation energy for the α′ process is comparable to the Ea = 83.1 kJ/mol value obtained from the thermogravimetric analysis [46]. Similar to the α′ process, the activation energies for the β and γ processes determined by assuming the Arrhenius temperature dependence are Ea,β = 32.0 kJ/mol and Ea,γ = 24.3 kJ/mol, respectively, which is the typical value of the activation energy.Figure 7
Reciprocal temperature dependence of the central frequency (f0) of the imaginary loss peak for the α′ (black circle), α (red circle), β (blue circle), γ (gray circle), and conductivity (navy circle) processes in NBR. The corresponding solid lines are the fitted results for the relaxation processes.On the other hand, the relaxation rate of theα relaxation process is normally fitted by the VFT law, as shown in equation (3), giving rise to the glass temperature, Tg = 242.4 K, as shown in Figure 7. The obtained glass transition temperature Tg is very similar to the 242.2 K measurement of the differential scanning calorimetry.As a result of fitting the first term of equation (2) for the conductivity contribution, the DC conductivity σdc is determined as a function of temperature; as shown in Figure 7σdc is found to increase with increasing temperature, which confirms the Arrhenius dependence as expected. The reciprocal temperature dependence of conductivity satisfies the Arrhenius temperature dependence of equation (6), which gives an activation energy of Ea,c = 38.7 kJ/mol.


## 3.2. EPDM

The complex dielectric permittivity for the EPDM specimen is measured in the 0.01 Hz to 1 MHz and 233 K to 404 K ranges. Similar to NBR, we also decompose the superimposed dielectric spectra by employing the dispersion analysis program. In Figure8(a), a screenshot for the complex permittivity spectrum for all temperatures with a 5 K interval is shown. Figure 8(b) shows an example of the decomposition of relaxation spectra at 272 K, which comprises three relaxation processes. The decomposition procedure is performed by the three HN functions, as shown in Figure 8(b).Figure 8
An example of the decomposition of relaxation spectra by the dispersion analysis program at 272 K in EPDM. (a) Real (red line) and imaginary (blue line) fitted permittivity line versus frequency for all temperatures and (b) the decomposition of three relaxation processes by the dispersion analysis program. In (b), the measured data of the real (red ×) and imaginary (blue ×) permittivity and the fitted real (red line) and imaginary (blue line) permittivity graphs at 272 K are shown, and their fitted parameters are listed in the bottom table.The five relaxation processes,β, α, αβ, MWS, and conductivity, are determined from imaginary dielectric spectra for all the temperatures investigated. The relaxation processes are observed following the temperature ranges.(i)
α relaxation process at temperatures ranging from 233 K to 302 K(ii)
β relaxation process at temperatures ranging from 233 K to 305 K(iii)
Theαβ merged relaxation at temperatures ranging from 305 K to 404 K originating from the β process(iv)
MWS relaxation process at temperatures ranging from 233 K to 404 K(v)
Conduction contribution at temperatures ranging from 308 K to 404 K.Typical dielectric loss spectra of EPDM at seven temperatures (404 K, 374 K, 344 K, 313 K, 281 K, 257 K, and 233 K) are displayed in Figure9, which shows the evolution of four relaxation peaks as the temperature increases. The black solid lines through the data (open squares) are the fitted sums of the MWS (orange dashed line), β (blue dashed line) processes by two HN fits, and the conductivity process (navy dashed line) by the fitted conductivity function (404 K, 374 K, 344 K, and 313 K) or the fitted sums of the MWS (orange dashed line), α (red dashed line), and β (blue dashed line) processes by three HN fits (281 K, 257 K, and 233 K). The αβ merged relaxation process from 305 K to 404 K originates from the β process. All the relaxation peaks shift toward a higher frequency with increasing temperature.Figure 9
Relaxation processes of the EPDM specimen at seven temperatures ranging from 233 K to 404 K. Four relaxation processes,α (red dashed line), β (blue dashed line), MWS (orange dashed line), and conductivity (navy dashed line), contribute to the imaginary permittivity data. (a) 404 K (b) 374 K (c) 344 K (d) 313 K (e) 281 K (f) 257 K, and (g) 233 K.
(a)(b)(c)(d)(e)(f)(g)Theα dielectric spectra of the EPDM specimen are observed below room temperature. Each dielectric loss spectrum versus frequency displays a wide loss peak with a symmetrical shape, which is attributed to the relaxation of local segmental motion of the main chain. The β relaxation spectra detected in the entire temperature range investigated display a narrower loss peak than that of the α process with a symmetrical shape. The β process is attributed to the relaxation of the rotational motion of the side groups. The α process is responsible for the rotation of the main chain (C-C bond) in the ethylene, propylene, and ENB units, as shown in Figure 10. The β relaxation process is based on the reorientation of the CH3 bond (methyl) in the side chain of the propylene structure.Figure 10
Proposedα and β relaxation processes in the EPDM molecular structure.Theα relaxation spectra are well fitted by the Cole-Cole function with b = 1 in equation (2). Figure 11(a) shows the temperature dependence of a for the α process obtained from the dispersion analysis program. The a values for the α process are independent of temperature. The value for a is found to be 0.85 with b = 1. This finding implies that relaxation spectra have a symmetrical shape with a linewidth above 1015 Hz.Figure 11
Temperature dependence of the (a), (b) shape parameter and (c) relaxation strength for theα and β processes.
(a)(b)(c)Theβ relaxation spectra detected in the entire temperature range investigated display a narrower loss peak compared to the α process with a symmetrical shape, fitted by the Cole-Cole equation in the case of b = 1 in equation (2). Figure 11(b) shows the temperature dependence of the shape parameter obtained from the dispersion simulation with equation (2). The a value for the β process remains constant for temperature variations within the measurement uncertainty. The values for a are found to be 0.50∼0.65, corresponding to a linewidth of 1012 Hz∼1014 Hz with a symmetrical shape.As a result of a simulation from the dispersion analysis program, theΔε values for both the α and β processes as a function of temperature are also determined, as shown in Figure 11(c). The relaxation strength shows Δεα>Δεβ for all the temperatures investigated. This is generally observed for most amorphous polymers having a dipole moment attached to the main chain. For the α process, Δεα decreases with increasing temperature. This temperature dependence is predicted by the Debye theory of dielectric relaxation. However, the Δεβ value for the β relaxation process increases with temperature, which is caused by an increase in the fluctuation angle of the dipole vector or by the number of moving dipoles increasing the temperature. These are also general behaviors for the β process. The peaks of the α and β processes, which are separated below 305 K (indicated by arrows), merge above this temperature, as shown in Figure 11(c). The merged α and β processes are called the αβ process. However, Figure 11(c) shows that the temperature behavior for the merged peak may follow that for the β process.MWS spectra are detected in the entire temperature range investigated. Each dielectric loss spectrum versus frequency displays a clear loss peak with a symmetrical shape above 290 K (Figure9), which may be attributed to the interfacial interaction between the filler and polymer matrix. This MWS relaxation may be due to the accumulation of charged ions at the interfaces of inhomogeneous media. This relaxation arises from the free ingredients added to EPDM, which are present at the stage of processing and are now immobilized in the specimen, which leads to conductivity, and consequently, the charges can migrate under the applied electric field. MWS relaxation is usually confounded with ionic conduction, as these two phenomena originate from charged carrier mobility.The MWS relaxation spectra are fitted with the HN function in equation (2). Figure 12(a) shows the temperature dependence of a and b obtained from the dispersion analysis program. The values for a for the MWS relaxation spectra range from 0∼0.4. The relaxation spectra have narrow linewidths in the range from 104 Hz∼106 Hz. The linewidth increases with increasing temperature. The shape parameter b is found to be 1 at all the temperatures investigated, implying that the imaginary permittivity for MWS relaxation is symmetrical with respect to the frequency. This finding again indicates that this behavior of the relaxation spectra also satisfies the Cole-Cole equation. Moreover, Figure 12(b) shows that ΔεMWS for the MWS relaxation process decreases with temperature, similar to that of the α process.Figure 12
Temperature dependence of the (a) shape parameter and (b) relaxation strength for the MWS relaxation spectra.
(a)(b)Conductivity spectra are detected above room temperature, which is shown in Figure9 at the four representative temperatures (404 K, 374 K, 344 K, and 313 K). The conductivity spectra are attributed to the conductive carbon black filler contained in EPDM. It is normally observed that the imaginary permittivity increases with decreasing frequency and with increasing temperature. Two relaxation processes corresponding to both the MWS polarization and conductivity are normally observed in the composite containing conducting filler. The conductivity (navy solid lines) through the data (open square) is the fitting result by the first term of equation (2), and the DC conductivity σdc and exponent N are determined as a function of temperature. The determined N value is nearly 0.85∼0.95, implying that σdc is almost constant and independent of the frequency.Thef0 values of the imaginary loss peak for the corresponding α, β and MWS relaxation processes are plotted as a function of the reciprocal temperature, as shown in Figure 13. The fitting results (solid line) of these relaxations are also displayed. According to equation (4), the Arrhenius temperature dependence of f0 in the relaxation rates for the β and MWS relaxation processes results in Ea,β = 16.4 kJ/mol and Ea, MWS = 68.0 kJ/mol, respectively. By assuming a non-Arrhenius temperature dependence of the α relaxation rate, the relaxation rate can be fitted by the VFT law, as shown in equation (3); this approach provides glass transition temperature values of Tg = 211.4 K (Figure 13).Figure 13
Reciprocal temperature dependence of the central frequency (f0) of the imaginary loss peak for α (red circles), β (blue circles), MSW (orange circles), and conductivity (navy circles) processes in EPDM. The corresponding solid lines are the fitted results for the relaxation processes.The navy solid lines through the data (navy open squares) are the fitting results by the first term of equation (2), and the DC conductivity σdc is determined as a function of the temperature, as shown in Figure 13σdc is found to increase with increasing temperature, which follows the Arrhenius dependence as expected. The reciprocal temperature dependence of the conductivity data fit by equation (6) gives activation energy of Ea = 10.9 kJ/mol.Theαβ peaks that merge in the temperature range from 305 K to 404 K are separated below approximately 305 K (Figure 11(c)). The intersection at 305 K in the relaxation strength for the α and β processes, as shown in Figure 11(c), is the same phenomenon as the relaxation peak indicated by the arrow in Figure 13. The merging and splitting of both the α and β relaxations are observed for the first time in our work. There are three ways to interpret the splitting of the α and β processes [8]. Among them, the merged α and β relaxations originate from the β relaxation because the merged relaxation is described by the Arrhenius dependence of the β relaxation (Figure 13), which is also found in the temperature dependence of the relaxation strength for the α and β relaxations, as shown in Figure 11(c).


## 4. Conclusions

Dielectric relaxation spectroscopy was substantially employed to evaluate the dielectric response of a material. Polymer and polymer composites were subjected to dielectric relaxation spectroscopy for a wide range of times and temperatures to study the relaxation-related phenomena occurring between them. This dielectric spectroscopy was applied to two rubbery polymers to elucidate the relaxation processes with a self-developed dispersion analysis program. The analysis program decomposed the superposed dielectric spectra in the frequency domain, determined the relaxation processes, and made it possible to assign their peaks, thereby resulting in the characterization of entire relaxation processes.In the NBR investigation results, five relaxation processes were identified, namely,α, α′, β, γ, and conductivity relaxations, which could be explained in relation to chain structure and filler. The dielectric strength and shape parameters for the relaxation process were investigated as functions of temperature. The temperature behavior on the conductivity followed the Arrhenius equation at temperatures above the glass transition temperature. The activation energies for the α′, β, γ, and conductivity relaxations were obtained by applying the Arrhenius equation to the position of the imaginary permittivity peak versus the temperature. The VFT fits of the observed α relaxation were used to determine the glass transition temperature, and these values were compared with those obtained in a different study. The molecular chains responsible for the relaxation processes, α′ and α, were also assigned by calculating the corresponding dipole moment through molecular modeling.In the EPDM investigation results, five relaxation processes were also identified:α, β, αβ, MWS, and conductivity relaxations. The relaxation strength and shape parameter for these relaxation processes were obtained as a function of temperature, which presents information on relaxation. The temperature behavior on the conductivity relaxation followed an Arrhenius dependence. The activation energies for both the β, MWS and conduction relaxations were obtained by assuming Arrhenius behavior. The VFT fits of the observed α relaxation were used to determine the glass transition temperature. We discovered all the relaxations in the polymers by employing our developed dispersion analysis program. The merging and splitting of the α and β relaxation peaks in the EPDM samples were observed for the first time in this work. We know that the αβ relaxation originates from the β process from the temperature dependence of the relaxation strength and relaxation peak. Therefore, dielectric relaxation spectroscopy can characterize the relaxation process and the related molecular dynamics in various polymers.



---

*Source: 8815492-2020-12-07.xml*