Open access peer-reviewed chapter

# Advances in Thermodynamics of Ferroelectric Phase Transitions

By Shu-Tao Ai

Submitted: April 5th 2012Reviewed: August 4th 2012Published: November 19th 2012

DOI: 10.5772/52089

## 1. Introduction

The thermodynamics of ferroelectric phase transitions is an important constituent part of the phenomenological theories of them, as well as the interface dynamics of them. In particular, if we confine the thermodynamics to the equilibrium range, we can say that the Landau-Devonshire theory is a milestone in the process of the development of ferroelectric phase transition theories. This can be found in many classical books such as [1,2]. Many studies centering on it, especially the size-effects and surface-effects of ferroelectric phase transitions, have been carried out. For the reason of simplicity, we just cite a few [3,4]. But we think these are a kind of technical but not fundamental progress.

Why we think so is based on that the Landau-Devonshire theory is confined to the equilibrium range in essene so it can’t deal with the outstanding irreversible phenomonon of first-order ferroelectric phase transitions strictly, which is the „thermal hysteresis“. The Landau-Devonshire theory attributes the phenomenon to a series of metastable states existing around the Curie temperature TC. In principle, the metastable states are not the equilibrium ones and can not be processed by using the equilibrium thermodynamics. Therefore, we believe the Landau-Devonshire theory is problematic though it is successful in mathematics. The real processes of phase transition were distorted. In this contribution, the Landau-Devonshire theory will be reviewed critically, then the latest phenomenological theory of ferroelectric phase transitions will be established on the basis of non-equilibrium or irreversible thermodynamics.

This contribution are organized as the follows. In Section 2, we will show the unpleasant consequence caused by the metastable states hypothesis, and the evidence for the non-existence of metastable states, i.e. the logical conflict. Then in Section 3 and 4, we will give the non-equilibrium ( or irreversible ) thermodynamic description of ferroelectric phase transitions, which eliminates the unpleasant consequence caused by the metastable states hypothesis. In Section 5, we will give the non-equilibrium thermodynamic explanation of the irreversibility of ferroelectric phase transitions, i.e. the thermal hysteresis and the domain occurrences in ferroelectrics. At last, in Section 6 we will make some concluding remarks and look forward to some possible developments.

## 2. Limitations of Landau-Devonshire theory and demonstration of new approach

The most outstanding merit of Landau-Devonshire theory is that the Curie temperature and the spontaneous polarization at Curie temperature can be determined simply. However, in the Landau-Devonshire theory, the path of a first-order ferroelectric phase transition is believed to consist of a series of metastable states existing around the Curie temperature. This is too difficult to believe because of the difficulties encounted ( just see the follows) TC.

### 2.1. Unpleasant consequence caused by metastable states hypothesis

Basing on the Landau-Devonshire theory, we make the following inference. Because of the thermal hysteresis, a first-order ferroelectric phase transition must occur at another temperature, which is different from the Curie temperature [5]. The state corresponding to the mentioned temperature ( i.e. actural phase transition temperature ) is a metastable one. Since the unified temperature and spontanous polarization can be said about the metastable state, we neglect the heterogeneity of system actually. In other words, every part of the system, i.e. either the surface or the inner part, is of equal value physically. When the phase transition occurs at the certain temperature, every part of the system absorbs or releases the latent heat simultaneously by a kind of action at a distance. ( The concept arose in the electromagnetism first. Here it maybe a kind of heat transfer. ) Otherwise, the heat transfer in system, with a finite rate, must destroy the homogeneity of system and lead to a non-equilibrium thermodynamic approach. The unpleasant consequence, i.e. the action at a distance should be eliminated and the lifeforce should be bestowed on the non-equilibrium thermodynamic approach.

In fact, a first-order phase transition process is always accompanied with the fundamental characteristics, called the co-existence of phases and the moving interface ( i.e. phase boundary ). The fact reveals that the phase transition at various sites can not occur at the same time. Yet, the phase transition is induced by the external actions ( i.e. absorption or release of latent heat ). It conflicts sharply with the action at a distance.

### 2.2. Evidence for non-existence of metastable states: logical conflict

In the Landau-Devonshire theory, if we neglect the influence of stress, the elastic Gibbs energy G1can be expressed with a binary function of variables, namely the temperature Tand the electric displacement D( As G1is independent of the orientation of D, here we are interested in the magnitude of Donly )

G1=g1(T,D)E1

The long-standing, close correlation between analytical dynamics and thermodynamics implies that Equation (1) can be taken as a scleronomic constraint equation

f1(G1,T,D)=G1g1(T,D)=0E2

where G1,T,Dare the generalized displacements. The possible displacements dG1,dTand dDsatisfy the following equation

dG1g1TdTg1DdD=0E3

