A Domain Wall Theory for Ferro electric Hysteresis
Ralph C. Smith
Center for Research in Scienti c Computation
Department of Mathematics
North Carolina State University
Raleigh, NC 27695-8205
Craig L. Hom
Advanced Technology Center
Lo ckheed Martin Missiles & Space
Palo Alto, CA 94304
Abstract
This pap er addresses the mo deling of hysteresis in ferro electric materials through consideration
of domain wall b ending and translation. The development is considered in two steps. In the rst
step, dielectric constitutive relations are obtained through consideration of Langevin, Ising spin and
preferred orientation theory with domain interactions incorp orated through mean eld relations. This
yields a mo del for the anhysteretic p olarization that o ccurs in the absence of domain wall pinning.
In the second step, hysteresis is incorp orated through the consideration of domain wall dynamics and
the quanti cation of energy losses due to inherent inclusions or pinning sites within the material. This
yields a mo del analogous to that develop ed by Jiles and Atherton for ferromagnetic materials. The
viability of the mo del is illustrated through comparison with exp erimental data from a PMN-PT-BT
actuator op erating at a temp erature within the ferro electric regime. i
1 Intro duction
An inherent prop erty of ferro electric materials is the presence of hysteresis in the relation b etween
an applied eld and the resulting p olarization. This phenomenon is asso ciated with the domain
structure of the materials and must be accommo dated in applications which employ ferro electric
comp onents. In certain materials, such as relaxor ferro electrics, hysteresis can be minimized in a
di use transition region through the choice of stoichiometry and op erating conditions. For example,
a signi cant advantage of PMN-PT actuators lies in the fact that they pro duce large strains with
minimal hysteresis when employed near their Curie p oint. Similarly the hysteresis exhibited by PZT
actuators at low to mo derate drive levels is suciently small to p ermit the use of linear design and
control metho dologies. However, at the high drive levels and general temp eratures encountered in
many current applications, the level of hysteresis in b oth PMN-PT and PZT elements is signi cant
and must be accommo dated to attain the full p otential of these materials in high p erformance
actuator design. At a fundamental level, the presence and nature of hysteresis lo ops represents a
basic prop ertycharacterizing ferro electric materials. At the system level, hysteresis generates internal
heat in the actuator which can cause thermal runaway in insulated systems. Highly hysteretic devices
also require closed lo op control for precision displacement actuation. To b etter understand and utilize
these actuator materials, the fundamental mechanisms underlying the hysteresis must b e ascertained
and eventually mo deled.
In an e ort to provide a mo del which incorp orates asp ects of the underlying physical mecha-
nisms and is appropriate for high p erformance actuator and control design, we consider a theory for
ferro electric hysteresis which is based on domain wall dynamics and energy losses due to inherent in-
clusions in the material. Although various typ es of hysteresis curves can b e generated in ferro electric
materials, sigmoid curves such as that shown in Figure 1 are common in applications, and the form
of hysteresis that we shall address here. This theory is the ferro electric analogue of that develop ed
by Jiles and Atherton for ferromagnetic materials [27 ].
Ferro electric Domain Theory
It has long b een recognized that domain and domain wall mechanisms within ferro electric materi-
als result in phenomena suchashysteresis, aging and fatigue. Initial investigations fo cussed primarily
on mechanisms leading to domain formation, nucleation and growth in single crystal BaTiO . Cas-
3
pari and Merz provided an early analysis of the mechanisms leading to linear and quadratic strains
in single crystal ferro electric crystals [5 ]. They concluded that the sp ontaneous p olarization and
domain e ects on the electromechanical resp onse of the material are signi cant and demonstrated
that a linear piezo e ect results when domain e ects were minimized through the application of an
external eld which aligned domains. Caspari and Merz p ostulated that hysteresis is due to the
presence of parallel and antiparallel domains which switchatvarious eld strengths.
Fundamental investigations regarding the nucleation and growth of domains in BaTiO were
3
provided by Merz [34 ] and Little [32 ], and were extended by Drougard [11 ] and Miller and Weinreich
[35 , 37 ]. In concert, these investigations established that b oth nucleation and domain wall motion
contributed to p olarization reversal in BaTiO . The dep endence of p olarization reversal on domain
3
wall motion illustrated that in certain asp ects, this e ect is analogous to magnetization reversal in
ferromagnetic materials. Subsequent exp eriments have elucidated further details regarding domain
nucleation and growth e.g., see [40 , 41] and numerous mo dels have b een prop osed to quantify
domain phenomena and their e ect on hysteresis and fatigue [20 , 23 , 24 , 43 , 45 ].
An imp ortant comp onent in domain and domain wall theories is the e ect due to impurities or
inclusions in the materials. Miller and Savage observed that domain wall motion is dep endentupon 1
the impurity levels of the crystal and considered a mo del based on a surface layer adjacent to the
domain wall [36 ]. Laikhtman extended this idea to develop a mo del for the b ending of domain walls
pinned by inclusions or defects in the material [30 ]. This mo del was develop ed under the premise
that the application of an external eld pro duces a force which acts as a pressure on the wall and
causes it to b end in a manner analogous to a pinned membrane. The resulting mo del was derived
using elasticity theory and accounts for the inertia of the domain wall as well as certain anisotropy
and temp erature e ects. Finally, Laikhtman observed from exp erimental evidence in [14 , 19] that
domain walls exhibit an \e ective viscosity" which dep ends on the concentration of defects. This
provides an irreversible comp onent to the p olarization which Laikhtman incorp orated through the
inclusion of factors mo deling the dissipation of heat. Finally, it is noted that this theory breaks down
if the mo deled wall comes in contact with remote defects. We note that an overview of the \friction"
e ects and dielectric losses due to domain wall motion can b e found in [6] and [50 , pp. 402-427].
The ideas underlying the in uence of pinning sites on domain wall dynamics were extended by
Fedosov and Sidorkin [12 ]. They concluded that for low displacement levels on the order of the
lattice constant, domain walls can be considered as moving as a whole. For higher displacement
levels, however, pinning e ects b ecome signi cant and lo cal b ending b etween pinning sites should b e
considered. A comprehensive analysis of eld, frequency and thermal e ects on domain wall dynamics
in relaxor ferro electrics has b een provided by Chen and Wang [7 ]. These authors note that at low
eld levels with high frequencies, the dielectric b ehavior of the material can b e primarily attributed
to the b ending of domain walls pinned at inclusions. For large applied or internal elds, the lo cal
energy barriers asso ciated with pinning sites are overcome and the domain walls move for extended
distances. We note that these lo calized energy barriers provide a thermo dynamic mechanism for
describing the \e ective viscosity" and \frictional" domain wall e ects rep orted in [30 , 50 ].
The combination of the b ending mechanisms describ ed in [7 , 12 , 30 ], and the energy barriers
detailed in [7 ] indicates that b oth reversible and irreversible mechanisms contribute to ferro electric
hysteresis. Furthermore, the analysis and exp eriments rep orted in these investigations indicate that
these mechanisms are due to the pinning of domain walls at defects in the material and the energy
input required to move domain walls across pinning sites. The discussion in the subsequent section
summarizes certain previous approaches used to quantify these e ects and provides a background
for the mo del presented here.
0.25
0.2
0.15
0.1 ) 2 0.05
0
−0.05 Polarization (C/m
−0.1
−0.15
−0.2
−0.25
−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1
Electric Field (MV/m)
o
Figure 1. Hysteresis lo op measured at a temp erature of T = 20 C in a PMN-PT-BT sample
o
whose Curie temp erature is T =45 C.
c 2
Hysteresis Mo dels
The nature of currenthysteresis mo dels for ferro electric materials can b e roughly categorized as
follows: 1 microscopic theories based on quantum mechanics, classical elasticity and electromagnetic
relations, or thermo dynamic laws, 2 macroscopic theories based up on phenomenological principles,
or 3 semi-macroscopic theories which are derived using a combination of approaches.
Microscopic theories for ferro electric materials generally address the previously describ ed p olar-
ization pro cesses at the lattice or domain level. Mo del development at the lattice level is illustrated
by the theory of Omura et al [42 ]. In their approach, the free energy for a one-dimensional ferro-
electric lattice system is represented in terms of the summed moments and coupling co ecients at
the atoms in the lattice. Hysteresis is incorp orated through the inclusion of viscosity co ecients at
each lattice site. Simulation results demonstrate that the mo del yields reasonable families of hys-
teresis curves but requires a large number of degrees of freedom to accommo date the atoms in the
one-dimensional lattice the examples were generated with N = 250.
For simple stoichiometries, microscopic mo dels of this typ e are theoretically plausible. For the
crystalline conditions encountered in present actuator designs, however, it is dicult to construct
tractable microscopic mo dels which accommo date attributes such as grain b oundaries and intergran-
ular stresses as well as the p olycrystalline nature of the materials. Moreover, microscopic mo dels at
the atomic or lattice level are often dicult to employ in control design due to the large number of
required states and parameters.
The domain wall mo dels of Fedosov and Sidorkin [12 ] and Laikhtman [30 ] are also lo cal in
nature and hence would require either macroscopic averaging or a large numb er of parameters when
employed in actuator characterization or design. However, the concepts underlying these mo dels are
in accordance with the approach employed in this pap er for mo deling the reversible comp onent of the
p olarization with the primary di erence in philosophy residing in our use of macroscopic parameters
which quantify \average" prop erties of the material.
Numerous macroscopic mo dels for hysteresis pro cesses in ferro electric materials have b een de-
rived using empirical or phenomenological principles. In some cases, these mo dels are motivated by
physical domain pro cesses while other mo dels circumvent p o orly understo o d physics through purely
mathematical characterizations. Mo dels in the latter category are advantageous when the underlying
physics is dicult to quantify. Due to their nature, however, it is dicult to employ known physics
or physical measurements to directly construct such purely phenomenological mo dels or up date them
to accommo date changing op erating conditions.
A large class of hysteresis mo dels have b een derived through ts obtained with hyp erb olic tangent
and cotangent functions. As detailed in Section 2 of this pap er, Langevin functions involving the
hyp erb olic cotangent and Ising spin mo dels yielding the hyp erb olic tangent arise when Boltzmann
statistics are used to derive the anhysteretic hysteresis-free resp onse in ferro electric and ferromag-
netic materials. These representations provide natural envelops for phenomenological quanti cation
of hysteresis families.
Hom and Shankar have demonstrated that certain hysteresis curves can be generated through
the appropriate choice of parameters in the Ising spin mo del [22 ]. As illustrated by their results,
however, this metho d pro duces hysteresis curves in which the transition from one saturation region
to the other is much steep er than that observed in typical ferro electric applications.
More gradual transitions can b e obtained by employing Langevin or hyp erb olic cotangent func-
tions which are shifted to coincide with the remanence p olarizations. This technique was employed
by Zhang and Rogers when mo deling hysteresis in piezo ceramic materials [54 ] and Miller et al for
characterization hysteresis in ferro electric capacitors [38 ]. While this technique often provides ac-
curate mo del ts with a small number of required parameters in some cases two, its capabilities 3
for general applications involving symmetric and asymmetric minor lo ops app ears limited. In the
context of the underlying physics, the shifting of the anhysteretic curves provides a phenomenological
mechanism for mo deling the \viscosity" or energy losses which o ccur during domain growth through
domain wall movement.
Another class of phenomenological hysteresis mo dels are based up on the determination of empir-
ical laws which quantify the rate of dip ole or p olarization switching. Based on observations by Merz
[34 ], Landauer et al [31 ] prop osed a mo del which quanti es the rate of p olarization switching in terms
of the current and saturation p olarization and applied eld. Applying a similar philosophy, Chen
and Montgomery [6] develop ed a mo del in which the rate at which dip oles align with the applied
eld is sp eci ed as a nonlinear function of the eld. Both approaches provide phenomenological
characterizations of the rate at which domain walls move.
Preisach and generalized Preisach mo dels have also b een widely used to characterize hysteresis in
ferro electric materials. In their classical form, Preisach op erators are constructed from linear com-
binations of multivalued kernels [33 ]. The co ecients for each kernel represent the input level and
magnitude when switching o ccurs b etween two saturation states which de ne the kernel. When used
to quantify magnetic hysteresis, the kernel has b een interpreted as representing diametrically opp o-
site dip ole con gurations with the co ecients interpreted as quantifying the eld inputs required for
dip ole switching and degree to which they switch at the p oints. While a similar interpretation can
be made for ferro electric pro cesses, the techniques are more commonly considered as phenomeno-
logical and hence provide a purely mathematical mo del for the pro cess. Generalized Preisach or
Krasnoselskii-Pokrovskii mo dels di er through the choice of kernel. As detailed in [2 ], the classical
Preisach op erators are discontinuous with resp ect to b oth the parameters and time. This is allevi-
ated in Krasnoselskii-Pokrovskii op erators through the use of smo oth ridge functions which provide
an envelop of admissible families. We note that the construction of Krasnoselskii-Pokrovskii kernels
through translations of ridge functions is quite similar in philosophy to the use of shifted anhysteretic
curves employed in the mo dels of Miller et al [38 ] and Zhang and Rogers [54 ].
The mo di cations necessary for employing the classical Preisach techniques for ferro electric hys-
teresis are illustrated by the mo dels develop ed by Ge and Jouaneh [17 , 18 ]. As illustrated by b oth nu-
merical and exp erimental examples, this technique has b een used accurately characterize b oth ma jor
and minor hysteresis lo ops in piezo ceramic actuators. The development of Krasnoselskii-Pokrovskii
hysteresis mo dels for piezo ceramic actuators has b een considered by Galinaitis and Rogers [15 , 16].
References [16 ] also provides an overview concerning the construction of inverse op erators for control
metho ds based on output linearization.
Finally, semi-macroscopic mo dels are typically derived using a combination of thermo dynamic
and phenomenological principles. The goal in these approaches is to employ energy-based relations to
the extent p ossible so that known physics can b e used to facilitate mo del construction and up dating.
Phenomenological principles contribute to asp ects of the constitutive equations and the development
of macroscopic averages of microscopic prop erties such as the intergranular stresses and domain wall
pinning site energies.
An illustration of this typ e of approach can be found in the hysteresis mo del of Yoo and Desu
for ferro electric thin lms [53 ]. The hysteresis is characterized through consideration of three stages
in the p olarization pro cess: domain nucleation and growth, domain merging, and domain shrink-
age. The dynamics in these stages are develop ed through the combination of physical mechanisms
underlying p olarization reversal and macroscopic averages of certain domain prop erties. The do-
main wall mo del of Chen and Wang [7 ] for relaxor ferro electric materials can also b e considered as
semi-macroscopic. In their development, Chen and Wang employ energy relations in combination
with phenomenological relations to obtain microscopic p olarization rate equations and corresp ond-
ing macroscopic constitutive equations. Asp ects of the mo del are formulated in the framework of 4
dislo cation theory for plasticity and anelasticity which leads to a fairly general mo del for domain
wall dynamics. While low temp erature hysteresis lo ops are used to provide evidence of domain wall
motion, this investigation do es not address p er se the mo deling of hysteresis in relaxor ferro electric
materials.
This pap er describ es a hysteresis mo del which is semi-macroscopic in nature and is based up on
the quanti cation of the domain wall energy in the material. For an ideal material which is free
from inclusions or second-phase materials, domain wall movement is unimp eded and the relation
between applied stresses, electric elds and temp eratures, and the resulting hysteresis-free or anhys-
teretic p olarization is mo deled by nonlinear constitutive equations. As detailed in [22 ], appropriate
constitutive relations are derived through thermo dynamic theory combined with phenomenological
relations such as Ho oke's law. In actual materials, however, domain wall movement is imp eded by
inherent inclusions within the material. The mo del presented here quanti es the resulting hysteresis
through consideration of the energy required to translate domain walls and break pinning sites. An
`average' measure of this energy is employed to obtain a mo del which requires the identi cation of
only ve parameters. Due to the small numb er and physical nature of certain parameters e.g., the
saturation p olarization P , the mo del is easily constructed and up dated. We note that this mo del
s
is the ferro electric analogue of that develop ed by Jiles and Atherton for ferromagnetic materials
[27 ]. While asp ects of the Jiles/Atherton theory were employed in [51 ] to characterize hysteresis in
stress/strain relations for LiCsSO crystals, this is the rst application of this approach to quantify
4
hysteresis in the E -P relation for ferro electric materials. We also p oint out that while certain ferro-
electric actuator materials provide the motivation, the mo del is general in nature and can b e applied
toavariety of ferro electric materials including capacitor and electro-optic materials.
Certain asp ects of the constitutive theory of Hom and Shankar are summarized in Section 2.
This discussion includes the development of the classical Langevin and Ising spin mo dels as well
as a third mo del which combines attributes of the two. This section also provides a discussion
concerning the incorp oration of interdomain coupling to provide the e ective eld observed at the
domain level and outlines the extension of the mo dels from the microscopic to macroscopic levels.
The domain wall theory contributing to b oth reversible and irreversible changes in p olarization is
presented in Section 3. To simplify the development, the nal hysteresis mo del is formulated through
scalar-valued relations between the applied eld E and resulting p olarization P . This formulation
is suitable for a variety of current applications and materials. Atvarious p oints in the development,
however, intermediate vector-valued relations are required and the notation E and P are used to
3
denote the resp ective eld and p olarization in lR . This also provides a framework for developing the
corresp onding vector-valued hysteresis mo del if required by the application. In the nal section, the
viability of the mo del is illustrated through comparison with exp erimental data from a PMN-PT-BT
actuator.
2 Constitutive Theory for Ferro electric Materials
In this section we consider the development of constitutive equations which describ e the anhysteretic
hysteresis-free mechanical and dielectric b ehavior of certain ferro electric materials. Hysteresis
e ects are included in the next section through quanti cation of energy loss due to pinning sites in
the material.
Following classical developments for electrostrictive materials, the constitutive relations are as-
sumed to b e of the form
2
= Y e Q P
33
1
E = 2Q P + FT; P :
33 5
Here denotes the axial stress, e denotes the longitudinal strain, Y is the closed-circuit elastic mo du-
lus, T is the temp erature, P denotes the p olarization and Q denotes the longitudinal electrostrictive
33
coupling co ecient. We note that these constitutive relations can b e mo di ed for piezo electric ma-
terials by replacing the quadratic p olarization term in the mechanical relation by a linear term.
In the remainder of the section, we summarize three techniques for mo deling the function F
which quanti es the dep endence of the electric eld E on the temp erature and p olarization. All three
mo dels are based up on the application of Boltzmann's lawtovarious dip ole orientations within the
ferro electric crystal. This follows the approach employed by Hom and Shankar [21 , 22 ] and in one case
provides a Langevin mo del analogous to that typically employed in ferromagnetic applications. We
note that other techniques can b e employed for mo deling the anhysteretic constitutive relations for
certain ferro electric materials and refer the reader to [8 , 44 ] for discussion concerning two alternative
mo dels.
2.1 Noninteracting Cells
We view a ferro electric crystal as a lattice of cells, where each cell p ossesses a p ermanent dip ole
moment with an asso ciated direction. Normally, the cell has a nite numb er of p ossible orientations
eight for a tetragonal system. When free of electric eld and ab ove the Curie temp erature, thermal
agitation creates a random distribution of orientations and the ceramic has an average p olarization of
zero. The application of an electric eld forces some dip ole moments to switch to an orientation closer
to the direction of the eld. The resulting distribution is then no longer random and a macroscopic
p olarization develops. We p oint out that collections of neighb oring cells having the same p olarization
form the domain structure of the crystal. The degree of switching increases with increasing eld until
the global p olarization eventually saturates when all the cells have optimally oriented their dip ole
moments relative to the eld. For the development in this section, we will assume that the cells do
not interact and use statistical mechanics to derive three nonlinear mo dels for the dielectric resp onse
of this system. The mo dels di er only in the assumptions made ab out the cells' p ossible orientation.
The e ects of cell interaction are then considered in Section 2.2.
For a dip ole moment p in an electric eld E, the p otential energy is given by
0
E = p E = p E cos 2
0 0
where p = jp j;E = jEj. For noninteracting cells and moments, classical Boltzmann statistics can
0 0
be used to express the probability of dip oles o ccupying certain energy states. Letting k denote
B
Boltzmann's constant, the thermal energy is k T and the probability that a dip ole o ccupies the
B
energy state E is