Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures

A new definition of a dipole glass parameter for ferro-antiferroelectric solid mixtures is suggested. It is constructed on the nearest neighbours pair correlation functions for interacting dipole momenta. The behaviour of the dipole glass parameter is calculated and discussed.

Gespeichert in:
Bibliographische Detailangaben
Datum:2013
Hauptverfasser: Korynevskii, M.A., Solovyan, V.B.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут фізики конденсованих систем НАН України 2013
Schriftenreihe:Condensed Matter Physics
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/121071
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2013. — Т. 16, № 1. — С. 13703:1–8. — Бібліогр.: 35 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-121071
record_format dspace
spelling irk-123456789-1210712017-06-14T03:02:36Z Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures Korynevskii, M.A. Solovyan, V.B. A new definition of a dipole glass parameter for ferro-antiferroelectric solid mixtures is suggested. It is constructed on the nearest neighbours pair correlation functions for interacting dipole momenta. The behaviour of the dipole glass parameter is calculated and discussed. Запропоновано нове означення параметра дипольного скла для сегнето-антисегнетоелектричних твердих сумiшей. Цей параметр будується на парних кореляцiйних функцiях взаємодiючих дипольних моментiв частинок найближчих сусiдiв. Розраховано та дослiджено поведiнку параметра дипольного скла, а також фiзичнi властивостi фази дипольного скла. 2013 Article Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2013. — Т. 16, № 1. — С. 13703:1–8. — Бібліогр.: 35 назв. — англ. 1607-324X PACS: 76.30.K, 77.80Bh DOI:10.5488/CMP.16.13703 arXiv:1303.5582 http://dspace.nbuv.gov.ua/handle/123456789/121071 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description A new definition of a dipole glass parameter for ferro-antiferroelectric solid mixtures is suggested. It is constructed on the nearest neighbours pair correlation functions for interacting dipole momenta. The behaviour of the dipole glass parameter is calculated and discussed.
format Article
author Korynevskii, M.A.
Solovyan, V.B.
spellingShingle Korynevskii, M.A.
Solovyan, V.B.
Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
Condensed Matter Physics
author_facet Korynevskii, M.A.
Solovyan, V.B.
author_sort Korynevskii, M.A.
title Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
title_short Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
title_full Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
title_fullStr Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
title_full_unstemmed Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
title_sort dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures
publisher Інститут фізики конденсованих систем НАН України
publishDate 2013
url http://dspace.nbuv.gov.ua/handle/123456789/121071
citation_txt Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2013. — Т. 16, № 1. — С. 13703:1–8. — Бібліогр.: 35 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT korynevskiima dipoleglassparameterbehaviourforferroantiferroelectricsolidmixtures
AT solovyanvb dipoleglassparameterbehaviourforferroantiferroelectricsolidmixtures
first_indexed 2025-07-08T19:08:41Z
last_indexed 2025-07-08T19:08:41Z
_version_ 1837106959623389184
fulltext Condensed Matter Physics, 2013, Vol. 16, No 1, 13703: 1–8 DOI: 10.5488/CMP.16.13703 http://www.icmp.lviv.ua/journal Dipole glass parameter behaviour for ferro-antiferroelectric solid mixtures M.A. Korynevskii1,2,3, V.B. Solovyan1 1 Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii St., 79011 Lviv, Ukraine 2 Lviv Polytechnic National University, 12 Bandera St., 79013 Lviv, Ukraine 3 Institute of Physics of the University of Szczecin, 15 Wielkopolska St., 70451 Szczecin, Poland Received May 25, 2012, in final form September 24, 2012 A new definition of a dipole glass parameter for ferro-antiferroelectric solid mixtures is suggested. It is con- structed on the nearest neighbours pair correlation functions for interacting dipole momenta. The behaviour of the dipole glass parameter is calculated and discussed. Key words: ferroelectrics, antiferroelectrics, mixed systems, correlation functions, dipole glass PACS: 76.30.K, 77.80.Bh 1. Introduction The problem of a dipole glass state in frustrated ferroelectrics is difficult both in experimental detec- tion and for theoretical interpretation. It still remains actual, because this state of non-ordering system at low temperatures is widely observed in different materials, for example in dielectrics with non-central ions, relaxors, mixed ferro-antiferroelectric compounds [1–4]. Dipole glass state is a characteristic one for structurally complex systems and corresponds to non-equilibrium thermodynamic phases in them. The most important property of the system and a simple criterion for detection of the process of dipole glass transition is a dispersion of dielectric characteristics. Therefore, in experimental investigations, the dynamical methods for the study of real and imaginary parts of dielectric susceptibility behaviour under high and low frequencies are used. Naturally, this state does not depend on the external field effect and must be observed in a static case. On the other hand, due to a competition between different types of interparticle interactions in a com- plex crystalline matrix there is a possibility of different local states coexisting with their non-equilibrium occupation. Those states correspond to a local minima of the free energy of the system with high poten- tial barriers between them. In the thermodynamic limit (N →∞, V →∞, N /V = const), the relaxation time for transition from some state to another one tends to infinity. As far as the system is close to such a local minimum, it remains there for a long time. These systems are not ergodic and their theoretical de- scription is quite complicated. The dispersion of dielectric characteristics of the systems in a dipole glass state is caused by the stimulated transitions from a certain energy state to another one. So, the dynamic behaviour of the investigated system can be illustrated by their microscopic structure, or in other words, by their static correlation functions. The experimental study of dipole glass phases in frustrated ferroelectrics is presented mostly in papers devoted to mixed ferro-antiferroelectric compounds based on the KDP-type crystals [Kn(NH4)1−nH2PO4 , Rbn(NH4)1−nH2PO4 , Rbn(NH4)1−nD2PO4 , Rbn(NH4)1−nH2AsO4]. The most fun- damental of them, in our opinion, are dielectric investigations [5–11]. In all of them, the temperature behaviour of real (ǫ′) and imaginary (ǫ′′) parts of dielectric permittivity at low and high frequencies looks similar for concentrations 0.2 < n < 0.8. Namely, the existence of smooth (ǫ′) or sharp (ǫ′′) peaks at tem- peratures about 50 K indicate some ordering processes in the system. Since no long-range ordered phases © M.A. Korynevskii, V.B. Solovyan, 2013 13703-1 http://dx.doi.org/10.5488/CMP.16.13703 http://www.icmp.lviv.ua/journal M.A. Korynevskii, V.B. Solovyan in this range of concentrations were observed, only an abnormal increase of short-range correlations is possible there. This state was called a dipole glass phase. The EPR [12, 13], NMR [14] and X-ray [15] structural investigations have shown the existence of signif- icant local fields in ferro-antiferroelectric mixtures at low temperatures. Those fields are connected with strong short-range correlations between charged particles (protons). The theoretical study of mixed ferro-antiferroelectric systems is based on the use of a replica method for calculation of the basic thermodynamic functions. This method was first suggested in the theory of diluted magnetic or mixed ferro-antiferromagnetic systems [16–20]. As for mixed hydrogen bonded ferro-antiferroelectric compounds, the Ising model with a random transverse and longitudinal fields was explored. The Gaussian distribution for random interactions and the corresponding fields was used in [21, 22]. A cluster theory ofmixed KDP-type ferroelectrics, including the effect of piezoelectric coupling, was proposed in [23–25]. The phase diagrams, the behaviour of thermodynamic functions, long-range and dipole glass order parameters were obtained and discussed. Most of those results are in good agreement with the experiment. In our previous papers [26–28], a new approach to the study of hydrogen bonded ferro-antiferro- electric solid mixtures of Rbn(NH4)1−nH2AsO4 (RADA) type was suggested. It is based on the concept of a decisive role of microscopic structure of a system (single and pair correlation functions) in the formation of a dipole glass phase. The cluster mean field approximation with replica method was used for their calculation. So, the “concentration-phase transition temperature” phase diagram was obtained in good agreement with the experimental one. The areas of the existence of dipole glass phases were interpreted as multiplication of pair correlation functions at a certain concentration for one of the components at relatively low temperature. Since the pair correlation functions characterize a short-range order in the system, the dipole glass state is determined by different possibilities of short-range arrangement of inter- acting particles, i.e., by their stable and metastable states all over the sample. As for a strict definition of a dipole glass parameter, there appear to be some problems. Usually [16] it is defined as amean value of the dipole momentum square. However, in this case, the multi-stage charac- ter of dipole glass phase and the real interparticle correlations are not taken into account. Moreover, the temperature behaviour of dipole glass parameter does not react adequately upon passing the “freezing” point (the point of the dipole glass phase appearance). This parameter remains non-equal to zero even in the paraelectric phase where no glass state can be observed. A more appropriate definition of this parameter and the static interpretation of a dipole glass phase will be presented in the present paper. 2. Model. Analytical results It is well known that [29] the ferroelectric RbH2AsO4 and antiferroelectric NH4H2AsO4 single crystals are isomorphic with tetragonal lattice D12 2d space group symmetry in paraelectric phase. Below 110 K in the first crystal and below 213 K in the second one, a ferroelectric and antiferroelectric ordering are correspondingly observed. The axis of the orientation of dipole moments in RbH2AsO4 is directed along 0Z , but in NH4H2AsO4 it is directed along 0X (0Y ). Solid mixtures of both crystals due to frustration of ferro- and antiferroelectric types of interactions demonstrate a possibility of a new “ordered” dipole glass phase appearance. In this phase, no type of long-range ordering is preferential and only short-range correlations between particles remain. Dipole glass phase is realized for intermediate concentration (0.50 < n < 0.80) and at relatively low tempera- tures (below 100 K). The physical reason for the existence of these phases are numerous arrangements of nearest neighbours (particles) in a disordered ferro-antiferroelectric system. We will regard a regular KDP-type crystalline lattice, each site of which is randomly occupied by rubidium atoms or ammonium groups, that is by Rb or NH4. The effective Hamiltonian of this mixed system is taken in the form [28, 30]: H = − 1 2 ∑ i , j ∑ mm′ [ V mm′ i j n̂i n̂ j Sz im Sz j m′ +U mm′ i j (1− n̂i )(1− n̂ j )Sim Sx j m′ ] + 1 2 ∑ i ∑ m [ E cosΘn̂i Sz im +E sinΘ(1− n̂i )Sx im ] , (2.1) 13703-2 Dipole glass parameter behaviour where both V mm′ i j and U mm′ i j are the intensities of interaction between m-th dipole particle in i -th site with m′-th dipole particle in j -th site when dipole moments are oriented along z-axis and along x-axis, correspondingly. Sz and Sx are projections of the unit classical vector ~S, the site occupation operators n̂i have the following eigenvalues: ni = { 1, when i -th site is occupied by Rb atom, 0, when i -th site is occupied by NH4 group. (2.2) Here, E is the external electric field, θ is the angle between z-axis and the direction of this field. To take into account the random distribution of dipoles oriented along z-axis (Rb-AsO4 groups) and dipoles oriented along x-axis (NH4-AsO4 groups) we have used the replica method [18] for the configu- ration averaging of thermodynamic functions. Since there are two types of dipoles having ferroelectric and antiferroelectric orientation, we must correspondingly introduce two types of replica variables σk im for Sz im and ξk im for Sx im . The following property of a binomial random variable χc will be used. If the probability of χc is Pr { χc = r } = ( c r ) pr qc−r , then the expectation value of exp { aχc } is equal to ( pea +q )c , which behaves like 1+ c ln ( pea +q ) , when the analytic continuation for c → 0 is performed. In the con- text of the investigated system, the p and q are as follows: p = n, q = 1−n, n is a concentration (density) of Rb atoms. Since the unit cell for KDP-type crystal contains two formula units (there are also two sublattices below the phase transition point [29]), the numbers m, m′ run from 1 to 2. So, we consider the next form for the interactions between all particles in the crystalline lattice: the strict two-particle interaction for particles belonging to the same crystalline lattice site and the self-consistent mean field approximation for particles from different sites, the so-called two-particle cluster approximation [31]. With the accuracy up to the second order of symmetric replica expansion and using two-particle cluster approximation, the following expression for the free energy of the system investigated has been obtained: F = − 1 β ln { 4 [ e−A1 +eA1 cosh ( 2B1 +βE cosΘ )] × [ eA2 cosh ( βE sinΘ ) +e−A2 cosh(2B2) ]} + 1 β C . (2.3) Here, A1 = J1 + J2 { Ṽ 2 f − Ṽ 2 a f + [ 1+ ( Ṽ f − Ṽa f )2 ] gz −Ũ gx } , A2 = ( tanhβV − J1 ) Ũ + J2 { Ũ 2 f −Ũ 2 a f + [ Ũ 2 + ( Ũ f −Ũa f )2 ] gx −Ũ gz } , B1 = [ J1Ṽ f +4J2 ( Ṽ 2 f P 2 −3Ṽ f Ũa f q2 )] P , B2 = [ ( tanhβV − J1 ) Ũa f +4J2 ( Ũ 2 a f q2 −3Ṽ f Ũa f P 2 )] q , (2.4) C = J1Ṽ f P 2 + ( tanhβV − J1 ) Ûa f q2 + J2 2 [ −1−Ũ + ( gz −Ũ gx )2 + ( Ṽ f + Ṽa f )2 + ( Ṽ f − Ṽa f )2 g 2 z + ( Ũ f +Ũa f )2 + ( Ũ f +Ũa f )2 g 2 x +12 ( Ṽ f P 2 −Ũa f q2 )2 ] , Ṽ f = V f V , Ṽa f = Va f V , Ũ f = U f V , Ũa f = Ua f V , V f = 1 2 N ∑ j=1 ( V 11 i j +V 12 i j ) , Va f = 1 2 N ∑ j=1 ( V 11 i j −V 12 i j ) , U f = 1 2 N ∑ j=1 ( U 11 i j +U 12 i j ) , Ua f = 1 2 N ∑ j=1 ( U 11 i j −U 12 i j ) , V =V 12 i i , U =U 12 i i . (2.5) 13703-3 M.A. Korynevskii, V.B. Solovyan J1 = ∞ ∑ t=1 (−1)t 2ln ( 1−n+ne−2βV t ) , J2 = ∞ ∑ t=1 (−1)t 4t ln ( 1−n+ne−2βV t ) , (2.6) β= (kT )−1, k is the Boltzmann constant, T is the absolute temperature. The form of J1(n,T ) and J2(n,T ) functions is determined by the binomial random distribution of different types of “particles” (Rb or NH4) and the replica method of configuration averaging. J1(n,T ) is proportional to the concentration n, but J2(n,T ) reflects the correlation between both subsystems and naturally tends to zero in the limits n → 0 and n → 1. Ferroelectric P , antiferroelectric q order parameters and gz , gx nearest neighbours pair correlation functions are defined in the following way: P =<σi1 >, q =< ξi1 >, gz =<σi1σi2 >, gx =< ξi1ξi2 >, (2.7) where < . . . > denotes thermal averaging for a given distribution of z- and x-oriented dipole momenta of Rb-AsO4 and NH4-AsO4 groups, respectively, whereas a bar denotes a stochastic averaging over different random distributions. From the principle of stationarity for the free energy (2.3) under order parameters P , q and both pair correlation functions gz and gx (they determine a structure of nearest neighbours particles) the set of equations is obtained: P = sinh(2B1 +βE cosΘ) cosh(2B1 +βE cosΘ)+e−2A1 , q = sinh(2B2) cosh(2B2)+e2A2 cosh(βE sinΘ) , gz = −1+e2A1 cosh(2B1 +βE cosΘ) 1+e2A1 cosh(2B1 +βE cosΘ) , gx = cosh(βE sinΘ)−e−2A2 cosh(2B2) cosh(βE sinΘ)+e−2A2 cosh(2B2) . (2.8) Spontaneous values (in zero external field E ) of P, q, gz , gx parameters are limited by±1. So, naturally, they may be interpreted as a set of order parameters of the system under investigation. 3. Phase diagram. Dipole glass parameter Computing the equations (2.8) one can find the dependencies of P , q , gz , gx parameters on tem- perature T and concentration n. Those dependencies form a phase diagram of the investigated ferro- antiferroelectric mixed system in the plane n −T . The principal outline of a phase diagram for RADA systemwas obtained in the paper [28], but due to its considerable role in the analysis of the mixed system behaviour, we reproduced it here (figure 1). The main features of the phase diagram are: 1) the existence of two non-crossing areas for ferro- electric (F) and antiferroelectric (AF) phases; 2) the existence of pure dipole glass phase at intermediate concentration and relatively small temperature (this phase is divided into several (G3-G9) subphases ac- cording to the order of the chaotic state); 3) the possibility of mixed ferroelectric-dipole glass (F-G) and antiferroelectric-dipole glass (AF-G) phases existence. It should be noted that the upper dipole glass phase G3 is rather a precursor of the real dipole glass phase appearance. This conclusion is in good agreement with the experiment [32, 33]. The difference between RADA phase diagram in [28] and the one presented here is the line divid- ing the paraelectric phase and the antiferroelectric phase (P-AF). In figure 1 the upper line represents the points of stable non-zero root arising for antiferroelectric order parameter q while in [28] the cor- responding line indicates the points of its bifurcation (lower line in figure 1). It should be noted that the 13703-4 Dipole glass parameter behaviour Figure 1. Extended phase diagram of RADA type compound built for the model Hamiltonian parameters: V f = 70 K, U f =−40 K, Va f =−80 K, Ua f = 120 K, V = 65 K, U =−210 K. only two fitting points were used to build figure 1. The first one is the antiferroelectric phase transition point temperature for pure NH4H2AsO4 crystal, and the second one is the ferroelectric phase transition point temperature for pure RbH2AsO4 crystal. The dipole glass phase is interpreted as a possibility of different arrangements for the particles of nearest neighbours. This fact is reflected by numerous solutions for pair correlation functions gz and gx at different values of temperature and concentration. When the temperature decreases, the dipole glass state becomes more complicated. In this context, the behaviour of correlation functions gz , gx is characteristic of the dipole glass struc- ture. However, there are a some peculiarities connected with their relation to dipole glass state parame- ter. Usually [34], this parameter is introduced to ferroelectrics theory by analogy with spin glass param- eter in ferromagnetic systems [16], named Edwards and Anderson parameter. Edward and Anderson parameter is defined as a second moment of magnetic order parameter. This moment is maximum for low temperature and vanishes at relatively large temperature Tg named glassy temperature. As for ferro- electric systems, the situation is not that simple. The recent investigations [25] of mixed hydrogen bonded ferro-antiferroelectric systems have shown some restriction to using the Edwards and Anderson param- eter for a complete definition of a dipole glass state. The second moment of ferroelectric order parameter does not vanish even for very high temperatures when no glass state can exist. Pair correlation functions gz and gx , which are analogous to the second moment of order parameter, are nonzero even in paraelectric phase P (see figure 1) and tend to zero for very high temperatures (theoretically limT→∞ gz = limT→∞ gx = 0). According to our concept, the dipole glass phase is a state with numerous types of short-range correlations between particles, or numerous numbers of solutions for a set of equations (2.8). The intensity of correlations (values of gz and gx ) is important but not decisive for dipole glass phase characterization. Dipole glass phase arises only at the moment when the number of different solutions for gz or gx changes from 1 to a bigger number. Such a situation is observed on the upper line of G3 area (figure 1). As we have calculated (the phase diagram has been built on this basis), there is only one solution for (gz , gx ) pair in the paraelectric phase (P ), three solutions for (gz , gx ) pair in the G3 area, five solutions for (gz , gx ) pair in the G5 area and so on. With an increase of the number of different solutions for gz , gx , a bigger number of different types of correlations of the nearest particles are observed. As a result, a more complicated arrangement of dipoles all over the crystal takes place. This situation is accompanied by a more chaotic distribution of correlations and by the formation of a dipole glass phase. When temperature decreases, a higher stage of dipole glass phase is realized. Thus, dipole glass parameter should take into account two factors: the number of solutions for gz , gx parameters and their absolute values. It is obvious that the unique (trivial) solution for gz , gx , which remains in the 13703-5 M.A. Korynevskii, V.B. Solovyan Figure 2. Temperature dependencies of the “partial” dipole glass parameters Gz and Gx for different concentrations. paraelectric phase must be eliminated. Therefore, since both gz and gx functions form a dipole glass parameter G , we propose the next expression for its calculation: G = 1 2M0 m−1 ∑ i=1 [ n|g (i) z |+ (1−n)|g (i) x | ] , (3.1) where g (i) z , g (i) x are the roots of a set of equations (2.8), m is the number of these roots for any temperature and concentration, N0 is the maximum number of roots, which are taken into account (in the presented calculations M0 = 9), n is the concentration. The “partial” dipole glass parameters Gz and Gx : Gz = 1 M0 m−1 ∑ i=1 |g (i) z | , Gx = 1 M0 m−1 ∑ i=1 |g (i) x | (3.2) describe the power of dipole glass state intensity according to z − z and x − x correlation of Rb-AsO4 and NH4-AsO4 dipoles, correspondingly. Figure 3. Temperature dependencies of the dipole glass parameters G for different concentrations. 13703-6 Dipole glass parameter behaviour The temperature dependencies for Gz and Gx at different concentrations are presented in figure 2, and for G they are presented in figure 3. The most characteristic feature in the behaviour of Gz , Gx and G are their step-like dependence. The multigraded character of a dipole glass parameter is caused by the complex nature of dipole glass phase to exactly change the number of roots for gz , gx correlation func- tions (see figure 1). The position of “jumps” for Gz , Gx , G is determined by the lines dividing the areas of P-G3, G3-G5, . . . for “pure” dipole glass phases and F-G3, G3-G4, . . . ; AF-G3, G3-G5, . . . for the mixed or- dered and dipole glass phases. The existence of those lines is determined by the finite number of different roots of equations (2.8) in every area of temperature and concentration. Theoretically, this number is in- finite [32, 35], so the step-like behaviour of dipole glass parameter in these limits may transform into a quasi-continuous curve. 4. Conclusions A new approach to the calculation of the dipole glass parameter in mixed ferro-antiferroelectric sys- tems has been developed. It is based on the conception of a decisive role of the pair correlations of the nearest neighbours particles in a crystalline lattice for the dipole glass phase formation. The cor- responding correlation function has been calculated in the second order of replica expansion and using the cluster approximation. This microscopic approach leads to the basic physical statement: the dipole glass phase is a polystate arrangement of the nearest neighbours particles with the possibility of non- equilibrium occupation of them. In such a way, the proposed definition differs from the well known Edwards and Anderson spin glass parameter in the quenched ferromagnetic system theory. The calculated dipole glass parameter takes into account not only the intensity of interparticle correla- tions but also the multiplicity of a different type of those correlations. The latter corresponds to different possible types of short-range order in the mixed system and usually describes a metastable states. The to- tal dipole glass parameter is formed both by ferro-ferro and antiferro-antiferro types of correlations. The step-like behaviour of the dipole glass parameter reflects the non-uniform structure of the dipole glass phase and the dependence of the intensity of the disorder on temperature and concentration. When the number of different roots for pair correlation functions increases, the temperature dependence of this parameter becomes more smooth. The obtained phase diagram “temperature-concentration” demonstrates the possibility of the exis- tence of both pure ferroelectric, antiferroelectric, dipole glass phases and mixed ferroelectric-dipole glass, antiferroelectric-dipole glass phases. This fact was recently confirmed by numerous experimental investigations. The precursors of the dipole glass phase which appears at rather high temperatures are interpreted as the state with the finite number of different types of interparticle correlations, in counter- balance to paraelectric state with only one type of those correlations. References 1. Vugmeister B.E., Glinchuk M.D., Usp. Fiz. Nauk, 1985, 146, 459 (in Russian); doi:10.3367/UFNr.0146.198507c.0459 [Sov. Phys. Uspekhi, 1985, 28, 589; doi:10.1070/PU1985v028n07ABEH003872]. 2. Korotkov L.N., Shuvalov L.A., Crystallogr. Rep., 2004, 49, 832; doi:10.1134/1.1803315. 3. Rout D., Subramanian V., Hariharan K., Sivasubramanian V., J. Phys. Chem. Solids, 2006, 67, 1629; doi:10.1016/j.jpcs.2006.02.009. 4. Trybula Z., Stankowski J., Condens. Matter Phys., 1998, 1, 311. 5. Courtens E., J. Physique Lett., 1982, 43, L199; doi:10.1051/jphyslet:01982004306019900. 6. Samara G.A., Phys. Rev. Lett., 1987, 59, 347; doi:10.1103/PhysRevLett.59.347. 7. Kim S., Kwun S., Phys. Rev. B, 1990, 42, 638; doi:10.1103/PhysRevB.42.638. 8. Trybula Z., Schmidt V.H., Drumheller J.E., He D., Li Z., Phys. Rev. B, 1989, 40, 5289; doi:10.1103/PhysRevB.40.5289. 9. Kutnjak Z., Pire R., Levstik A., Levstik I., Filipic C., Blinc R., Phys. Rev. B, 1994, 50, 12421; doi:10.1103/PhysRevB.50.12421. 10. Gridnev S.A., Korotkov L.N., Shuvalov L.A., Rogova S.P., Fedosyuk R.M., Ferroelectrics Lett., 1991, 13, 67; doi:10.1080/07315179108203323. 11. Kwan O.J., Kim J.J., Phys. Rev. B, 1993, 48, 6639; doi:10.1103/PhysRevB.48.6639. 13703-7 http://dx.doi.org/10.3367/UFNr.0146.198507c.0459 http://dx.doi.org/10.1070/PU1985v028n07ABEH003872 http://dx.doi.org/10.1134/1.1803315 http://dx.doi.org/10.1016/j.jpcs.2006.02.009 http://dx.doi.org/10.1051/jphyslet:01982004306019900 http://dx.doi.org/10.1103/PhysRevLett.59.347 http://dx.doi.org/10.1103/PhysRevB.42.638 http://dx.doi.org/10.1103/PhysRevB.40.5289 http://dx.doi.org/10.1103/PhysRevB.50.12421 http://dx.doi.org/10.1080/07315179108203323 http://dx.doi.org/10.1103/PhysRevB.48.6639 M.A. Korynevskii, V.B. Solovyan 12. Cevc P., Zalar B., Blinc R., Solid State Commun., 1989, 70, 461; doi:10.1016/0038-1098(89)91080-6. 13. Grinberg E.S., Izotov V.V., Efimov V.N., Ferroelectrics Lett., 1993, 15, 61; doi:10.1080/07315179308204240. 14. Slak J., Kind R., Blinc R., Courtens E., Zumer S., Phys. Rev. B, 1984, 30, 85; doi:10.1103/PhysRevB.30.85. 15. Hayase S., Sakashita H., Teranchi H., Ferroelectrics, 1987, 72, 245; doi:10.1080/00150198708017949. 16. Edwards S.F., Anderson P.W., J. Phys. E, 1975, 5, 965; doi:10.1088/0305-4608/5/5/017. 17. Thouless D.J., Anderson P.W., Palmer R.G., Philos. Mag., 1977, 35, 593; doi:10.1080/14786437708235992. 18. Parisi C., Phys. Rev. Lett., 1979, 43, 1754; doi:10.1103/PhysRevLett.43.1754. 19. Kirkpatrick S., Sherrington D., Phys. Rev. B, 1978, B17, 4384; doi:10.1103/PhysRevB.17.4384. 20. Bray A.J., Moore M.A., J. Phys. C, 1980, 13, L469; doi:10.1088/0022-3719/13/19/002. 21. Blinc R., Dolinsek J., Zalar B., Z. Phys. B, 1997, 104, 629; doi:10.1007/s002570050499. 22. Kim D.H., Kim J.J., Ferroelectrics, 2002, 268, 263; doi:10.1080/713715957. 23. Levitskii R.R., Sorokov S.I., Stankowski J., Trybula Z., Condens. Matter Phys., 2008, 11, 523. 24. Sorokov S.I., Levitskii R.R., Vdovych A.S., Korotkov L.N., Ferroelectrics, 2010, 397, 43; doi:10.1080/00150193.2010.484717. 25. Levitskii R.R., Sorokov S.I., Zachek I.R., Vdovych A.S., Moina A.P., Korotkov L.M., Bocharov A.I., Preprint of the Institute for Condensed Matter Physics, ICMP-11-13E, Lviv, 2011. 26. Korynevskii N.A., Solovyan V.B., Phase Transitions, 2007, 80, 55; doi:10.1080/01411590601092712. 27. Korynevskii N.A., Solovyan V.B., Condens. Matter Phys., 2009, 12, 267; doi:10.5488/CMP.12.2.267. 28. Korynevskii N.A., Solovyan V.B., Physica B, 2010, 405, 2232; doi:10.1016/j.physb.2010.02.017. 29. Lines M.E., Glass A.M., Principles and Applications of Ferroelectrics and Related Materials, Clarendon, Oxford, 1979. 30. Korynevskii N.A., Baran O.R., Ferroelectrics, 2004, 300, 151; doi:10.1080/00150190490446458. 31. Strukov B.A., Levaniuk A.P., Physical Principles of Ferroelectric Phenomena in Crystals, Nauka, Moskow, 1983 (in Russian). 32. Buixaderas E., Kamba S., Petzelt J., Ferroelectrics, 2004, 308, 131; doi:10.1080/00150190490508909. 33. Trybula Z., Kaszynska K., Trybula M., Levitskii R.R., Abstracts of IX Polish-Ukrainian Meeting and XXIX Interna- tional School on Ferroelectric Physics (Krakow, Poland, 2008), Pub. Krakow Acad. Pedagogical, 2008, p. 24. 34. Bidaux R., Carton J.P., Sarma J., J. Phys. A: Math. Gen., 1976, 9, L87; doi:10.1088/0305-4470/9/8/002. 35. Moskalenko V.A., Vladimir M.I., Dogotar’ L.A., Kozhukar’ S.P., In: “Problems of Modern Statistic Physics”, Naukova Dumka, Kyiv, 1985 (in Russian). Поведiнка параметра дипольного скла для сегнето-антисегнетоелектричних твердих сумiшей М.А.Кориневський1,2,3, В.Б.Солов’ян1 1 Iнститут фiзики конденсованих систем НАН України, вул. I. Свєнцiцького, 1, 79011 Львiв, Україна 2 Нацiональний унiверситет “Львiвська Полiтехнiка”, вул. С. Бандери 12, 79013 Львiв, Україна 3 Iнститут фiзики, Щецiнський унiверситет, вул. Вєлькопольска 15, 70451 Щецiн, Польща Запропоновано нове означення параметра дипольного скла для сегнето-антисегнетоелектричних твер- дих сумiшей. Цей параметр будується на парних кореляцiйних функцiях взаємодiючих дипольних мо- ментiв частинок найближчих сусiдiв. Розраховано та дослiджено поведiнку параметра дипольного скла, а також фiзичнi властивостi фази дипольного скла. Ключовi слова: сегнетоелектрики, антисегнетоелектрики, змiшанi системи, кореляцiйнi функцiї, дипольне скло 13703-8 http://dx.doi.org/10.1016/0038-1098(89)91080-6 http://dx.doi.org/10.1080/07315179308204240 http://dx.doi.org/10.1103/PhysRevB.30.85 http://dx.doi.org/10.1080/00150198708017949 http://dx.doi.org/10.1088/0305-4608/5/5/017 http://dx.doi.org/10.1080/14786437708235992 http://dx.doi.org/10.1103/PhysRevLett.43.1754 http://dx.doi.org/10.1103/PhysRevB.17.4384 http://dx.doi.org/10.1088/0022-3719/13/19/002 http://dx.doi.org/10.1007/s002570050499 http://dx.doi.org/10.1080/713715957 http://dx.doi.org/10.1080/00150193.2010.484717 http://dx.doi.org/10.1080/01411590601092712 http://dx.doi.org/10.5488/CMP.12.2.267 http://dx.doi.org/10.1016/j.physb.2010.02.017 http://dx.doi.org/10.1080/00150190490446458 http://dx.doi.org/10.1080/00150190490508909 http://dx.doi.org/10.1088/0305-4470/9/8/002 Introduction Model. Analytical results Phase diagram. Dipole glass parameter Conclusions