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

Authors: Vijay Singh Bist; Narayan Singh Panwar
Journal: Physics Research International
 (2016)

---

## 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 H2PO4 group in KDP-type ferroelectrics. By fitting model values of physical quantities, the temperature dependence of Sqz and Sqx 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.


---

## Body

## 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 (KH2PO4), 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 (KH2PO4) 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 x(yx)y configuration. Since then B2(z) soft mode which is connected to the susceptibility along the crystalline c-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 a-b plane is coupled to other ion modes, bearing an electrical dipole moment along c-axis is considered to play an essential role for the ferroelectric transition, and therefore little attention has been paid to the modes other than B2(z) 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 B2(z), a doubly degenerate E(x,y), and A2 modes [18]. E(x,y) mode, which is both infrared and Raman active, provides valuable information, because this mode and B2(z) mode reflect the nature of collective proton motion which triggers the phase transition, and, moreover, in contrast to B2(z), it is directly related to the polarization along the crystalline a-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 siteq (Sqz,Sqx) have been derived and discussed in KDP-type crystals. By fitting model values of physical quantities, the temperature dependence of Sqz and Sqx 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(1)H=-2Ω∑iSix-12∑ijJijSizSjz-14∑ijklJijkl′SizSjzSkzSlz+14ωkAk+Ak+Bk+Bk-∑iV-ikSixAk+∑k→1,k→2,k→3V3k→1k→2k→3Ak→1,Ak→2,Ak→3+∑k→1,k→2,k→3,k→4V4k→1k→2k→3k→4Ak→1,Ak→2,Ak→3,Ak→4,where the first two terms constitute the original pseudospin model Hamiltonian and the third is the quadrupole contribution (the four-body interaction). Six is the tunneling operator which measures the tunneling power of the proton between the hydrogen double well, Ω the tunneling frequency, and Siz the half of the difference of the occupation probabilities for the proton in the two equilibrium positions of hydrogen bond. Jij is the two-body coupling coefficient and is the same for every pair of protons in KDP and the four-body coupling coefficient, and Jijkl′ refers to the four hydrogen bonds in the PO4 group in KDP. In the last fourth terms ωk is bare phonon frequency, Ak and Bk are displacement and momentum operators, V-ik is proton-lattice interaction term, and V3 and V4 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:(2)Gqq′zzt-t′=Sqzt;Sq′zt′=-jθt-t′Sqzt,Sq′zt′,where the angular brackets denote the average over the large canonical ensemble and θ(t) is the Heaviside step function having properties(3)θt=1fort>1,θt=0fort<1,j=-11/2.Differentiating (2) twice with respect to time “t,” using Hamiltonian (1) taking Fourier transformation, one obtains (4)ω2-Ω2Gqq′zz=ΩSqxδqq′π+ΩπEω,where E(ω) is higher order Green’s functions:(5)Eω=Fqt;Sq′z,with(6)
F
q
t
=
π
-
2
V
-
q
S
q
x
A
q
-
J
q
S
q
z
S
q
′
x
+
S
q
x
S
q
′
z
-
1
2
·
J
q
′
S
q
x
S
q
′
z
S
-
q
z
S
-
q
′
z
+
S
q
′
z
S
q
x
S
-
q
z
S
-
q
′
z
+
S
-
q
z
S
q
z
S
q
′
x
S
-
q
′
z
+
S
-
q
′
z
S
q
z
S
q
′
z
S
-
q
x
.E(ω) is calculated by differentiating (5) twice with respect to “t” using Hamiltonian (1) and then taking Fourier transformation; one obtains(7)Eω=πω2-Ω2Gqq′zzjΩ-E-ωj,and higher order Green’s functions(8)E-ω=Fqt;Fq′t′with(9)
F
q
′
t
′
=
π
-
2
V
-
q
S
q
x
A
q
-
J
q
S
q
x
S
q
′
z
+
S
q
x
S
q
′
z
-
1
2
·
J
q
′
S
q
x
S
q
′
z
S
-
q
z
S
-
q
′
z
+
S
q
′
z
S
q
x
S
-
q
z
S
-
q
′
z
+
S
-
q
z
S
q
z
S
q
′
x
S
-
q
′
z
+
S
-
q
′
z
S
q
z
S
q
′
z
S
-
q
x
.Substituting the value of E-(ω) from (7) into (4) and using Dyson’s equation, one obtains(10)Gqq′zzω+jε=ΩSqxδqq′πω2-Ω~-Ω/2πSqxE--ω,where the renormalized frequency is (11)Ω~2=a2+b2-bc,with11aa=J0Sqz+J0′Sqz3,11bb=2Ω,11cc=J0Sqx+3J0′SqxSqz2,11dJ0=∑qJq,J0′=∑q′Jq′′.Higher order Green’s functions E--ω are calculated using symmetrical decoupling scheme, the cross combinations are not considered because they do not contribute significantly, and (10) can be written as(12)Limitε→0⁡Gqq′zzω+jε=ΩSqxδqq′πω2-Ω~~2+jΓsq,ω,where Ω~~ is the proton renormalized frequency of the coupled system, which on solving self-consistently takes the form(13)Ω~~2=Ω~2+2ΩΔsq,ω.The real and imaginary parts of (12) are obtained by using the formula(14)limitε→0⁡1x±jε=1x∓jπδx.Δs(q,ω) and Γs(q,ω) represent collective proton mode frequency shift and width given as(15)Δsq,ω=2V-q2ωqSqxω2-ω~q2Ωω2-ω~q22+4ωq2Γp2+bc2ω2-Ω~2+a2Ω^2ω2-Ω^2,with(16)Ω^2=a2+nqV-q2,(17)Γsq,ω=-4πV-q2ωq2Sqxδqq′ΓpΩω2-ω~~q22+4ωq2Γp2+πbc22Ω~δω-Ω~-δω+Ω~+πa2Ω^2bδω-Ω^-δω+Ω^.


