About this Journal Submit a Manuscript Table of Contents
Journal of Ceramics
Volume 2013 (2013), Article ID 795827, 7 pages
http://dx.doi.org/10.1155/2013/795827
Research Article

Mean-Field Approach to Dielectric Relaxation in Giant Dielectric Constant Perovskite Ceramics

1College of Materials Science and Engineering, Shenzhen University, Shenzhen 518060, China
2Shenzhen Key Laboratory of Special Functional Materials, Shenzhen 518060, China
3Department of Applied Physics and Materials Research Center, The Hong Kong Polytechnic University, Hong Kong
4School of Materials Science and Engineering, Northwestern Polytechnical University, Xi’an 710072, China

Received 25 June 2013; Accepted 25 August 2013

Academic Editor: Matjaz Valant

Copyright © 2013 Shanming Ke et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The dielectric properties of CaCu3Ti4O3 (CCTO) and A2FeNbO6 (AFN, A = Ba, Sr, and Ca) giant dielectric constant ceramics were investigated in the frequency range from 1 Hz to 10 MHz. The relaxation properties can be perfectly described by a polaron model, indicating that the dielectric relaxation is intimately related to the hopping motion caused by localized charge carriers.

1. Introduction

As driven by the impetus of smaller and smaller feature size of devices in microelectronics, researchers are looking for the so-called high- materials. Perovskites such as lead zirconate titanate (PZT) usually possess high dielectric constant of a few hundreds at room temperature [1, 2] and are widely used as capacitive components. The perovskite-based oxide CaCu3Ti4O12 (CCTO, space group Im3) was reported to have a colossal dielectric constant (CDC) in the order of 104 which is nearly constant from 400 K down to 100 K but drops rapidly to less than 100 below 100 K [3]. A huge amount of work [37] has thereafter been carried out in attempts to understand the origin of the remarkable dielectric properties. It has been found that the temperature at which the step-like decrease in dielectric constant takes place strongly depends on the measuring frequency and roughly follows an Arrhenius behavior. Similar phenomenon has also been reported in a number of materials, such as A2FeBO6 (A = Ba, Sr, and Ca; B = Nb, Ta, etc.) [8], (LCMO) [9], Pr0.7Ca0.3MnO3 (PCMO) [10], LaCuLiO4 [11], LaSrNiO4 [12], TbMnO3 [13], and Li/Ti doped NiO [14]. The remarkable low temperature dielectric relaxation in many manganites, cuprates, and nickelates [912, 1517] has been attributed to localized hopping of polarons between lattice sites with a characteristic timescale. Furthermore, this relaxation behavior clearly displays a freezing temperature following a glass-like process, suggesting that an electronic glass state is realized. It seems plausible to think of the influence of the electric state on the dielectric response of a solid.

Polaronic relaxation usually involves a variable range hopping (VRH) or a nearest-neighbor hopping conduction process. Polarons as the main sources contributing to the conduction of CCTO ceramics at the low temperature relaxation region has been confirmed [18]. Tselév et al. [19] found a power-law frequency dependent dielectric response in CCTO epitaxial thin films suggesting that the dielectric contribution came from hopping carriers. Based on these results, the low temperature dielectric response in CCTO has been attributed to polaron relaxation related localized charge carriers [20]. It should be mentioned that Maxwell-Wagner (MW) origin, such as the internal barrier layer capacitance (IBLC) should not be excluded in the explanation of the colossal dielectric constants, especially at high temperatures [4, 6, 7, 21]. Although, the MW type mechanisms which were resulted from the surface layer [21], IBLC [4], and/or contact electrode [6] do have an influence on the dielectric constant values of CDC materials, they nearly have no influence on the relaxation process at low temperature.

In this paper, we present a systematic study on the low temperature dielectric relaxation process of CCTO and A2FeNbO6 (AFN, A = Ba, Sr, and Ca) systems. A polaron relaxation model was proposed as a possible explanation to the low temperature dielectric relaxation in CCTO and AFN ceramics. It should be mentioned that MW model and polaron model are both related to the moving of space charge carriers. The difference is that MW model emphasizes the interface, but polaron model considers the space electrons/ions exist in the whole bulk materials. Although the polaron concept has been used to explain the dipolar effects induced by charge carrier hopping motions inside the CCTO grains [20, 22], a clear theory formula to describe polaron relaxation is absence. Based on polaron theory, Jonscher’s law (universal dielectric response (UDR)) could be used to fit the frequency dependence of the dielectric permittivity. By using a mean-field approach, we will present a simplicity expression on the temperature dependence of the dielectric permittivity driven from polaron relaxation.

