Introduction
The dielectric properties of ice, though extensively investigated in recent years, are still worth a thorough examination in order to be completely determined.
Apart from some discrepancies among the results of certain authors, further investigations are justified by the acknowledgement that even in undoped ice other mechanisms are present (Von Hippel and others, 1972; Ruepp, 1973), beyond the already known Debye spectrum, and that other quantities such as for example the dielectric strength, must be investigated in more detail. Furthermore the temperature range in which the ice spectrum is known is not wide and its extension to lower temperatures is therefore desirable.
The experimental technique of thermally stimulated depolarization has been proposed by some authors (Bishop and Glen, 1969; Jeneveau and others, 1972; Johari and Jones, 1975) as a complement to dielectric constant measurements in the lower temperature range.
We wished to investigate the behaviour of samples of monocrystalline “pure” ice Ih by means of both these techniques to verify the coincidence of the results, and, if this is the case, to extend towards lower temperatures the range in which the dielectric relaxation of ice is known.
I. Admittance Measurements
Experimental
For both admittance and depolarization measurements the electrodes were circular parallel plates, the upper electrode having a diameter of 18 mm with a guard-ring with inner diameter 20 mm and outer diameter 35 mm mounted on a teflon disk; the samples were monocrystalline disks (24 mm in diameter, 2 to 5 mm in thickness) with the c-axis normal or parallel to the electrodes, which were pressed against the sample surfaces by the long steel rods on which they were mounted, somewhat like sugar tongs.
The samples were cut from cylindrical monocrystals grown in the laboratory in lucite tubes, at the rate of 0.05 μm s–1. We used a modified Bridgman method, starting from out-gassed water with a conductivity of 10–5 Ω–1 m–1 (at 25°C and 2 kHz).
To avoid accidental sample fractures caused by mechanical stresses induced by the electrodes during the thermal cycles, the interfaces with the ice were never frozen. Some measurements were also made with mica layers between the ice and the electrodes. The sample surfaces were flattened by lapping with alumina paper and silk cloth, after the removal of an ice layer of about 1 mm contiguous to the cut surfaces. A pressure of about 1 bar was applied to the electrodes to avoid air gaps at the interfaces.
Diagram (not to scale) of the apparatus used for the measurements of the complex admittance at frequencies higher than 0.1 Hz.
To measure the admittance at frequencies higher than 0.1 Hz, we made use of the phase-sensitive system sketched in Figure 1. For lower frequencies, in the range 2 to 10−3 Hz, a circuit quite similar to the other (Fig. 2) measured the integral of the current across the sample during a half-period, starting at the time the applied voltage reached zero or its peak value.
Diagram (not to scale) of the apparatus used for the measurements at frequencies between 10−3 Hz and 2 Hz.
Results and observations
The relative dielectric constant k′ and relative loss k″ were calculated from the experimental values of the equivalent parallel conductance and capacitance, whose absolute values were obtained by comparison with standard resistors and capacitors. As the plots of k′ and k″ as functions of the frequency clearly show the composite nature of the spectrum, we analysed it with the sum of between one and three parallel mechanisms each with Debye relaxation, plus a d.c. conductivity:
where Δki is the dielectric strength of the ith dispersion, τi = τ o exp (Ei/kBT) is its relaxation time, Ei its activation energy, ω the angular frequency, ε 0 the permittivity of free space, and k∞′ the high-frequency relative permittivity of ice, taken as 3.2. The data were fitted to these equations and the parameters obtained.
In Figure 3 we present the values of the relaxation times for one of our samples obtained from the admittance data, plotted against the reciprocal temperature. As stated above, at constant temperature the existence of superimposed relaxation mechanisms is suggested by the many steps in the curves of k′ and k″. When the temperature is changed, the continuity of the relaxation time and of the dielectric strength Δki convinces us that the same mechanism is being observed. In the case of the spectrum usually referred to as “Debye spectrum” (S2 in Fig. 3) there are specific reasons for dropping the continuity requirement.
We observed four spectra, indicated as S1, S2, S3 and S4 in order of increasing relaxation frequency. We did not investigate in detail the spectrum S1 but we simply ascertained its presence, always at the lowest frequencies of measurement. We attributed this spectrum to space charge at the electrodes.
Considering spectrum S2, we observe two distinct behaviours as a function of the temperature. In the range T > T c = 190 K the plot of ln τ 2 has a slope equivalent to an activation energy of 0.64 eV, with a value of τo = 7.7 × 10−16 s. As may be expected a certain scatter in the results is found from sample to sample, a few of them showing a somewhat lower slope (down to 0.60 eV). The value of 0.64 eV coincides with that found by Camplin and Glen (1973) for the Debye spectrum, and is the highest among those reported in literature. Since the pre-exponential factor is also of the right order of magnitude, we identify the spectrum S2 with that due to the reorientation of the water dipoles caused by intrinsic Bjerrum defects.
In Figure 3 it can be seen that at temperatures lower than T c the slope of ln τ2 decreases and then increases again, as in the spectrum reported by Ruepp (1973).
Values of the relaxation times obtained from measurements of admittance (open symbols) and from TSD measurements (full symbols) as a function of the reciprocal temperature. The figures below the designations of the spectra indicate the energies in electron volts deduced from the best-fitting straight lines.
In Figure 4 the dielectric strength of the spectrum S2 is plotted as a function of the reciprocal temperature. As can be seen, at a temperature of around T c Δk has a pronounced decrease. The existence of a minimum in Δk and of a related change in the activation energy of τ is also found for crystals doped with HF (e.g. by Camplin and Glen, 1973). In the treatment given by Jaccard (1959, 1964) this behaviour is expected when the numerical current of Bjerrum defects equals that of ions. These latter defects are replaced by Bjerrumion aggregates in the model proposed by Bilgram and Gränicher (1974). An increase of the minimum temperature for Δk with increasing HF concentration seems now established. A second minimum in Δk appears in samples with an HF concentration higher than 10−5 mol 1−1; this second minimum occurs at higher temperatures than the first one, and seems not to be associated with anomalies of τ. The lowest measured temperature of a Δk minimum in HF doped samples, as far as we know, is about 200 K, reported by Von Hippel and others (1972); the dopant concentration in their crystal was less than 10−8 mol 1−1. The curve of Δk 2 of Figure 4 is similar to a curve of HF-doped ice with a minimum of the first type. The temperature at which our minimum occurs is lower than that which would occur with the doping concentration of about 10−8 mol 1−1 reported above. The authors just mentioned also report that for crystals which are chemically pure but with some amount of crystalline imperfections, the Δk tends to decrease at low temperatures, together with a bending of the ln τ graph. Unfortunately even for the samples that do not show the bending, the reported measurements stop at 190 K. Therefore in the common temperature range of measurement, our results are in substantial agreement with theirs.
Values of dielectric strength Δk obtained from measurements of admittance (open symbols) and from TSD measurements (full symbols) as a function of the reciprocal temperature. The same symbols as in Figure 3 have been used for each spectrum.
However as regards the results of Ruepp and Käss (1969), Ruepp (1973), Gough and Davidson (1970), and Camplin and Glen (1973), the agreement is restricted to the simple bending of the ln τ graph. In the first two papers a change in the activation energy is reported for the intrinsic spectrum of “pure” ice between 230 K and 210 K, while only Ruepp reports that the dielectric strength tends to decrease, at much lower temperatures, about 120 K. Camplin and Glen, for “pure” ice in the range 270–220 K, and Gough and Davidson, for reconverted ice (I → II → I), measured activation energies similar to that we found for τ 2 at the higher temperatures, but did not find any decrease in Δk down to the lowest temperatures of measurement.
At temperatures a little lower than that at which ln τ 2 changes its slope a new spectrum S3 appears, whose dielectric strength, reported in Figure 4, increases rapidly. This spectrum also has a Debye behaviour, with an activation energy of 0.33 eV and a value of τ 0of 4.0 × 10−11 s.
A further spectrum S4, with a dielectric strength almost constant with the temperature (Fig. 4) is observed at the higher frequencies of measurement, down to a temperature of 155 K, after which it was no more observable. Apart from the larger values of the activation energy we found, 0.37 eV, it behaves like that referred to as the high-frequency spectrum by Ruepp (1973). The appropriate value of τ 0 for S4 is 2.5 × 10−14 s.
II. TSD Measurements
Experimental
We investigated the behaviour of the ice samples in the low-temperature range, where the extrapolated values of the previously reported relaxation times are of the order of tens of seconds, by thermally stimulated depolarization current measurements, without changing the experimental set-up.
The TSD measurements were made by cooling the ice samples under an applied static electric field and by subsequently heating them, after the field had been removed. The constant heating rate dT/dt = α was of the order of 0.1 K s−1, and the current was measured as potential difference across a calibrated resistance of 108 to 1010 Ω, by means of an electrometer with 1015 Ω input impedance and 10−3 s rise time.
TSD spectrum obtained by polarizing the sample between 190 K and 77 K, heating rate 0.11 K s−1, polarizing field Fa = 3 × 104 V m−1. (?) experimental points, (——) best fitting curves obtained as described in the text.
Results and observations
The TSD spectrum of “pure” ice crystals has been found to be quite complex from the experiments of Bishop and Glen (1969), Jenevau and others (1972), and Johari and Jones (1975); that of the same sample to which Figure 3 refers is shown in Figure 5. It was obtained by applying a polarizing field F a = 30 kV m–1 from 190 K to 77 K; by so doing it was possible to suppress the whole high-temperature portion of the spectrum, which clearly appears to derive from a charge storage in front of the electrodes. As it is sufficient to apply the field at temperatures lower than 190 K to avoid this effect, we may think we are dealing with mobile charges probably trapped at the low temperatures in sufficiently deep traps, in agreement with the results of Gelin and Stubbs (1965).
The result shown in Figure 5 depends neither on the time during which F a is applied nor on the use of partially blocking (i.e. gold or steel) or totally blocking (with two mica layers) electrodes. Two current peaks, P2 and P3, can be clearly seen, while a third smaller peak is almost hidden in the first part of the curve; this peak can be shown up with respect to the other two, by applying a field of higher intensity than the previous one, in a short temperature range just below the temperature at which its current maximum is expected to occur. In general by the use of this procedure and also by partial discharges of the polarization, one can isolate the contiguous peaks of a TSD spectrum, provided their maximum temperatures are sufficiently distant. We succeeded in doing so for peaks P2 and P3 but not in the case of P4, because of the much higher efficiency of mechanism 3 with respect to 4. Therefore we think that these three peaks are mutually independent, that is their mechanisms can be activated separately.
In the analysis of a complete spectrum or of a peak isolated in the way described above, we made the following assumptions:
-
a. the TSD spectrum is given by the sum of the currents Fi produced by distinct and independent mechanisms;
-
b. each mechanism in isothermal conditions is governed by a relaxation equation of the form
(2)
where P is the polarization and τ the relaxation time;
-
c. for each mechanism τ has a temperature dependence of the form:
(3)with τ 0 and E constant.
Following Kauzmann (1942), τ 0 should be proportional to T –1, but for our purposes this further temperature dependence can, to a first approximation, be neglected with respect to the exponential one.
Equation (2) with the use of Equation (3) and d T/dt = α can be integrated to give the current density J for the relaxation kinetics:
where P(o) is the polarization at the temperature T 0 at which the heating started.
According to Equation (2) the ratio P/J gives the relaxation time of the mechanism at each temperature, whatever its temperature dependence, and if hypothesis (c) also holds, a plot of In P/J versus 1/T would be a straight line. The value of P being found by integrating the values of J since 
these plots were drawn for each of the isolated peaks. We excluded all points with J values lower than 10% of the maximum, but we found that only that plot which referred to P4 was straight in a sufficiently wide temperature range; on the contrary the two plots referring to P2 and P3 showed humps. Furthermore in the plot of P4 the points at the highest temperatures bent upwards. This bending was due to the P3 peak contributing to the J values at the higher temperature, it being impossible to suppress the contribution completely. Therefore to be able to determine the value of the activation energy by a different method, we used that proposed by Land (1969), which makes use of the condition to nullify the second derivative of the current with respect to the temperature, in either of the two curvatures.
The straight line corresponding to the values of E and of τ 0 obtained by this method (0.37 eV and 1.7 × 10–14 s) fits quite well the points of the low-temperature portion of the plot of ln(P/J) as shown in Figure 3; the values are in good agreement with those calculated at the higher temperatures for spectrum S4.
The irregular shape of the plots of ln(P/J) for peaks 3 and 4 may be due to one or both of the following reasons:
-
Equation (2) is not the correct one to describe the relaxation;
-
the observed peak is the result of a continuous distribution of relaxation times or the result of a number of inseparable mechanisms.
We succeeded in maintaining hypothesis (c) by dropping hypothesis (b) and assuming that the isothermal relaxation is described by second order kinetics; that is:
where P* is a constant with the same dimensions as the polarization. With this hypothesis P 2/J gives the product P*τ and the logarithm of this quantity was actually linear as a function of 1/T for the peaks 2 and 3. Also in this case from the slopes of the two best-fitting straight lines we were able to obtain the values of the activation energies of mechanisms 2 and 3, which were 0.38 eV and 0.35 eV respectively; these values are very similar to those obtained for the dielectric spectra S2 and S3 at higher temperatures, but the τ 0 values are still unknown.
Assuming P* to be equal to P(o) (the polarization of the sample at time zero) we found that the pre-exponential factors τ 0 for both P2 and P3 (2.8 × 10–11 s and 1.4 × 10–11 s respectively) were of the right order of magnitude to give relaxation times comparable with those extrapolated from the higher temperatures values. The calculated τ are plotted in Figure 3 together with the best-fitting straight lines. These results were confirmed by means of Land’s method, which in this case refers to the equation:
The total current density .J tot = J 2?+ J 3?+ J 4 calculated with the values of E and of τ 0 reported above, fits the experimental data quite well, as shown by the solid line in Figure 5.
If we wish to investigate also the correlation between total initial polarization (stored charge) and applied field, we cannot have recourse to the experimentally isolated peaks, as the techniques for obtaining them reduce the value of P(o) caused by the applied field. Therefore we determined its true value from the fit of the complete spectrum. Each value so obtained can be thought of as a contribution to the sample polarization that freezes-in when the temperature is lowered while an electric field is being applied. We estimate that this happens when the relaxation time of the mechanism is larger than 103 s. Therefore the values of P(o) i obtained from the fit of the complete spectrum, are associated with the temperatures at which the respective τ have such a value.
The contributions Δki to the dielectric strength were calculated by using the expression:
and are plotted in Figure 4.
III. Discussion
The main purpose of our investigation was the measurement of the single relaxation times of the spectrum of “pure” ice in the widest possible temperature range.
By the use of the same criterion of continuity used in plotting the dielectric spectra, we can conclude formally that the TSD relaxation times refer to the same spectra as are observed by the admittance measurements. To achieve this result we had to use a second-order differential equation to describe the depolarization of mechanisms 3 and 4.
Equation (5), which is formally similar to that describing the recombination kinetics of objects of opposite nature with equal non-equilibrium concentrations, implies a depolarization kinetics different from that valid for a uniform distribution of dipoles given by Equation (2). Therefore one should admit that the probability of each dipole changing its reorientation in unit time, 1/τ, is proportional to the polarization.
In ice, where the reorientation is thought to occur by the migration of defects, we should admit that their velocity or their out-of-equilibrium concentration is proportional to the sample polarization; these same ideas were also suggested by Bishop and Glen (1969). It is not easy to specify them in more detail in microscopic terms, but we think that the movement of space charge might be the dominant part of mechanisms 2 and 3.
Considering more closely the Debye spectrum of ice, S2, we must note first of all, that the slope we found indicates an activation energy for the formation of a pair of L and D defects of 0.80 eV, which is higher than the value usually reported, but is in agreement with the theoretical calculations (Hobbs, 1974, p. 131–40).
If in the ice crystal only the “classic” traffic of Bjerrum and ionic defects (or of Bjerrumion aggregates) occurs, and if only these latter can sustain a d.c. conductivity σsD, then:
when ω ≫ 1/τ D at each temperature.
The current interpretation of the dielectric relaxation in ice is that at higher temperatures σ?D is determined by the majority Bjerrum defects, while σsD is determined by the minority ionic defects. On lowering the temperature a cross-over condition can be reached, where the numerical currents of the two kinds of defects are equal. Then
at each ω.
At temperatures lower than that of the cross-over, an exchange in the roles occurs, with σ?D determined by ions and σsD by the Bjerrum defects.
In measurements like ours, in which more than one mechanism is present, we have:
for w such that ω τ 2 ≫ 1, and therefore also ω τ 1 ≫ 1.
We could also obtain from our measurements σ?2 and σdc separately, but we cannot identify σdc with σ0D, the minority carrier conductivity in the Jaccard model, because from the TSD measurements we also identified, at temperatures higher than 190 K, extrinsic mobile charges.
However we can observe that, far from the cross-over, where σ?D ≫ σsD in Equation (8), τ 2 and σ?2, have the same activation energy. The values we found for τ 2 above and below T c = 190 K, are 0.64 eV and 0.38 eV, corresponding to 14.7 kcal/mol and 8.7 kcal/mol respectively, and are in good agreement with those obtained by Kahane (1969) for σ?D (g kcal/mol and 14 kcal/mol) by the use of ion-exchange membranes. Also the temperature at which the temperature dependence of σsD changes, in Figure 3 of his paper, is very similar to our T c.
Therefore we can reasonably think that we are dealing with a defect cross-over at 190 K.
The efficiency of the ions is suggested by the energy found, 0.38 eV, comparable to that reported for their formation (Bullemer and others, 1969; Engelhardt and others, 1969).
We also found it interesting to compare our results with the model proposed by Bilgram and Gränicher (1974) proposing the presence of Bjerrum-ion aggregates. To do this we calculated σ?D by Equation (8) from our values of τ 2 and Δk 2 measured at high temperatures, and obtained the extrapolated value of σ?2 at T c.
This value was of the order of 10–9 Ω–1 m–1 in good agreement with that at 190 K deduced from the plot of the cross-over conductivity σc reported by those authors.
The meaning of so low a cross-over temperature can be found in the observation that asp must also be lower in our case. Since we do not find any variation of the high-temperature portion of spectrum S2, determined by Bjerrum defects, we must suppose that the concentration of ions available to form aggregates is, in our crystals, appreciably lower than those giving the results reported by Bilgram and Gränicher and referred to in their Figure 1.
Spectrum S4 is typical of a relaxation of dipoles uniformly distributed in the sample. Because of its dielectric strength and of the high relaxation frequency, it is reasonable to attribute it to water molecules loosely bonded to the lattice. Ruepp (1973) proposed that their origin could be found in the microstructures described by Truby (1955).
Although this is certainly possible, we are also attracted by the possibility suggested by Johari and Jones (1976), that inert gas molecules, like those that can form clathrate hydrates, could locally modify the reorientation rate of water dipoles and the corresponding activation energy.
As regards spectrum S3 although it only arises below 170 K and its Δk rapidly attains high values, and could therefore be related to the cross-over, we observe that because of its values of activation energy and relaxation times, it could have an explanation similar to the second one advanced for spectrum 4.
Gough and others (1968), in fact obtained for water molecules in nitrogen clathrate a dielectric relaxation with E = 7.9 kcal/mol and τ (0°C) = 1.5 ? 10–5 s, where we find for spectrum S3 E = 7.6 kcal/mol and τ (0°C) = 4.9 ? 10–5 s.
IV. Conclusions
The apparent purity of some of our samples (if not a compensation) has enabled us to ascertain that in certain conditions Bjerrum defects may remain the majority carriers in governing the dielectric relaxation in ice, at least down to 190 K, and that at lower temperatures this role is played either by ions or by Bjerrum-ion aggregates.
The TSD technique seems to give results similar to those obtained by dielectric relaxation measurements, provided that, in certain cases, one assumes a depolarization kinetics different from that shown to be valid at the higher temperatures, in frequency-dependent measurements.
Further experimental investigations and furthermore a careful analysis of the thermally stimulated depolarization kinetics in real ice crystals will be necessary to clarify this discrepancy.
Acknowledgements
The authors wish to thank Drs S. Mantovani and S. Valeri for many helpful discussions. Thanks are also due to Mr G. Prampolini for valuable technical assistance. This research was supported by the Consiglio Nazionale delle Ricerche under Contract No. 76.00100-02-115.1489.
Discussion
M. ELDRUP : Could you possibly associate one of your relaxation times with the migration of vacancies?
U. DEL PENNINO : In Mogensen and Eldrup (1978) you report values of the vacancy concentration and of their energies of formation and migration that could be compatible with the results that refer to the high-temperature portion of spectrum 2, but the idea of substituting vacancies for Bjerrum defects as the origin of this spectrum would involve a general revision of the theory. As regards spectra S3 and S4, with relaxation times lower than S2 and that can be independently activated, their contributions to the dielectric strength seem to be too high to be associated with the reorientation of molecules adjacent to vacancies or vacancy clusters.
R. TAUBENBERGER : Have you also measured the σ0(1/T) of the so-called Debye dispersion, in other words measured the entire complex ε*(ω)?
DEL PENNING: As already said in the full paper, we have not yet investigated in detail the temperature dependence of σ0.
G. P. JOHARI : Have you looked at the complex-plane plots of the real and imaginary part of the complex permittivity?
DEL PENNING: Yes, in some cases we also considered the plot of k? versus k?. It could be described as a “sum” of semicircular curves.
JOHARI : I would like to point out that there is a serious problem in thermally stimulated depolarization measurement, namely, one of identifying the orientation polarization peak from space-charge polarization peaks. The only method of distinguishing between them that I know of is to measure the total charge associated with the peak as a function of the magnitude of the electric field applied and also as a function of either the thickness and/or the area of the ice sample. If the total charge under a peak increases linearly with the field and also linearly with the geometric capacitance of the cell, the peak is due to orientation polarization. If not, it is due to space-charge polarization. I wonder if your data satisfy the criteria for the identification of the peak as due to the orientation polarization ?
DEL PENNINO: I agree with you on the criteria; furthermore in a case like this in which the temperature during the polarization changes linearly with the time, the maximum temperature of a peak T m should be independent of the way in which the sample is polarized if the peak is due to orientational polarization. In the case of space-charge, T m should depend on the procedure used to polarize the sample, because in this case the initial charge distribution changes accordingly. During the procedures to isolate the peaks experimentally we observed that T m is independent of them. As regards the relationship between field intensity and the charge associated with the peak, we refer, for the reasons stated above, to the values of P(0) i used to fit the complete TSD spectrum. For the same sample these values change linearly with the field up to 30 kV/m, but we intend to verify this behaviour at higher field intensities to have conclusive results. To change the width of the sample we had to use, of course, three samples of different thickness. In this case we have some scatter in the results, but they are still compatible with a linear relationship between charge and width.