### 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:(18)Gt-t′=Aq;Aq′+.Differentiating (18) twice with respect to time “t,” using Hamiltonian (1) taking Fourier transformation, one has (19)Gpq,ω=ωqδqq′πω2-ω~~q2+2jωqΓpq,ω,where(20)ω~~q2=ω~q2+2ωqΔpq,ω=ω~q2+2ωq∑i=13Gip′q,ω,20aω~q2=ωq2+8ωq2V3+V4coth⁡βωq′2.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(21)Gpω=ωqδqq′πω2-ω~~q2-Pq,ω.The real part of P(q,ω) is obtained by using (14) and the collective mode frequency shift is obtained as(22)Δpq,ω=RePq,ω=∑i=13Gip′q,ω,with22
a
G
1
p
′
q
,
ω
=
2
V
-
q
2
Ω
S
q
x
ω
2
-
Ω
~
~
2
δ
q
q
′
ω
q
ω
2
-
Ω
~
~
2
2
+
4
Ω
2
Γ
s
2
q
,
ω
,
22
b
G
2
p
′
q
,
ω
=
24
∑
q
V
3
2
q
q
′
ω
q
n
q
ω
~
ω
2
-
4
ω
~
q
2
,
22
c
G
3
p
′
q
,
ω
=
24
∑
q
V
4
2
q
q
′
ω
q
ω
q
′
ω
~
q
2
ω
~
q
′
1
+
2
n
q
n
q
′
+
n
q
2
2
ω
~
q
+
ω
~
q
′
ω
2
-
2
ω
~
q
+
ω
~
q
′
2
+
n
q
2
-
1
2
ω
~
q
+
ω
~
q
′
ω
2
-
2
ω
~
q
-
ω
~
q
′
2
+
2
ω
~
q
′
n
q
2
-
1
ω
2
-
ω
q
′
2
+
72
∑
q
V
4
2
q
q
′
ω
q
2
ω
~
q
2
1
+
3
n
q
2
ω
2
-
9
ω
~
q
2
+
n
q
2
-
1
ω
2
-
ω
q
′
2
,where nq=ωq/ω~qcoth⁡(βω~q/2) is occupation number and β=(kBT)-1, kB being Boltzmann’s constant and T being the absolute temperature. Calculating (20) self-consistently and approximating, the collective mode frequency is given by(23)ω~~q±2=12ω~q2+Ω~~2±12ω~q2-Ω~~22+16V-q2ωqΩSx1/2.Equation 22a is obtained using model Hamiltonian (1), and 22b and 22c are obtained without decoupling and using the renormalized Hamiltonian (24)Hren=-2Ω∑qSqx-12∑qq′JqSqzSq′z-14∑qq′-q-q′Jq′SqzSq′zS-qzS-q′z+14∑qωq2ωqAk+Ak+Bk+Bk.The imaginary part of P(q,ω) is obtained by using (19) and collective phonon half width is obtained as(25)Γpq,ω=πIm⁡Pq,ω=π∑i=33Gip′′q,ω.


