Advances in Condensed Matter Physics

Volume 2012, Article ID 259712, 5 pages

http://dx.doi.org/10.1155/2012/259712

## Calculation of the Spontaneous Polarization and the Dielectric Constant as a Function of Temperature for

Department of Physics, Middle East Technical University, 06531 Ankara, Turkey

Received 13 March 2012; Revised 23 August 2012; Accepted 28 August 2012

Academic Editor: Biljana Stojanovic

Copyright © 2012 Hamit Yurtseven and Sema Şen. 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 temperature dependence of the spontaneous polarization *P* is calculated in the ferroelectric phase of KH_{2}PO_{4} (KDP) at atmospheric pressure (*T*_{C} = 122 K). Also, the dielectric constant *ε* is calculated at various temperatures in the paraelectric phase of KDP at atmospheric pressure. For this calculation of *P* and *ε*, by fitting the observed Raman frequencies of the soft mode, the microscopic parameters of the pure tunnelling model are obtained. In this model, the proton-lattice interaction is not considered and the collective proton mode is identified with the soft-mode response of the system. Our calculations show that the spontaneous polarization decreases continuously in the ferroelectric phase as approaching the transition temperature *T*_{C}. Also, the dielectric constant decreases with increasing temperature and it diverges in the vicinity of the transition temperature (*T*_{C} = 122 K) for KDP according to the Curie-Weiss law.

#### 1. Introduction

Potassium dihydrogen phosphate (KH_{2}PO_{4} or KDP) undergoes a phase transition at = 122 K at atmospheric pressure [1]. In the paraelectric phase it has the space group , whereas in the ferroelectric phase its space group is *Fdd2 *. The paraelectric phase has the tetragonal structure with the four K-PO_{4} groups in the unit cell and the tetrahedral PO_{4} groups are connected by O–HO bonds. In this phase the protons are randomly distributed between the two equivalent minima of the double-well potential. Below in the ferroelectric phase, the protons occupy in one minima of this double-well potential so that the transition occurs by this ordering of the protons. Thus, the spontaneous polarization, which grows as the temperature decreases below in the ferroelectric phase, is parallel to the *c*-axis and it is mainly due to the motion of the protons in KDP. This occurs when the protons tunnel from one minima of the double-well potential to the other, which is considered in the pure tunnelling model [1–5]. However, it has been considered that the proton motion couples with an optic mode of the lattice in KDP, as studied by Kobayashi [6] in the coupled proton-optic-mode model. In the pure tunnelling model, it has been shown that the collective proton-proton tunnelling mode represents a soft mode [4], whereas in the coupled-mode model the soft mode is the lower frequency branch of the coupled system [6]. The soft mode has been observed experimentally in the Raman spectra of KDP [7–9]. At high pressures, the soft mode has also been observed experimentally in the ferroelectric phase of KDP [10]. The wavevector dependence of the soft mode has been detected at the pressures above 40 kbar [11]. Recently, we have calculated the Raman frequency and the damping constant in KH_{2}PO_{4} [12].

In this work, we calculate the spontaneous polarization from the observed Raman frequencies of the collective proton mode which behaves as a soft-mode response of the pure tunnelling model. By extrapolating the Raman frequencies of this mode to zero pressure, which were measured at various temperatures as a function of pressure [10] in the ferroelectric phase of KDP, the temperature dependence of the spontaneous polarization is calculated. This calculation is performed in the absence of the proton-lattice interaction for KDP. We also calculate here the temperature dependence of the dielectric constant in the paraelectric phase of KDP using the pure tunnelling model. For this calculation, the proton-lattice interaction is also not considered and the critical behaviour of the dielectric constant is predicted close to in the paraelectric phase of KDP.

Below, in Section 2, we give a theoretical outline for our calculations of the spontaneous polarization and the dielectric constant. In Section 3 we give our results. Sections 4 and 5 give our discussion and conclusions, respectively.

#### 2. Theory