In the Landau-Devonshire theory, the scleronomic constraint equation, i.e. Equation (1) is expressed in the form of the power series of D( For simplicity, only the powers whose orders are not more than six are considered )

G1=G10+12αD2+14βD4+16γD6E4

where α,β,γare the functions of T, and G10is the elastic Gibbs energy of paraelectric phase. The relation between G1and Dat various temperature, which belongs to first-order phase transition ferroelectrics is represented graphically in Figure 1. The electric displacements which correspond to the bilateral minima of G1are identified as ±D*, and the electric displacement which corresponds to the middle minimum of G1equals zero. The possible electric displacements should be the above ones which correspond to the minima of G1.

Equivalently, imposed on the generalized displacements G1,T,Dis a constraint, which is

G1D=0E5

So, the possible displacement dDshould be the follows dD*,dD*,0(±D*)=D*,±D*0=±D*,0. After all, if our discussion are limited in the equilibrium thermodynamics strictly, there must be the third constraint, i.e. the equilibrium Dand Tshould satisfy

h(D,T)=0E6

where his a binary function of the variables Dand T. It can be determined by the principle of minimum energy

G1=minE7

for certain T. Then, the metastable states are excluded. Thus, the thermal hysteresis does not come into being. The corollary conflicts with the fact sharply. This reveals that the first-order ferroelectric phase transition processes must not be reversible at all so as not to be dealt with by using the equilibrium thermodynamics.

How can this difficulty be overcome? An expedient measure adopted by Devonshire is that the metastable states are considered. However, do they really exist?

Because the metastable states are not the equilibrium ones, the relevant thermodynamic variables or functions should be dependent on the time t. In addition, the metastable states are close to equilibrium, so the heterogeneity of system can be neglected. Here, the elastic Gibbs energy G1should be

G1=g2(T,D,t)E8

For the same reason as was mentioned above, Equation (8) can be regarded as a rheonomic constraint on the generalized displacements G1,T,D

f2(G1,T,D,t)=G1g2(T,D,t)=0E9

In this case, the possible displacements dG1,dTand dDsatisfy the following equation

dG1g2TdTg2DdDg2tdt=0E10

Comparing Equation (3) with Equation (10), we may find that the possible displacements here are not the same as those in the former case which characterize the metastable states for they satisfy the different constraint equations, respectively. ( In the latter case, the possible displacements are time-dependent, whereas in the former case they are not. ) Yet, the integral of possible displacement dDis the possible electric displacement in every case. The possible electric displacements which characterize one certain metastable state vary with the cases. A self-contradiction arises. So the metastable states can not come into being.

What are the real states among a phase transition process? In fact, both the evolution with time and the spatial heterogeneity need to be considered when the system is out of equilibrium [6-9]. Just as what will be shown in Section 2.3, the real states should be the stationary ones, which do not vary with the time but may be not metastable.

### 2.3. Real path: Existence of stationary states

The real path of a first-order ferroelectric phase transition is believed by us to consist of a series of stationary states. At first, this was conjectured according to the experimental results, then was demonstrated reliable with the aid of non-equilibrium variational principles.

Because in the experiments the ferroelectric phase transitions are often achieved by the quasi-static heating or cooling, we conjetured that they are stationary states processes [8]. The results on the motion of interface in ferroelectrics and antiferroelectrics support our opinion [10-12]. From Figure 2, we may find that the motion of interface is jerky especially when the average velocity vais small. A sequence of segments of time corresponding to the states of rest may be found. The experimental results about other materials such as PbTiO3 are alike [11]. This reveals that in these segments of time ( i.e. characteristic time of phase transition ) the stationary distributions of temperature, heat flux, stress, etc. may be established. Otherwise, if the motion of interface is continuous and smooth, with the unceasing moving of interface ( where the temperature is TC) to the inner part, the local temperature of outer part must change to keep the temperature gradient Tof this region unchanged for it is determined by ±lρv=Jqdiff=κT, where lis the latent heat ( per unit mass ), ρis the mass density, vis the velocity of interface ( where the phase transition is occurring ), Jqdiffis the diffusion of heat, i.e. heat conduction, κis the thermal conductivity ( and maybe a tensor. ) Then, the states are not stationary.

The non-equilibrium variational principles are just the analogue and generalization of the variational principles in analytical dynamics. The principle of least dissipation of energy, the Gauss’s principle of least constraint and the Hamiltonian principle etc., in non-equilibrium thermodynamics play the fundamental roles as those in analytical dynamics. They describe the characteristics of stationary states or determine the real path of a non-equilibrium process.

For the basic characteristics of non-equilibrium processes is the dissipation of energy, the dissipation function φis defined as

φ=σsπE11