### 2.3. Order Parameter Values ofSqx, Sqy, and Sqz

The expectation values of the proton collective mode component at site “q” have been obtained by Blinc and Zeks [26] as (26)Sqx=ΩΩ~tanh⁡βΩ~2,Sqy=0,Sqz=a2Ω~tanh⁡βΩ~2.In PE phase (T>Tc), (26) represent a system of 3N 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 Sqz=Sqy=0, and so(27)Sqx=Ω4Ω2-2ΩJ0Sqx1/2·tanh⁡4Ω2-2ΩJ0Sqx1/22kBT,in the ferroelectric phase (T<Tc) [27]. Consider(28)Sqx=2ΩJ0+J0′Sqz2,Sqy=0,Sqz=12tanh⁡a2kBT.


### 2.4. Dielectric Constant and Tangent Loss

Following Kubo [28] and Zubarev [25] the real part of dielectric constant is given by (29)ε′ω-1=-8π2Nμ2G′ω,where μ is the effective dipole moment per unit cell and N is the number of unit cells in the sample.The dielectric loss(tan⁡δ) is defined as the ratio of imaginary and real parts of dielectric constant and can be written as(30)tan⁡δ=G′′ωG′ω.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(31)ε′ω-1=-8πNμ2ω~ω2-ω~~2ω2-ω~~22+4ω2Γp2ω.For the experimental range of frequencies ω≪ω~~ (as well as ωτp≪1 for KDP crystals), (31) can be reduced to(32)ε′ω≅8πNμ2ω~ω~~2,asε′ω≫1,where ω~~2 is given by (20) and ω~ is given by 20a.The tangent loss is given by(33)tan⁡δ=ε′′ε=-ωΓpω2-ω~~2,where Γp is given by (25). For ω≪ω~~, (33) becomes(34)tan⁡δ=-ωΓpω~~2=ωα+βT+γT2T-Tc,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(35)αpq,ω=Γsq,ωCμq,where damping constant is given by (17) and Cμ(q) is the ultrasonic velocity. For small q limit, we obtain the collective proton mode frequency width, for ω≪ω~~, from (17) expressed as(36)Γs=2V-q2Sxω2τpω~~q2Ω1+ω2τp2+πbc22Ω~δω-Ω~-δω+Ω~+πa2Ω^2bδω-Ω^-δω+Ω^,where τp is polarization relaxation time and is given by Litov and Garland [29] as (37)τp=2Γpω~~-2.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 T=Tc, since ω~~-2 is proportional to (Tc-T) for T≤Tc. the effect of factor 1+ω2τp2-1 is to smooth out the discontinuity.The ultrasonic attenuation peaks forωτp≅1 and ωq (i.e., ω~~-) goes to zero as T→Tc. For temperature such as ωτp≪1, the attenuation will have the form(38)αq=Λω2τpTc-T+α0,where(39)α0=πbc22Ω~Cμδω-Ω~-δω+Ω~+πa2Ω^2bCμδω-Ω^-δω+Ω^,Λ=2V-q2SxΩCμγ.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 q=0, q=qk, or at a general point in the Brillouin zone. Though the linear dependence of (ω~~-) on (T-Tc) has been assumed for T≅Tc, this dependence presumably breaks down sufficiently close to Tc.


### 2.6. Transition Temperature

