Physics Research International
Volume 2016 (2016), Article ID 9475740, 10 pages
http://dx.doi.org/10.1155/2016/9475740

Research Article

Study of Dielectric Properties and Ultrasonic Attenuation in KDP-Type Ferroelectrics

Vijay Singh Bist and Narayan Singh Panwar

USIC, HNB Garhwal University, Srinagar, Garhwal, Uttrakhand 24 6174, India

Received 12 October 2015; Accepted 19 January 2016

Academic Editor: Sergey B. Mirov

Copyright © 2016 Vijay Singh Bist and Narayan Singh Panwar. 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 soft mode dynamical model has been used to study dielectric properties and ultrasonic attenuation in KDP-type ferroelectric crystals. The model Hamiltonian proposed by Blinc and Zeks has been modified by considering lattice anharmonicity up to fourth-order. The correlations appearing in the dynamical equation have been evaluated using double-time thermal retarded Green’s functions method and Dyson’s equation. Without any decoupling, the higher order correlations, appearing in the dynamical equation, have been evaluated using the renormalized Hamiltonian. The expressions for collective frequencies, width, dielectric constant, ultrasonic attenuation, and tangent loss have been calculated. The dielectric properties and ultrasonic attenuation strongly depend on the relaxational mode behavior of stochastic motion of H₂PO₄ group in KDP-type ferroelectrics. By fitting model values of physical quantities, the temperature dependence of and for different value of four-body coupling coefficient and dielectric constant and loss tangent has been calculated. The calculated and observed results have been found in good agreement.

1. Introduction

The oscillations of atoms in solids are responsible for different characteristics, such as specific heat, optical, dielectric, and electrical properties. The anharmonicity in solids is responsible for existence of thermal expansion, temperature variation of elastic constants, lattice thermal conductivity, deviation of specific heat from Dulong-Petit law at high temperature, existence of ferroelectricity in certain materials, and so forth. Many attempts have been made theoretically and experimentally to find the explanation of these phenomena in terms of anharmonicity. Extensive reviews [1, 2] are available discussing the contribution of anharmonicity to various properties of crystals.

In order-disorder type ferroelectrics, as Potassium Dihydrogen Phosphate (KH₂PO₄), the transition is associated with the tunneling of proton through a barrier between two positions of minimum potential energy in double well potential in the hydrogen bond at the transition temperature [3]. Busch [4] was the first to show that KDP (KH₂PO₄) exhibits a phase transition at low temperature. KDP is prototype of a family of crystals with bridging hydrogen bonds and its physical properties have been extensively studied [3, 5–10].

Kaminow and Damen [11] first observed the soft mode associated with the ferroelectric phase transition of the KDP-type crystal at 122.3 K by measuring the low frequency Raman scattering in configuration. Since then soft mode which is connected to the susceptibility along the crystalline -axis through the Lyddane-Sachs-Teller relation [12] has been extensively studied by Scott [13] and interpreted using the pseudospin model [14–16] and its modifications in [3, 17]. In these theories, such a particular mode of proton motions along hydrogen bonds in plane is coupled to other ion modes, bearing an electrical dipole moment along -axis is considered to play an essential role for the ferroelectric transition, and therefore little attention has been paid to the modes other than () soft mode. In KDP crystal, however, there are four tunneling protons in a primitive unit cell of the paraelectric (hereafter referred to as PE) phase and consequently four normal modes are belonging to (), a doubly degenerate , and modes [18]. () mode, which is both infrared and Raman active, provides valuable information, because this mode and () mode reflect the nature of collective proton motion which triggers the phase transition, and, moreover, in contrast to (), it is directly related to the polarization along the crystalline -axis.

At first, Pak [19] employed Green’s function methods in the order-disorder type ferroelectrics and, however, did not consider the anharmonic interactions. The phonon anharmonic interactions have been found very important in explaining dielectric, thermal, and scattering properties of solids by many authors [6–9, 20, 21] in the past. Pak’s theory was further developed by Ramakrishnan and Tanaka [22], who calculated the excitation spectrum of the system but did not consider the anharmonic interactions. Their attempt, however, established the superiority of Green’s function method over the other methods. Ganguli et al. [23] modified Ramakrishnan and Tanaka theory by considering anharmonic interaction. Their treatment explains many features of order-disorder ferroelectrics. However, due to insufficient treatment of anharmonic interactions, they could not explain quantitatively good results and could not describe some very interesting properties, like dielectric properties, ultrasonic attenuation, relaxation rate, and so forth.