where σsis the rate of local entropy production and πis the external power supply ( per unit volume and temperature ). After the rather lengthy deducing and utilizing the thermodynamic Gauss’s principle of least constraint which makes the system choose a real path [13], the evolution with time tof the deviation from a given non-equilibrium stationary state, ξ(t)={ξi(t)}, was obtained in two cases. If no external power supply,

ξi(t)=ξi(0)e2SiRitE12

where ξi(0)is the initial value of ξi(t), Siis the coefficient between the linear variation in the thermodynamic force χTi(1)(ξ)and the deviation of the extensive pseudo-thermodynamic variable from a given non-equilibrium stationary state ξi, Riis the coefficient between the linear variation in the dissipative force χDi(1)(ξ˙)and the time-derivative of the deviation of the extensive pseudo-thermodynamic variable from a given non-equilibrium stationary state ξ˙i. Equation (12) define the real path with the addition that Rishould be a suitable value Ri*. If the external power supply exists, similarly the evolution of deviation ξican be obtained

ξi(t)=ξi(0)e2(SiVi)RitE13

where Viis the coefficient between the linear variation in the force related to the external power χEi(1)and the deviation of the extensive pseudo-thermodynamic variable from a given non-equilibrium stationary state ξi. If the coefficients Vi,Riassume the suitable values Vi**,Ri**, the system choose a real path.

Both the real paths in the two cases reveal that the deviations decrease exponentially when the system regresses to the stationary states. Stationary states are a king of attractors to non-equilibrium states. The decreases are steep. So the regressions are quick. It should be noted that we are interested in calculating the change in the generalized displacements during a macroscopically small time interval. In other words, we are concerned with the determination of the path of an irreversible process which is described in terms of a finite difference equation. In the limit as the time interval is allowed to approach zero, we obtain the variational equation of thermodynamic path.

So, if the irreversible process is not quick enough, it can be regarded as the one that consists of a series of stationary states. The ferroelectric phase transitions are usually achieved by the quasi-static heating or cooling in the experiments. So, the processes are not quick enough to make the states deviate from the corresponding stationary states in all the time. In Figure 3, three types of regions and their interfaces are marked with I, II, III, 1, 2 respectively. The region III where the phase transition will occur is in equilibrium and has no dissipation. In the region I where the phase transition has occurred, there is no external power supply, and in the region II ( i.e. the paraelectric-ferroelectric interface as a region with finite thickness instead of a geometrical plane ) where the phase transition is occurring, there exists the external power supply, i.e. the latent heat ( per unit volume and temperature ). According to the former analysis in two cases, we may conclude that they are in stationary states except for the very narrow intervals of time after the sudden lose of phase stability.

## 3. Thermo-electric coupling

In the paraelectric-ferroelectric interface dynamics induced by the latent heat transfer [6,7], the normal velocity of interface vnwas obtained

vn=1lρ[kfer(T)ferkpar(T)par]nE14

where lis the latent heat ( per unit mass ), ρis the density of metastable phase (paraelectric phase ), kferis the thermal conductivity coefficient of ferroelectric phase, kparis the thermal conductivity coefficient of paraelectric phase, (T)feris the temperature gradient in ferroelectric phase part, (T)paris the temperature gradient in paraelectric phase part, nis the unit vector in normal direction and directs from the ferroelectric phase part to the paraelectric phase part. The temperature gradients can be studied from the point of view that a ferroelectric phase transition is a stationary, thermo-electric coupled transport process [8].

### 3.1. Local entropy production

In the thermo-electric coupling case, the Gibbs equation was given as the following [8]

Tds=duEdDiμidniE15

where T,E,Dis the temperature, the electric field intensity and the electric displacement within a random small volume, respectively; s,u,μi,niis the entropy density, the internal energy density, the chemical potential and the molar quantity density in the small volume, respectively. And there, it was assumed that the crystal system is mechanically-free ( i.e. no force is exerted on it ). Differentiating Equation (15) and using the following relations

ut+Ju=0E16
Dt=(ε0E+P)t=JPE17
nit+Jni=0E18

we have

st=1TJuETJP+iμiTJniE19

where Ju,JP,Jniis the energy flux, the polarization current and the matter flux; Pis the polarization. Jushould consist of three parts: the energy flux caused by the heat conduction, the energy flux caused by the charge transport and the energy flux caused by the matter transport

Ju=Jq+φJP+iμiJniE20

where Jqis the heat flux, φis the electrical potential and satisfies

E=φE21

Then we deduce the following

st=(Jq+φJPT)+Jq(1T)+JP(φT)iJniTμiE22

If we define a entropy flux Jsand a rate of local entropy production σsas

Js=Jq+φJPTE23
σs=Jq(1T)+JP(φT)iJniTμiE24

Equation (22) can be written as

st+Js=σsE25

