Piezoelectric materials have wide band gap and no inversion symmetry. Only the orthorhombic phase of lead metaniobate (PbNb2O6) can be ferroelectric and piezoelectric below Curie temperature, but not the rhombohedral phase. High temperature piezoelectric applications in current decades have revived international interest in orthorhombic PbNb2O6, synthesis of which in pure form is difficult and not well documented. Second problem is that its impedance spectroscopy (IS) data analysis is still incomplete. Present work attempts to fill up these two gaps. Presently found synthesis parameters yield purely orthorhombic PbNb2O6, as checked by X-ray Rietveld analysis and TEM. Present 20 Hz to 5.5 MHz IS from room temperature to 700°C shows its ferroelectric Curie temperature to be one of the highest reported, >574°C for 0.5 kHz and >580°C for 5.5 MHz. Dielectric characteristics and electrical properties (like capacitance, resistance and relaxation time of the equivalent CR circuit, AC and DC conductivities, and related activation energies), as derived here from a complete analysis of the IS data, are more extensive than what has yet been reported in the literature. All the properties show sharp changes across the Curie temperature. The temperature dependence of activation energies corresponding to AC and DC conductivities has been reexamined.
1. Introduction
A generator of pressure or movement or ultrasonic wave and their sensor/detector can be fabricated utilizing inverse and direct piezoelectric effect, respectively. These possibilities opened up a huge array [1] of medical, industrial, and other applications of piezoelectric materials with significant commercial implications. Commercial piezoelectric materials have mostly been barium titanate (BT), lead zirconate titanate (PZT), or materials based on BT or PZT. But Curie temperature (Tc), the upper limit for piezoelectricity, is at best 130°C [2] for BT and 365°C for a modified PZT. So, higher Tc materials are being developed worldwide for high temperature (HT) applications. Lead metaniobate (PbNb2O6), shortened here as PNO, with a higher Curie temperature (517 to 570°C in different papers), is one of the present candidate materials [2–14]. It was discovered [3] in 1953. However, synthesis of its piezoelectric phase, the metastable orthorhombic structure, in a pure state is difficult [4]. During the preparation of the orthorhombic PNO by quenching from a temperature above 1250°C, a few competing compounds and phases (like the rhombohedral PNO phase) tend to form. The rhombohedral PNO is not piezoelectric. In fact, the difficulty [4] of forming single-phase and high-density orthorhombic PNO, together with the success of BT, PZT, and related piezomaterials for near room temperature applications, somewhat halted R & D on PNO after the initial years of pioneering work and publications [5–9]. In the meantime, a need has developed for high temperature (HT) piezoelectric materials (i) for HT ultrasonic imaging of nuclear fuel rods through molten metal coolants [15] in Fast Breeder Reactors and (ii) for HT applications like sensors in car exhaust systems. So, since the late nineties, work on PNO has been revived by a few groups (Table 1) in search of higher Curie temperature materials. The stable forms of PbNb2O6 are rhombohedral (at low temperature) and tetragonal (at high temperature). Quenching from the tetragonal form, as mentioned earlier, leads to metastable orthorhombic form at room temperature. Roth showed by X-ray diffraction [5, 6] that the orthorhombic PNO (ferroelectric) changes to a tetragonal structure (paraelectric) above ~570°C.
Comparison of recently reported Curie temperatures of orthorhombic PbNb2O6 samples.
Goodman [3] had earlier concluded the T(Curie) of quenched PNO to be 570°C, from their dilatometric and dielectric characterizations. Although this temperature is regarded even in recent papers [10] as the Curie temperature of PbNb2O6, almost all the recent electrical measurements have obtained much lower values of its Curie temperature (Table 1). Reference [10] quoted T(Curie) of 500°C for Ferro-perm (Denmark) PNO samples (1995) and 495°C for Morgan Matroc (UK) PNO samples (1997). Lower T(Curie) in these commercial samples may or may not be due to efforts to optimize other piezoelectric parameters. Table 1 has been compiled from recent publications. It shows T(Curie) to be in the range 517°C to 540°C only, except in present samples. So, there is a need to optimize preparation conditions of quenched PbNb2O6 to get highest possible Curie temperature, high purity and high density. A pointer in this direction is our observation on quenched PNO [11] of an endothermic signal at 570.6°C in Differential Scanning Calorimetry [12].
Here, our preparation of high quality orthorhombic PNO and results of our impedance spectroscopy up to 700°C and over a wide frequency range from 20 Hz to 5.5 MHz, with a detailed analysis, have been presented. Intrinsic cationic defects in niobates like LiNbO3 and their role are well known [16]. Likely lattice defects in a quench-prepared PbNb2O6 sample should similarly be regarded as an intrinsic part of the sample. The structure of the samples has, therefore, been characterized fully before other studies. Although some work on PNO (Table 1) and a few on various targeted chemical substitutions in PNO have been reported in the literature, not all possible dielectric and electrical parameters for PNO itself have yet been derived from IS data. This has been attempted here. The dielectric response in such an electroceramic is known to be complex with many contributions that need to be known and understood. So, desirable observations like the highest T(Curie) as well as all uncommon observations of the present work are reproduced. All possible data analysis has been carried out, yielding important parameters like bulk resistance, relaxation times, activation energies, AC and DC conductivity, and phase angle.
2. Materials and Methods
Aldrich lead(II) oxide (PbO, 99.999%, yellowish orange and Melting Point = 886°C) and niobium pentoxide (99.9% GR) of Loba Chemie, weighed in proper proportion to give Pb : Nb = 1 : 2, plus 2% of excess PbO, have been well ground under isopropyl alcohol, dried, and pelletized. Additions of excess PbO were introduced by some earlier workers to compensate possible evaporative loss of PbO. To get best PbNb2O6, present firings have been done in stages at 1050°C, 1290°C and 1270°C in a box furnace with PID temperature controller. Extra PbO (3%) has been added before the 3rd firing. All our firings have been done in palletized form to provide intimate contact and reduce Pb loss. But grinding and repalletizing at each step have also been undertaken, to improve mixing and the sample quality. After the last firing, the samples have been quenched. The samples have been characterized by XRD with Rietveld analysis [11] and Transmission Electron Microscopy (TEM).
The impedance spectroscopy has been carried out using Solatron “SI 1260 Impedance Gain Phase Analyzer” with a high temperature attachment. Our PNO pellets (diameter ~ 12 mm and 2.2 to 2.4 mm thick) have not been poled before measurement. Silver paint coat on the flat faces of the pellet under study served as the electrodes for IS measurements. This measurement of ε′ and ε′′ (the real and imaginary components, respectively, of the relative dielectric constant i.e., permittivity) has been done for different measuring frequencies at different values the sample temperature (T) while heating from room temperature to 700°C.
3. Results and Discussion
XRD pattern and Rietveld analysis [11], Figure 1(a), confirmed the orthorhombic structure and absence of any other phase. The smallness of the difference of the fitted graph and raw data proved our Rietveld fit to be good. It yields a=17.6468(7) Å, b=17.9512(7) Å, and c=3.8704(1) Å. The samples showed about 80% of the theoretical density on average. TEM diffractograms [6] with the electron beam on individual grains of the ceramic sample always gave diffraction patterns that matched these lattice parameters from XRD result. One TEM diffraction pattern is shown in Figure 1(b). This way, all of the nearly forty grains examined have been found to be in orthorhombic phase.
(a) Rietveld fit (red graph) to XRD data points (black circles) for orthorhombic PbNb2O6 with their difference shown by the continuous curve in the smaller-sized graph. The vertical blue bars below the main graph are the calculated peak positions. (b) Electron diffraction pattern for HR Transmission Electron Microscopy on a selected grain (microcrystal) of quenched lead metaniobate (PNO).
Figures 2–4 show ε′ versus frequency, ε′′ versus frequency, and tanδ=ε′′/ε′ versus frequency results (over 20 Hz to 5.5 MHz range) for 12 temperatures from 50°C to 700°C. In Figure 2, real part (ε′) of dielectric constant decreases with increase of measuring frequency up to about ln f=12 to 14 region, implying ~162 kHz to ~1203 kHz region, and then increases sharply with increase of frequency as reported earlier [13] for all temperatures except at low temperatures like 200°C. Moreover, ε′ is low and less dependent on frequency at low temperatures up to 300°C, leading to overlap of several graphs. It increases with increase of temperature, up to 573°C. Graph of ε′ shows significantly increased values in the 573°C graphs with respect to the 570°C graph. This temperature is the region of ferroelectric to paraelectric phase transition, as is clearer in graph on temperature variation (Figure 5). One must note the interesting decrease of ε′ with further increase of temperature above 573°C. This, in fact, shows up in ε′ versus T plot as peaking at around 573°C. Here, it was studied up to 700°C. The values at different temperatures of dielectric dispersion with respect to frequency, which will be proportional to dε′/d(lnf), will be discussed from the low frequency, specifically lnf (Hz) <7, part of Figure 2. At a low temperature of 100°C, dε′/d(lnf)=2.483±0.739 only, increasing to 41.158±7.805 at 300°C, whereas at the high temperature of 573°C, dε′/d(lnf)=-4121.578±226.784, a high value.
Frequency dependence of the real part of dielectric constant (ε′) in the frequency range 20 Hz (lnf=2.995) to 5.5 MHz (lnf=15.52), compiled from heating runs at different temperatures for quenched PNO sample, the orthorhombic phase of PbNb2O6. The inset is an enlarged view of the variation of ε′ in the higher frequency range.
Frequency dependence of the imaginary part (ε′′) of dielectric constant in the frequency range 20 Hz to 5.5 MHz, compiled from heating runs at different temperatures, for quenched PNO sample, the orthorhombic phase of PbNb2O6. The inset is an enlarged view of the variation of ε′′ in the higher frequency range.
Frequency dependence of the dielectric loss (tanδ) in the frequency range 20 Hz to 5.5 MHz, compiled from heating runs at different temperatures, for quenched PNO sample, the orthorhombic phase of PbNb2O6.
Temperature dependence of the real part (ε′) of dielectric constant in the temperature range of 50°C to 700°C, at different frequencies, for quenched PNO sample, the orthorhombic phase of PbNb2O6.
In Figure 3, frequency dependence of ε′′ or the imaginary part of relative permittivity is more at low frequencies and higher temperatures, the ε′′ versus frequency graph for 700°C being at top in the lower frequency region. Here, ε′′ decreases steeply as frequency increases from 20 Hz to about lnf=8 implying f=2.98 kHz. A weak peaking of ε′′ at lnf=13.8, implying f=984.6 kHz, is visible for graphs for low temperatures like 50°C, 300°C, 500°C, and 550°C in the enlarged view. Our higher temperature graphs show no peak in this temperature region but only a slow increase of ε′′ to allow peaking at higher frequencies. This fall and peaking of ε′′ has been shown and discussed by earlier workers [13], although the details differ. Dielectric loss (tanδ) in Figure 4 approximately follows ε′′ with respect to frequency and temperature dependence.
Figures 5 and 6 are plots with respect to temperature, of ε′ and tanδ=ε′′/ε′, respectively. The peak value of ε′ is denoted by εmax′ and its position on the temperature axis by Tm. Data for eight frequencies, selected out of 20 Hz to 5.5 MHz data in Figures 2–4, have been plotted in Figures 5 and 6 to avoid overlaps and provide a clearer view of the graphs.
Temperature dependence of the dielectric loss (tanδ) in the range 50°C to 700°C at different frequencies for quenched PNO sample, the orthorhombic phase of PbNb2O6. Two insets magnify two important features.
Fall of ε′ and, hence, of the dielectric susceptibility, χ′= (ε′-1), with increasing T, a characteristic of the paraelectric phase and exhibited here for T>Tm region, is due to opposition of the increasing thermal agitation, to dipole alignment. The rise of ε′ with T in the T<Tm region is believed to be achieved by the increase in the number of dipoles available for alignment [21]. Sharpness of the peak at Tm indicates near absence of diffused scattering in present samples, as discussed elaborately later. Let a minor peculiarity be discussed first. A wide minimum (Figure 5), appearing at 300°C, for frequencies 0.5 kHz, 1 kHz, and 5 kHz in ε′ versus T, as well as in tanδ versus T graph (Figure 6), cannot be readily explained. However, the more important feature of this work is the large permittivity peak at Tm with Tm ranging from 574°C to 576°C as the measurement frequency is increased from 0.5 kHz to 4 MHz, as mentioned earlier. In Figure 7, the peak value, εmax′, of ε′ shows a sharp fall initially and then a slow rise with increase of measuring frequency. The value of εmax′ decreases from 11623 to 1450 as the frequency (f) increases from 20 Hz to 1 MHz. Then, for f>1 MHz, the value of εmax′ increases slowly, being about 5465 at 5.5 MHz. The prominent peak of the imaginary part ε′′ in the ε′′ versus T graph can be seen to be around Tm as expected. The dielectric loss (tanδ) also peaks at around Tm.
The maximum value (εmax′) of the real part of dielectric constant in the in the temperature range of 50°C to 700°C, as a function of frequency (Hz), for quenched PNO sample, the orthorhombic phase of PbNb2O6.
Above Tm, ε′ decreases with increase of T, and, hence, the material is obviously paraelectric above Tm. It is ferroelectric below Tm. This transition temperature gives [22] the Curie temperature. So, present impedance spectroscopy measurements indicate a high Curie temperature that depends on frequency, 576°C for 4.0 MHz and 580°C for 5.5 MHz, for example.
The plot of Z′′ (imaginary part of the impedance) versus Z′ (real part of the impedance), called Nyquist diagram [23], is based on the well-known Cole-Cole plot of ε′′ versus ε′. The Z′′ versus Z′ plot for our orthorhombic PNO is shown in Figures 8(a) and 8(b). The straight line segments joining of the data points in Figures 8(a) and 8(b) are drawn only as a guide to the eye. Single semicircular nature of the plots is clearer in the higher temperature plots, as the lower temperature plots, like that for 200°C in the insert of Figure 8(a), show a near linear rise, possible for a large semicircle. An important observation is that the Nyquist diagram here shows only one semicircular arc [23] over the studied frequency range 20 Hz to 5.5 MHz. One semicircle implies only one contribution, here from the sample grains and not from any second source (e.g., grain boundaries or electrode effect). A sample with one semicircle Nyquist diagram or Cole-Cole plot is equivalent [23] to a bulk resistance R and a capacitance C in parallel. These elements give rise to a time constant τ=CR, the dielectric relaxation time of the basic material. The frequency ωm of the peak of the Nyquist diagram can find the relaxation time τ as (ωm)(τ)=1 for this peak point [23]. At this point, Z′′ has its maximum value. So, many authors [24, 25] use the frequency of peak (i.e., maximum) of Z′′ versus frequency plot to find the relaxation time τ.
(a) Plot of the imaginary part, Zimaginary in ohm, of impedance versus its real part, Zreal in ohm, for quenched PNO sample, PbNb2O6 in orthorhombic phase, at temperatures: 500°C, 550°C, 570°C, 573°C, 600°C, 650°C, and 700°C. Since the complex impedance Z*=Zreal-jZimaginary, quantity (-Zimaginary) has actually been plotted in the figure. (b) Plot of the imaginary part, Zimaginary in ohm, of impedance versus its real part, Zreal in ohm, for the quenched PNO sample, PbNb2O6 in orthorhombic phase, at temperatures: 400°C, 440°C, 460°C, and 480°C. Since the complex impedance Z*=Zreal-jZimaginary, quantity (-Zimaginary) has actually been plotted in the figure.
Present work additionally suggests an alternative method (based [23] on Z′′=R/2=Z′ criterion) of determining relaxation time (τ). We later compare the two methods, with the finding that they match (Figure 14) fairly well, match being better at higher temperatures. Often the centre of the semicircular arc of the Nyquist diagram is not on the Z′-axis, so that the semicircle is seen to be depressed by an angle (θ) below the Z′-axis. This is the case, if there is a distribution of relaxation times. This θ is determined by the width of the distribution of relaxation time.
If the semicircle starts near the origin at Z′=R∞ and intercepts the Z′-axis again at R0, R=(R0-R∞) gives [23] the bulk resistance. Fitting the impedance spectroscopy (Z′′ versus Z′) data, usually forming an arch, to the equation of a circle, its radius r and θ and, hence, R0 can be determined. Bulk resistance R=(R0-R∞)=2r·cosθ has thus been calculated from our data in Figure 8 and plotted in Figure 9. It is found that the semicircle of data points for a particular temperature in Figures 8(a) and 8(b) reduces in size with increase of temperature, implying smaller bulk resistance of the sample at higher temperatures as documented in Figure 9. Real part ε′ of permittivity in present samples in the paraelectric region has been fitted to the modified Curie-Weiss law [26], (1/ε′-1/εmax′)=(T-Tm)γ/Cw, where Cw is Curie-Weiss constant and γ is the diffusivity parameter with 1<γ<2. For a ferroelectric relaxor, γ is at its highest value of 2.
Bulk resistance R in MΩ and capacitance C in nF of the equivalent RC circuit of quenched PbNb2O6 sample, plotted against temperature, T. Here, R and C have been calculated from Nyquist diagram of Z′ and Z′′ in 20 Hz to 5.5 MHz range and the corresponding time constant τ=CR, related to the dielectric relaxation time.
The previous equation reduces to Curie-Weiss law for γ=1. The full equation can be rewritten as ln(1/ε′-1/εmax′)=γln(T-Tm)-lnCw. So, plotting of ln(1/ε′-1/εmax′) versus ln(T-Tm) for different frequencies, as shown in Figures 10(a) and 10(b), gives the fit parameters (Table 2). Present finding that γ is nearly 1 in all cases proves that the relaxor property is minimum for these samples. This γ=1 leads to the observed sharpness of the peak at Tm. Relaxor piezoelectric materials, on the other hand, show a broadened peak.
Result of fitting ε versus T data to the modified Curie-Weiss equation, as detailed in the text.
Frequency (Hz)
Tm (°C)
Cw×105 (°C)
γ
εmax′
5.5 M
579.5 ± 4
2.187
1.141
5472.6
3.0 M
575.2 ± 2
1.739
1.136
1876.1
2.0 M
574.8 ± 1
1.471
1.110
1574.2
50 k
574.9 ± 1
1.236
1.035
1752.2
1.0 k
573.8 ± 1
5.237
1.283
3613.9
0.5 k
573.0 ± 1
2.635
1.117
3632.5
0.1 k
572.5 ± 1
0.877
0.744
7103.5
(a) Plot of ln(1/ε′-1/εmax′) versus ln(T-Tm) for quenched PNO sample, the orthorhombic phase of PbNb2O6, at various frequencies (3 MHz, 2 MHz, 1 MHz, 50 kHz, and 1 kHz). It examines the relationship, expected for the paraelectric phase of a normal ferroelectric. (b) A specific example of finding the parameters (γ=1.0346 and C=123569.04°C) of the modified Curie-Weiss law from the 50 kHz plot of ln(1/ε′-1/εmax′) versus ln(T-Tm) for the paraelectric phase of quenched PbNb2O6, which is orthorhombic.
Next, AC conductivity (σac) has been calculated from the dielectric loss (tanδ) data using the empirical formula, σac=ωε0ε′tanδ=2πfε0ε′′, where ε0=107/4πc2=8.845×10-12 F·m−1 = vacuum permittivity and c = speed of light in vacuum. Figure 11 presents this result on AC conductivity as a function of 1000/T. Now the maxima in different graphs of Figure 11 (prominent only in high frequency graphs) are at 1000/T~1.18 K−1, that is, at T~573°C that corresponds to the ferroelectric-to-paraelectric transition. AC activation energy (Eac) was calculated using the empirical Arrhenius relation σac=σ0 exp (-Eac/kT).
Frequency dependence of AC conductivity in the temperature range 50°C to 700°C, at different frequencies, for quenched PNO sample, the orthorhombic phase of PbNb2O6.
The aforementioned maximum, known to occur in the AC conductivity graph at around the ferroelectric-paraelectric phase transition temperature (Tc), is due to domain reorientation, domain wall motion, and the dipolar behavior [27]. In our data in Figure 11, there are three distinct temperature regions of Arrhenius type conduction, (a) ~50 to ~300°C section, corresponding to [28] electron or hole hopping, (b) ~400 to ~550°C part, due to conduction by small polarons and oxygen vacancies, and (c) ~580 to ~700°C high temperature region due to intrinsic ionic conduction [29]. Data in these three regions have been separately fitted to the Arrhenius relation as shown in Figure 12. The slope of ln(σac) versus 1/T plot gives AC activation energy Eac.
AC conductivity activation energy, Ea, as calculated from the slope of the Arrhenius plot of ln(σac) versus 1/T, resulting in three Ea values for three temperature ranges. This figure depicts the variation of Ea with measuring frequency.
The previous fitting avoided the region between 300 and 400°C, as it showed an unexplained peak. However, the rise in some graphs at ~300°C to ~322°C region of Figure 11 reminds one of the anomaly at ~300°C in ε′′ versus T as well as in tanδ versus T graphs in Figures 5 and 6. It is noteworthy that temperature dependence of σac is radically different below and above this temperature zone, indicating different AC activation energies for AC conduction. However, somewhat wavy nature of the data at low frequencies and low temperatures, in this and other work [30], and suggestion of the rise by only one point in high frequency graphs prevent definite conclusions on the rise. These data at low frequencies (here 0.5, 1, and 10 kHz in Figure 11) and low temperature (T<401°C) have not been considered for present AC activation energy calculation in next paragraph. In the high temperature graph in Figure 12, Eac falls rapidly with increase of measuring frequency up to 500 kHz and then becomes weakly dependent on frequency at higher frequencies. The low temperature graph of Eac is practically independent of frequency over the range (500 kHz to 5 MHz) studied here. As seen in the graphs, the AC activation energy is higher at lower frequencies, 0.5 kHz data giving 0.954 eV (with σ0=57.56 Ω·m−1) for 400 to 550°C, and 1.17 eV (with σ0=1.508×103 Ω·m−1) for 580 to 700°C. These values fall to 0.523 eV (with σ0=24.62 Ω·m−1) for 400 to 550°C and −1.994 eV (with σ0=6.253×10-13 Ω·m−1) for 580 to 700°C, that is, in paraelectric region, respectively. Figure 13 gives the variation of Zreal=Z′ and Zimaginary=Z′′ against the measuring frequency. Here, Z′ and Z′′ graphs are seen to merge in the high frequency region approaching zero value. The single semicircular Cole-Cole plot suggests an equivalent circuit [23] of resistance R and capacitance C in parallel giving the well-known complex impedance Z*=R/(1+jωRC)=Z′-jZ′′. Applying the relation τ=CR and the earlier-mentioned condition: ωmτ=1, one gets Z′′=R/2 and Z′=R/2. So, relaxation time related frequency is the frequency at which Z′ versus frequency and Z′′ versus frequency graphs cross each other making Z′=Z′′. We calculate this relaxation time and present it as the relaxation time τcross in Figure 14. We determine relaxation time also from the frequency at the Z′′ peak of Z′′ versus frequency graphs (Figure 13) as done by others [30], calling it τpeak to compare it in Figure 14 with τcross. Figure 14 shows that τcross and τpeak behave similarly proving correctness of each approach. But these two graphs show different values with more relative scatter in the τpeak graph. This point shows up again in DC conductivity discussion in next paragraphs.
Variation of Zreal (real part of impedance) and Zim (imaginary part of impedance) with frequency (plotted in log scale) at different temperatures. Frequency at the intersection of the two graphs gives a value of the relaxation time of the system, to be called τcross. Relaxation time estimated from the frequency) at the peaking in Z′′ is to be called τpeak.
Variation of the relaxation time (τpeak and τcross) with temperature (T) has been shown in one graph, relaxation times τpeak and τcross being estimated by two different techniques as defined in Figure 13.
The decrease of relaxation time with temperature has been adequately discussed in the literature [31]. Relaxation time for this material is seen to fall with increasing temperature, except for a prominent peaking at the ferroelectric to paraelectric transition temperature. R=(R0-R∞)~R0, since R∞≪R0. This R has already been calculated from Figure 8. So, presently determined CR=τ value gives us C. These two equivalent circuit parameters C and R have been presented in Figure 9. Here, there is an impressive peaking of capacitance C at the Curie temperature. Another related parameter, the angle (θ), has been plotted against temperature in Figure 15. Clearly, it may be positive or negative, depending on the position of the centre of the semicircular arc. Variation of the impedance phase angle, δ=tan-1-(Z′/Z′′), with frequency, plotted on a logarithmic scale, is called Bode plot [23, 26]. Bode plots at different temperatures have been shown in Figures 16(a) and 16(b), graphs for 570°C, 573°C, and 600°C being given in Figure 16(b) to avoid overlap. At low frequencies and high temperatures (e.g., the 700°C graph), the constant phase part of the phase angle can be seen to be predominating. For a pure capacitor, δ=-90°, while for a resistor it is zero. For a capacitor, the voltage, Vc, lags behind the response current by 90°. In the phase angle plot, an approach to pure capacitive behavior at low frequency values can usually be identified by δ→-90°. For an RC circuit, it will have an intermediate value and depends on ω as seen in the graph. Figure 16(b) contains data sets at Tm region temperatures and shows a prominent increase of the phase angle between 1 MHz and 5.5 MHz. The frequency ranges for near constancy of the phase angle (e.g., over 20 Hz to ~1 kHz for 700°C) can be observed from the figures. This constant phase region appears at low frequencies for the high temperature runs and at high frequencies for low temperature runs.
Diagrammatic definition of the angle (θ) of Nyquist Plot (in the inset), as detailed in the text, and its variation with temperature.
(a) Variation of the impedance phase angle, δ=tan-1-(Zreal/Zim), with frequency at different temperatures (255°C, 340°C, 400°C, 450°C, 500°C, and 700°C). (b) Variation of the impedance phase angle, δ=tan-1-(Zreal/Zim), with frequency at a few important temperatures (570°C, 573°C, and 600°C).
DC conductivity has been calculated using the relation [32], σdc=(ε0ε∞)/τ from experimentally found τ and ε∞(=ε′ value as f→∞). Using τcross and τpeak data, as found in Figure 14, temperature dependence of σdc(cross) and σdc(peak), and, hence, that of two sets of DC activation energy have been calculated for the present PNO samples. Result from each technique shows expected increase of σdc with temperature (T), as shown in Figure 17. DC activation energy (Edc) due to DC conductivity [14] was calculated using the empirical Arrhenius relation σdc=σ0exp(-Edc/kT), where k is Boltzmann constant, σ0 is preexponential factor and Edc is the DC activation energy for the conduction process [33–36], by plotting ln(σdc) versus 1/T. Figure 19 gives this DC activation energy for different temperature ranges, indicated by the horizontal bars, as found independently from τcross and τpeak data, for two- and three-region fits. Almost all earlier authors determined DC activation energy, Edc, from relaxation time via DC conductivity. The 2-segment fit has been utilized in the literature for the usual τpeak-based evaluation of Edc below and above T(Curie).
Variation of DC conductivity, σdc(cross/peak), with temperature in 400°C to 700°C range, for quenched PNO sample. Here, two values of DC conductivity, σdc(cross)=(ε0ε∞)/τcross in the lower graph and σdc(peak)=(ε0ε∞)/τpeak in the upper graph have been calculated from τcross and τpeak given in Figure 14.
Gonzalez et al [36] calculated relaxation times (τ=τCC,τT,τ*) and other parameters for PNO sample from their 450°C to 590°C measurement of the temperature dependence of the imaginary dielectric modulus, separately by frequency domain analysis and by time domain analysis. Their τ*, relaxation time due to thermal effect, was calculated in time domain using the Kohlrausch-Williams-Watts function, f(t)=e-(t/τ*)β. Finally 〈τT〉, the average relaxation time, was calculated using the relation, 〈τT〉=(Γ(1/β)/β)τ*, where Γ is the Euler gamma function. Next, σdc was calculated using the relation σdc=ε∞/〈τT〉. They made separate liner fits to their ln(σdc) versus 1/T data for data below T(Curie) and for data above T(Curie), finding, however, similar values Edc=1.08 eV below T(Curie) and Edc=1.11 eV above T(Curie). Finding practically same values apparently implies a single conduction mechanism [36]. However, the fact that more than one conduction mechanism is expected [36] from the observed (Figure 2) “departure of the real permittivity data from a low frequency plateau-like behaviour towards the lowest frequencies and high temperature” has also been pointed out by these authors, with the comment that “further measurements and data treatment to better elucidate this additional low-frequency effect are in progress.” This motivated us (i) to use τcross to find σdc and Edc, and also (ii) to try three segment fitting (Figure 18) as needed for the data resulting from the τcross-based analysis. The data actually demands a middle segment (here 570°C to 600°C) fit.
Variation of ln σdc(cross) with 1/T(K-1) in temperature 400°C to 700°C range, and least square fitting of the liner equation in three different temperature regions that yield three different values of DC activation energy.
Activation energy (Edc) for DC conduction in different temperature regions. P in symbol name indicates the use of peak technique of finding relaxation time, and C indicates the cross technique case. Here, 2 stands for fit to 2 temperature regions and 3 for 3 temperature regions fit, the chosen temperature region of fit being indicated by corresponding horizontal bar through the symbol.
Here, one notes in Figure 17 that σdc(cross) versus T data clearly fall into three linear segments:
400°C to 565°C giving (Figure 18) Ea=1.0724±0.043 eV,
570°C to 600°C giving (Figure 18) Ea=1.595±0.120 eV,
610°C to 700°C giving (Figure 18) Ea=1.263±0.156 eV,
which implies a peaking of the activation energy Ea=Edc (cross) in the Curie temperature (Tc) region and practically the same value for low and high temperature regions (“C3” data points in Figure 19). This discovers the expected [36] existence of different Edc and, hence, different conduction mechanisms, retaining same or similar values of Edc below and above the Tc region.
A 2-segment fit to σdc(cross) versus T data of Figure 17 yields (“C2” points in Figure 19)
Ea=1.072±0.043 eV over 400°C to 565°C,
Ea=1.026±0.086 eV over 570°C to 700°C,
missing the significant rise and fall within a narrow Tc region (570 to 600°C) to wrongly conclude no significant change across Tc (the Curie temperature). This is in line with the limited knowledge, currently available in the literature as already discussed.
σdc(peak) versus T data of Figure 17 can be seen to have more scatter than σdc(cross) versus T data with the former not falling well either on a set of straight lines or even on a smooth graph, as should be the case from already discussed higher scatter of τpeak data. We feel that τcross should be used as relaxation time if there is a choice to choose from τpeak and τcross data. Still we tried linear fits to σdc(peak) versus T data. The 3-segment fit could tackle this scattered data set relatively better as follows:
400°C to 480°C Ea=1.028±0.101 eV,
500°C to 590°C Ea=1.523±0.130 eV,
600°C to 700°C Ea=1.387±0.168 eV,
as shown in Figure 19 as “P3” data points.
So, the three-region fits of τcross and τpeak suggest a peaking in Edc versus T variation, at or below (resp.) the Curie temperature. This better match of the Edc peak with the Curie temperature for τcross data further supports the better reliability of τcross data. Near the transition temperature, DC activation energy is highest (1.595±0.120 eV), a new observation that needs further verification and explanation.
4. Conclusion
Presently prepared PbNb2O6 samples revealed pure orthorhombic phase by X-ray Rietveld analysis, with confirmation down to granular level, by High Resolution Transmission Electron Microscopy. These structural characterizations as well as impedance spectroscopy over 20 Hz to 5.5 MHz and up to 700°C established the high quality of these samples.
Present evaluations of DC and AC conductivities, capacitance (C) and resistance (R) of the equivalent CR circuit, relaxation times and related activation energies from the analysis of our wide-ranging impedance spectroscopy data on orthorhombic PNO, and discussions have been either more extensive than in the literature or done for the first time. This is the first reporting of the Nyquist diagram, the plot of Z′′ versus Z′, for orthorhombic lead metaniobate. An impressive peak (~760 nF) of the equivalent CR circuit capacitance (C) in the C versus T graph has been obtained at Tm as shown. AC conductivity (σac) is seen to peak at Tm, more prominently in the higher frequency plots. Relaxation time registers a clear peak at Tm, while falling slowly with increasing temperature in other temperature regions.
Measured DC electrical conductivity appears to show three regions of temperature dependence (Figure 17) rather than the usually tried ferroelectric and paraelectric regions. A suitably chosen transition region in the plot of ln(σdc) versus 1/T in Figure 18 forms the third region, interestingly giving the highest activation energy (Figure 19). A 3-region fit should be carefully examined for different systems, since the third region may be missed due to data points, inadequate with respect to either the temperature range or the density of data points or both. Moreover, a forced two-region fit (C2 points in Figure 19) to data needing a three-region fit will mask the large up and down (C3 points) of the activation energy value in the Curie temperature region and indicate a small change of activation energy across the Curie temperature, similar to the result of 2-segment fit to time domain analysis [35] result. It appears that the DC activation energy can be estimated better from τcross and 3-segment fit (Figure 19), unless the result in a particular experiment is against 3-segment fit. Presently observed peaking of the activation energy Edc(cross) in the Curie temperature (Tc) region is a new result.
The peak at Tm in the ε′ versus T graph is sharp, leading to a value of ~1 for the diffusivity parameter and implying near absence of any relaxor property. Impedance spectroscopy showed (Table 1) highest record of ferroelectric-paraelectric transition temperature (Tm), Tm>580°C for 5.5 MHz, and Tm>574°C for 0.5 kHz. This implies potential use of this material for piezoelectric sensors and actuators up to temperatures much higher than the limits of BT and PZT based piezomaterials.
Acknowledgments
The dielectric measurement set-up at the Deptartment of Ceramic Engineering, NIT, Rourkela 769008, India, has been used with the active encouragement of Professor Swadesh K Pratihar. Authors also thank Shri Pulak Ray and his team at SINP for help in TEM operation. U. De thanks Alexander von Humboldt Foundation (Germany) for partial support that led to preliminary presentation of a part of this work in international conferences and discussions.
IharaI.Ultrasonic sensing: fundamentals and its applications to nondestructive evaluation20082128730510.1007/978-3-540-69033-7_14KazysR.VoleisisA.VoleisieneB.High temperature ultrasonic transducers: review2008632717GoodmanG.Ferroelectric properties of lead metaniobate1953361136837210.1111/j.1151-2916.1953.tb12820.xNesterovA. A.KaryukovE. V.MasurenkovK. S.Synthesis of phases of composition PbNb2O6 and BaNb2O6 with the use of active precursors20098233703732-s2.0-6544912392210.1134/S1070427209030057RathR. S.Unit-cell data of the lead niobate PbNb206195710part 643743710.1107/S0365110X57001437DeU.TEM, HT XRD, electrical, optical and thermal characterizations of materials for high temperature piezoelectric applicationsProceedings of the IUMRS-ICA2012Busan, KoreaBEXCOAbstract no. 21_1078http://www.iumrs-ica2012.org/program.phpFrancombeM. H.LewisB.Structural, dielectric and optical properties of ferroelectric leadmetaniobate195811part 10696670310.1107/S0365110X58001882SubbaraoE. C.ShiraneG.Nonstoichiometry and ferroelectric properties of PbNb2O6-type compounds1960326184618512-s2.0-0006014215SubbaraoE. C.X-ray study of phase transitions in ferroelectric PbNb2O6 and related materials196043943944210. 1111/j. 1151-2916. 1960.tb13694.xRayS.GünterE.Ritzhaupt-KleisslH. J.Manufacturing and characterization of piezoceramic lead metaniobate PbNb2O62000352462216224ChakrabortyK. R.SahuK. R.DeA.DeU.Structural characterization of orthorhombic and rhombohedral lead meta-niobate samples201012011021132-s2.0-7865087452410.1080/10584587.2010.503804SahuK. R.DeU.Thermal characterization of piezoelectric and non-piezoelectric lead meta-niobate20094901-2757710.1016/j.tca.2008.12.028GuerreroF.LeyetY.VenetM.de LosJ.EirasJ. A.Dielectric behavior of the PbNb2O6 ferroelectric ceramic in the frequency range of 20 Hz to 2 GHz20072713–15404140442-s2.0-3434734057610.1016/j.jeurceramsoc.2007.02.091BarbarS. K.RoyM.Synthesis, structural, electrical, and thermal studies of Pb1−xBaxNb2O6 (x = 0.0 and 0.4) ferroelectric ceramics20122012571017310.5402/2012/710173GriffinJ. W.PetersT. J.PosakonyG. J.ChienH. T.BondL. J.2009Richland, Wash, USAU.S. Department of Energy, Pacific Northwest National LaboratoryPNNL-18292KitykI. V.Makowska-JanusikM.FontanaM. D.AillerieM.FahmiA.Influence of non-stoichiometric defects on optical properties in LiNbO32001366577589VenetM.VendraminiA.ZabottoF. L.GarciaF. D.EirasJ. A.Piezoelectric properties of undoped and titanium or barium-doped lead metaniobate ceramics200525122443244610.1016/j.jeurceramsoc.2005.03.079LeeH. S.KimuraT.Effects of microstructure on the dielectric and piezoelectric properties of lead metaniobate1998811232283236SoejimaJ.SatoK.NagataK.Preparation and characteristics of ultrasonic transducers for high temperature using PbNb2O6200039530833085LiY. M.ChengL.GuX. Y.ZhangY. P.LiaoR. H.Piezoelectric and dielectric properties of PbNb2O6-based piezoelectric ceramics with high Curie temperature20081971–3170173MantasP. Q.Dielectric response of materials: extension to the Debye model199919122079208610. 1016/S0955-2219(98)00273-8ZhangS.RandallC. A.ShroutT. R.High Curie temperature piezocrystals in the BiScO3-PbTiO3 perovskite system20038315315031522-s2.0-024230444710.1063/1.1619207MacdonaldJ. R.JohnsonW. B.BarsoukovE.MacdonaldJ. R.20052ndHoboken, NJ, USAJohn Wiley & SonsSenS.ChoudharyR. N. P.Impedance studies of Sr modified BaZr0.5Ti0.95O3 ceramics2004872-32562632-s2.0-424305063510.1016/j.matchemphys.2004.03.005ShuklaA.ChoudharyR. N. P.ThakurA. K.Thermal, structural and complex impedance analysis of Mn4+ modified BaTiO3 electroceramic200970111401140710.1016/j.jpcs.2009.08.015UchinoK.NomuraS.Critical exponents of the dielectric constants in diffused phase transition crystals19824411556110.1016/j.jpcs.2009.08.015CrossL. E.Relaxor ferroelectrics198776124126710.1080/00150198708016945RaymondO.FontR.Suárez-AlmodovarN.PortellesJ.SiqueirosJ. M.Frequency-temperature response of ferroelectromagnetic Pb (Fe1/2Nb1/2)O3 ceramics obtained by different precursors—part I: structural and thermo-electrical characterization200597882-s2.0-2144445680710.1063/1.1870099084107BeheraB.ArajoE. B.ReisR. N.GuerraJ. D. S.AC conductivity and impedance properties of 0.65Pb(Mg1/3Nb2/3)O3-0.35 PbTiO3 ceramics2009200962-s2.0-7654912923910.1155/2009/361080361080ChoudharyR. N. P.PradhanD. K.BonillaG. E.KatiyarR. S.Effect of La-substitution on structural and dielectric properties of Bi(Sc1/2Fe1/2)O3 ceramics20074371-22202242-s2.0-3424724281810.1016/j.jallcom.2006.07.077RaoK. S.PrasadD. M.KrishnaP. M.LeeJ. H.Synthesis, electrical and electromechanical properties of tungsten-bronze ceramic oxide: Pb0.68K0.64Nb2O6200840320792087LuZ.BonnetJ. P.RavezJ.HagenmullerP.Correlation between low frequency dielectric dispersion (LFDD) and impedance relaxation in ferroelectric ceramic Pb2KNb4TaO151992573-423524410.1016/0167-2738(92)90153-GSahooP. S.PanigrahiA.PatriS. K.ChoudharyR. N. P.Impedance spectroscopy of Ba3Sr2DyTi3V7O30 ceramic20103321291342-s2.0-7795384000110.1007/s12034-010-0018-8RachnaS.GuptaS. M.BhattacharyyaS.Impedance analysis of Bi3.25La0.75Ti3O12 ferroelectric ceramic20087135996102-s2.0-5264914375110.1007/s12043-008-0135-8DeviS.JhaA. K.Dielectric and complex impedance studies of BaTi0.85W0.15O3+δ ferroelectric ceramics2010336683690GonzalezR. L.LeyetY.GuerreroF.Relaxation dynamics of the conductive process for PbNb2O6 ferroelectric ceramics in the frequency and time domain2007191313621810.1088/0953-8984/19/13/136218