In the PE phaseSqz=0, and the stability limit of PE phase is determine by the temperature where ω~~q- approaches zero as T→Tc. Consider(40)ω~~q-2=∂ω~~q-2∂TT=TcT-Tc,(41)ω~~q-2≅γT-Tc,using (13), where (42)γ=Ω2J^kBTc2cosh⁡2Ω/kBTc,with effective exchange coupling constant(43)J^=J0+2V-q2ωqω~qT=Tc,as well as transition temperature(44)Tc=Ωkβtanh-1⁡4ΩJ^-1.


## 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:(2)Gqq′zzt-t′=Sqzt;Sq′zt′=-jθt-t′Sqzt,Sq′zt′,where the angular brackets denote the average over the large canonical ensemble and θ(t) is the Heaviside step function having properties(3)θt=1fort>1,θt=0fort<1,j=-11/2.Differentiating (2) twice with respect to time “t,” using Hamiltonian (1) taking Fourier transformation, one obtains (4)ω2-Ω2Gqq′zz=ΩSqxδqq′π+ΩπEω,where E(ω) is higher order Green’s functions:(5)Eω=Fqt;Sq′z,with(6)
F
q
t
=
π
-
2
V
-
q
S
q
x
A
q
-
J
q
S
q
z
S
q
′
x
+
S
q
x
S
q
′
z
-
1
2
·
J
q
′
S
q
x
S
q
′
z
S
-
q
z
S
-
q
′
z
+
S
q
′
z
S
q
x
S
-
q
z
S
-
q
′
z
+
S
-
q
z
S
q
z
S
q
′
x
S
-
q
′
z
+
S
-
q
′
z
S
q
z
S
q
′
z
S
-
q
x
.E(ω) is calculated by differentiating (5) twice with respect to “t” using Hamiltonian (1) and then taking Fourier transformation; one obtains(7)Eω=πω2-Ω2Gqq′zzjΩ-E-ωj,and higher order Green’s functions(8)E-ω=Fqt;Fq′t′with(9)
F
q
′
t
′
=
π
-
2
V
-
q
S
q
x
A
q
-
J
q
S
q
x
S
q
′
z
+
S
q
x
S
q
′
z
-
1
2
·
J
q
′
S
q
x
S
q
′
z
S
-
q
z
S
-
q
′
z
+
S
q
′
z
S
q
x
S
-
q
z
S
-
q
′
z
+
S
-
q
z
S
q
z
S
q
′
x
S
-
q
′
z
+
S
-
q
′
z
S
q
z
S
q
′
z
S
-
q
x
.Substituting the value of E-(ω) from (7) into (4) and using Dyson’s equation, one obtains(10)Gqq′zzω+jε=ΩSqxδqq′πω2-Ω~-Ω/2πSqxE--ω,where the renormalized frequency is (11)Ω~2=a2+b2-bc,with11aa=J0Sqz+J0′Sqz3,11bb=2Ω,11cc=J0Sqx+3J0′SqxSqz2,11dJ0=∑qJq,J0′=∑q′Jq′′.Higher order Green’s functions E--ω are calculated using symmetrical decoupling scheme, the cross combinations are not considered because they do not contribute significantly, and (10) can be written as(12)Limitε→0⁡Gqq′zzω+jε=ΩSqxδqq′πω2-Ω~~2+jΓsq,ω,where Ω~~ is the proton renormalized frequency of the coupled system, which on solving self-consistently takes the form(13)Ω~~2=Ω~2+2ΩΔsq,ω.The real and imaginary parts of (12) are obtained by using the formula(14)limitε→0⁡1x±jε=1x∓jπδx.Δs(q,ω) and Γs(q,ω) represent collective proton mode frequency shift and width given as(15)Δsq,ω=2V-q2ωqSqxω2-ω~q2Ωω2-ω~q22+4ωq2Γp2+bc2ω2-Ω~2+a2Ω^2ω2-Ω^2,with(16)Ω^2=a2+nqV-q2,(17)Γsq,ω=-4πV-q2ωq2Sqxδqq′ΓpΩω2-ω~~q22+4ωq2Γp2+πbc22Ω~δω-Ω~-δω+Ω~+πa2Ω^2bδω-Ω^-δω+Ω^.