This is the local entropy balance equation. We know, the system is in the crystalline states before and after a phase transition so that there is no diffusion of any kind of particles in the system. So, Jni=0. The local entropy production can be reduced as

σs=Jq(1T)+JP(φT)E26

We know the existence of ferroics is due to the molecular field. It is an internal field. So we must take it into account. Here, the electric field should be the sum of the external electric field Eeand the internal electric field Ei

E=Ee+EiE27

Correspondingly, there are the external electrical potential φeand the internal electrical potential φiand they satisfy

Ee=φeE28
Ei=φiE29

If the external electric field is not applied, φecan be a random constant. There is no harm in letting the constant equal zero. Then the entropy production equals

σs=Jq(1T)+JP(φiT)E30

According to the crystal structures of ferroelectrics [2], we know the polarization current JPoriginates from the displacement or ordering of ions in ferroelectrics. We may consider it as the transport of charges influenced by the internal electric field.

### 3.2. Description of phase transitions and verification of interface dynamics

Assume the external electric field is not applied. Here are the thermodynamic fluxex Jq,JPand the corresponding thermodynamic forces Xq,XP

Xq=(1T)E31
XP=(φiT)E32
Jican be expanded linearly with Xj(i,j=q,P)
Jq=LqqXq+LqPXPE33
JP=LPqXq+LPPXPE34

where Lqq,LqP,LPqand LPPare the transport coefficients, which are four second-order tensors. They should satisfy the generalized Onsager relations [14]

Lqq=LqqT, LqP=LPqT,
LPP=LPPTE35

Because to a first-order ferroelectric phase transition the electric displacement changes suddenly and so does the internal electrical potential, the force XPof the region where the phase transition is occurring can be regarded as a large constant roughly in the characteristic times of phase transition ( i.e. the times in which the interface keeps rest ). For the ferroelectric phase transition may be regarded as a stationary states process, the principle of minimum entropy production must be satisfied [15].

According to Equations (30)-(35), we have

σs=Lqq:XqXq+2LqP:XqXP+LPP:XPXPE36

If there is no any restriction on Xqand XP, according to the conditions on which the entropy production is a minimum

(σsXq)XP=2LqqXq+2LqPXP=2Jq=0E37
(σsXP)Xq=2LPqXq+2LPPXP=2JP=0E38

We know the stationary states are equilibrium ones actually. If we let Xq( or XP) be a constant, according to Equation (38) ( or (37) ) we know JP( or Jq) which is corresponded to another force XP( or Xq) should be zero.

Then, a first-order ferroelectric phase transition can be described by the second paradigm. Since the force XPof the region where the phase transition is occurring is a large constant, the flux Jqof this region should be zero ( but JP0). This states clearly that the pure heat conduction and the heat conduction induced by the thermo-electric coupling cancel out each other so as to release or absorb the latent heat. It is certain that the latent heat passes through the region where the phase transition has occurred ( at the outside of the region where the phase transition is occurring ) and exchanges itself with the thermal bath. Accompanied with the change of the surface’s temperature and the unceasing jerky moving of the region where the phase transition is occurring, a constant temperature gradient is kept in the region where the phase transition has occurred, i.e. the force Xqis a constant. So, the flux JP=0( but Jq0). This states clearly that the electric displacement of the region where the phase transition has occurred will not change but keep the value at Curie temperature or zero until the phase transition finishes. Differently, the region where the phase transition will occur should be described by the first paradigm for there is no restriction on the two forces XP,Xq. The states of this region are equilibrium ones. So the temperature gradient Tshould be zero.

Considering that (T)par0for the region where the phase transition will occur ( i.e. the paraelectric phase part ) can be regarded as an equilibrium system, we modify Equation (14) as

vn=kfer(T)fernlρE39

In order to compare it with the experiments, we make use of the following values which are about PbTiO3 crystal: ρ=7.1g/cm3 [16], l=900cal/mol [17], kfer=8.8 ×105erg/cm﹒s﹒K [18]. The value of the velocity of the interface’s fast motion, which has been measured in the experiments, is 0.5mm/s [11]. According to Equation (39), we calculate the corresponding temperature gradient to be 57.35K/cm. However, in [19] it is reported that the experimental temperature gradient varies from 1.5 to 3.5K/mm while the experimental velocity of the interface’s motion varies from 732 to 843μm/s. Considering the model is rather rough, we may conclude that the theory coincides with the experiments.

### 3.3. Relation between latent heat and spontaneous polarization

In the experiments, the latent heat and the spotaneous polarization are measured often for first-order ferroelectric phase transitions. So in the follows, we will establish the relation between latent heat and spontaneous polarization in the realm of non-equilibrium thermodynamics.

All the quantities of the region where the phase transition has occurred are marked with the superscript “I”; all the quantities of the region where the phase transition is occurring are marked with the superscript “II”; and all the quantities of the region where the phase transition will occur are marked with the superscript “III”. Let’s consider the heating processes of phase transition firstly. In the region where the phase transition has occurred,

