What can we learn about ferroelectrics using methods of nonlinear dynamics?

The nonlinear series resonance circuit with ferroelectric capacitor has been extensively investigated. If the ferroelectric within the capacitor is in its polar phase many of the features well known from model systems of nonlinear dynamics may be observed. These characteristics are the shift of r...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2003
Автор: Diestelhorst, M.
Формат: Стаття
Мова:English
Опубліковано: Інститут фізики конденсованих систем НАН України 2003
Назва видання:Condensed Matter Physics
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/120757
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:What can we learn about ferroelectrics using methods of nonlinear dynamics? / M. Diestelhorst // Condensed Matter Physics. — 2003. — Т. 6, № 2(34). — С. 189-196. — Бібліогр.: 13 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-120757
record_format dspace
spelling irk-123456789-1207572017-06-13T03:06:09Z What can we learn about ferroelectrics using methods of nonlinear dynamics? Diestelhorst, M. The nonlinear series resonance circuit with ferroelectric capacitor has been extensively investigated. If the ferroelectric within the capacitor is in its polar phase many of the features well known from model systems of nonlinear dynamics may be observed. These characteristics are the shift of resonance frequency with increasing driving voltage, bifurcations and chaotic behaviour. Considerations of the system on the basis of simple Landautheory suggest to describe the resonance circuit by means of a Duffingequation. Because of the switching process during the nonlinear vibrations the real situation is more complicated as can be shown by nonlinear time series analysis [1]. Об’єктом дослідження у даній роботі є нелінійний послідовний резонансний контур з сегнетоелекричним конденсатором. Якщо сегнетоелектрик у конденсаторі знаходиться у полярному стані, спостерігаються властивості, відомі для модельних систем у нелінійній динаміці. Це, зокрема, зсув резонансної частоти зі збільшенням прикладеної напруги, біфуркації та хаотична поведінка. Розгляд системи в рамках простої теорії Ландау веде до опису резонансного контура за допомогою рівняння Дуффінга. Насправді ж, як показано з допомогою аналізу нелінійних часових послідовностей [1], процеси перемикання підчас нелінійних коливань значно ускладнюють ситуацію. 2003 Article What can we learn about ferroelectrics using methods of nonlinear dynamics? / M. Diestelhorst // Condensed Matter Physics. — 2003. — Т. 6, № 2(34). — С. 189-196. — Бібліогр.: 13 назв. — англ. 1607-324X PACS: 05.45.a, 05.45.Tp, 77.80.-e, 77.84.Fa DOI:10.5488/CMP.6.2.189 http://dspace.nbuv.gov.ua/handle/123456789/120757 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description The nonlinear series resonance circuit with ferroelectric capacitor has been extensively investigated. If the ferroelectric within the capacitor is in its polar phase many of the features well known from model systems of nonlinear dynamics may be observed. These characteristics are the shift of resonance frequency with increasing driving voltage, bifurcations and chaotic behaviour. Considerations of the system on the basis of simple Landautheory suggest to describe the resonance circuit by means of a Duffingequation. Because of the switching process during the nonlinear vibrations the real situation is more complicated as can be shown by nonlinear time series analysis [1].
format Article
author Diestelhorst, M.
spellingShingle Diestelhorst, M.
What can we learn about ferroelectrics using methods of nonlinear dynamics?
Condensed Matter Physics
author_facet Diestelhorst, M.
author_sort Diestelhorst, M.
title What can we learn about ferroelectrics using methods of nonlinear dynamics?
title_short What can we learn about ferroelectrics using methods of nonlinear dynamics?
title_full What can we learn about ferroelectrics using methods of nonlinear dynamics?
title_fullStr What can we learn about ferroelectrics using methods of nonlinear dynamics?
title_full_unstemmed What can we learn about ferroelectrics using methods of nonlinear dynamics?
title_sort what can we learn about ferroelectrics using methods of nonlinear dynamics?
publisher Інститут фізики конденсованих систем НАН України
publishDate 2003
url http://dspace.nbuv.gov.ua/handle/123456789/120757
citation_txt What can we learn about ferroelectrics using methods of nonlinear dynamics? / M. Diestelhorst // Condensed Matter Physics. — 2003. — Т. 6, № 2(34). — С. 189-196. — Бібліогр.: 13 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT diestelhorstm whatcanwelearnaboutferroelectricsusingmethodsofnonlineardynamics
first_indexed 2025-07-08T18:31:44Z
last_indexed 2025-07-08T18:31:44Z
_version_ 1837104633345998848
fulltext Condensed Matter Physics, 2003, Vol. 6, No. 2(34), pp. 189–196 What can we learn about ferroelectrics using methods of nonlinear dynamics? M.Diestelhorst Martin-Luther-University, Department of Physics, Friedemann-Bach-Platz 6, D-06114 Halle, Germany Received September 3, 2002 The nonlinear series resonance circuit with ferroelectric capacitor has been extensively investigated. If the ferroelectric within the capacitor is in its polar phase many of the features well known from model systems of nonlinear dynamics may be observed. These characteristics are the shift of reso- nance frequency with increasing driving voltage, bifurcations and chaotic behaviour. Considerations of the system on the basis of simple Landau- theory suggest to describe the resonance circuit by means of a Duffing- equation. Because of the switching process during the nonlinear vibrations the real situation is more complicated as can be shown by nonlinear time series analysis [1]. Key words: ferroelectricity, domains, TGS, nonlinear dynamics, Duffing-equation, chaos PACS: 05.45.a, 05.45.Tp, 77.80.-e, 77.84.Fa 1. Introduction Switching processes in ferroelectrics are caused by electric fields of magnitudes above the coercive field strength EC. One of the characteristics of switching is the occurrence of hysteresis-loops which can be observed below the so-called Curie- temperature TC. They are accompanied by large dielectric nonlinearities [2]. The description of these switching processes is one of the problems to be solved in the physics of ferroelectricity and in the other fields of condensed matter physics. On the other hand, great efforts have been made during the recent years to analyse nonlinear dynamic systems [3]. So, the motivation of our work is twofold. Using ferroelectric materials like TGS in its polar phase as a capacitor in a series resonance circuit we have a nonlinear dynamic system that makes it possible to test the predictions of the results of the theory of nonlinear dynamics, e.g. bifurcations and chaotic behaviour. Another problem is the very switching process of the ferroelectrics underlying these bifurcations. Therefore, we are interested in answering the question whether it is possible to draw some conclusions concerning the ferroelectric switching by analysing the experimentally recorded time series of our nonlinear resonance circuit. c© M.Diestelhorst 189 M.Diestelhorst Figure 1. The nonlinear series resonance circuit with ferroelectric capacitor. 2. Experimental setup The nonlinear series resonance circuit is shown in figure 1. It consists of a linear inductance L, and a nonlinear capacitance CNL, which is a ferroelectric crystal. In our special case we used triglycine sulfate (TGS) as ferroelectric. The driving voltage Uext = U0 cos(ωt + φ) may be varied both in amplitude and frequency to force the circuit in order to exhibit bifurcations and chaotic behaviour. The linear capacitor C and the resistor R have been added to derive signals proportional to the charge of the TGS-capacitor and the current of the circuit, respectively. To minimize the influence of these two elements it is necessary to fulfill the condition R � 1/C � 1/CNL which reminds of the Sawyer-Tower circuit. So, two signals UC ∼ D and UR ∼ Ḋ = j may be derived [1,4], where D and j are the dielectric displacement and the current density, respectively. These signals are recorded by a digital storage oscilloscope and can be used for direct observation of the phase portrait as demonstrated in figure 2. Alternatively, it is possible to record only the signal UC and to calculate the temporal derivatives numerically [1,5]. According to Kirchhoff’s laws and neglecting the small voltages UC and UR, the behaviour of the circuit may be written by the following equation: Uext = UL + URV + UNL , (1) where UL and UNL are the voltage drops in the inductance and the ferroelectric, respectively, and URV is a voltage drop both across the ferroelectric and the induc- tance. From equation (1) one gets a second order differential equation [1,5,6], that is the equation of motion in the sense of the theory of dynamic systems: D̈ = − h LS ENL(D, Ḋ) − RV L Ḋ + Uext(t) LS . (2) Here, D, Ḋ, D̈ are the dielectric displacement and its temporal derivatives, ENL is the electric field strength in the ferroelectric which is assumed to depend on both D 190 Learn ferroelectrics using methods of nonlinear dynamics Figure 2. Experimentally observed phase portrait constructed from the two time series D(t) and j(t) (both in arbitrary units). and Ḋ in general case. S and h are the electroded surface and the thickness of the ferroelectric sample, respectively. For a further analysis, time series of D and Ḋ were recorded at different tem- peratures of the ferroelectric (see figure 3). The ferroelectric was a rectangular TGS sample with thickness h = 0.28 mm and the electroded surface S = 7.31×4.66 mm2. The other parameters were as follows: inductance L = 100 mH, amplitude of the si- nusoidal driving voltage U0 = 15.5 V and frequency f = 2369 Hz. The temperatures were Ta = 311.65 K, Tb = 308.67 K and Tc = 306.67 K. 3. Analysis of the experimental data As mentioned above, this relatively simple circuit is interesting both from the viewpoint of nonlinear dynamics and ferroelectricity. In our case we can combine these fields and hope to get additional information for the characterization of the ferroelectric switching process using the methods proposed in order to extract the parameters of the equation of motion. For a review of such methods see [3] and references therein. As is well known, many properties of ferroelectrics in its polar phases are dom- inated by the effect of the domain structure (see e.g. [2,7,8]). Nevertheless, already the simple Landau-theory of second order phase transitions motivates the consider- ation of the resonance circuit with ferroelectric TGS-capacitor as an experimental realization of a nonlinear dynamic system. Using the relation ENL = α0(T − TC)D + γD3 equation (2) yields directly: D̈ = − h LS (α0(T − TC)D + γD3) − RV L Ḋ + Uext(t) LS . (3) This equation is the so-called Duffing-equation which is one of the model equations of nonlinear dynamics [9,10]. An interesting feature of equation (3) is that it represents a single-well oscillator above TC and a double-well oscillator if the temperature of the 191 M.Diestelhorst Figure 3. Three phase portraits which were recorded experimentally at differ- ent temperatures of the TGS-crystal [5]. Figure 4. One-step prediction errors χn for the investigated phase portraits [5]. 192 Learn ferroelectrics using methods of nonlinear dynamics ferroelectric is below TC. A wide variety of theoretical results for that model could be confirmed experimentally in our circuit. So, it was possible to measure the shift of resonance frequency with the increasing amplitude of the driving voltage [11] as well as period-doubling bifurcations and chaotic behaviour by tuning the parameters of the circuit such as temperature of the ferroelectric, amplitude and frequency of the driving voltage [4,6]. But more detailed investigations revealed that the situation is more complicated [12]. So, it is not possible to describe the behaviour of the circuit with TGS in its ferroelectric phase by parameters α0 and γ determined from the shift of resonance frequency at temperatures of the same TGS capacitor above TC. The reason is that linear and nonlinear dielectric properties in the paraelectric phase are not affected by domain switching processes. This is a confirmation of the well-known fact that Landau-theory is not capable of describing the switching [2]. The task is now to modify the model such that it becomes capable of characteriz- ing the nonlinear resonance circuit in the polar region of the ferroelectric. Therefore, we tested the possibilities of time series analysis of chaotic signals from the view- point of an experimentalist which may be a little questionable at first glance. But as we have checked even in the regime of chaotic vibrations, the behaviour of the circuit is reproducible. It means that bifurcations, chaos and periodic windows [13] occur at the same control parameters (temperature of the ferroelectric, amplitude and frequency of the driving voltage) even when the ferroelectric is heated above the phase-transition temperature and cooled back to the initial temperature. A first attempt was made to get a better agreement between the experimental data and a model equation of extended Duffing-type: D̈ = n∑ i=0 aiD i − RV L Ḋ ︸ ︷︷ ︸ F (D,Ḋ) + Uext(t) LS , (4) where the degree of the polynomial in D was increased up to n = 9 [1]. What follows is a short description of the fitting procedure. For a detailed outline see [1,5]. Having recorded a long time series of UC ∼ D (129907 data points Dn per time series in our experiments) it is possible to estimate the temporal derivatives Ḋn and D̈n from adjacent data points of the time series Dn according to the following formulas: Ḋn = 1 60τ (Dn+3 − 9Dn+2 + 45Dn+1 −45Dn−1 + 9Dn−2 − Dn−3) + O(τ 5), D̈n = 1 180τ 2 (2Dn+3 − 27Dn+2 + 270Dn+1 −490Dn + 270Dn−1 − 27Dn−2 + 2Dn−3) + O(τ 5). (5) Having calculated the derivatives for each point of time tn = nτ we have quadruples (Dn, Ḋn, D̈n, n) which are used to determine the parameters in the ansatz equa- 193 M.Diestelhorst tion (4) by a least-squares fit which minimizes: χ2 = N∑ n=1 ( D̈n − F (Dn, Ḋn) − Uext(nτ) LS )2 = N∑ n=1 χ2 n . (6) The ansatz equation (3) yielded a relative error of the fit εN of 30%–27% for N = 5−9, where εN is calculated as εN = 100 √ χ2 √ ∑N n=1 D̈2 n . (7) A further improvement of the fit was achieved with a more general ansatz for the function F (D, Ḋ) (see underbrace in equation (4)). Taking into account possible couplings between D and Ḋ which may be physically motivated by various attempts of modelling the nucleation processes [2] we made an ansatz of a complete bivariate polynomial with the terms up to the power k in D and l in Ḋ, including all mixed terms: F (D, Ḋ) = k∑ i=0 l∑ j=0 αijD iḊj. (8) It could be shown ([5]) that the relative error εN saturates at approximately 8.5% for k > 12 and l > 5. A further increase of the degree of the polynomial F (D, Ḋ) to k = 15 and l = 7 gave no significant reduction of the error. But the model contains a lot of parameters (128 αij for the basis functions DiḊj plus 2 parameters for amplitude and phase of the driving voltage Uext). So, the question arises whether it is possible to reduce the number of coefficients without significantly increasing the relative error εNv omitting the terms which give only a small contribution to the model. The first step towards such a simplified model results from the knowledge of the symmetry of the system. From the symmetry of the driving voltage Uext(t) = −Uext(t + T/2), where T is the period of the driving, it can be concluded that the attractor in the phase space X = (D, Ḋ) has the same symmetry or there exist two attractors related to each other by the relation XA1 (t) = −XA2 (t + T/2). Taking into account this symmetry we can neglect all the basis functions DiḊj where i + j equals an even number except α00 to describe an asymmetry (e.g. inner field) in the ferroelectric. The result of this reduction was a small increase of the relative errors εNodd ≈ εNall +0.5% for the three different cases. Further reduction of the number of the remaining 65 αij was achieved by a modified backward elimination algorithm (for details see [5]). Nevertheless, it was necessary to keep about 20–30 basis functions to avoid a dramatic increase of the error. As a result of the procedure, we present the numerically calculated phase portraits for the “optimal” models for the three experimental time series (figure 4). Obviously, the main features of the experimental data are correctly reproduced. Additionally, figure 4 contains the absolute values of the one-step prediction errors χn (see equation (6)) in color coding over the phase portrait. It is not surprising that the largest errors occur at the normalised values |D| > 0.5, that is in the vicinity of the values of the spontaneous polarization, where nucleation processes cannot be neglected. 194 Learn ferroelectrics using methods of nonlinear dynamics The most surprising result was the fact that we can extract a nonlinear relation ENL(D, Ḋ) for the ferroelectic TGS from our model with no significant indication of a double-well potential [5]. Nevertheless, we showed that it is possible to fit even chaotic time series of the nonlinear resonance circuit with high accuracy with only few assumptions concerning the ferroelectric. Further steps towards an improvement of the model will follow including more physically motivated models of ferroelectric switching processes. References 1. Hegger R., Kantz H., Schmueser F., Diestelhorst M., Kapsch R.-P., Beige H. Dy- namical properties of a ferroelectric capacitor observed through nonlinear time series analysis. // Chaos, 1998, vol. 8, No. 3, p. 727–736. 2. Lines M.E., Glass A.M. Principles and Applications of Ferroelectrics and Related Materials. Oxford, University Press, 1979. 3. Kantz H., Schreiber T. Nonlinear Time Series Analysis. Cambridge, University Press, 1997. 4. Diestelhorst M., Beige H. Nonlinear behaviour of TGS-crystals. // Ferroelectrics, 1988, vol. 81, p. 15–18. 5. Kapsch R.-P., Kantz H., Hegger R., Diestelhorst M. Determination of the dynamical properties of ferroelectrics using nonlinear time series analysis. // Int. J. Bifurcation and Chaos, 2001, vol. 11, No. 3, p. 1019–1034. 6. Diestelhorst M. Nichtlineares Verhalten im ferroelektrischen TGS. – In: Vortraege der 15. Fruehjahrsschule “Ferroelektrizitaet”, Fincken, 6–10 April 1987, p. 68–75 (in German). 7. Strukov B.A., Levanyuk A.P. Physical Foundations of Ferroelectric Phenomena in Crystals. Moskva, Fizmatlit, 1995 (in Russian). 8. Sidorkin A.S. Domain Structure in Ferroelectrics and Related Materials. Moskva, Fiz- matlit, 2000 (in Russian). 9. Guckenheimer J., Holmes P. Nonlinear Oscillations, Dynamical Systems, and Bifur- cations of Vector Fields. New York, Springer, 1990. 10. Duffing G. Erzwungene Schwingungen bei veraenderlicher Eigenfrequenz. Braun- schweig, Vieweg, 1918 (in German). 11. Diestelhorst M., Hofmann R., Beige H. Nonlinear properties of TGS- and BaTiO3- crystals with defects. // Japan J. Appl. Phys., 1985, vol. 24, Suppl. 25–2, p. 1019–1021. 12. Beige H., Diestelhorst M., Forster R., Krietsch T. Chaos near structural phase tran- sitions. // Phase Transitions, 1992, vol. 37, p. 213–238. 13. Stewart H.B., Thompson J.M.T., Ueda Y., Lansbury A.N. Optimal escape from po- tential wells – patterns of regular and chaotic bifurcation. // Physica D, 1995, vol. 85, p. 259–295. 195 M.Diestelhorst Що можна дізнатися про сегнетоелектрики методами нелінійної динаміки? М.Дістельхорст Університет ім. Мартіна Лютера, фізичний факультет, Німеччина, D-06114 Галле, пл. Фрідемана-Баха, 6 Отримано 3 вересня 2002 р. Об’єктом дослідження у даній роботі є нелінійний послідовний резо- нансний контур з сегнетоелекричним конденсатором. Якщо сегне- тоелектрик у конденсаторі знаходиться у полярному стані, спостері- гаються властивості, відомі для модельних систем у нелінійній дина- міці. Це, зокрема, зсув резонансної частоти зі збільшенням прикла- деної напруги, біфуркації та хаотична поведінка. Розгляд системи в рамках простої теорії Ландау веде до опису резонансного контура за допомогою рівняння Дуффінга. Насправді ж, як показано з допо- могою аналізу нелінійних часових послідовностей [1], процеси пере- микання підчас нелінійних коливань значно ускладнюють ситуацію. Ключові слова: сегнетоелектрики, домени, ТГС, нелінійна динаміка, рівняння Дуффінга, хаос PACS: 05.45.a, 05.45.Tp, 77.80.-e, 77.84.Fa 196