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

The soft mode dynamical model has been used to study dielectric properties and ultrasonic attenuation in KDP-type ferroelectric crystals.Themodel Hamiltonian proposed by Blinc and Zeks has beenmodified by considering lattice anharmonicity up to fourthorder. 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 2 PO 4 group in KDP-type ferroelectrics. By fitting model values of physical quantities, the temperature dependence of ⟨Sz q ⟩ and ⟨Sx q ⟩ 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.


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 2 PO 4 ), 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 2 PO 4 ) 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][6][7][8][9][10].
Kaminow and Damen [11] first observed the soft mode associated with the ferroelectric phase transition of the KDPtype crystal at 122.3 K by measuring the low frequency Raman scattering in () configuration.Since then  2 () 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][15][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  2 () 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  2 (), a doubly degenerate (, ), and  2 modes [18].(, ) mode, which is both infrared and Raman active, provides valuable information, because this mode and  2 () mode 2 Physics Research International reflect the nature of collective proton motion which triggers the phase transition, and, moreover, in contrast to  2 (), 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 fourthorder 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.

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 fourbody coupling coefficient, and    refers to the four hydrogen bonds in the PO 4 group in KDP.In the last fourth terms   is bare phonon frequency,   and   are displacement and momentum operators,   is proton-lattice interaction term, and  3 and  4 are the third-and fourth-order anharmonic coefficients.

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

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 Physics Research International with where   = (  /ω  ) coth(ω  /2) is occupation number and  = (  ) −1 ,   being Boltzmann's constant and  being the absolute temperature.Calculating (20) self-consistently and approximating, the collective mode frequency is given by Equation (22a) is obtained using model Hamiltonian (1), and (22b) and (22c) 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 3 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 ⟨   ⟩ = ⟨   ⟩ = 0, and so in the ferroelectric phase ( <   ) [27].Consider ) . (28) 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 (tan ) 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   ≪ 1 for KDP crystals), (31) can be reduced to where ω2 is given by (20) and ω is given by (20a).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.

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 ω2 − is proportional to (  − ) for  ≤   .the effect of factor (1+ 2  2  ) −1 is to smooth out the discontinuity.The ultrasonic attenuation peaks for   ≅ 1 and   (i.e., ω− ) goes to zero as  →   .For temperature such as   ≪ 1, 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  = 0,  =   , 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 ⟨   ⟩ = 0, 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 ( ω− ), transverse dielectric constant {  ()}, observed dielectric constant {  ()}, and tangent loss along -axis (tan   ) and -axis (tan   ) for KDP-type crystals are listed in Tables 2  and 3. Their variations with temperature are shown in Figures 1-4.

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   = /3 < 4/3, curve (b) is   = 4/3 < 7/3, and curve (c) is   = 7/3 > 4/3.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.
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].corresponds to transverse (, ) mode, which is responsible for the observed transverse dielectric properties of KDP.In the simplest approximation ω+ can be written ( ω+ =  1 +  2 ), where  1 and  2 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 ω+ {(, ) mode} is closely related to the macroscopic dielectric constant    .This also suggests that Physics Research International the -mode Raman 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 ω− { 2 ()} 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 ( 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 -mode Raman spectrum of KH 2 PO 4 [27].
The observed dielectric constant (   ) of KH 2 PO 4 along -axis is shown in Figure 3.The ω− { 2 ()} mode may be assigned for the observed temperature dependence of    .As from (41) ω2  ∝ (−  ), 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 2 PO 4 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    → 0 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  ̸ = 0 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.

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 KDPtype 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 fourparticle cluster Hamiltonian.The dielectric properties and ultrasonic attenuation strongly depend on the relaxational
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

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