## 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:(18)Gt-t′=Aq;Aq′+.Differentiating (18) twice with respect to time “t,” using Hamiltonian (1) taking Fourier transformation, one has (19)Gpq,ω=ωqδqq′πω2-ω~~q2+2jωqΓpq,ω,where(20)ω~~q2=ω~q2+2ωqΔpq,ω=ω~q2+2ωq∑i=13Gip′q,ω,20aω~q2=ωq2+8ωq2V3+V4coth⁡βωq′2.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(21)Gpω=ωqδqq′πω2-ω~~q2-Pq,ω.The real part of P(q,ω) is obtained by using (14) and the collective mode frequency shift is obtained as(22)Δpq,ω=RePq,ω=∑i=13Gip′q,ω,with22
a
G
1
p
′
q
,
ω
=
2
V
-
q
2
Ω
S
q
x
ω
2
-
Ω
~
~
2
δ
q
q
′
ω
q
ω
2
-
Ω
~
~
2
2
+
4
Ω
2
Γ
s
2
q
,
ω
,
22
b
G
2
p
′
q
,
ω
=
24
∑
q
V
3
2
q
q
′
ω
q
n
q
ω
~
ω
2
-
4
ω
~
q
2
,
22
c
G
3
p
′
q
,
ω
=
24
∑
q
V
4
2
q
q
′
ω
q
ω
q
′
ω
~
q
2
ω
~
q
′
1
+
2
n
q
n
q
′
+
n
q
2
2
ω
~
q
+
ω
~
q
′
ω
2
-
2
ω
~
q
+
ω
~
q
′
2
+
n
q
2
-
1
2
ω
~
q
+
ω
~
q
′
ω
2
-
2
ω
~
q
-
ω
~
q
′
2
+
2
ω
~
q
′
n
q
2
-
1
ω
2
-
ω
q
′
2
+
72
∑
q
V
4
2
q
q
′
ω
q
2
ω
~
q
2
1
+
3
n
q
2
ω
2
-
9
ω
~
q
2
+
n
q
2
-
1
ω
2
-
ω
q
′
2
,where nq=ωq/ω~qcoth⁡(βω~q/2) is occupation number and β=(kBT)-1, kB being Boltzmann’s constant and T being the absolute temperature. Calculating (20) self-consistently and approximating, the collective mode frequency is given by(23)ω~~q±2=12ω~q2+Ω~~2±12ω~q2-Ω~~22+16V-q2ωqΩSx1/2.Equation 22a is obtained using model Hamiltonian (1), and 22b and 22c are obtained without decoupling and using the renormalized Hamiltonian (24)Hren=-2Ω∑qSqx-12∑qq′JqSqzSq′z-14∑qq′-q-q′Jq′SqzSq′zS-qzS-q′z+14∑qωq2ωqAk+Ak+Bk+Bk.The imaginary part of P(q,ω) is obtained by using (19) and collective phonon half width is obtained as(25)Γpq,ω=πIm⁡Pq,ω=π∑i=33Gip′′q,ω.


## 2.3. Order Parameter Values ofSqx, Sqy, and Sqz

The expectation values of the proton collective mode component at site “q” have been obtained by Blinc and Zeks [26] as (26)Sqx=ΩΩ~tanh⁡βΩ~2,Sqy=0,Sqz=a2Ω~tanh⁡βΩ~2.In PE phase (T>Tc), (26) represent a system of 3N 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 Sqz=Sqy=0, and so(27)Sqx=Ω4Ω2-2ΩJ0Sqx1/2·tanh⁡4Ω2-2ΩJ0Sqx1/22kBT,in the ferroelectric phase (T<Tc) [27]. Consider(28)Sqx=2ΩJ0+J0′Sqz2,Sqy=0,Sqz=12tanh⁡a2kBT.


## 2.4. Dielectric Constant and Tangent Loss