Using the microscopic theory, the mechanism of the phase transition in KDP has been explained mainly by the two models, namely, the tunnelling model of the proton and the coupled proton-optic-mode model. The ordering of the protons, which tunnels between the two minima of a double-well potential on the O–HO bonds, as studied by Slater [1], causes collective proton excitations. This collective proton-proton tunnelling mode has been considered as a soft mode when the protons do not interact with the lattice [4]. Due to the fact that the spontaneous polarization (order parameter) does not have large values in this pure tunnelling model, the collective proton tunnelling mode interacts with the transverse optic mode (proton-lattice interaction) in the coupled proton-optic-mode model, as treated by Kobayashi [6]. In a previous study, the soft mode and coupled modes in the paraelectric and ferroelectric phases of KDP have been studied experimentally and the experimental results have been analyzed using both tunnelling and coupled proton-optic mode models [10]. In this study, we calculate the temperature dependence of the spontaneous polarization (order parameter) and the dielectric constant within the pure tunnelling model using the observed Raman frequencies [10] of the soft mode in KDP.

##### 2.1. Tunnelling Model

In the tunnelling model, the proton can occupy two equivalent potential minima assigned by the pseudospin values of ±1/2. Protons can tunnel between the two minima and the energy levels of the double-well potential are separated by Δ*E* = 2*ħ*Ω where Ω is the tunnelling frequency.

The proton Hamiltonian can be assumed to consist of a single-particle contribution and a contribution due to proton-proton interaction [3], expressed as
where (*α* = *x*, *y*, *z*) represents the *α* component of the pseudospin and is the interaction parameter between the nearest-neighbour protons. The proton Hamiltonian can be solved in the mean-field approximation above the transition temperature in terms of the collective proton excitation frequency Ω_{0}, as given in a previous study [10],
In (2) *J* is the Fourier transform of as the interaction parameter between the nearest-neighbour protons and is defined as
As approaching the transition temperature, when there is no proton-lattice coupling, the collective excitation frequency can be assumed to disappear (Ω_{0}→0), which results in [10]

using (2) and (3) with the transition temperature →.

According to the tunnelling model, the dielectric constant can be obtained as a function of temperature [13]
where is the lattice dielectric constant independent of the temperature. The temperature dependence of the dielectric constant is governed by the second term in (5), which is the soft-mode contribution to *ε*(*T*). In (5), *N* is the number of dipoles per unit volume and *μ* is the electric dipole moment.

The temperature dependence of the dielectric constant *ε*(*T*) near the transition temperature in the paraelectric phase can be obtained by expanding (5) about , which gives the Curie-Weiss law expressed as
where the Curie constant is defined as
In relation to the dielectric constant *ε*(), the temperature dependence of the order parameter (spontaneous polarization) can also be obtained in the tunnelling model, as given in the previous study [10],
where represents the pseudospins directed along the *z*-axis. Since the spontaneous polarization (order parameter) is defined in the ferroelectric phase only, in the paraelectric phase . At very low temperatures , the saturated spontaneous polarization is obtained as
for the pure tunnelling model [10]. In order to obtain the temperature dependence of the spontaneous polarization (order parameter), we substitute (3) into (8) since near the transition temperature can be replaced by for the pseudospins orientated along the *z*-axis. This then gives the temperature dependence of the spontaneous polarization as

The temperature dependence of the collective proton excitation frequency Ω_{0} can also be obtained by identifying it with the soft-mode response frequency *w*_ [4, 10]. Experimentally, the pressure dependence of *w*_ at constant temperatures and its temperature dependence at 6.54 kbar have been observed for KDP [10]. On the assumption that there is no proton-lattice interaction, by substituting (3) and (4) into (2) the temperature dependence of Ω_{0} or the soft-mode response frequency *w*_ can be expressed as
where is the transition temperature, as before.

#### 3. Calculations and Results

The temperature dependence of the spontaneous polarization and the dielectric constant for KDP were calculated here. For this calculation, the interaction parameter *J* and the tunnelling frequency Ω were determined by fitting (11) to the observed data for the soft-mode response frequency *w*_ at various temperatures at zero pressure [10]. The extrapolated values of the observed *w*_ for various temperatures at zero pressure as given in Table 1 for the ferroelectric phase of KDP [10], were used to determine the fitting parameters Ω and in (11). The fitting parameters Ω and determined are given in Table 2. In order to determine the interaction parameter *J* between the nearest-neighbour protons at the transition temperature, Ω = 54.5 cm^{−1} and K at zero pressure were used in (4) for KDP. The Curie constant was also determined using the *J* value at and the Ω value according to (7) which is given as a ratio in Table 2. Figure 1 gives the ratio of the spontaneous polarization with respect to the saturated polarization which was calculated according to (9) and (10), where cm^{−1} was used, as a function of temperature below the transition temperature ( K) in the ferroelectric phase of KDP. Finally, the dielectric constant was calculated as a function of temperature above the transition temperature K in the paraelectric phase according to (6), as plotted in Figure 2. Here, the saturation value of the spontaneous polarization was calculated as a ratio of using the Ω and *J* values according to (9), as given in Table 2.