JqIpara=LqqIparaXqIpara=lρvaE40

where we have ignored the difference between the mass density of ferroelectric phase and that of paraelectric phase ( almost the same ) and denote them as ρ, vais the average velocity of interface. In the region where the phase transition is occurring,

JqII=LqqIIXqII+LqPIIXPII=0E41
JPII=LPqIIXqII+LPPIIXPIIE42

The heat which is transferred to the region where the phase transition is occurring is absorbed as the latent heat because the pure heat conduction and the heat conduction induced by the thermo-electric coupling cancel out each other. So,

lρva=LqqIIXqIIE43

According to Eqations (41)-(43), we work out

JPII=lρ[LPqII(LqqII)1LPPII(LqPII)1]vaE44

where the superscript “-1” means reverse. While

JPIIdt=Dspon=ε0φiIIIPsponIII=PsponIIIE45

where we utilized the boundary condition of Dand considered the region where the phase transition will occur is in equilibrium, the superscript “spon” means spontaneous. So,

Pspon=PsponIII=JPIIdtE46

The relation between latent heat and spontaneous polarization are obtained. In the cooling processes of phase transition,

lρva=LqqIferrXqIferr+LqPIferrXPIferr=LqqIIXqIIE47
JPIIdt=DsponI=ε0φiI+PsponIE48

Repeating the above steps, we obtain

PsponI=JPIIdt+ε0φiIE49
With
JPII=lρ[LPqII(LqqII)1LPPII(LqPII)1]va=JPIIE50

Then we find that PsponIII( serves as the equilibrium polarization ) is not equal PsponI( serves as the non-equilibrium polarization ). For the region where the phase transition will occur is stressed [9], there is some difference between Pspon(=PsponIII)and the equilibrium spontaneous polarization without the affects of stress Psponbecause of the piezoelectric effect.

## 4. Thermo-electro-mechanical coupling

The comprehensive thermo-electro-mechanical coupling may be found in ferroelectric phase transition processes. Because there exists not only the change of polarization but also the changes of the system’s volume and shape when a ferroelectric phase transition occurs in it, the mechanics can not be ignored even if it is mechanically-free, i.e. no external force is exerted on it. To a first-order ferroelectric phase transition, it occurs at the surface layer of system firstly, then in the inner part. So, the stress may be found in the system.

Since one aspect of the nature of ferroelectric phase transitions is the thermo-electro-mechanical coupling, we take the mechanics into account on the basis of Section 3, where only the thermo-electric coupling has been considered. This may lead to a complete description in the sense of continuum physics.

### 4.1. Deformation mechanics

For a continuum, the momentum equation in differential form can be written as

σ+ρf=ρaE51

where σ,f,a,ρis the stress, the volume force exerted on unit mass, the acceleration and the mass density, respectively. Let k=12v2be the local kinetic energy density (per unit mass), with vis the velocity. Then

dkdt=vaE52

where tis the time. In terms of

(vσ)=v(σ)+(v):σE53
σ:(v)=σ:dE54

we can deduce the following balance equation of mechanical energy basing on Equations (51) (52)

ρdkdt(vσ)=ρvfσ:dE55

which is in differential form and in Lagrangian form. Or in Eulerian form

(ρk)t+(ρkvvσ)=ρvfσ:dE56

where d=12(v+v)=12[v+(v)T]is the rate of deformation or strain rate ( the superscript “T” means transposition ).

To a ferroelectric phase transition, fand σmay be the nominal volume force and stress, which are the embodiments of the actions of thermo-electro-mechanical coupling and are two internal fields. Generally, they are the sums of real and nominal volume force or stress

f=freal+fnomE57
σ=σreal+σnomE58

The nominal volume force and stress are not zero until the eigen ( or free ) deformation of system finishes in phase transitions. If they are zero, the eigen ( or free ) deformation finishes.

### 4.2. Local entropy production and description of phase transitions

The Gibbs equation was given as the following [9]

Tds+iμidni+1ρEdD+fdr+1ρ(vσ)dt=du+dkE59

where T,E,Dis the local temperature, the local electric field intensity and the local electric displacement, respectively; s,u,ni,μiis the local entropy density ( per unit mass ), the local internal energy density ( per unit mass ), the local molar quantity density ( per unit mass ) and the chemical potential, respectively; ris the displacement vector. If the external electric field is not applied, the quantity Eis the internal electric field Eionly [8].

Make the material derivative of Equation (59) with t, then obtain

ρdsdt+ρTiμidnidt+1TEidDdt+ρTfdrdt+1T(vσ)=ρTdedtE60