In the present study, the four-particle cluster model Hamiltonian with the phonon anharmonicity up to fourth-order has been taken to theoretical study of dielectric properties and ultrasonic attenuation in KDP-type crystals, using double time Green’s functions method and Dyson’s equation. Proton Green’s function and phonon Green’s function have been solved for the collective system. Expressions for collective mode frequency shifts, widths, transition (Curie) temperature, and the expectation value of the proton collective mode components at site (,) have been derived and discussed in KDP-type crystals. By fitting model values of physical quantities, the temperature dependence of and for different value of four-body coupling, dielectric constant, and loss tangent has been calculated. The calculated and observed results have been found in good agreement.

2. The Model Hamiltonian and Green’s Function

For KDP crystal, the four-particle cluster model Hamiltonian [24] along with third- and fourth-order phonon anharmonic interaction terms is expressed as

where the first two terms constitute the original pseudospin model Hamiltonian and the third is the quadrupole contribution (the four-body interaction). is the tunneling operator which measures the tunneling power of the proton between the hydrogen double well, the tunneling frequency, and the half of the difference of the occupation probabilities for the proton in the two equilibrium positions of hydrogen bond. is the two-body coupling coefficient and is the same for every pair of protons in KDP and the four-body coupling coefficient, and refers to the four hydrogen bonds in the PO₄ group in KDP. In the last fourth terms is bare phonon frequency, and are displacement and momentum operators, is proton-lattice interaction term, and and are the third- and fourth-order anharmonic coefficients.

2.1. Collective Proton Wave Width and Shift

The correlations appearing in the proton response function can be evaluated using double time thermal retarded Green’s function [25] using the symmetrical decoupling scheme, after applying Dyson’s treatment:

where the angular brackets denote the average over the large canonical ensemble and is the Heaviside step function having properties

Differentiating (2) twice with respect to time “,” using Hamiltonian (1) taking Fourier transformation, one obtains

where is higher order Green’s functions:

with

is calculated by differentiating (5) twice with respect to “” using Hamiltonian (1) and then taking Fourier transformation; one obtains

and higher order Green’s functions

with

Substituting the value of from (7) into (4) and using Dyson’s equation, one obtains

where the renormalized frequency is

with

Higher order Green’s functions are calculated using symmetrical decoupling scheme, the cross combinations are not considered because they do not contribute significantly, and (10) can be written as

where is the proton renormalized frequency of the coupled system, which on solving self-consistently takes the form

The real and imaginary parts of (12) are obtained by using the formula

and represent collective proton mode frequency shift and width given as

with

2.2. Collective Phonon Half Width and Mode Frequency Shift

The acoustic phonon frequency width and shift are obtained analogously from acoustic phonon Green’s function:

Differentiating (18) twice with respect to time “,” using Hamiltonian (1) taking Fourier transformation, one has

where

Higher order Green’s functions are evaluated without any decoupling and using renormalized Hamiltonian. Putting the evaluated value of higher order Green’s function in (19), one gets

The real part of is obtained by using (14) and the collective mode frequency shift is obtained as

with

where is occupation number and , being Boltzmann’s constant and being the absolute temperature. Calculating (20) self-consistently and approximating, the collective mode frequency is given by

Equation is obtained using model Hamiltonian (1), and and are obtained without decoupling and using the renormalized Hamiltonian

The imaginary part of is obtained by using (19) and collective phonon half width is obtained as

2.3. Order Parameter Values of , , and

The expectation values of the proton collective mode component at site “” have been obtained by Blinc and Zeks [26] as

In PE phase (), (26) represent a system of equations for the average value of the collective mode components. The solution of this system will, however, be stable only if they minimize the free energy, that is, if , and so

