Influence of the glass transition on the secondary relaxation of an ...

44 downloads 0 Views 140KB Size Report
Dielectric spectra reveal the existence of two relaxation processes: the ... relaxation strength and of the low-frequency slope of the secondary relaxation.
PHYSICAL REVIEW B

VOLUME 56, NUMBER 6

1 AUGUST 1997-II

Influence of the glass transition on the secondary relaxation of an epoxy resin R. Casalini INFM and Department of Physics, University of Pisa, Piazza Torricelli 2, 56126 Pisa, Italy

D. Fioretto INFM and Department of Physics, University of Perugia, via Pascoli, 06100 Perugia, Italy

A. Livi, M. Lucchesi, and P. A. Rolla INFM and Department of Physics, University of Pisa, Piazza Torricelli 2, 56126 Pisa, Italy ~Received 1 July 1996! The dynamics of diglycidyl ether of bisphenol-A ~DGEBA! was studied by employing broadband dielectric spectroscopy over a wide temperature range extending from below to above the glass-transition temperature. Dielectric spectra reveal the existence of two relaxation processes: the structural relaxation, slowing down for decreasing temperature and freezing at T g , and a secondary process present in both liquid and glassy phase. Above the glass transition the temperature behavior of the structural relaxation time is properly described by a Vogel-Fulcher-Tamman equation with the singularity at T 0 5234 K. The analysis of the variation of the relaxation strength versus temperature of both processes clearly shows the onset of the structural relaxation at T on'350 K. Below this temperature the two relaxations progressively separate their time scales and change their shapes. The occurrence of the glass-transition phenomenon is markedly evidenced by the changes of the relaxation strength and of the low-frequency slope of the secondary relaxation. @S0163-1829~97!08830-9#

INTRODUCTION

The glass transition is no longer considered a distinguishing feature of a few systems but rather a physical phenomenon which can occur in every material provided that the melt is cooled at sufficiently high rate. In fact, as a melt is cooled, the relaxation processes progressively slow down and, if the cooling rate becomes faster than the relaxation rate, nonequilibrium structures can be eventually frozen and a glassy solid obtained. In recent years many efforts have been devoted to find an unified description of the complex behavior of supercooled liquids going through the glass transition. Much of the theoretical approaches, such as Adam and Gibbs’ free volume, or more recent percolation models1–3 and the scaling form discussed by Dixon,4 deal with the dynamics of the structural relaxation only, though it is experimentally well known that secondary relaxations and splitting of the single hightemperature process into at least two processes ~a slow cooperative process and a faster secondary one! is a very common scenario in both glass-forming simple liquids5,6 and polymeric systems.7 Such subglass relaxations are generally characterized by an Arrhenius temperature behavior of the relaxation time and have been associated with localized motions of different origin.8–11 A special interest in such phenomena was induced by the mode coupling theory12,13 ~MCT! which, since its early formulation, predicted the bimodal relaxation as a main feature of the glass transition. Moreover, a bifurcation of two relaxation processes was experimentally found in some systems close to the critical temperature, T c , of the MCT.14 More generally, recent dielectric investigations have shown that the decoupling of structural and secondary relaxations is a crucial feature of the 0163-1829/97/56~6!/3016~6!/$10.00

56

relaxation of both glass-forming simple liquids and polymeric systems for temperatures higher than the glass transition.7,15 In this paper we analyze the dielectric behavior of structural and secondary relaxation of dyglycidyl ether of bisphenol-A ~DGEBA! epoxy resin, a system which shows a secondary relaxation strong enough to be easily detected in the whole region from the decoupling of the structural relaxation down to temperatures below the glass transition, in the glassy state. Dielectric measurements were made in a wide frequency and temperature interval and the temperature behavior of the whole set of dielectric parameters, i.e., relaxation times, relaxation strengths, and shape parameters, is analyzed to have an insight into the splitting phenomenon and into the characteristics of the dynamics of the system near the glass transition.

EXPERIMENT

The epoxy resin used in this study was a commercial sample ~EPON828 by Shell Co.! of a low molecular weight liquid diglycidyl ether of bisphenol-A ~DGEBA! with an epoxy equivalent weight of about 190. The complex dielectric constant of the sample was measured over a frequency range spanning from 103 to 1010Hz, by using two different techniques. In the low-frequency interval ~ 103 – 107 Hz! the impedance analyzer HP4194A was employed. The apparatus and the measurement procedure were previously described.16 In the high-frequency interval ~ 107 – 1010 HZ! measurements were carried out via Hewlett-Packard Network Analysers ~HP8753A and HP8720C!. Dielectric spectra were collected for temperatures ranging from 163 to 363 K. 3016