where eis the total energy with e=u+k. We know, the system is in the crystalline states before and after a phase transition so that there is no diffusion of any kind of particles in the system. So, nias the local molar quantity ( per unit mass ) does not change with t, i.e. dnidt=0. dDdtstands for the polarization current JP, while drdtstands for the velocity v, Ei=φi( φistands for the internal electric potential ).

After the lengthy and troublesome deduction [9], the local entropy balance equation in Lagrangian form can be obtained

ρdsdt=Js+σsE61

with the entropy flux Js

Js=Jqdiff+φiJPvσTE62

and the rate of local entropy production σs

σs=Jqdiff(1T)+JP(φiT)σ:[v(1T)]ρfvTE63

where Jqdiffis the diffusion of heat, i.e. heat conduction.

Here are the thermodynamic fluxes Jqdiff(=Jq),JP,σ(=Jσ),ρf(=Jf)and the corresponding thermodynamic forces Xq,XP,Xσ,Xf

Xq=(1T)E64
XP=(φiT)E65
Xσ=v(1T)E66
Xf=vTE67
Ji(i=q,P,σ,f)can be expanded linearly with Xj(j=q,P,σ,f)
Jq=LqqXq+LqPXP+Lqσ:Xσ+LqfXfE68
JP=LPqXq+LPPXP+LPσ:Xσ+LPfXfE69
Jσ=LσqXq+LσPXP+Lσσ:Xσ+LσfXfE70
Jf=LfqXq+LfPXP+Lfσ:Xσ+LffXfE71

where Lqq,LPP,Lff,LqP,LPq,Lqf,Lfq,LPf,LfPare nine second-order tensors, Lqσ,Lσq,LPσ,LσP,Lσf,Lfσare six third-order tensors, Lσσis a fourth-order tensor. They should satisfy the generalized Onsager relations [14]

Lαα=LααT
(α=q, P, σ, f)E72
Lαβ=LβαT
(α, β=q, P, σ,f, and αβ)E73

So the rate of local entropy production can be written as

σs=Lqq:XqXq+LqPXPXq+LqσXσTXq+LqfXfXq     +LPqXqXP+LPP:XPXP+LPσXσTXP+LPfXfXP     +LσqXqXσT+LσPXPXσT+LσσXσTXσT+LσfXfXσT     +LfqXqXf+LfPXPXf+LfσXσTXf+Lff:XfXfE74

According to the condition on which the local entropy production is a minimum, from Equation (74) we can deduce the following

(σsXq)XP, Xσ, Xf=2LqqXq+LqPXP+Lqσ:Xσ+LqfXf             +LPqTXP+LσqT:Xσ+LfqTXf             =2(LqqXq+LqPXP+Lqσ:Xσ+LqfXf)             =2Jq             =0E75
(σsXP)Xq, Xσ, Xf=LqPTXq+LPqXq+2LPPXP+LPσ:Xσ              +LPfXf+LσPT:Xσ+LfPTXf              =2(LPqXq+LPPXP+LPσ:Xσ+LPfXf)              =2JP              =0E76
(σsXσ)Xq, XP, Xf=LqσTXq+LPσTXP+LσqXq+LσPXP              +2Lσσ:Xσ+LσfXf+LfσTXf              =2(LσqXq+LσPXP+Lσσ:Xσ+LσfXf)              =2Jσ              =0E77
(σsXf)Xq, XP, Xσ=LqfTXq+LPfTXP+LσfT:Xσ+LfqXq             +LfPXP+Lfσ:Xσ+2LffXf             =2(LfqXq+LfPXP+Lfσ:Xσ+LffXf)             =2Jf             =0E78

This reveals that if the kforces among those are kept constant, i.e. Xi=const(i=1,2,,k,k<4), the fluxes corresponding to the left 4kforces are zero. Of course, if there are no restrictions on Xi(i=1,2,3,4), all the flues are zero. ( For convenience, we have modified the superscripts q,P,σ,fto be 1,2,3,4 ).

We may describe a ferroelectric phase transition by using the two paradigms above similarly as we have done in Section 3. To a first-order ferroelectric phase transition, the forces XP,Xσ,Xfof the region where the phase transition is occurring can be regarded as three large constants roughly in the characteristic times of phase transition ( i.e. the times in which the interface keeps rest ) because the electric displacement, the volume and the shape change suddenly. So, the flux Jqdiffof the region should be zero (but JP0,σ0,ρf0). This states clearly that the pure heat conduction and the heat conduction induced by the thermo-electric coupling and the thermo-mechanical coupling cancel out each other so as to release or absorb the latent heat. The phase transition occurs at the surface layer firstly, which is mechanically-free. So, when the phase transition occurs in this region, the flux σmaybe the nominal stress σnomonly, which does work to realize the transformation from internal energy to kinetic energy. When the phase transition occurs in the inner part, the flux σshould be the sum of σrealand σnombecause the sudden changes of the inner part’s volume and shape have to overcome the bound of outer part then σrealarises. The region where the phase transition is occurring, i.e. the phase boundary is accompanied with the real stress σrealusually, which does work to realize the transformation from kinetic energy to internal energy. This has been predicted and described with a propagating stress wave [20].