#### 4. Discussion

Using the pure tunnelling model, the temperature dependence of the spontaneous polarization in the ferroelectric phase (Figure 1) and of the dielectric function in the paraelectric phase (Figure 2) at zero pressure was calculated here. For this calculation, the observed [10] Raman frequencies of the soft-mode, which decrease with increasing temperature at zero pressure (Table 1), were used.

The spontaneous polarization decreases continuously as the temperature is increased towards the transition temperature at zero pressure (Figure 1). This is the soft-mode behaviour where the frequency decreases abruptly as the temperature approaches . In fact, it has been obtained that for the pure tunnelling model the soft mode frequency is directly related to the spontaneous polarization according to the relation
in the ferroelectric phase of KDP [13, 14]. This indicates that the mechanism of the phase transition in KDP which is associated with the soft-mode frequency *w*_, is due to the ordering of the protons in one minima of the double-well potential in this ferroelectric material.

The soft-mode response frequency *w*_ as it was first observed experimentally in KDP [7], is a heavily damped response centered at *w* = 0 in the Raman spectrum, as also pointed out previously [10]. The collective proton-proton tunnelling mode Ω_{0} was taken here as the soft-mode response frequency *w*_ since displays a soft mode character in the absence of proton-lattice interactions [4, 10]. The temperature dependence of the soft mode frequency *w*_ can be described by
where *T*_ is different from the in the presence of the optic-acoustic mode interaction [10]. It has been concluded [10] that the proton motion is strongly coupled to an optic mode of the lattice so that when the protons order, the proton-lattice interaction leads to the distortion of the K-PO_{4} sublattices to produce the spontaneous polarization in the modified tunnelling model or the coupled proton-optic mode model [6]. In regard to the proton-proton tunneling model, (13) can also be used for the proton-proton tunnelling model in the case that the proton tunnelling frequency Ω is taken as the temperature dependent.

The soft-mode response also displays a temperature dependence of the Raman frequency *w*_ near the transition at at a constant pressure of 6.54 kbar , as observed experimentally in KDP [10]. It has been indicated that the difference between and is primarily due to the optic-acoustic-mode interaction in KDP. It has also been indicated that this transition occurs at ~86 K when *w*_ vanishes for the coupled proton-optic-mode model, whereas for the pure tunnelling model in the absence of proton-lattice interactions the proton system undergoes a transition at ~58 K [10]. So, the temperature difference of nearly 30 K emerges due to proton-lattice interactions at 6.54 kbar in KDP. At atmospheric pressure, our fitting parameter of (Table 2) is also about 30 K lower than K determined previously [9]. In fact, the microscopic parameter , which we obtained by fitting (11) to the observed Raman frequencies of the soft mode *w*_ (or collective proton excitation frequency ), was not well determined since a large ratio of was obtained in our analysis. This gives rise to tanh which means that the fit is relatively insensitive to , as also pointed out previously for the analysis of KDP at 6.54 kbar [10]. So, the macroscopic properties of the ferroelectric materials are not well described by the microscopic parameters of the pure tunnelling model and also of the coupled proton-optic mode model, which both use the Raman data [10]. On this basis, in (4) we used the value of = 122 K instead of the fitted parameter in the pure tunnelling model, which gave the ratio of = 0.64 or tanh = 0.57 for KDP. Thus, the interaction parameter was determined as *J *= 384.6 cm^{−1} at zero pressure in this study, which can be compared with the value of 345.1 cm^{−1} at 6.54 kbar [10].