© 1997 The American Physical Society

56

INFLUENCE OF THE GLASS TRANSITION ON THE . . .

3017

introduced to describe the dielectric behavior of amorphous polymers, and successively recognized also suitable to fit the relaxation behavior of low molar mass glass-forming liquids. The HN equation is a generalization of both the Cole-Cole and Cole-Davidson equations which account for the dielectric behavior of many simple liquids. Concerning the claimed universality of the KWW function, consider the work of Jonscher21 who analyzed 100 systems and found no experimental evidence of the universality of the KWW function at low frequencies, and the analysis of Havriliak and Havriliak22 for nearly 1000 compounds who came to the same conclusion using a different method. On the other hand, Alvarez et al.23 compared the relaxation time distributions of KWW and HN functions and demonstrated that each value of the stretching parameter in KWW corresponds to a specific pair of shape parameters in the HN function. This result was confirmed more recently by Havriliak and Havriliak24 using a direct fitting procedure. It has been shown that for representing the dielectric relaxation data the HN function is more general than the KWW, being the additional parameter necessary to take into account for the behavior of both of low- and high-frequency slopes of the relaxation.22 About the matter, a model which related the high- and low-frequency slopes of the relaxation with intraand intermolecular interactions ~large- and small-scale modes!, respectively, was proposed2,25 and confirmed by the experimental analysis of the behavior of low molar mass and polymeric systems.2,26 In a system where a bimodal relaxation is found, the most suitable fitting function is a superposition of two relaxation functions. In the case of our system the experimental data were conveniently fitted by a superposition of two Havriliak-Negami functions: FIG. 1. Dielectric spectra of «8 and «9 vs frequency for the epoxy resin at different temperatures. Solid lines from fitting equation. RESULTS AND DISCUSSION

Some representative dielectric spectra of the real ~e8! and imaginary ~e9! parts of the dielectric constant for DGEBA at different temperatures are shown in Figs. 1~a! and 1~b!, respectively. They confirm the existence of two relaxations: the structural relaxation ~or a! and one secondary process which, in this system, is referred as a g relaxation. In fact, in DGEBA systems the b relaxation was assigned to the molecular group including the hydroxyl group, which in EPON828 gives a negligible contribution being present one in a ten molecules only. The discussion on these assignments has been presented elsewhere.17,18 As the temperature increases, the a relaxation shifts toward higher frequencies, where it merges into the g peak whose position slightly depends on temperature. The spectra of our system show a non-Debye and asymmetric behavior of both main and secondary relaxations as it is expected for the relaxations in complex systems. In fact, both low and high molar mass supercooled systems shows a non-Debye behavior and there is no evidence of the existence of an universal law described by a common fitting function. Among the several phenomenological equations proposed to describe the observed behaviors, the Kohlrausch-Williams-Watts ~KWW!,19 in the time domain, and the Havriliak-Negami ~HN!,20 in the frequency domain, are the most widely employed. Both were

« ~ v ! 2« 2 5

~ « 0 2« 1 ! @ 11 ~ i v t 1

! 12 a 1 # b 1

1

~ « 1 2« 2 ! @ 11 ~ i v t 2 ! 12 a 2 # b 2

, ~1!

where e~v! is the complex dielectric constant; « 0 is the completely relaxed dielectric constant; « 2 is the unrelaxed dielectric constant of the g relaxation; « 1 is the unrelaxed ~relaxed! dielectric constant with respect to the a~g! relaxation; a 1 , b 1 , t 1 and a 2 , b 2 , t 2 are the shape parameters and the relaxation times of the a and g relaxation, respectively. The fitting equation was suggested by the behavior of the relaxations when their time scales were sufficiently apart, i.e., just above the glass-transition temperature; then Eq. ~1! was maintained up to the high temperature except for the two highest temperature measurement points where a single HN function was employed. The use of a simple superposition of two relaxations could not be satisfactory especially in the splitting region, where the two relaxations can couple, and a different model equation could be considered. A suitable model equation was very recently proposed by Arbe et al.27 who described the dielectric response in the splitting region of polybutadiene on the basis of the low temperature main and secondary relaxation processes and recovered the agreement between viscosity and the characteristic time of the main relaxation around the splitting region. This interesting approach leads to results appreciably different from those obtainable from Eq. ~1! only when the time scales of the two relaxations are closer than two decades. In our case this con-