It is certain that the latent heat passes through the region where the phase transition has occurred ( at the outside of the region where the phase transition is occurring ) and exchange itself with the thermal bath. For ±lρva=Jqdiff=κT, a constant temperature gradient Tis kept in the region where the phase transition has occurrd, i.e. the force Xqat every site is a constant ( which does not change with the time but may vary with the position ). So, the fluxes JP=σ=ρf=0( but the flux Jqdiff0). This states clearly that the electric displacement Dwill not change but keep the value at Curie temperature or zero until the phase transition finishes and σreal=σnomin this region. Because the electric displacement Dand the strain ( or deformation ) are all determined by the crystal structure of system, JP=0reveals that Dof this region does not change so does not the crystal structure then does not the strain ( or deformation ). According to [20], we know the region where the phase transition has occurred is unstressed, i.e. σreal=0, then σnom=0. This reveals that the eigen ( or free ) strain ( or deformation ) of system induced by the thermo-electro-mechanical coupling of phase transition is complete and the change of it terminates before the phase transition finishes. The two deductions coincide with each other. σrealmay relaxes via the free surface.

The region where the phase transition will occur should be in equilibrium because there are no restrictions on the forces Xq,XP,Xσ,Xf. Whereas, according to [20], the region is stressed, i.e. σreal0. To the heating process of phase transition, this may lead to a change of the spontaneous polarization of this region because of the electro-mechanical coupling ( piezoelectric effect ).

An immediate result of the above irreversible thermodynamic description is that the action at a distance, which is the kind of heat transfer at phase transitions, is removed absolutely. The latent heat is transferred within a finite time so the occurrence of phase transition in the inner part is delayed. ( Of course, another cause is the stress, just see Section 5 ) In other words, the various parts absorb or release the latent heat at the various times. The action at a distance does not affect the phase transition necessarily.

## 5. Irreversibility: Thermal hysteresis and occurrences of domain structure

### 5.1. Thermal hysteresis

The “thermal hysteresis” of first-order ferroelectric phase transitions is an irreversible phenomenon obviously. But it was treated by using the equilibrium thermodynamics for ferroelectric phase transitions, the well-known Landau-Devonshire theory [2]. So, there is an inherent contradiction in this case. The system in which a first-order ferroelectric phase transition occurs is heterogeneous. The occurrences of phase transition in different parts are not at the same time. The phase transition occurs at the surface layer then in the inner part of system. According to the description above, we know a constant temperature gradient is kept in the region where the phase transition has occurred. The temperature of surface layer, which is usually regarded as the temperature of the whole system in experiments, must be higher ( or lower ) than the Curie temperature. This may lead to the thermal hysteresis.

No doubt that the shape and the area of surface can greatly affect the above processes. We may conclude that the thermal hysteresis can be reduced if the system has a larger specific surface and, the thermal hysteresis can be neglected if a finite system has an extremely-large specific surface. So, the thermal hysteresis is not an intrinsic property of the system.

The region where the phase transition will occur can be regarded as an equilibrium system for there are no restrictions on the forces Xq,XP,Xσ,Xf. In other words, the forces and the corresponding fluxes are zero in this region. To a system where a second-order ferroelectric phase transition occurs, the case is somewhat like that of the region where a first-order ferroelectric phase transition will occur. The spontaneous polarization, the volume and the shape of system are continuous at the Curie temperature and change with the infinitesimal magnitudes. This means Xq,XP,Xσ,Xfand Jqdiff,JP,σ,ρfcan be arbitrary infinitesimal magnitudes. The second-order phase transition occurs in every part of the system simultaneously, i.e. there is no the co-existence of two phases ( ferroelectric and paraelectric ). So, there is no the latent heat and stress. The thermal hysteresis disappears.

The region where a first-order ferroelectric phase transition will occur is stressed. This reveals that the occurrences of phase transition in the inner part have to overcome the bound of outer part, where the phase transition occurs earlier. This may lead to the delay of phase transition in the inner part.

### 5.2. Occurrences of domain structure

Though the rationalization of the existence of domain structures can be explained by the equilibrium thermodynamics, the evolving characteristics of domain occurrences in ferroelectrics can not be explained by it, but can be explained by the non-equilibrium thermodynamics.

In the region where the phase transition is occurring, the thermodynamic forces XP(= (φiT)), Xσ(= v(1T)), Xf(= vT) can be regarded as three large constants in the characteristic times of transition and the thermodynamic flux Jqdiff=0(but JP0, σ0, ρf0). The local entropy production (cf. Equation (63)) reduces to