2. Experimental

Single phase CCTO and AFN (A = Ba, Sr, and Ca) ceramics were prepared through a conventional mixed oxide route, and the detailed processing parameters can be found elsewhere [23, 24]. The Wolframite method, using FeNbO4 as the -site precursor, has been used to synthesize AFN with no secondary iron-oxide phases.

X-ray diffraction (XRD) was conducted on sintered ceramics samples of CCTO and AFN. Data were collected on an automated diffractometer (X’Pert PRO MPD, Philips) with Cu  radiation. The XRD results confirm that the sintered ceramics are single phase. The fracture surfaces of the ceramic pellets were examined by scanning electron microscopy (SEM, JSM-6335F, and JEOL). Micro-Raman spectral measurements were performed on a JY HR800 Raman spectrometer under backscattering geometry. An Argon ion laser was used as the excitation source with an output power of 15 mw at 488 nm. For dielectric measurement, top and bottom electrodes were made by coating silver paint on both sides of the sintered disks and followed by a firing at 650°C for 20 minutes. The temperature and the frequency dependences of the dielectric permittivity of the samples were measured by a frequency response analyzer (Novocontrol Alpha-analyzer) over a broad frequency range (1 Hz–10 MHz).

3. Results and Discussion

3.1. Crystal Structure and Microstructure

Figure 1 shows the room temperature XRD pattern of sintered CCTO ceramics. All the diffraction peaks can be indexed according to a cubic cell of space group Im3 [3]. Standard group theory analysis predicts that the Raman active modes are distributed among the irreducible representations as [25]. Inset of Figure 1 shows the Raman spectrum at room temperature. The Raman lines at 448 and 512 cm−1 have exact symmetry, and a line of symmetry is clearly pronounced at 575 cm−1, which is assigned to the Ti–O–Ti antistretching mode of the oxygen octahedra [25]. The Raman spectrum also confirms a single cubic phase of our CCTO sample.

795827.fig.001
Figure 1: XRD pattern of CCTO ceramics. The inset shows the micro-Raman scattering spectra.

The XRD patterns of the Ba2FeNbO6 (BFN), Sr2FeNbO6 (SFN), and Ca2FeNbO6 (CFN) are given in Figure 2. The diffraction peaks of BFN samples can be indexed according to a monoclinic structure [26], while SFN and CFN samples show an orthorhombic structure with the space group Pnma(62) [27, 28]. Figure 3 displays typical SEM photographs of the fracture surface of CCTO and SFN ceramics. The average grain size of CCTO is found to be in the range of 2~8 μm, while SFN exhibits a morphology with the bimodal distribution of grain size. One has the grain size of several micrometers, and the other has the smaller grain size in the order of several hundreds of nanometers. It should be noted that BFN and Ca2FeNbO6 (CFN) display a similar morphology to that of SFN, while BFN has a slightly larger average grain size than SFN and CFN.

795827.fig.002
Figure 2: XRD patterns of BFN, SFN, and CFN ceramics.
fig3
Figure 3: SEM images of the fracture surface of (a) CCTO and (b) SFN ceramics. The arrows denote nanograins in SFN.
3.2. Dielectric Behavior

Figure 4 shows the frequency dependences of the real [] and imaginary [] parts of the dielectric constant of CCTO, BFN, CFN, and SFN under different temperatures. The temperatures were chosen because they reveal the typical relaxation features of the mentioned ceramics. The high frequency region of CCTO and BFN shows the presence of a Debye-like dipolar relaxation process, and SFN and CFN show a generalized dipolar loss peak possessing a more serious departure from Debye, as summarized by Jonscher [29]. The peak frequency of decreases with decreasing temperatures for all the samples, indicating the freezing in dipolar moment. The width of the peak of AFN (A = Ba, Sr, Ca) is much larger than that of CCTO (as shown in Figure 4), which implies a broader distribution of relaxation time. Above the peak frequencies, the high-frequency branches of CCTO and AFN follow the fractional power law [29]: where is the angular frequency () and is the exponential constant. The calculated exponent from Figure 4 is 0.11, 0.12, 0.25, and 0.38 for CCTO, BFN, CFN, and SFN, respectively. It should be noted that is nearly independent of temperature for all these CCTO and AFN samples. Larger means a much lower loss () and flatter frequency dependence than that in a Debye relaxation.