CASALINI, FIORETTO, LIVI, LUCCHESI, AND ROLLA

3018

56

TABLE I. The dielectric parameters and their standard deviations calculated by fitting the experimental data at different temperatures T by Eq. ~1!. T ~K! 163 173 183 193 203 213 233 243 253 263 268 273 278 283 288 293 298 303 313 323 333 343 353

«0

11.960.5 11.660.4 11.3060.05 11.0560.01 10.7560.01 10.5460.02 10.2360.01 9.9860.01 9.7460.01 9.2460.01 8.8460.02 8.4560.01 8.0860.03

«1

«2

2.4560.01 3.4960.01 3.5260.01 3.5660.01 3.6060.01 3.6460.01 3.7260.01 3.7960.02 4.0060.02 4.6660.05 5.0760.03 5.3060.03 5.5060.04 5.7060.07 5.9760.13 6.0760.21 6.2860.29 6.3060.27 6.5560.38 6.9060.58 7.2060.45 7.4960.34 7.7160.20

2.7360.01 2.7460.02 2.7560.03 2.7760.03 2.7760.05 2.7760.08 2.7660.08 2.6860.03 2.7160.03 2.6860.03 2.6560.03 2.6560.02 2.6760.02 2.6560.02 2.6960.02 2.6560.02 2.5960.03 2.6060.03 2.5660.05 2.6060.05 2.5060.19 2.4060.01 2.4060.06

1-a 1

0.8360.01 0.8460.01 0.8760.02 0.8760.02 0.8660.02 0.8860.02 0.8860.01 0.9260.03 0.9860.04 0.9960.05 1.060.1

1-a 2 0.3060.02 0.3360.01 0.3660.01 0.460.01 0.4060.01 0.4360.01 0.4660.01 0.4960.01 0.5160.02 0.4960.01 0.5260.01 0.5360.01 0.5560.01 0.5760.02 0.5660.02 0.5760.03 0.5960.03 0.6060.02 0.6960.04 0.7660.03 0.8060.03 0.8160.01 0.8360.01

dition is approached for a few experimental points only, for the temperatures higher than 313 K; consequently the temperature behavior of the parameters would be, in any case, slightly affected by that treatment. On this basis, the fitting procedure was simultaneously performed on both «8 and «9 data by Eq. ~1!. Since the dielectric relaxation function depends nonlinearly on the parameters, we used a fitting procedure based on the Levenberg-Marquardt method.28 The values of the dielectric parameters obtained by the fitting procedure are listed in Table I. The dielectric measurements were mainly affected by the calibration procedure so that the errors are not normally distributed and an estimate of the accuracy of the dielectric parameters ~errors of the parameters! could be obtained by means of the bootstrap Monte Carlo method.28 The relaxation times, t 1 and t 2 of the a and g relaxation, respectively, have been reported in Fig. 2 as a function of the reciprocal temperature. The temperature dependence of the average structural relaxation time t 1 is non-Arrhenius, with an apparent activation energy that increases with decreasing temperature. This behavior has been described by the empirical Vogel-Fulcher-Tamman equation, that can be written in the following form:29

t 1 5 t 01 exp

S D

DT 0 , T2T 0

b1

b2

0.4860.01 0.4960.02 0.4560.02 0.4660.05 0.4760.04 0.4760.04 0.4860.05 0.4960.09 0.4760.14 0.5360.18 0.5360.25

t 1 ~s!

0.9260.15 0.8860.13 0.8260.14 0.7660.13 0.7760.15 0.7460.17 0.7560.15 0.6760.06 0.6760.07 0.7460.06 0.7260.05 0.7060.04 0.7060.04 0.6860.05 0.7060.05 0.6960.06 0.6760.07 0.7060.06 0.6060.12 0.5360.06 0.5160.07 0.5160.03 0.4860.02

t2 ~s!