σs=JP(φiT)σ:[v(1T)]ρfvTE79

Now, we are facing a set of complicated fields of T,  vand φi, respectively. Assume that the phase transition front is denoted by S. The points included in Sstand for the locations where the transition is occurring. Because the transition occurs along all directions from the outer part to the inner part, we may infer that the orientations of (φiT), or of v(1T)and or of vTvary continuously such that they are differently oriented at different locations.

There are always several ( at least two ) symmetry equivalent orientations in the prototype phase ( in most cases it is the high temperature phase ), which are the possible orientations for spontaneous polarization ( or spontaneous deformation or spontaneous displacement ). Therefore, the spontaneous polarization, the spontaneous deformation and the spontaneous displacement must take an appropriate orientation respectively to ensure σsis a positive minimum when the system transforms from the prototype ( paraelectric ) phase to the ferroelectric ( low temperature ) phase. The underlying reasons are that

JP=dDdt=d(ε0Ei+P)dtE80
σ=L:εE81
ρf=ρdvdtσE82

where L,  εand ε0are the modulus of rigidity, the strain and the permittivity of vacuum, respectively. Therefore, P,  εat different locations will be differently oriented. The domain structures in ferroelectrics thus occur.

It seems that the picture of domain occurrences for first-order ferroelectric phase transition systems should disappear when we face second-order ferroelectric phase transition systems. This is true if the transition processes proceed infinitely slowly as expounded by the equilibrium thermodynamics. But any actual process proceeds with finite rate, so it is irreversible. Then the above picture revives.

In [21], the domain occurrences in ferromagnetics can be described parallelly by analogy. And the case of ferroelastic domain occurrences is a reduced, simpler one compared with that of ferroelectrics or ferromagnetics.

It is well known that the Landau theory or the Curie principle tells us how to determine the symmetry change at a phase transition. A concise statement is as follows [22]: for a crystal undergoing a phase transition with a space-group symmetry reduction from G0to G, whereas Gdetermines the symmetry of transition parameter ( or vice versa ), it is the symmetry operations lost in going from G0to Gthat determine the domain structure in the low-symmetry phase. The ferroic phase transitions are the ones accompanied by a change of point group symmetry [23]. Therefore, the substitution of “point group” for the “space group” in the above statement will be adequate for ferroic phase transitions. From the above statement, the domain structure is a manifestation of the symmetry operations lost at the phase transition. In our treatment of the domain occurrences in ferroics, we took into account the finiteness of system ( i.e. existence of surface ) and the irreversibility of process ( asymmetry of time ). The finiteness of system make the thermodynamic forces such as (φiT),  v(1T),  vThave infinite space symmetry. The infinite space symmetry, combined with the asymmetry of time, reproduces the symmetry operations lost at the phase transition in the ferroic phase. It can be viewed as an embodiment of time-space symmetry.

After all, for the domain structures can exist in equilibrium systems, they are the equilibrium structures but not the dissipative ones, for the latter can only exist in systems far from equilibrium [24].

## 6. Concluding remarks

In order to overcome the shortcoming of Landau-Devonshire theory, the non-equilibrium thermodynamics was applied to study the ferroelectric phase transitions. The essence of transitions is the thermo-electro-mechanical coupling. Moreover, the irreversibility, namely thermal hysteresis and domain occurrences can be explained well in the realm of non-equilibrium thermodynamics.

The non-equilibrium thermodynamic approach utilized here is the linear thermodynamic one actully. In order to get the more adequate approaches, we should pay attention to the new developments of non-equilibrium thermodynamics. The thermodynamics with internal variables [25] and the extended ( irreversible ) thermodynamics [26] are two current ones. They all expand the fundamental variables spaces to describe the irreversible processes more adequately. Whereas, the relevant theoretical processing must be more complicated undoubtly. This situation needs very much effort.

## Acknowledgments

The work is supported by the Natural Science Foundation Program of Shandong Province, China ( Grant No. ZR2011AM019 ).

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2012 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Shu-Tao Ai (November 19th 2012). Advances in Thermodynamics of Ferroelectric Phase Transitions, Advances in Ferroelectrics, Aimé Peláiz Barranco, IntechOpen, DOI: 10.5772/52089. Available from:

### chapter statistics

1Crossref citations

### Related Content

Next chapter

#### Pb(Mg1/3Nb2/3)O3 (PMN) Relaxor: Dipole Glass or Nano-Domain Ferroelectric?

By Desheng Fu, Hiroki Taniguchi, Mitsuru Itoh and Shigeo Mori

First chapter

#### Raman Spectra of Soft Modes in Ferroelectric Crystals

By Takeshi Shigenari

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books