fig4
Figure 4: Frequency dependence of the real and imaginary parts of dielectric permittivity for (a) CCTO, (b) BFN, (c) CFN, and (d) SFN. The straight gray line in (a) is a guide to the eye.

The peak frequency of the above relaxation process obeys the following Arrhenius law: where is a prefactor and is the activation energy. The Arrhenius fits for CCTO, BFN, SFN, and CFN ceramics are summarized in Figure 5, where = 0.1, 0.18, 0.23, and 0.34 eV for CCTO, BFN, SFN, and CFN, respectively. The obtained values are very close to the previous reported ones for the same materials [47, 30].

795827.fig.005
Figure 5: Arrhenius fit of the frequency dependence of the dynamic transition temperatures for the dielectric relaxation in CCTO, BFN, SFN, and CFN.

Although, the Maxwell-Wagner type mechanism is always responsible for the dielectric permittivity of ceramics as an extrinsic source due to the existence of grain boundaries. This model cannot explain why the one with very similar microstructure and electrical resistivity (BFN, CFN, and SFN) shows so different relaxation regions and activation energies. In addition, in terms of the Maxwell-Wagner model, it is very difficult to explain the similar dielectric relaxation behavior appearing in the same temperature region for single crystals, polycrystalline ceramics, and epitaxial thin films whose microstructure is certainly different. In Figure 6, the imaginary part of dielectric permittivity of CCTO and BFN is shown as a function of the real part , namely, the Cole-Cole plot. From Figure 6, all the data at different temperatures collapse on an arc described by the Cole-Cole equation, showing universal scaling behavior in the dielectric response. We note that in the lower frequency region there is a relaxation relating to the extrinsic effect. This relaxation is probably due to a Schottky-type barrier existing in this ceramic [31].

fig6
Figure 6: The imaginary part of dielectric permittivity of CCTO and BFN is plotted against the real part .
3.3. Mean-Field Approach

In a previous study, It has been confirmed that localized charge carriers are the main sources contributing to the conduction of CCTO [18] and AFN [32] ceramics at low temperature relaxation region. Under an external electric field, the hopping motion of localized carriers like polarons gives rise to dipolar effect and sizable polarization contribution to the dielectric permittivity. Therefore a quantitative description of the polaron contribution to the dielectric response is highly demanding.

Ramirez et al. [33] have proposed a defect model as a possible explanation of the low temperature dielectric response of CCTO. In this model, isolated defects such as Cu vacancies produce a local disruption of the ideal cubic structure and then the defective regions relax between alternate equivalent configurations. However, high resolution experiments [33] have revealed that Cu vacancies most probably exist in the grain boundary regions, and thus they cannot be shown in the bulk response. Therefore, polarons are mainly responsible for the bulk dielectric relaxation. Essentially, it is reasonable to presume that polarons are isolated; then the cooperative effects are absent at low densities of polarons. The relaxation of polarons occurs at the rate determined by the energy barrier between alternate equivalent configurations and the temperature: where depends on the effective mass of the defect and is in the order of 1010–1012 Hz. The local polarizability of a polaron is where is the static polarizability of a polaron. It should be noted that (3) and (4) could be obtained directly from [33]. As suggested by Ramirez et al. [33], a mean-field approximation can be adopted and the dielectric response can be calculated by the Clausius-Mossotti equation: where is the density of the polarons, and the details can be found in [34]. It is worth noting that the Clausius-Mossotti equation works best for gases and is only approximately true for liquids or solids, particularly for polar materials, or the dielectric constant is large. In addition, the limitations of the approximation in (5), which does not consider the variations around each polaron properly, should also be kept in mind.

Let ; then we get The real and imaginary parts of the dielectric permittivity are