The temperature dependence of the spontaneous polarization was calculated here in the absence of the proton-lattice interactions for KDP, as stated previously. Since the spontaneous polarization does not have large values in the pure tunnelling model, the collective proton tunnelling mode can interact with the transverse optic mode (proton-lattice interaction) in the coupled proton-optic mode model introduced by Kobayashi [6], as pointed out previously [10]. The amplitude of the soft-optic mode is maximum at the phase transition temperature and there is no coupling at . Below , proton-lattice interaction can be considered to enhance the temperature dependence of the spontaneous polarization in KDP. However, the coupled mode features can disappear below in the ferroelectric phase in KDP, was also observed experimentally at approximately 30 K below in BaTiO_{3} [15]. Thus, in the absence of the proton-lattice interaction, for the temperature dependence of the soft-mode frequency *w*_ (or collective proton excitation frequency Ω_{0}) can be described below in the ferroelectric phase of KDP.

In regard to the collective proton excitation frequency Ω_{0} which was taken as the soft mode response frequency *w*_ for KDP as studied here, Ω_{0} can also be considered as a soft mode for those ferroelectrics exhibiting displacive phase transition such as BaTiO_{3}. Soft mode in BaTiO_{3} which consists of the vibration of Ti ions against oxygen [16] is the transverse optic mode of long wavelength at *k* = 0. Its frequency goes to zero ( → 0) at the critical temperature and the lattice displacements associated with this mode cause the crystal distorting to the ferroelectric phase [17]. Very recently, we have investigated the temperature dependence of this soft-optic mode in ferroelectric barium titanate [18].

Using the pure tunnelling model, we also calculated the temperature dependence of the dielectric constant in the paraelectric phase of KDP, as shown in Figure 2. It decreases as the temperature increases below and it diverges in the vicinity of the transition temperature ( = 122 K) according to Curie-Weiss law (see (6)). This diverging behaviour of the dielectric constant is due to the soft-mode response frequency *w*_ or the collective proton-tunnelling frequency Ω_{0} which vanishes as in the pure tunnelling model of KDP. Thus, the temperature or pressure dependence of the soft-mode response frequency from the Raman spectra determines the temperature or pressure dependence of the dielectric properties of those ferroelectric materials such as KDP. Also, from the pressure dependence of the transition temperature , the phase diagram calculated using proton- or-deuteron-tunnelling model [19] can be tested by the experimental Raman spectra for KDP and DKDP. Alternatively, a temperature-uniaxial pressure phase diagram of KH_{2}PO_{4} which was obtained experimentally [20] can be calculated using the proton tunnelling model.

#### 5. Conclusions

The spontaneous polarization *P* and the dielectric constant *ε* were calculated at various temperatures in the ferroelectric and paraelectric phases, respectively, at atmospheric pressure in KDP. The observed soft-mode frequencies from the Raman data were used to evaluate the temperature dependence of *P* and *ε* by the pure tunnelling model close to the transition temperature ( = 122 K) in KDP.

The critical behaviour of the spontaneous polarization *P* and the dielectric constant *ε* was predicted in both the ferroelectric and paraelectric phases and they can be compared with the experimental measurements. The critical behaviour of the soft-mode response frequency *w*_, spontaneous polarization *P*, and the dielectric constant *ε* can also be investigated at higher pressures in KDP.

#### References