~1.1060.35!31021 ~3.6260.93!31023 ~1.7860.05!31024 ~2.0760.06!31025 ~3.6760.10!31026 ~8.8460.31!31027 ~3.0360.15!31027 ~1.1060.06!31027 ~4.4560.31!31028 ~1.1360.16!31028 ~4.7360.68!31029 ~2.0960.49!31029 ~9.9965.10!310210

~3.162.1!31026 ~1.1060.57!31026 ~5.1260.90!31027 ~2.6160.90!31027 ~1.0460.40!31027 ~4.7461.90!31028 ~1.460.46!31028 ~7.363.0!31029 ~5.662.0!31029 ~3.5860.69!31029 ~3.2060.49!31029 ~2.9560.30!31029 ~2.6060.23!31029 ~2.2060.22!31029 ~1.7060.18!31029 ~1.5860.18!31029 ~1.0160.17!31029 ~8.061.4!310210 ~6.161.9!310210 ~4.260.9!310210 ~3.060.8!3 10210 ~2.1760.40!310210 ~1.4060.08!310210

the fragility parameter which controls the departure from the Arrhenius law and is useful to distinguish between strong ~high D! and fragile ~low D! glasses. The values of the parameters obtained by fitting Eq. ~2! to the relaxation time t 1 are

t 015 ~ 1.260.3! 310212 s;

D5 ~ 3.160.1! ;

T 0 5 ~ 234.260.6! K.

~2!

where t 01 is the relaxation time in the high temperature limit; T 0 is the temperature of the structural blocking. Some authors identify T 0 with the Kauzmann temperature.30,31 D is

FIG. 2. Relaxation times, t 1 and t 2 , of the a and g relaxation, respectively, vs reciprocal temperature 1/T. Solid lines from fitting equations. Error bars smaller than the symbol size are omitted.

INFLUENCE OF THE GLASS TRANSITION ON THE . . .

56

3019

The glass-transition temperature T g can be conventionally defined as the temperature at which t 1 becomes 102 s, i.e., comparable to the maximum reasonable measurement time in a relaxation experiment;32 accordingly the value T g 5257 K was calculated by Eq. ~2!. Such T g value is about 23 K greater than T 0 and is consistent with previously reported calorimetric measurements.33,34 The characteristic time of the g relaxation below and above T g were fitted by the Arrhenius equation:

t 2 ~ T ! 5 t 02 exp

S D

Ea , RT

~3!

where E a is the activation energy per mole and R is the gas constant; t 02 is the relaxation time in the high temperature limit. The following values of the parameters were obtained:

t 025 ~ 7.062.7! 310214 s;

E a 5 ~ 5.760.2! K cal/mol.

The data of secondary relaxation in Fig. 2 show a perceptible change of the slope at the glass-transition temperature but the accuracy of the experimental data is not sufficient to support any specific conclusion on this point. We limit ourselves to point out that, in some systems, a change of the activation energy of secondary relaxations at the glass transition was observed35 while, in other systems, the activation energy was found stable.8 More evident changes have been observed for the temperature behavior of « 0 , « 1 , and « 2 near the glass transition as described in the following. The temperature dependence of the relaxed dielectric constant « 0 is described by the Kirkwood-Fro¨hlich equation, which, under the assumption that orientational effects linearly superimpose, can be written36 ~ « 0 2« ` !~ 2« 0 1« ` !

« 0 ~ « ` 12 !

2

5

4p 9KT

( N i ^ m 2i & .

~4!

This equation exploits the Onsager internal electric field and relates the relaxed, « 0 , and completely unrelaxed, « ` , dielectric constants to the volume concentration of dipoles, N i , and to the mean square dipole moment of each dipolar species ^ m 2i & . Equation ~4! predicts that the quantity on the left-hand side, indicated in the following as f (« 0 ), depends linearly on the reciprocal temperature and approaches zero at very high temperatures. The experimental results plotted in Fig. 3 confirm the theoretical predictions of Eq. ~4! except for the intercept for f (« 0 )50 which occurs at T'700 K. Here the « ` value was calculated from the optical refractive index n51.5. The zeroing of the overall strength at a finite temperature was also observed in other systems and it was generally ascribed to the dipole-dipole interactions not completely taken into account by the model.37,38 The a-relaxation strength, D« 1 5(« 0 2« 1 ), behaves as the overall relaxation strength, (« 0 2« 2 )~Fig. 4!, but with a different slope. It zeros at the temperature T on'350 K, where the onset of a relaxation occurs. The existence of an onset for the a-relaxation strength has been recently reported for the poly~n-alkylmethacrylate!s7,35 and for polybutadiene,27 the linear increase of the a-relaxation strength with decreasing temperature has been also verified by simulations with a modified Fredrickson model.15