in the ferroelectric phase () [27]. Consider

2.4. Dielectric Constant and Tangent Loss

Following Kubo [28] and Zubarev [25] the real part of dielectric constant is given by

where is the effective dipole moment per unit cell and is the number of unit cells in the sample.

The dielectric loss is defined as the ratio of imaginary and real parts of dielectric constant and can be written as

Thus retarded phonon Green’s function is enough to determine the dielectric constant and loss tangent. Using (29) and (22), the real part of dielectric constant can be written as

For the experimental range of frequencies (as well as for KDP crystals), (31) can be reduced to

where is given by (20) and is given by .

The tangent loss is given by

where is given by (25). For , (33) becomes

where is harmonic and defect contribution, and are due to three and four phonon anharmonic interaction terms of the lattice.

2.5. Ultrasonic Attenuation

The expression for ultrasonic attenuation constant is given by

where damping constant is given by (17) and is the ultrasonic velocity. For small limit, we obtain the collective proton mode frequency width, for , from (17) expressed as

where is polarization relaxation time and is given by Litov and Garland [29] as

Similar expressions for frequency shift and ultrasonic attenuation have been derived by Zurek [30], using the Landau-Khalatnikov theory. If the collective mode damping is neglected, a discontinuity in sound velocity is obtained at , since is proportional to () for the effect of factor is to smooth out the discontinuity.

The ultrasonic attenuation peaks for and (i.e., ) goes to zero as For temperature such as , the attenuation will have the form

where

The dependence of ultrasonic attenuation on applied frequency and the collective mode frequency () may be expected to apply also to the displacive phase transition [31], for which the coupling is linear to the strain and bilinear to the soft optical normal mode coordinates, independent of whether the phonon instability occurs at , , or at a general point in the Brillouin zone. Though the linear dependence of () on has been assumed for , this dependence presumably breaks down sufficiently close to .

2.6. Transition Temperature

In the PE phase , and the stability limit of PE phase is determine by the temperature where approaches zero as Consider

using (13), where

with effective exchange coupling constant

as well as transition temperature

3. Comparison with Experiments and Discussion

3.1. Numerical Calculations

The parameters in our calculation are listed in Table 1. The calculated values of and for KDP-type crystals for different values of four-body coupling coefficient , collective phonon mode frequency , transverse dielectric constant , observed dielectric constant , and tangent loss along -axis and -axis for KDP-type crystals are listed in Tables 2 and 3. Their variations with temperature are shown in Figures 1–4.

Table 1: Blinc-de Gennes model parameters for KDP-type crystals as given by Ganguli et al. [23].

(cm⁻¹)
(cm⁻¹)
(cm⁻¹)
(cm⁻¹)
(cm⁻¹)
(cgs)
(K)

82
334
440
25.56
153
1.8
123
0.299

Table 2: Calculated values of and for KDP crystal for different values of four-body coupling coefficient .

Temperature (K)

20
0.500
0.499
0.499
0.454
0.369
0.311

40
0.498
0.499
0.499
0.455
0.369
0.311

60
0.459
0.481
0.498
0.460
0.376
0.312

80
0.257
0.278
0.487
0.482
0.446
0.317

100
0.137
0.143
0.454
0.489
0.479
0.332

120
0.102
0.104
0.409
0.490
0.485
0.354

122
0.100
0.100
0.403
0.491
0.486
0.357

123
0.000
0.000
0.400
0.492
0.492
0.358

125
0.000
0.000
0.398
0.484
0.484
0.449

130
0.000
0.000
0.389
0.468
0.468
0.457

135
0.000
0.000
0.373
0.446
0.446
0.443

140
0.000
0.000
0.352
0.426
0.426
0.436

145
0.000
0.000
0.320
0.408
0.408
0.410

150
0.000
0.000
0.300
0.396
0.396
0.392

Table 3: Calculated values of , and for KDP-type crystals.

Temperature (K)
Collective phonon mode frequency (cm⁻¹) (present study)
Transverse dielectric constant (present study)
Observed dielectric constant (present study)
Loss tangent along the -axis (present study)
Loss tangent along the -axis (present study)