The measured temperature dependent and (see Figures 7, 8, and 9) for CCTO, BFN, and CFN are similar to those reported in previous work [37, 30]. The real part dielectric permittivity reveals a pronounced step-like increase accompanied by a loss peak. The dielectric permittivity of CCTO, BFN, and CFN is fitted according to (7) by a least-square method. For CCTO measured at 100 kHz, the best fitted result was obtained with , , and  K. It is very interesting to note that the obtained barrier energy (obtained from ) is 0.08 eV, a value quite close to the characteristic activation energy obtained from (2) (Section 3.2,  eV for CCTO). The same fitting parameters were used to calculate the temperature dependence of and at different frequencies. The solid lines in Figure 7 are calculated from the previous model, which reproduce the experimental data quite well. We have also performed the least square fitting process on BFN and CFN. The best result was obtained with , , and  K for BFN, and , , and  K for CFN, respectively. The calculated barrier energy is 0.17 eV for BFN, which is also quite close to the result obtained from (2) (see Figure 5). However, the calculated barrier energy for CFN is only 0.24 eV, which is much smaller than the one obtained from (2) (0.34 eV). Presumably, this may be attributed to an overlap between the polaron relaxation and a high temperature relaxation process as shown in Figure 9.

795827.fig.007
Figure 7: The temperature dependence of dielectric constant of CCTO. The solid lines are the best fits from polaron model (see text).
795827.fig.008
Figure 8: The dielectric constant of BFN is shown against temperature. The solid lines are the beat fits of polaron model.
795827.fig.009
Figure 9: The dielectric constant of CFN is shown as a function of temperature. The solid lines are the beat fits of polaron model.

4. Conclusions

In summary, we have prepared single phase CCTO and AFN (BFN, CFN, and SFN) ceramics by conventional solid state reaction. Based on the dielectric study, it can be concluded that the low temperature step-like dielectric relaxation in CCTO and AFN can be well explained by a polaron model involving a dipolar-type relaxation. By using a mean-field approach, we get a dielectric function which describes the dielectric relaxation in CCTO and BFN quite well.

Acknowledgments

This work was supported by the Hong Kong Polytechnic University (Projects nos. A-PL54 and 1-BD08), the National Natural Science Foundations of China (Nos. 51302172 and 51272161), and Shenzhen Innovation and Technology Commission under the strategic emerging industries development project (Contract no. ZDSY20120612094418467).