FIG. 3. Dielectric function f (« 0 ) plotted vs reciprocal temperature 1/ T. Solid line from linear fit equation.

Theg-relaxation strength, D« 2 , was found to increase with temperature both below and above the glass transition ~Fig. 4!; its slope markedly changes through the glass transition. The D« 2 data well above and well below T g can be fitted by two different linear equations which cross together at the temperature T x '238 K, which lies in between the glass transition and the Kauzmann temperature. This result clearly points out that not only the structural relaxation but also the secondary one is markedly sensitive to the glass transition. Near the glass-transition temperature, secondary processes are located in a time scale more easily accessible to spectroscopic investigations than to the structural relaxation whose characteristic time attains values too long to be experimentally detected. This means, in particular, that during the measurement the secondary relaxation strength can be considered as an equilibrium quantity, the related relaxation time being much shorter than the measurement time. Consequently, the observed change of the slope of the secondary relaxation strength could be useful to characterize the glass-transition temperature itself. The relaxation strength of the secondary process, D« 2 5(« 1 2« 2 ), equals the overall relaxation

FIG. 4. Dielectric strengths plotted vs reciprocal temperature 1/ T. Solid lines from linear fit equations.

CASALINI, FIORETTO, LIVI, LUCCHESI, AND ROLLA

3020

strength for temperatures higher than T on . As the two relaxations have merged, i.e., above the onset temperature, it could be meaningless to distinguish between a and g relaxation; in fact, in principle the relaxation is no more the same as below the merging temperature, even if it could be very close to either the main or secondary relaxation, depending on which is the bigger one. Accordingly, on the basis of our experimental data we conclude that the relaxation at temperatures above the onset temperature seems to be mainly related to the previous g relaxation. To the best of our knowledge, at the present time there are no models of the glass transition which involve secondary relaxation and their temperature behavior. Moreover, none of the numerous models that have been developed to treat relaxation phenomena is able to completely model the dynamics of the systems, i.e., positions and intensities of all of their relaxations as functions of the temperature and frequency, starting from the molecular structure.39–42 Concerning the strength of dielectric subglass relaxations for temperatures below T g , a molecular model has been developed to account for the increase of D« 2 for increasing temperatures in polymers,43 by assuming that below T g an isolated-chain model is adequate to calculate the relaxation strength. Such a model describes the subglass motional modes by two states of unequal energy; consequently, the temperature dependent population of the states leads to a relaxation strength which increases with temperature. At higher temperatures, the increase of the strength of secondary relaxations above T g was observed in many other experiments.6,35,44–46 However, we found no indications concerning how the total relaxation strength is apportioned between the structural and the secondary relaxations. It is reasonable to think that local motions can couple with diffusive motions dominating at temperatures higher than T g and continue with even larger frequencies; it is expected that the effect of the coupling can vary greatly from system to system so that, in principle, the strength of the secondary relaxation can both increase ~as in the case here discussed! and decrease. To this respect, experimental data describing how the total relaxation strength above T g is divided between the a and subglass relaxations have to be considered a critical test for detailed molecular models. Further significant information can be obtained from the analysis of the shape parameters: m 1 512a 1 , n 1 5(1 2 a 1 ) b 1 and m 2 512 a 2 , n 2 (12 a 2 ) b 2 ; they represent the low- and the high-frequency tails of their respective relaxation spectra:47,48 « 9 } v m 1,2

for v ! v 01,2 ,

« 9 } v 2n 1,2

for v @ v 01,2 ,

where v 01,2 are the angular frequencies at the maxima of «9 for the a and g relaxation, respectively. The observed temperature dependence of these parameters @Figs. 5~a! and 5~b!# means that the time-temperature scaling law does not apply to this system in the whole range analyzed; in particular m 1 and m 2 display a stronger temperature dependence than n 1 and n 2 . The failure of this scaling law proves that the simplified MCT does not apply to our system; an extended version of the MCT theory,13 though specifically developed

56

FIG. 5. Shape parameters, m and n, for ~a! the main relaxation and ~b! the secondary relaxation. Dashed lines are guides for the eyes.