125
45.65
63
35741
0.004
0.067

130
57.04
62
6144
0.00398
0.035

135
58.69
61
2286
0.00397
0.0261

140
63.04
60
874
0.00396
0.0255

145
65.91
59
486
0.00365
0.0250

Figure 1: Temperature dependence of and for KDP-type crystal for different values of . (a) , (b) , and (c) (present study).

Figure 2: Temperature dependence of collective phonon mode frequency () in PE phase for KDP-type crystals. Present calculation (shown by ) and solid line representing experimental results of Baumgartner [32] and Choi and Lockwood [33].

Figure 3: Temperature dependence of dielectric constant for KDP type crystals obtained by Busch [4] (shown by +), Deguchi and Nakamura [37] (shown by ×), and Kaminow and Harding [35] (shown by ), and obtained from Raman intensity [5] (shown by ) and solid line represents the theoretical results of present study.

Figure 4: Temperature dependence of tangent loss for KDP-type crystals at 9.2 GHz for fields along the -axis , for field along the -axis . Present calculation is represented by ‘’, and solid line represents experimental results of Kaminow and Harding [35].

3.2. Temperature Dependence of and for KDP-Type Crystal for Different Values of

Using Blinc-de Gennes model parameter values for KDP-type ferroelectrics crystals as given by Ganguli et al. [23], putting these values into (27) and (28), we have calculated temperature dependence of and for KDP-type crystal for different values of , and variation is shown in Figure 1.

In Figure 1, the curve (a) is the case of , curve (b) is , and curve (c) is In curves (a) and (b) the value of increases to the saturated value 0.5 from zero, when temperature decreases from transition temperature. That is the case of second-order phase transition. But in curve (c) the change of with the temperature starts from a nonzero value at point “A”; that is to say, when temperature decreases increases to the saturation value from the finite value of . This is the case of first-order phase transition. The temperature at point “A” is transition temperature, and the value of at “A” is the discontinuity of . The value of decreases when temperature decreases in the ferroelectric phase. On the other hand, in PE phase the value of decreases when temperature increases from transition temperature.

3.3. Temperature Dependence of Collective Phonon Mode Frequency in PE Phase

Using (23) and (41) and model parameters values from Table 1, the temperature variation of collective phonon mode frequency in PE phase is shown in Figure 2. In the PE phase the temperature dependence of normalized collective phonon frequency enables one to calculate the transition temperature (44), as well as effective exchange coupling constant (43), which increases due to proton-phonon coupling and decreases due to anharmonic interactions. By fitting model values of physical quantities, temperature dependence of collective phonon mode frequency has been calculated which compares well with experimental results of Baumgartner [32] and Choi and Lockwood [33].

3.4. Temperature Dependence of Dielectric Constant

Putting calculated values for different temperature into (32) we obtain dielectric constant for KDP-type crystals. mode corresponds to transverse () mode, which is responsible for the observed transverse dielectric properties of KDP. In the simplest approximation can be written (), where and are temperature independent parameters. The results for transverse dielectric constant obtained from the integrated intensity of Raman spectroscopy [34] and those measured by Busch [4] and Kaminow and Harding [35] are shown in Figure 3, together with the theoretical results of Havlin et al. [36]. This indicates that low frequency is closely related to the macroscopic dielectric constant . This also suggests that the -mode Raman spectrum originates neither from the second-order Raman scattering nor from the density of states due to the local disorder above but from one of the collective modes at the centre of the Brillouin zone. It should be mentioned here that the low frequency -mode continuum appears also in a deuterated KDP (DKDP), although the intensity is about one-third of that of KDP, which indicates the possibility that the spectrum is due to the hydrogen collective motion. Using (23) for mode, it can be seen that the -mode collective hydrogen motion has characteristics damping factor which slowly increases as the temperature approaches , while for that of soft mode the damping factor slowly decreases down to a finite value, which agrees with the observations of Kaminow and Damen [11].

The present results agree with the behaviour of the observed -mode Raman spectrum in the following aspects: (i) does not change appreciably as in PE phase, (ii) is weakly dependent on temperature, and (iii) because of the factor () in the numerator of (31), the susceptibility derived changes the corresponding spectrum from a simple overdamped form to a more flat one, like the -mode Raman spectrum of KH₂PO₄ [27].