Following Kubo [28] and Zubarev [25] the real part of dielectric constant is given by (29)ε′ω-1=-8π2Nμ2G′ω,where μ is the effective dipole moment per unit cell and N is the number of unit cells in the sample.The dielectric loss(tan⁡δ) is defined as the ratio of imaginary and real parts of dielectric constant and can be written as(30)tan⁡δ=G′′ωG′ω.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(31)ε′ω-1=-8πNμ2ω~ω2-ω~~2ω2-ω~~22+4ω2Γp2ω.For the experimental range of frequencies ω≪ω~~ (as well as ωτp≪1 for KDP crystals), (31) can be reduced to(32)ε′ω≅8πNμ2ω~ω~~2,asε′ω≫1,where ω~~2 is given by (20) and ω~ is given by 20a.The tangent loss is given by(33)tan⁡δ=ε′′ε=-ωΓpω2-ω~~2,where Γp is given by (25). For ω≪ω~~, (33) becomes(34)tan⁡δ=-ωΓpω~~2=ωα+βT+γT2T-Tc,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(35)αpq,ω=Γsq,ωCμq,where damping constant is given by (17) and Cμ(q) is the ultrasonic velocity. For small q limit, we obtain the collective proton mode frequency width, for ω≪ω~~, from (17) expressed as(36)Γs=2V-q2Sxω2τpω~~q2Ω1+ω2τp2+πbc22Ω~δω-Ω~-δω+Ω~+πa2Ω^2bδω-Ω^-δω+Ω^,where τp is polarization relaxation time and is given by Litov and Garland [29] as (37)τp=2Γpω~~-2.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 T=Tc, since ω~~-2 is proportional to (Tc-T) for T≤Tc. the effect of factor 1+ω2τp2-1 is to smooth out the discontinuity.The ultrasonic attenuation peaks forωτp≅1 and ωq (i.e., ω~~-) goes to zero as T→Tc. For temperature such as ωτp≪1, the attenuation will have the form(38)αq=Λω2τpTc-T+α0,where(39)α0=πbc22Ω~Cμδω-Ω~-δω+Ω~+πa2Ω^2bCμδω-Ω^-δω+Ω^,Λ=2V-q2SxΩCμγ.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 q=0, q=qk, or at a general point in the Brillouin zone. Though the linear dependence of (ω~~-) on (T-Tc) has been assumed for T≅Tc, this dependence presumably breaks down sufficiently close to Tc.


## 2.6. Transition Temperature

In the PE phaseSqz=0, and the stability limit of PE phase is determine by the temperature where ω~~q- approaches zero as T→Tc. Consider(40)ω~~q-2=∂ω~~q-2∂TT=TcT-Tc,(41)ω~~q-2≅γT-Tc,using (13), where (42)γ=Ω2J^kBTc2cosh⁡2Ω/kBTc,with effective exchange coupling constant(43)J^=J0+2V-q2ωqω~qT=Tc,as well as transition temperature(44)Tc=Ωkβtanh-1⁡4ΩJ^-1.


## 3. Comparison with Experiments and Discussion

### 3.1. Numerical Calculations