- J. C. Slater, “Theory of the transition in KH
_{2}PO_{4},”*The Journal of Chemical Physics*, vol. 9, no. 1, pp. 16–33, 1941. View at Google Scholar · View at Scopus - R. Blinc, “On the isotopic effects in the ferroelectric behaviour of crystals with short hydrogen bonds,”
*Journal of Physics and Chemistry of Solids*, vol. 13, no. 3-4, pp. 204–211, 1960. View at Google Scholar · View at Scopus - P. G. de Gennes, “Collective motions of hydrogen bonds,”
*Solid State Communications*, vol. 1, no. 6, pp. 132–137, 1963. View at Google Scholar · View at Scopus - R. Brout, K. A. Müller, and H. Thomas, “Tunnelling and collective excitations in a microscopic model of ferroelectricity,”
*Solid State Communications*, vol. 4, no. 10, pp. 507–510, 1966. View at Google Scholar · View at Scopus - M. Tokunaga and T. Matsubara, “Theory of ferroelectric phase transition in KH
_{2}PO_{4}type crystals I,”*Progress of Theoretical Physics*, vol. 35, pp. 581–599, 1966. View at Google Scholar - K. K. Kobayashi, “Dynamical theory of the phase transition in KH
_{2}PO_{4}-type ferroelectric crystals,”*Journal of the Physical Society of Japan*, vol. 24, no. 3, pp. 497–508, 1968. View at Google Scholar · View at Scopus - I. P. Kaminov and T. C. Damen, “Temperature dependence of the ferroelectric mode in KH
_{2}PO_{4},”*Physical Review Letters*, vol. 20, no. 20, pp. 1105–1108, 1968. View at Publisher · View at Google Scholar · View at Scopus - J. F. Scott and C. M. Wilson, “Dynamical theory of the phase transition in KH
_{2}PO_{4}-type ferroelectric crystals,”*Solid State Communications*, vol. 10, pp. 597–600, 1972. View at Google Scholar - C. Y. She, T. W. Broberg, L. S. Wall, and D. F. Edwards, “Effect of proton-phonon coupling on the ferroelectric mode in KH
_{2}PO_{4},”*Physical Review B*, vol. 6, no. 5, pp. 1847–1850, 1972. View at Publisher · View at Google Scholar · View at Scopus - P. S. Peercy, “Measurement of the “soft” mode and coupled modes in the paraelectric and ferroelectric phases of KH
_{2}PO_{4}at high pressure,”*Physical Review B*, vol. 12, no. 7, pp. 2725–2740, 1975. View at Publisher · View at Google Scholar · View at Scopus - T. Yagi, A. Sakai, and M. Arima, “Raman scattering study of KDP under high pressure with a diamond anvill cell,”
*Phase Transitions*, vol. 28, pp. 63–78, 1990. View at Google Scholar - H. Karacali, A. Kiraci, and H. Yurtseven, “Calculation of the raman frequency and the damping constant of a coupled mode in the ferroelectric and paraelectric phases in KH
_{2}PO_{4},”*Physica Status Solidi B*, vol. 247, no. 4, pp. 927–936, 2010. View at Publisher · View at Google Scholar · View at Scopus - R. Blinc and B. Žekš, “Dynamics of order-disorder-type ferroelectrics and anti-ferroelectrics,”
*Advances in Physics*, vol. 21, pp. 693–757, 1972. View at Google Scholar - R. Blinc, B. Lavrencic, I. Levstek, V. Smolej, and B. Zeks, “Underdamped soft tunnelling modes in hydrogen-bonded ferroelectrics below T
_{C},”*Physica Status Solidi B*, vol. 60, no. 1, pp. 255–259, 1973. View at Google Scholar · View at Scopus - M. Yamaguchi, M. Watanabe, K. Inoue, Y. Akishige, and T. Yagi, “Light scattering study of the coupled soft-optic and acoustic mode in hexagonal barium titanate,”
*Physical Review Letters*, vol. 75, no. 7, pp. 1399–1402, 1995. View at Publisher · View at Google Scholar · View at Scopus - Y. Luspin, J. L. Servoin, and F. Gervais, “Soft mode spectroscopy in barium titanate,”
*Journal of Physics C*, vol. 13, no. 19, pp. 3761–3773, 1980. View at Publisher · View at Google Scholar · View at Scopus - G. A. Samara, “Pressure and temperature dependences of the dielectric properties of the perovskites BaTiO
_{3}and SrTiO_{3},”*Physical Review*, vol. 151, no. 2, pp. 378–386, 1966. View at Publisher · View at Google Scholar · View at Scopus - A. Kiraci and H. Yurtseven, “Temperature dependence of the Raman frequency, damping constant and the activation 3 energy of a soft-optic mode in ferro-electric barium titanate,”
*Ferroelectrics*, vol. 432, no. 1, pp. 14–21, 2012. View at Google Scholar - H. Mashiyama, “Proton potential and a tunneling model for KDP,”
*Journal of the Korean Physical Society*, vol. 46, no. 1, pp. 63–68, 2005. View at Google Scholar · View at Scopus - F. E. A. Melo, S. G. C. Moreira, A. S. Chaves, I. Guedes, P. T. C. Freire, and J. Mendes-Filho, “Temperature-uniaxial pressure phase diagram of KH
_{2}PO_{4},”*Physical Review B*, vol. 59, no. 5, pp. 3276–3279, 1999. View at Google Scholar · View at Scopus