The observed dielectric constant of KH₂PO₄ along -axis is shown in Figure 3. The mode may be assigned for the observed temperature dependence of . As from (41) , the real part of the dielectric constant associated with this mode, from (32), can be expressed as

which explains the Curie-Weiss behaviour of dielectric constant along the -axis of KH₂PO₄ crystal in the PE phase observed by Deguchi and Nakamura [37], Busch [4], and Kaminow and Harding [35], shown in Figure 3. For temperature , tends to maximum value, which is consistent with the theory of Hill and Ichiki [38] for TGS and KDP crystals, while Mason monodisperse theory [39] gives as . The origin of this difference in the temperature dependence of is easily traced back in monodisperse theory, and the critical slowing down of the relaxation time has a dominant effect over the Curie-Weiss law of static dielectric constant, while the Hill-Ichiki theory of distribution function of relaxation time makes contribution to finite to more dominant.

There are actually, however, many cases in which takes a minimum of finite value at being neither zero as in Mason’s theory nor maximum as in Hill and Ichiki theory.

3.5. Temperature Dependence of Tangent Loss

The tangent loss is associated with damping parameter (25). Damping can be understood as the creation of a virtual polarization mode excited by the transverse electromagnetic radiation and its subsequent decay into the real phonons by scattering from crystal defects, higher order phonon anharmonicities, and so forth. At higher temperature the loss deviates from the Curie-Weiss type behaviour and increases linearly with temperature. This behaviour suggests that at higher temperatures the phonon anharmonicity contributes significantly to the observed loss.

The calculated values of tangent loss from (33) and (34) and experimental result of Kaminow and Harding [35] are shown in Figure 4 and given in Table 3. The loss was calculated at 9.2 GHz because experimental data is available only at that frequency range. The theoretical results of the present study are in good agreement with experimental results of Kaminow and Harding [35]. The temperature dependence of loss does not appear to be exponential. Thus, third- and fourth-order anharmonicity may be responsible for the observed behaviour of loss tangent. In the microwave frequency rage, an increase in frequency is followed by an increase in transverse and longitudinal dielectric loss tangent. The loss decreases with increase in temperature for KDP-type crystals, in their PE phase. This shows Curie-Weiss type behaviour of the dielectric loss tangent.

In this paper the four-particle cluster model Hamiltonian with the anharmonic contributions up to fourth-order has been taken into consideration in study of dielectric properties and ultrasonic attenuation in KDP-type ferroelectric crystals. Using double time thermal Green’s function method and Dyson’s equation the collective mode frequencies and widths have been calculated. These parameters lead to the expressions for the dielectric constant and loss tangent. The observed dielectric properties have been explained in terms of present study. The expressions for dielectric constant and loss tangent have been derived and compared with the experimental results. Using Blinc-de Gennes model parameter values given by Ganguli et al. [23] we have calculated temperature variations of these quantities for KDP- type crystals. It is observed that these results are in good agreement with each other and with the results obtained by other methods. The present results reduce to the results of others [22, 23, 40] if the width and shift are neglected. Only Ganguli et al. [23] modified the Ramakrishnan and Tanaka theory by considering anharmonic interaction. Their treatment explains many features of order-disorder ferroelectrics. However, due to insufficient treatment of anharmonic interactions, they could not obtain quantitatively good results and could not describe some very interesting properties, like dielectric properties, acoustic attenuation, and so forth.

Thus, from the present study, it is concluded that the consideration of four-cluster Hamiltonian along with the third- and fourth-order anharmonicities for the KDP-type ferroelectrics leads to the renormalization and stabilization of the relaxational soft mode and the renormalization of the pseudospin exchange interaction constant. The decoupling of the correlations appearing in the dynamical equation after applying Dyson’s equation results in shift in frequency and facilitates the calculation of damping parameter, which is related to the loss tangent. The anomalous behaviour in order-disorder KDP-type ferroelectrics finds explanation by the consideration of collective proton-phonon interaction and third- and fourth-order anharmonicities in the four-particle cluster Hamiltonian. The dielectric properties and ultrasonic attenuation strongly depend on the relaxational mode behaviour of stochastic motion of H₂PO₄ group in KDP-type ferroelectrics.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