for polymers, might be considered and this analysis will be probably carried out in future investigation. The monotonic increase of m 1 and n 1 with increasing temperature agrees with the predictions of the model relating the high-frequency part of the dielectric loss curve to the local-chain dynamics and the low-frequency one to the correlation of segments of different molecular chains.2 In fact, this model predicts that the increase of m 1 with the temperature is produced by the corresponding reduction of the intermolecular interactions and the hindrance of the local chains. The corresponding increase of n 1 is originated by the increase of the mobility of local chains. In the system here considered, the parameter m 2 of the g relaxation @Fig. 5~b!# features a linear dependence on the reciprocal temperature below T g and shows a marked change of the slope near the glass-transition temperature. This result demonstrates that, besides the strength, also the shape of the secondary relaxation is modified by the glass transition. Moreover, it has to be noticed that m 2 and n 2 are getting the same value as temperature decreases; this means that the secondary relaxation at low temperatures is approaching the Cole-Cole behavior as observed in many other glass-forming systems.44 For the system here considered, the coupling between motions associated to the main and secondary relaxations can give a rationale of the observed be-

56

INFLUENCE OF THE GLASS TRANSITION ON THE . . .

havior. However, further investigations on different systems are needed to clarify this point. CONCLUSIONS

In this paper broadband dielectric spectroscopy has been exploited for studying the dynamics of diglycidyl ether of bisphenol-A ~DGEBA!. Two relaxation processes have been evidenced: the structural relaxation freezing at T g , and a secondary process present in both liquid and glassy phase. The fitting procedure of the experimental data was based on the superposition of two Havriliak–Negami functions and provided all the parameters describing the relaxations. The full analysis of the temperature behavior of relaxation times,

G. Adams and J. H. Gibbs, J. Chem. Phys. 43, 139 ~1965!. A. Scho¨nhals and E. Schlosser, Colloid Polym. Sci. 267, 125 ~1989!. 3 J. A. Campbell, J. M. Flesselles, R. Jullien, and R. Botet, Phys. Rev. B 37, 3825 ~1988!. 4 P. K. Dixon, L. Wu, S. R. Nagel, B. D. Williams, and J. P. Carini, Phys. Rev. Lett. 65, 1108 ~1990!. 5 G. P. Johari and M. Goldstein, J. Chem. Phys. 53, 2372 ~1970!. 6 G. P. Johari and M. Goldstein, J. Chem. Phys. 55, 4245 ~1971!. 7 R. Garwe, A. Scho¨nhals, M. Beiner, K. Schro¨ter, and E. Donth, J. Phys., Condens. Matter. 6, 6941 ~1994!. 8 G. P. Johari, J. Chem. Phys. 89, 2074 ~1992!. 9 P. Lunkeheimer, A. Pimenov, M. Dressel, Y. G. Goncharov, R. Bohemer, and A. Loidl, Phys. Rev. Lett. 77, 318 ~1996!. 10 D. E. Buerger and R. H. Boyd, Macromolecules 22, 2694 ~1989!. 11 F. Garwe, A. Scho¨nhals, H. Lockwenz, M. Beiner, K. Schro¨ter, and E. Donth, Macromolecules 29, 247 ~1996!. 12 W. Go¨tze and L. Sjo¨gren, Rep. Prog. Phys. 55, 241 ~1992!. 13 K. S. Schweizer, J. Chem. Phys. 91, 5802 ~1989!. 14 E. Ro¨ssler, Phys. Rev. Lett. 65, 1595 ~1990!. 15 M. Schulz and E. Donth, J. Non-Cryst. Solids 168, 186 ~1994!. 16 D. Fioretto, A. Livi, P. A. Rolla, G. Socino, and L. Verdini, J. Phys. Condens. Matter. 6, 5295 ~1994!. 17 E. Butta, A. Livi, G. Levita, and P. A. Rolla, J. Polym. Sci. Part B: Polym. Phys. 33, 2253 ~1995!. 18 J. M. Pochan, R. J. Gruber, and D. F. Pochan, J. Polym. Sci. Polym. Phys. Ed. 19, 143 ~1981!. 19 G. Williams and D. C. Watts, Trans. Faraday Soc. 66, 80 ~1970!. 20 S. Havriliak, Jr. and S. Negami, J. Polym. Sci. C 14, 99 ~1966!. 21 A. K. Jonscher, Dielectric Relaxation in Solids ~Chelsea Dielectric Press, London, 1983!, p. 199. 22 S. Havriliak, Jr. and S. J. Havriliak. J. Non-Cryst. Solids 172-174, 297 ~1994!. 23 F. Alvarez, A. Alegria, and J. Colmenero, Phys. Rev. B 44, 7306 ~1991!. 24 S. Havriliak, Jr. and S. J. Havriliak, Polymer 37, 4107 ~1996!. 25 L. A. Dissado and R. M. Hill, Proc. R. Soc. London, Ser. A 390, 131 ~1983!. 26 A. Scho¨nhals, F. Kremer, and E. Schlosser, Phys. Rev. Lett. 67, 999 ~1991!. 1 2