References

  1. L. Kerkache, A. Layadi, E. Dogheche, and D. Remiens, “Structural, ferroelectric and dielectric properties of In2O3:Sn (ITO) on PbZr0.53Ti0.47O3 (PZT)/Pt and annealing effect,” Journal of Alloys and Compounds, vol. 509, no. 20, pp. 6072–6076, 2011. View at Publisher · View at Google Scholar · View at Scopus
  2. N. Vittayakorn, G. Rujijanagul, and D. P. Cann, “Investigation of the influence of thermal treatment on the morphologies, dielectric and ferroelectric properties of PZT-based ceramics,” Journal of Alloys and Compounds, vol. 440, no. 1-2, pp. 259–264, 2007. View at Publisher · View at Google Scholar · View at Scopus
  3. M. A. Subramanian, D. Li, N. Duan, B. A. Reisner, and A. W. Sleight, “High dielectric constant in ACu3Ti4O12 and ACu3Ti3FeO12 phases,” Journal of Solid State Chemistry, vol. 151, no. 2, pp. 323–325, 2000. View at Publisher · View at Google Scholar · View at Scopus
  4. D. C. Sinclair, T. B. Adams, F. D. Morrison, and A. R. West, “CaCu3Ti4O12: one-step internal barrier layer capacitor,” Applied Physics Letters, vol. 80, no. 12, pp. 2153–2155, 2002. View at Publisher · View at Google Scholar · View at Scopus
  5. A. P. Ramirez, M. A. Subramanian, M. Gardel et al., “Giant dielectric constant response in a copper-titanate,” Solid State Communications, vol. 115, no. 5, pp. 217–220, 2000. View at Publisher · View at Google Scholar · View at Scopus
  6. P. Lunkenheimer, R. Fichtl, S. G. Ebbinghaus, and A. Loidl, “Nonintrinsic origin of the colossal dielectric constants in CaCu3Ti4O12,” Physical Review B, vol. 70, no. 17, Article ID 172102, 4 pages, 2004. View at Publisher · View at Google Scholar · View at Scopus
  7. C. C. Homes, T. Vogt, S. M. Shapiro, S. Wakimoto, and A. P. Ramirez, “Optical response of high-dielectric-constant perovskite-related oxide,” Science, vol. 293, no. 5530, pp. 673–676, 2001. View at Scopus
  8. I. P. Raevski, S. A. Prosandeev, A. S. Bogatin, M. A. Malitskaya, and L. Jastrabik, “High dielectric permittivity in AFe1/2B1/2O3 nonferroelectric perovskite ceramics (A = Ba, Sr, Ca; B = Nb, Ta, Sb),” Journal of Applied Physics, vol. 93, no. 7, pp. 4130–4136, 2003. View at Publisher · View at Google Scholar · View at Scopus
  9. K. P. Neupane, J. L. Cohn, H. Terashita, and J. J. Neumeier, “Doping dependence of polaron hopping energies in La1-xCaxMnO3 (0x0.15),” Physical Review B, vol. 74, no. 14, Article ID 144428, 2006. View at Publisher · View at Google Scholar · View at Scopus
  10. R. S. Freitas, J. F. Mitchell, and P. Schiffer, “Magnetodielectric consequences of phase separation in the colossal magnetoresistance manganite Pr0.7Ca0.3MnO3,” Physical Review B, vol. 72, no. 14, Article ID 144429, 2005. View at Publisher · View at Google Scholar
  11. T. Park, Z. Nussinov, K. R. A. Hazzard et al., “Novel dielectric anomaly in the hole-doped La2Cu1-xLixO4 and La2-xSrxNiO4 insulators: signature of an electronic glassy state,” Physical Review Letters, vol. 94, no. 1, Article ID 017002, 2005. View at Publisher · View at Google Scholar · View at Scopus
  12. J. Rivas, B. Rivas-murias, A. Fondado, J. Mira, and M. A. Señarís-Rodríguez, “Dielectric response of the charge-ordered two-dimensional nickelate La1.5Sr0.5NiO4,” Applied Physics Letters, vol. 85, no. 25, Article ID 6224, 3 pages, 2004. View at Publisher · View at Google Scholar
  13. C. C. Wang, Y. M. Cui, and L. W. Zhang, “Dielectric properties of TbMnO3 ceramics,” Applied Physics Letters, vol. 90, no. 1, Article ID 012904, 3 pages, 2007. View at Publisher · View at Google Scholar
  14. J. B. Wu, C. W. Nan, Y. H. Lin, and Y. Deng, “Giant dielectric permittivity observed in Li and Ti doped NiO,” Physical Review Letters, vol. 89, no. 21, Article ID 217601, 2002. View at Publisher · View at Google Scholar
  15. J. L. Cohn, M. Peterca, and J. J. Neumeier, “Low-temperature permittivity of insulating perovskite manganites,” Physical Review B, vol. 70, no. 21, Article ID 214433, pp. 1–6, 2004. View at Publisher · View at Google Scholar · View at Scopus
  16. A. Seeger, P. Lunkenheimer, J. Hemberger et al., “Charge carrier localization in La1-xSrxMnO3 investigated by ac conductivity measurements,” Journal of Physics Condensed Matter, vol. 11, no. 16, pp. 3273–3290, 1999. View at Publisher · View at Google Scholar · View at Scopus
  17. J. L. Cohn, M. Peterca, and J. J. Neumeier, “Giant dielectric permittivity of electron-doped manganite thin films, Ca1-xLaxO3 (0x0.03),” Journal of Applied Physics, vol. 97, no. 3, Article ID 034102, 2005. View at Publisher · View at Google Scholar · View at Scopus
  18. L. Zhang and Z. J. Tang, “Polaron relaxation and variable-range-hopping conductivity in the giant-dielectric-constant material CaCu3Ti4O12,” Physical Review B, vol. 70, no. 17, Article ID 174306, 6 pages, 2004. View at Publisher · View at Google Scholar
  19. A. Tselév, C. M. Brooks, S. M. Anlage et al., “Evidence for power-law frequency dependence of intrinsic dielectric response in the CaCu3Ti4O12,” Physical Review B, vol. 70, no. 14, Article ID 144101, 2004. View at Publisher · View at Google Scholar · View at Scopus
  20. C. C. Wang and L. W. Zhang, “Polaron relaxation related to localized charge carriers in CaCu3Ti4O12,” Applied Physics Letters, vol. 90, no. 14, Article ID 142905, 3 pages, 2007. View at Publisher · View at Google Scholar
  21. L. Liu, C. C. Wang, X. H. Sun, G. J. Wang, C. M. Lei, and T. Li, “Oxygen-vacancy-related relaxations of Sr3CuNb2O9 at high temperatures,” Journal of Alloys and Compounds, vol. 552, pp. 279–282, 2013. View at Publisher · View at Google Scholar
  22. L. Fang, M. R. Shen, F. Zheng, Z. Li, and J. Yang, “Dielectric responses and multirelaxation behaviors of pure and doped CaCu3Ti4O12 ceramics,” Journal of Applied Physics 104, vol. 104, no. 6, Article ID 064110, 8 pages, 2008. View at Publisher · View at Google Scholar
  23. S. M. Ke, H. T. Huang, and H. Q. Fan, “Relaxor behavior in CaCu3Ti4O12 ceramics,” Applied Physics Letters, vol. 89, no. 18, Article ID 182904, 3 pages, 2006. View at Publisher · View at Google Scholar
  24. S. Ke, H. Huang, H. Fan, H. L. W. Chan, and L. M. Zhou, “Colossal dielectric response in barium iron niobate ceramics obtained by different precursors,” Ceramics International, vol. 34, no. 4, pp. 1059–1062, 2008. View at Publisher · View at Google Scholar · View at Scopus
  25. N. Kolev, C. L. Chen, M. Gospodinov et al., “Raman spectroscopy of CaRuO3,” Physical Review B, vol. 66, no. 1, Article ID 014101, pp. 141011–141014, 2002. View at Scopus
  26. S. Saha and T. P. Sinha, “Low-temperature scaling behavior of BaFe0.5Nb0.5O3,” Physical Review B, vol. 65, no. 13, Article ID 134103, 7 pages, 2002. View at Publisher · View at Google Scholar
  27. K. Tezuka, K. Henmi, Y. Hinatsu, and N. M. Masaki, “Magnetic susceptibilities and Mossbauer spectra of perovskites A2FeNbO6 (A = Sr, Ba),” Journal of Solid State Chemistry, vol. 154, no. 2, pp. 591–597, 2000. View at Publisher · View at Google Scholar · View at Scopus
  28. Y. Y. Liu, X. M. Chen, X. Q. Liu, and L. Li, “Dielectric relaxations in Ca(Fe1/2Nb1/2)O3 complex perovskite ceramics,” Applied Physics Letters, vol. 90, no. 26, Article ID 262904, 3 pages, 2007. View at Publisher · View at Google Scholar
  29. A. K. Jonscher, “The universal dielectric response and its physical significance,” IEEE Transactions on Electrical Insulation, vol. 27, no. 2, pp. 407–423, 1992. View at Scopus
  30. Z. Wang, X. M. Chen, L. Ni, and X. Q. Liu, “Dielectric abnormities of complex perovskite Ba(Fe1/2Nb1/2)O3 ceramics over broad temperature and frequency range,” Applied Physics Letters, vol. 90, no. 2, Article ID 022904, 3 pages, 2007. View at Publisher · View at Google Scholar
  31. P. R. Bueno, M. A. Ramirez, J. A. Varela, and E. Longo, “Dielectric spectroscopy analysis of CaCu3Ti4O12 polycrystalline systems,” Applied Physics Letters, vol. 89, no. 19, Article ID 191117, 3 pages, 2006. View at Publisher · View at Google Scholar
  32. S. Ke, H. Fan, and H. Huang, “Dielectric relaxation in A2FeNbO6 (A = Ba, Sr, and Ca) perovskite ceramics,” Journal of Electroceramics, vol. 22, no. 1–3, pp. 252–256, 2009. View at Publisher · View at Google Scholar · View at Scopus
  33. A. P. Ramirez, G. Lawes, V. Butko, M. A. Subramanian, and C. M. Varma, “Colossal dielectric constants in braced lattices with defects,” http://arxiv.org/abs/cond-mat/0209498.
  34. B. I. Halperin and C. M. Varma, “Defects and the central peak near structural phase transitions,” Physical Review B, vol. 14, no. 9, pp. 4030–4044, 1976. View at Publisher · View at Google Scholar · View at Scopus