B. S. Semwal and P. K. Sharma, “Thermal conductivity of an anharmonic crystal,” Physical Review B, vol. 5, no. 10, pp. 3909–3914, 1972.

B. S. Semwal and P. K. Sharma, “Thermal conductivity of an isotopically disordered harmonic crystal,” Journal of Mathematical Physics, vol. 15, no. 5, pp. 648–653, 1974.

K. K. Kobayashi, “Dynamical theory of the phase transition in KH2PO4-type ferroelectric crystals,” Journal Physics Society Japan, vol. 24, p. 297, 1968.

G. Busch, “Neue seignette-elektrika,” Helvetica Physica Acta, vol. 11, no. 4, p. 269, 1938.

M. Tokunaga and I. Tatsuzaki, “Light scattering spectra of polarization fluctuations and models of the phase transition in KDP type ferroelectrics,” Phase Transitions, vol. 4, no. 2, pp. 97–155, 1984.

S. C. Deorani, U. C. Naithani, and B. S. Semwal, “Four—spin contribution in DKDP,” Indian Journal of Pure & Applied Physics, vol. 37, pp. 215–218, 1999.

T. C. Upadhyay and B. S. Semwal, “Temperature dependence of dielectric tangent loss in KDP and DKDP crystal at microwave frequencies,” Indian Journal of Pure & Applied Physics, vol. 40, pp. 615–619, 2002.

T. C. Upadhyay and B. S. Semwal, “Microwave dielectric tangent losses in KDP and DKDP crystals,” Pramana Journal of Physics, vol. 60, no. 3, pp. 525–533, 2003.

T. C. Upadhyay, R. S. Bhandari, and B. S. Semwal, “Dielectric properties of KDP-type ferroelectric crystals in the presence of external electric field,” Pramana, vol. 67, no. 3, pp. 547–552, 2006.

H. Yurtseven and S. Şen, “Calculation of the spontaneous polarization and the dielectric constant as a function of temperature for KH₂PO₄,” Advances in Condensed Matter Physics, vol. 2012, Article ID 259712, 5 pages, 2012.

I. P. Kaminow and T. C. Damen, “Temperature dependence of the Ferroelectric mode in KH₂PO₄,” Physical Review Letters, vol. 20, no. 20, pp. 1105–1108, 1968.

R. H. Lyddane, R. G. Sachs, and E. Teller, “On the polar vibrations of alkali halides,” Physical Review, vol. 59, no. 8, pp. 673–676, 1941.

J. F. Scott, “Soft-mode spectroscopy: experimental studies of structural phase transitions,” Reviews of Modern Physics, vol. 46, no. 1, pp. 83–128, 1974.

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.

M. Tokunaga and T. Matsubara, “Theory of ferroelectric phase transition in KH₂PO₄ type crystals I,” Progress of Theoretical Physics, vol. 35, pp. 581–599, 1996.

P. G. de Gennes, “Collective motions of hydrogen bonds,” Solid State Communications, vol. 1, no. 6, pp. 132–137, 1963.

R. S. Katiyar, J. F. Ryan, and J. F. Scott, “Proton-phonon coupling in CsH₂ AsO₄ and KH₂AsO₄,” Physical Review B, vol. 4, no. 8, pp. 2635–2638, 1971.

B. Labrenčič, I. Levstek, B. Žekš, R. Blinc, and D. Hadži, “Laser Raman study of quasi-spin wave hydrogen tunneling modes in KH₂PO₄ and KD₂PO₄,” Chemical Physics Letters, vol. 5, no. 7, pp. 441–444, 1970.

K. N. Pak, “Strong coupling between proton and phonon in KDP,” Physica Status Solidi B, vol. 60, no. 1, pp. 233–239, 1973.