The parameters in our calculation are listed in Table1. The calculated values of Sqz and Sqx for KDP-type crystals for different values of four-body coupling coefficient J′, collective phonon mode frequency ω~~-, transverse dielectric constant εao, observed dielectric constant εco, and tangent loss along a-axis (tan⁡δa) and c-axis (tan⁡δc) 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−1)
J (cm−1)
J
′ (cm−1)
V
i
k (cm−1)
ω
k (cm−1)
μ
×
1
0
18 (cgs)
T
c
(K)
V
/
k
T
c
82
334
440
25.56
153
1.8
123
0.299Table 2
Calculated values ofSqz and Sqx for KDP crystal for different values of four-body coupling coefficient (J′).
S
q
z
S
q
x
Temperature (K)
J
′
=
J
/
3
J
′
=
4
J
/
3
J
′
=
7
J
/
3
J
′
=
J
/
3
J
′
=
4
J
/
3
J
′
=
7
J
/
3
P
r
e
s
e
n
t
s
t
u
d
y
P
r
e
s
e
n
t
s
t
u
d
y
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.392Table 3
Calculated values ofω~~-,εa(0),εc(0),tan⁡δa, and tan⁡δc for KDP-type crystals.
Temperature (K)
Collective phonon mode frequencyω~~- (cm−1) (present study)
Transverse dielectric constantεa(0) (present study)
Observed dielectric constantεc(0) (present study)
Loss tangent along thea-axis (tan⁡δa) (present study)
Loss tangent along thec-axis (tan⁡δc) (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.0250Figure 1
Temperature dependence ofSqz and Sqx for KDP-type crystal for different values of J′. (a) J′=J/3, (b) J′=4J/3, and (c) J′=7J/3 (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 thea-axis (tan⁡δa), for field along the c-axis (tan⁡δc). Present calculation is represented by ‘∘’, and solid line represents experimental results of Kaminow and Harding [35].


### 3.2. Temperature Dependence ofSqz and Sqx for KDP-Type Crystal for Different Values of J′

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 Sqz and Sqx for KDP-type crystal for different values of J′, and variation is shown in Figure 1.In Figure1, the curve (a) is the case of J′=J/3<4J/3, curve (b) is J′=4J/3<7J/3, and curve (c) is J′=7J/3>4J/3. In curves (a) and (b) the value of Sqz 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 Sqz with the temperature starts from a nonzero value Sqz at point “A”; that is to say, when temperature decreases Sqz increases to the saturation value from the finite value of Sqz. This is the case of first-order phase transition. The temperature at point “A” is transition temperature, and the value of Sqz at “A” is the discontinuity of Sqz. The value of Sqx decreases when temperature decreases in the ferroelectric phase. On the other hand, in PE phase the value of Sqx 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 E(x,y) mode, which is responsible for the observed transverse dielectric properties of KDP. In the simplest approximation ω~~+ can be written (ω~~+=K1+K2T), where K1 and K2 are temperature independent parameters. The results for transverse dielectric constant εa′ 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 ω~~+E(x,y)mode is closely related to the macroscopic dielectric constant εa′. This also suggests that the E-mode Raman spectrum originates neither from the second-order Raman scattering nor from the density of states due to the local disorder above Tc but from one of the collective modes at the centre of the Brillouin zone. It should be mentioned here that the low frequency E-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 E-mode collective hydrogen motion has characteristics damping factor which slowly increases as the temperature approaches Tc, while for that of ω~~-B2z 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 observedE-mode Raman spectrum in the following aspects: (i) ω~~+ does not change appreciably as T→Tc in PE phase, (ii) Γp(ω+) is weakly dependent on temperature, and (iii) because of the factor (ω2-ω~~2) in the numerator of (31), the susceptibility derived changes the corresponding spectrum from a simple overdamped form to a more flat one, like the E-mode Raman spectrum of KH2PO4 [27].The observed dielectric constant(εc′) of KH2PO4 along c-axis is shown in Figure 3. The ω~~-B2z mode may be assigned for the observed temperature dependence of εc′. As from (41) ω~~q2∝T-Tc, the real part of the dielectric constant associated with this mode, from (32), can be expressed as(45)εc′=CT-Tc,which explains the Curie-Weiss behaviour of dielectric constant along the c-axis of KH2PO4 crystal in the PE phase observed by Deguchi and Nakamura [37], Busch [4], and Kaminow and Harding [35], shown in Figure 3. For temperature T→Tc, εc′ 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 εc′→0 as T→Tc. The origin of this difference in the temperature dependence of εc′ 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 τ≠0 to εc′ more dominant.There are actually, however, many cases in whichεc′ takes a minimum of finite value at T=Tc 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 H2PO4 group in KDP-type ferroelectrics.


## 3.1. Numerical Calculations

The parameters in our calculation are listed in Table1. The calculated values of Sqz and Sqx for KDP-type crystals for different values of four-body coupling coefficient J′, collective phonon mode frequency ω~~-, transverse dielectric constant εao, observed dielectric constant εco, and tangent loss along a-axis (tan⁡δa) and c-axis (tan⁡δc) 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−1)
J (cm−1)
J
′ (cm−1)
V
i
k (cm−1)
ω
k (cm−1)
μ
×
1
0
18 (cgs)
T
c
(K)
V
/
k
T
c
82
334
440
25.56
153
1.8
123
0.299Table 2
Calculated values ofSqz and Sqx for KDP crystal for different values of four-body coupling coefficient (J′).
S
q
z
S
q
x
Temperature (K)
J
′
=
J
/
3
J
′
=
4
J
/
3
J
′
=
7
J
/
3
J
′
=
J
/
3
J
′
=
4
J
/
3
J
′
=
7
J
/
3
P
r
e
s
e
n
t
s
t
u
d
y
P
r
e
s
e
n
t
s
t
u
d
y
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.392Table 3
Calculated values ofω~~-,εa(0),εc(0),tan⁡δa, and tan⁡δc for KDP-type crystals.
Temperature (K)
Collective phonon mode frequencyω~~- (cm−1) (present study)
Transverse dielectric constantεa(0) (present study)
Observed dielectric constantεc(0) (present study)
Loss tangent along thea-axis (tan⁡δa) (present study)
Loss tangent along thec-axis (tan⁡δc) (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.0250Figure 1
Temperature dependence ofSqz and Sqx for KDP-type crystal for different values of J′. (a) J′=J/3, (b) J′=4J/3, and (c) J′=7J/3 (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 thea-axis (tan⁡δa), for field along the c-axis (tan⁡δc). Present calculation is represented by ‘∘’, and solid line represents experimental results of Kaminow and Harding [35].


## 3.2. Temperature Dependence ofSqz and Sqx for KDP-Type Crystal for Different Values of J′

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 Sqz and Sqx for KDP-type crystal for different values of J′, and variation is shown in Figure 1.In Figure1, the curve (a) is the case of J′=J/3<4J/3, curve (b) is J′=4J/3<7J/3, and curve (c) is J′=7J/3>4J/3. In curves (a) and (b) the value of Sqz 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 Sqz with the temperature starts from a nonzero value Sqz at point “A”; that is to say, when temperature decreases Sqz increases to the saturation value from the finite value of Sqz. This is the case of first-order phase transition. The temperature at point “A” is transition temperature, and the value of Sqz at “A” is the discontinuity of Sqz. The value of Sqx decreases when temperature decreases in the ferroelectric phase. On the other hand, in PE phase the value of Sqx 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 E(x,y) mode, which is responsible for the observed transverse dielectric properties of KDP. In the simplest approximation ω~~+ can be written (ω~~+=K1+K2T), where K1 and K2 are temperature independent parameters. The results for transverse dielectric constant εa′ 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 ω~~+E(x,y)mode is closely related to the macroscopic dielectric constant εa′. This also suggests that the E-mode Raman spectrum originates neither from the second-order Raman scattering nor from the density of states due to the local disorder above Tc but from one of the collective modes at the centre of the Brillouin zone. It should be mentioned here that the low frequency E-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 E-mode collective hydrogen motion has characteristics damping factor which slowly increases as the temperature approaches Tc, while for that of ω~~-B2z 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 observedE-mode Raman spectrum in the following aspects: (i) ω~~+ does not change appreciably as T→Tc in PE phase, (ii) Γp(ω+) is weakly dependent on temperature, and (iii) because of the factor (ω2-ω~~2) in the numerator of (31), the susceptibility derived changes the corresponding spectrum from a simple overdamped form to a more flat one, like the E-mode Raman spectrum of KH2PO4 [27].The observed dielectric constant(εc′) of KH2PO4 along c-axis is shown in Figure 3. The ω~~-B2z mode may be assigned for the observed temperature dependence of εc′. As from (41) ω~~q2∝T-Tc, the real part of the dielectric constant associated with this mode, from (32), can be expressed as(45)εc′=CT-Tc,which explains the Curie-Weiss behaviour of dielectric constant along the c-axis of KH2PO4 crystal in the PE phase observed by Deguchi and Nakamura [37], Busch [4], and Kaminow and Harding [35], shown in Figure 3. For temperature T→Tc, εc′ 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 εc′→0 as T→Tc. The origin of this difference in the temperature dependence of εc′ 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 τ≠0 to εc′ more dominant.There are actually, however, many cases in whichεc′ takes a minimum of finite value at T=Tc 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 H2PO4 group in KDP-type ferroelectrics.



---

*Source: 9475740-2016-03-21.xml*