1. Introduction
The spin probe technique has been widely used for a long time to characterize the segmental mobility of polymers [
1]. It is currently a routine method for polymer studies, and the novel approaches based on pulsed and high-field EPR as well as spatial EPR imaging are constantly being proposed [
2,
3,
4]. The spin probe technique is based on the introduction of stable paramagnetic molecules with strongly anisotropic spin-Hamiltonian parameters, more often, nitroxide radicals, into the polymer matrix. The molecular mobility of probes inside the polymer, which can be determined by analyzing the shape of their EPR spectra, reflects the mobility of macromolecules. Of course, the information obtained in this way depends on the size and shape of the probe. Radicals with different geometry can exhibit various mobility in the same polymer under the same conditions. On the one hand, this makes the results obtained unambiguous; on the other hand, it is possible to characterize various dynamic processes in the polymer matrix using several probes of different structures (see, for example, [
5]).
A common parameter used to characterize the molecular mobility of spin probes is the rotational correlation time (τ
c) which is related to the rotational diffusion coefficient as follows: τ
c=1/(6·D
rot). There are two methods to extract information about the rotational mobility of radicals from their EPR spectra. The first and the most common approach is to measure spectral features, such as intensities and widths of the spectral lines and to derive correlation. A general theory of EPR linewidth for paramagnetic species has been developed in the 1960s by Kivelson [
6]. According to the theory, the peak-to-peak Lorentzian linewidth of an individual hyperfine line of the nitroxide is given by:
where m
N is the magnetic quantum number equal to -1, 0, and 1 for
14N nucleus. The A, B, C coefficients depend on the rotational correlation times. Substituting the three values of m
N into Eq. (1), and replacing
by
for the measurement convenience (it is valid for other m
N values,
I1, I0, I-1 are intensities of the corresponding lines), one can derive several expressions for the rotational correlation time [
7,
8]:
where the coefficients κ
1, κ
2, κ
3, κ
4 are determined from the principal values of g- and A-tensors of the nitroxide (see
Appendix A). The expressions (2)-(5) have been derived for isotropic rotation of nitroxide in the fast motion regime (5‧10
-11 < τ
c < 10
-9). Another assumption of exp. (2)-(5) is that the individual line shape is Lorentzian. However, it was shown in [Santiago 2008] that these formulas are applicable in the case when the inhomogeneous Gaussian broadening contributes to the spectrum due to unresolved proton hyperfine structure.
Another approach to the rotational correlation times determination is the simulation of EPR spectra. The most widely used software are EasySpin [
9] and the program package developed by Freed et al. [
10]. They are based on the stochastic Liouville equation. The EPR spectra simulation have been broadly applied for various spin-probed or spin-labelled systems due to high reliability of this approach. It should be noted that in the case of polymer matrices there is a distribution of probe molecules by rotational mobility. Indeed, any polymer is a heterogeneous medium at the molecular level; it includes regions of tightly coiled and loosely arranged macromolecular fragments. The probes being localized in these regions would differ in their rotational mobility. When a polymer includes several regions that differ significantly in the packing density, more mobile and less mobile radicals contribute to the EPR spectrum as a set of characteristic lines (see, for example, [
3]). Such a spectrum can be simulated as the sum of several signals related to the probes with different mobility. However, if the polymer is heterogeneous at the nanometer level, spin probes would be distributed in regions with gradually changing mobility, and the EPR line shape would differ significantly from the line shape for a homogeneous sample. This is the most complicated case, when it is necessary to simulate the spectra considering the continuous distribution of the probes by their rotational correlation times. This approach was implemented in [
11,
12,
13] for simulation of CW EPR spectra and in [
14,
15] for describing 2D ELDOR spectra with the use of lognormal radical distribution by rotational correlation times. In [
11,
13] it was shown that the width of distribution decreases with increasing temperature that means more uniform radical mobility at higher temperatures.
Unfortunately, the simulation of EPR spectra using distribution of radicals by rotational mobility is a complicated procedure. For this reason, spectra simulation is rarely used to analyze the spin probes mobility in polymers. Most researchers determine rotational correlation times from EPR spectra using the approximated formulas (2)-(5). It should be noted that the formulas were proposed for analysis of isotropic rotation and, generally speaking, were not intended for studying heterogeneous systems. However, since it is desirable to have a simple way for assessment of polymer mobility without complicated spectral simulation, the applicability of the formulas for calculation of “effective” τc* values should be considered. In paper we showed several examples of satisfactory agreement between the rotational correlation time values obtained by spectra simulation and calculated by the formulas for TEMPO in polymer matrices.
It seems appropriate to establish the range of applicability of the approximated approach by comparing τc* values obtained using different formulas among themselves and with those obtained by spectra simulation. That was the goal of this work. We chose a small compact nitroxide radical TEMPO as a spin probe and isotactic polypropylene (PP) as a matrix because of absence of functional groups in this polymer that could interact specifically with the probe molecules. The EPR spectra of the probe inside the polymer were recorded in a wide temperature range and simulated taking into account the lognormal distribution of the probe molecules by rotational mobility. The rotational correlation times obtained as a result of EPR spectra simulation were compared with the values calculated by formulas (2-5). It was shown that different formulas return different results in the case of strongly anisotropic rotation of radicals. The best approximation in a wide range 2‧10-11<τc < 5‧10-9 is given by eq. (5).
2. Materials and Methods
2.1. Materials
Isotactic polypropylene was kindly provided to us by Prof. P.M. Nedorezova (N.N. Semenov Federal Research Center for Chemical Physics, Russia). The polymer characteristics (stereoregularity, melting point, and degree of crystallinity) are given in [
16].
The stable nitroxide radical (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) was purchased from Sigma-Aldrich (
Figure 1a). The optimized radical geometry (
Figure 1b) was calculated using the ORCA software package [
17]. The B3LYP – 6-31g(d,p) model was used for the calculation.
2.2. Sample Preparation
TEMPO was introduced into powdered PP through diffusion from the gas phase. The number of spins as well as their distribution in PP were monitored by EPR spectroscopy. The spectra at 295 K were recorded periodically for 1-3 weeks until the spectrum shape and the line width remain constant. The sample contained 2·1014 spins per 1 mg of PP. It should be noted that the probes were localized only in the amorphous phase of PP.
2.3. EPR Spectroscopy
EPR spectra were recorded with BrukerEMXplus spectrometer equipped with a high- sensitive resonator ER 4119 HS and the temperature control unit (temperature setting accuracy ±1 K). The EPR measurements were performed at microwave power values of 0.5-1mW and magnetic field modulation amplitude of 0.05-0.10 mT.
2.4. EPR Spectra Simulation
The spin-Hamiltonian parameters of TEMPO in PP were found by simulation of the EPR spectrum recorded at a temperature of 100K in the absence of paramagnetic molecules mobility (in rigid limit). The ODF3 program available at
https://www.chem.msu.ru/rus/lab/chemkin/ODF3/ and described in detail in [
18] was used for simulation. The simulation result is shown in
Figure 2. The magnetic resonance parameters of the spin probe in the polymer were found to be g
xx=2.0096, g
yy=2.0061, g
zz=2.0022, A
xx=0.67 mT, A
yy=0.73 mT, A
zz=3.34 mT.
The EPR spectra recorded in the temperature range 165-435 K were simulated taking into consideration the distribution of radicals by their rotational mobility [
13]:
Here ρ is the probability density function, Drot is the rotational diffusion coefficient, Drot0 and σ are the center and dispersion of the distribution, respectively. The software developed by Prof. A. Kh. Vorobiev and Dr. D.A. Chernova (Moscow State University, Russia) was used. The principal axes of the g- and hyperfine interaction tensors of unpaired electron with 14N nuclei were supposed to coincide. The shape of the individual resonance line was expressed by Voigt function. The contributions of Gauss and Lorentz functions were described by tensors of the second rank, the principal axes of the tensors coincided with the magnetic axes of the radical.
Generally speaking, the simulation of EPR spectra is an inverse ill-posed mathematical problem. The main difficulty in solving such a problem lies in the correlation between the varied parameters. Simulating the spectra of mobile radicals, one should take into account that
τc significantly correlates with the contribution of Lorentz function to the spectral line width [
19]. According to our experience, in the case of distribution of the radicals by rotational mobility, a correlation between the width of the distribution and rotational correlation times may take place. To avoid the correlations, we simulated the spectra by varying the least number of parameters that yields an adequate simulation result, namely, the description of positions, widths and intensities of the spectral lines. A perfect fitting of the spectra was not considered necessary. The optimal parameters should be physically meaningful, as well as independent on the initial approximation.
3. Results
Figure 3 shows the examples of simulation of EPR spectra recorded at different temperatures. More spectra are given in the
Supplementary Materials (
Figure S1). It is seen that the simulation quality is quite good. In the high-temperature range the spectra can be simulated quite perfectly. With temperature decrease the simulation quality slightly worsens, but the positions, intensities, and shapes of the spectral lines are reproduced satisfactorily. It should be noted that we did not strive for a perfect simulation, although that was possible through increasing the number of variable parameters, but such increasing could lead to uniqueness decrease.
In
Table 1 one can see the parameters obtained as a result of EPR spectra simulation. Parameters that were not being varied during the simulation are marked with asterisks. More detailed information, including the contribution of Gauss and Lorentz functions to the Voigt function describing the shape of an individual spectral line is given in
Table S1 (
Supplementary Materials). It should be noted that we simulated EPR spectra over a wide temperature range using almost the same distribution of TEMPO molecules by rotational mobility. It can be assumed that the mobility distribution of the probes reflects their distribution by localization in the polymer matrix. Another important result is that the rotational mobility of TEMPO around the molecular X-axis (bound N-O) (see
Figure 1b) was found to be strongly hindered. Within the temperature range of 395-435 K the rotational correlation time
τx exhibited an increase with decreasing of temperature. At temperatures below 395 K this value turned out to be outside the range of X-band EPR spectra sensitivity, so, we fixed it as 1.7·10
−7 s. The phenomenon of hindered rotational mobility around X-axis for piperidine-type nitroxide radical TEMPOL in polymer matrices was previously observed in [
20,
21]. The authors explained this observation by a weak complexation between the probe molecule and polymer fragments. However, in our work, the polymer does not contain any fragments other than fully saturated -CH
3 groups, thus, the polymer contains only sp3-hybridized carbon atoms that are incapable of any specific interaction. We suppose that the strongly hindered rotational mobility of TEMPO around X-axis can be explained by its molecular geometry. Indeed, in a confined space, four CH
3-substituents would inhibit the rotation of the molecule around X-axis in the higher extent comparing to the rotation around Y-axis. The results of spectra simulation confirm this assumption – the fastest rotation of TEMPO in PP occurs around the Y-axis.
It should also be noted that in the temperature range of 245-345 K the quasi-librations of the spin probes [
13] were taken into consideration, that improved the quality of spectra simulation. At temperatures above 345 K, the librations cannot be observed due to the intense rotational motion of the radicals. At temperatures below 243 K, the libration amplitudes become small and do not manifest themselves in the EPR spectra. The separation of rotational and librational mobility is an ill-posed problem, since these motions both lead to averaging of the anisotropic spin-Hamiltonian parameters of the radicals. For this reason, we use the isotropic librational motion. As it was previously estimated [
13], anisotropy of rotation and anisotropy of librations cannot be extracted simultaneously from EPR spectra.
4. Discussion
Simulation of the EPR spectra revealed highly anisotropic rotational mobility of TEMPO spin probes in PP. To compare the rotational correlation times obtained by simulation of the spectra and those calculated by approximated formulas, we calculated the average rotational correlation time of the radicals at each temperature as
<τc> = (τx + τy + τz)/3.
Figure 4 shows the temperature dependences of this value as well as the parameters τ
c*(1), τ
c*(2), τ
c*(3), and τ
c*(4), calculated by formulas (2), (3), (4) and (5), respectively. The spectra, representative for different temperature intervals are shown on the insets. One can see that the shapes of all temperature dependences are similar. In all cases, the crossover between two modes of the polymer motion can be seen in the temperature range 305-315 K. These temperatures are higher than
Tg for isotactic PP (~ 263 K). Similar patterns were observed previously in [
22] for cis-1,4-poly(isoprene) using the spin probe method with TEMPO radical. The effective activation energies for TEMPO rotational diffusion (E
rota), obtained from the
<τc> temperature dependence, were found to be 43.5 kJ/Mol and 7.3 kJ/Mol for high and low temperature regions, respectively. The second value is close to E
rota for TEMPO in various polymers [
23]. It should be noted, that
<τc> changes monotonically with temperature, while τ
c(1), τ
c(2), τ
c(3), and τ
c(4) demonstrate non-monotonic behavior below T=200 K (
<τc> >1‧10
-8 s). We can conclude that the approximated formulas are inapplicable at such a slow rotational mobility of the radical.
Figure 5 shows differences between the <τ
c> values determined via spectral simulation and τ
c*(1), τ
c*(2), τ
c*(3), and τ
c*(4). One can see that formula (2) and (4) return values the closest to <τ
c> at τ
c < 7‧10
-10 s, while formula (5) gives the best approximation to the simulation result at longer τ
c, up to 5‧10
-9 s. One can see the shapes of EPR spectra within the range of 7‧10
-10 < τ
c < 5‧10
-9 are far from the fast-motional spectra, nevertheless, formula (5) allows obtaining of reasonable results. Formula (3) was shown to give the results which significantly decline from <τ
c> almost in all mobility range, therefore, this formula is not suitable for analysis of EPR spectra of the radicals with strongly anisotropic rotation in the polymer matrix. It should be reminded that all the approximated formulas were proposed for analysis of EPR spectra of isotropically rotating radicals.
The results obtained showed that, in the majority of cases, when the researcher is not interested in fine details of polymer organization, but is looking for a way to determine the change in the mobility of macromolecules under an impact on the system, the reliable data can be obtained using the spin probe method without simulation of EPR spectra, which is complicated, time-consuming, and requires special qualification of the researcher. It should be underlined that in the present work we investigated the behavior of the spin probe which did not contain any functional groups except for the paramagnetic N-O fragment. The studied polymer also did not contain groups capable to complexation with the probe molecules. In cases of polymers and/or probes bearing functional groups, the possibility of using the approximated formulas should be tested additionally.
5. Conclusions
The applicability of different approximated formulas for characterization of the mobility of spin probe in polymer matrix was analyzed on the example of TEMPO spin probe and isotactic polypropylene. The calculated values of rotational correlation time were compared with those obtained by simulation of EPR spectra taking into account the rotational anisotropy and lognormal distribution of radicals by rotational mobility. It was shown that in the case of strongly anisotropic rotation of probe molecules in polymer matrix, different formulas give different results. In general, the approximated method can be used for analysis of EPR spectra up to rotational correlation times τc < 5‧10-9 s.
Supplementary Materials
The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: EPR spectra of TEMPO in PP (black lines) and the results of their simulation (red lines); Table S1: Rotational characteristics of TEMPO in PP.
Author Contributions
Conceptualization and methodology, N.Ch.; formal analysis, N.Ch.; investigation, T.Ya.; resources, A.K.; writing—original draft preparation, N.Ch.; writing—review and editing, T.Ya. and A.K. All authors have read and agreed to the published version of the manuscript.
Funding
This work was funded within the framework of the Program of Fundamental Research of the Russian Federation (Reg. № 122040500068-0). The authors acknowledge M.V. Lomonosov Moscow State University Program of Development for EMXplus-10/12 PX X-band EPR Bruker spectrometer.
Acknowledgments
The authors acknowledge Prof. A.L. Iordansky for fruitful discussion.
Conflicts of Interest
The authors declare no conflicts of interest. .
Appendix A
In the case of nitroxide probes in the fast motion regime, as described in the main text, the equations derived by Kivelson [Kivelson, D. J. Chem. Phys. 1960, 33, 1094] can be combined to obtain the rotational correlation time (see Eqs. (2)-(5)).
The rotation correlation times
where h is the Planck constant, μB is the Bohr magneton, and g0 is the isotropic g-factor:
The magnetic parameters of the spin probe in the polymer were found to be gxx=2.0096, gyy=2.0061, gzz=2.0022, Axx=18.8 mHz, Ayy=20.5 mHz, Azz=93.78 mHz and central field B0=0.3356 T. Thus, κ1= κ3 = 5.54·10-10 G/s, κ2=6.34·10-10 G/s, κ4=5.17·10-10 G/s.
References
- Rånby, B.; Rabek, J.F. “ESR Spectroscopy in Polymer Research” Berlin Heidelberg New York,Springer-Verlag; Springer-Verlag Berlin Heidelberg: New York, 1977; Vol. 1, ISBN 9783642666025. [Google Scholar]
- Naveed, K.-R.; Wang, L.; Yu, H.; Ullah, R.S.; Haroon, M.; Fahad, S.; Li, J.; Elshaarani, T.; Khan, R.U.; Nazir, A. Recent Progress in the Electron Paramagnetic Resonance Study of Polymers. Polym. Chem. 2018, 9, 3306–3335. [Google Scholar] [CrossRef]
- Veksli, Z. ESR Spectroscopy for the Study of Polymer Heterogeneity. Prog. Polym. Sci. 2000, 25, 949–986. [Google Scholar] [CrossRef]
- Schlick, S. Advanced ESR Methods in Polymer Research; SCHLICK, S., Ed.; John Wiley & Sons, Inc. Publication:: Hoboken, NJ, USA, 2006; ISBN 0471731897. [Google Scholar]
- Švajdlenková, H.; Račko, D.; Bartoš, J. Spin Probe Reorientation and Its Connections with Free Volume and Relaxation Dynamics: Diglycidyl-Ether of Bis-Phenol A. J. Non. Cryst. Solids 2008, 354, 1855–1861. [Google Scholar] [CrossRef]
- Kivelson, D. Theory of ESR Linewidths of Free Radicals. J. Chem. Phys. 1960, 33, 1094–1106. [Google Scholar] [CrossRef]
- Schreier, S.; Polnaszek, C.F.; Smith, I.C.P. Spin Labels in Membranes Problems in Practice. Biochim. Biophys. Acta - Rev. Biomembr. 1978, 515, 395–436. [Google Scholar] [CrossRef] [PubMed]
- Santangelo, M.G.; Levantino, M.; Cupane, A.; Jeschke, G. Solvation of a Probe Molecule by Fluid Supercooled Water in a Hydrogel at 200 K. J. Phys. Chem. B 2008, 112, 15546–15553. [Google Scholar] [CrossRef] [PubMed]
- Stoll, S.; Schweiger, A. EasySpin, a Comprehensive Software Package for Spectral Simulation and Analysis in EPR. J. Magn. Reson. 2006, 178, 42–55. [Google Scholar] [CrossRef] [PubMed]
- Budil, D.E.; Lee, S.; Saxena, S.; Freed, J.H. Nonlinear-Least-Squares Analysis of Slow-Motion EPR Spectra in One and Two Dimensions Using a Modified Levenberg–Marquardt Algorithm. J. Magn. Reson. Ser. A 1996, 120, 155–189. [Google Scholar] [CrossRef]
- Faetti, M.; Giordano, M.; Leporini, D.; Pardi, L. Scaling Analysis and Distribution of the Rotational Correlation Times of a Tracer in Rubbery and Glassy Poly(Vinyl Acetate): An Electron Spin Resonance Investigation. Macromolecules 1999, 32, 1876–1882. [Google Scholar] [CrossRef]
- Bercu, V.; Martinelli, M.; Massa, C.A.; Pardi, L.A.; Leporini, D. Signatures of the Fast Dynamics in Glassy Polystyrene: First Evidence by High-Field Electron Paramagnetic Resonance of Molecular Guests. J. Chem. Phys. 2005, 123. [Google Scholar] [CrossRef] [PubMed]
- Chernova, D.A.; Vorobiev, A.K. Molecular Mobility of Nitroxide Spin Probes in Glassy Polymers: Models of the Complex Motion of Spin Probes. J. Appl. Polym. Sci. 2011, 121, 102–110. [Google Scholar] [CrossRef]
- Dubinskii, A.A.; Maresch, G.G.; Spiess, H.-W. Two-Dimensional Electron Paramagnetic Resonance Spectroscopy of Nitroxides: Elucidation of Restricted Molecular Motions in Glassy Solids. J. Chem. Phys. 1994, 100, 2437–2448. [Google Scholar] [CrossRef]
- Saalmueller, J.W.; Long, H.W.; Volkmer, T.; Wiesner, U.; Maresch, G.G.; Spiess, H.W. Characterization of the Motion of Spin Probes and Spin Labels in Amorphous Polymers with Two-Dimensional Field-Step ELDOR. J. Polym. Sci. Part B Polym. Phys. 1996, 34, 1093–1104. [Google Scholar] [CrossRef]
- Yankova, T.S.; Chumakova, N.A.; Nedorezova, P.M.; Palaznik, O.M.; Kokorin, A.I. Behavior of Spin Probe TEMPO in Composites Based on Polypropylene with Different Content of Single-Wall Carbon Nanotubes. Polym. Sci. Ser. A 2024. [Google Scholar] [CrossRef]
- Neese, F. The ORCA Program System. WIREs Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
- Vorobiev, A.K.; Chumakova, N.A. Simulation of Rigid-Limit and Slow-Motion EPR Spectra for Extraction of Quantitative Dynamic and Orientational Information. 2012. In Nitroxides–Theory, Experiment and Applications; InTech: Rijeka, Croatia, 2012; pp. 57–112. [Google Scholar]
- Vorobiev, A.K.; Bogdanov, A. V.; Yankova, T.S.; Chumakova, N.A. Spin Probe Determination of Molecular Orientation Distribution and Rotational Mobility in Liquid Crystals: Model-Free Approach. J. Phys. Chem. B 2019, 123, 5875–5891. [Google Scholar] [CrossRef] [PubMed]
- Chernova, D.A.; Vorobiev, A.K. Molecular Mobility of Nitroxide Spin Probes in Glassy Polymers. Quasi-libration Model. J. Polym. Sci. Part B Polym. Phys. 2009, 47, 107–120. [Google Scholar] [CrossRef]
- Kh. Vorobiev, A.; Gurman, V.S.; Klimenko, T.A. Rotational Mobility of Guest Molecules Studied by Method of Oriented Spin Probe. Phys. Chem. Chem. Phys. 2000, 2, 379–385. [Google Scholar] [CrossRef]
- Švajdlenková, H.; Šauša, O.; Adichtchev, S. V.; Surovtsev, N. V.; Novikov, V.N.; Bartoš, J. On the Mutual Relationships between Molecular Probe Mobility and Free Volume and Polymer Dynamics in Organic Glass Formers: Cis-1,4-Poly(Isoprene). Polymers (Basel). 2021, 13, 294. [Google Scholar] [CrossRef] [PubMed]
- Vasserman, A.M.; Barashkova, I.I.; Yasina, L.L.; Pupov, V.S. Rotary and Progressive Diffusion of the Nitroxyl Radical in Amorphous Polymers. Polym. Sci. U.S.S.R. 1977, 19, 2389–2398. [Google Scholar] [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).