S. C. Deorani, P. Singh, and M. S. Yadav, “Investigation of deuteration effect on sound attenuation in hydrogen bonded ferroelectric crystals,” Material Science Research India, vol. 8, no. 1, pp. 137–141, 2011.

N. S. Panwar and B. S. Semwal, “Soft mode dynamics of ABO3—type trigonal ferroelectrics: application to lithium niobate,” Indian Journal of Pure & Applied Physics (IJPAP), vol. 28, p. 706, 1990.

V. Ramakrishnan and T. Tanaka, “Green's-function theory of the ferroelectric phase transition in potassium dihydrogen phosphate (KDP),” Physical Review B, vol. 16, no. 1, pp. 422–426, 1977.

S. Ganguli, D. Nath, and B. K. Chaudhuri, “Green's-function theory of phase transitions in hydrogen-bonded ferroelectric crystals with pseudo-spin-lattice coupled mode model,” Physical Review B, vol. 21, no. 7, pp. 2937–2945, 1980.

R. Blinc and B. Žekš, “Proton order-disorder in KH₂PO₄-type ferroelectrics: slater theory and ising model in a transverse tunneling field,” Ferroelectrics, vol. 72, no. 1, pp. 193–227, 1987.

D. N. Zubarev, “Double-time green functions in statistical physics,” Uspekhi Fizicheskikh Nauk, vol. 71, pp. 71–116, 1960.

R. Blinc and B. Zeks, “In soft modes in ferroelectrics and antiferroelectrics,” Advance in Physics, vol. 29, p. 701, 1972.

Q. Zikai, Z. Jinbo, and W. Chunlei, “Application of Green's function method to the study of order-disorder ferroelectrics,” Ferroelectrics, vol. 101, no. 1, pp. 159–169, 1990.

R. Kubo, “Statistical-mechanical theory of irreversible processes. I. general theory and simple applications to magnetic and conduction problems,” Journal of the Physical Society of Japan, vol. 12, no. 6, pp. 570–586, 1957.

E. Litov and C. W. Garland, “Ultrasonic experiments in KDP-type crystals,” Ferroelectrics, vol. 72, no. 1, pp. 19–44, 1987.

R. Zurek, Personal Communication, 1969.

E. Pytte, “Soft-mode damping and ultrasonic attenuation at a structural phase transition,” Physical Review B, vol. 1, no. 2, pp. 924–930, 1970.

H. Baumgartner, “Elektrische sättigungserscheinungen und elektrokalorischer effekt von kaliumphosphat KH2PO4,” Helvetica Physica Acta, vol. 24, no. 7, p. 651, 1950.

B. Choi and D. Lockwood, “Dielectric anomalies in KH₂PO₄,” Ferroelectrics, vol. 106, no. 1, pp. 15–20, 1990.

Y. Takagi and T. Shigenari, “Transverse susceptibility and E-mode Raman spectra of a KH₂PO₄ crystal,” Journal of the Physical Society of Japan, vol. 39, no. 2, pp. 440–447, 1975.

I. P. Kaminow and G. O. Harding, “Complex dielectric constant of KH₂PO₄ at 9.2 Gc/sec,” Physical Review, vol. 129, no. 4, pp. 1562–1566, 1963.

S. Havlin, E. Litov, and E. A. Uehling, “Transverse susceptibility in KH₂PO₄-type crystals,” Physical Review B, vol. 9, no. 3, article 1024, 1974.

K. Deguchi and E. Nakamura, “Deviation from the Curie-Weiss Law in KH₂PO₄,” Journal of the Physical Society of Japan, vol. 49, no. 5, pp. 1887–1891, 1980.

R. M. Hill and S. K. Ichiki, “High-frequency behavior of hydrogen-bonded ferroelectrics: triglycine sulphate and KD₂PO₄,” Physical Review, vol. 132, no. 4, pp. 1603–1608, 1963.

W. P. Mason, “Theory of the ferroelectric effect and clamped dielectric constant of Rochelle salt,” Physical Review, vol. 72, no. 9, pp. 854–865, 1947.

H. Konwent, “On the theory of hydrogen-bonded ferroelectrics,” Physica Status Solidi B, vol. 28, no. 1, pp. 39–44, 1968.