3021

relaxation strengths, and shape parameters of both processes was performed. In particular, the onset of the main relaxation occurs well above the glass transition ~ T on5350 K!; while, as the glass-transition temperature was crossed, a change of the shape parameters and the relaxation strength of the secondary relaxation was observed. The temperature behavior of all parameters was exploited to discuss a possible scenario describing the temperature evolution of the main and secondary relaxation and their relationships. The correlation of the fitting parameters with the physical behavior of the system is satisfactory; the analysis we adopted to describe the behavior of this glass forming system provided reliable results and a deeper insight in the secondary relaxation behavior near the glass transition.

27

A. Arbe, D. Richter, J. Colmenero, and B. Farago, Phys. Rev. B 54, 3853 ~1996!. 28 W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes ~Cambridge University Press, Cambridge, England, 1992!. 29 C. A. Angell, J. Non-Cryst. Solids 131–133, 13 ~1991!. 30 C. Angell and K. J. Rao, J. Chem. Phys. 57, 470 ~1972!. 31 G. Floudas, J. S. Higgins, F. Kremer, and E. W. Fischer, Macromolecules 25, 4955 ~1992!. 32 J. Ja¨ckle, Rep. Prog. Phys. 49, 171 ~1986!. 33 N. F. Sheppard, Jr. and S. D. Senturia, J. Polym. Sci., Part B: Polym. Phys. 27, 753 ~1989!. 34 P. Hedvig, in Applied Polymer Analysis and Characterization, edited by J. Mitchell, Jr. ~Hanser, Munich, 1992!, Vol. II, p. 191. 35 F. Garwe, A. Scho¨nhals, H. Lockwenz, M. Beiner, K. Schro¨ter, and E. Donth, Macromolecules 29, 247 ~1996!. 36 N. E. Hill, in Dielectric Properties and Molecular Behavior, edited by T. M. Sugden ~Van Nostrand, London, 1969!, Chap. 1, p. 41. 37 G. Williams, Trans. Faraday Soc. 59, 1397 ~1962!. 38 N. F. Sheppard, Jr. and S. D. Senturia, Polym. Eng. Sci. 26, 354 ~1986!. 39 M. H. Cohen and G. S. Grest, Ann. ~N.Y.! Acad. Sci. 371, 199 ~1981!; R. Simha, in Molecular Basis of Transitions and Relaxations, edited by D. J. Meier ~Gordon and Breach Science, London, 1978!, pp. 203–223. 40 K. L. Ngai, Comments Solid State Phys. 9, 127 ~1979!. 41 K. L. Ngai, Comments Solid State Phys. 9, 141 ~1980!. 42 A. K. Rajagopal, K. J. Ngai, R. W. Rendell, and S. Teitler, Physica A 149, 358 ~1988!. 43 G. D. Smith and R. H. Boyd, Macromolecules 24, 2731 ~1991!. 44 N. G. MacCrum, B. E. Read, and G. Williams, Anelastic and Dielectric Effects in Polymeric Solids ~Wiley, New York, 1967!. 45 H. Sasabe and S. Saito, J. Polym. Sci., Part A-2 6, 1401 ~1968!. 46 J. L. Gomez Ribelles and R. Diaz Calleja, J. Polym. Sci. Polym. Phys. Ed. 23, 1297 ~1985!. 47 R. Casalini, A. Livi, P. A. Rolla, G. Levita, and D. Fioretto, Phys. Rev. B 53, 564 ~1996!. 48 A. K. Jonscher, Universal Relaxation Law ~Chelsea Dielectric Press, London, 1996!.