iMechanica - ferroelectric ceramics - Comments

atomistic simulation of cracks in ferroelectrics

Arash_Yavari — Wed, 30 Apr 2008 19:35:25 -0400

Dear Chad:

The best known and most common interatomic potentials for
ferroelectrics are shell potentials (here "shell" is very different
from what we usually understand as a shell in mechanics). These
potentials were introduced in the sixties by Dick and Overhauser:

Dick, B. G. and A. W., Overhauser [1964], Theory of the dielectric
constants of alkali halide crystals, Physical Review 112: 90-103.

In these potentials, each atom has a core (nucleus and inner electrons)
and a massless shell (valance electrons). The shell is assumed to be
spherical and uniformly charged. What enters into the energy expression
is the position of the center of the shell. Total energy has the
following three parts:

1) Short-range energy: Only shells contribute to this part and for the
most part this is a repulsive energy that prevents shells to get too
close to one another (Pauli repulsion).

2) Core-shell energy: In a given atoms core and shell interact usually
by a nonlinear spring (a fourth order polynomial in terms of the
relative distance). This prevents collapse of core and shell in the
same atom.

3) Electrostatic energy: This is the classical Coloumb potential. For a
given atoms shell and core do not interact electrostatically. This is
the troublesome part as long range interactions have to be very
carefully treated as the lattice sums representing energy and force are
conditionally convergent (Ewald summation method is usually used for
periodic systems).

In this model, spontaneous polarization (for example in the tetragonal
phase of PbTiO3) is mainly due to relative shifts of cores and shells
in the direction of polarization.

These potentials have been used by several groups in understanding the structure of domain walls, free surfaces, etc.

I have done some calculations for domain walls (surfaces on which
polarization is discontinuous, i.e. boundaries between different
energetically-equivalent variants in a single crystal). My experience
with these potentials is that they are very sensitive and the resulting
stiffness matrices could be close to ill-conditioned.

There have been analytic/semi-analytic lattice calculations for cracks
in the past mostly for idealized 2D lattices (an interesting observed
phenomenon was lattice trapping in early seventies). The main
contributers are Slepyan and Marder. However, as far as I can see it,
their methods cannot be used for analysis of something complicated like
PbTiO3 with its fairly complicated interatomic potential. My guess is
that the best approach would be MD.

Regards,
Arash

Reply to Arash

Chad Landis — Tue, 29 Apr 2008 14:25:08 -0400

Arash,

No I am not. I would be keenly interested in seeing something like this. I think that some of the methods that you are working with might be able to look at this problem. Perhaps you could post a reply detailing some of the issues that your methods would face while trying to address this.

Chad

atmistic simulations

Arash_Yavari — Sun, 27 Apr 2008 14:52:50 -0400

Dear Chad:

Are you aware of any atomistic calculations of fracture in ferroelectric crystals?

Arash

Reply to Kaushik

Chad Landis — Mon, 21 Apr 2008 11:37:06 -0400

Kaushik,

The question to answer is "what does the macro constitutive law need". If we take for example macro single crystal constitutive models (which are analogous to continuum slip plasticity models) we need to know the critical level of "driving force" required to cause switching on a given "transformation system", i.e. 180 or 90 degree switching for tetragonal ferroelectrics. My paper with Yu Su studied how charge defects affect the coercive field for 180 and 90 degree switching.

http://www.ae.utexas.edu/~landis/Landis/Research_files/JMPS2007.pdf

So one element that the phase field models might be able to provide to the macro models is how the coercive field is dependent on the concentration of defects. These phase field calculations have also suggested that the critical driving force can also depend on the other non-driving stresses and electric fields. For example, due to changes in piezoelectric properties across a 180 degree domain wall, the applied driving force for motion increases with combined electric field and tensile stress along the wall, suggesting that the observed critical electric field should decrease. However, the phase field calculations indicate that the critical driving force increases faster due to the tensile stress than the applied driving force, and so the critical electric field increases. This type of calculation suggests that the macro models should have a more complicated formulation than most (I think all) of them do.

Another part of the physics that is swept under the rug with the macro models is the distribution of residual stresses and electric fields associated with incompatible domain states. For example we see pictures of needle-like domain tips all of the time, and there are complicated fields associated with these. I have not thought about how to systematically study this (some type of lamination theory a la Jiangyu Li and Kaushik Bhattacharya is probably best suited for this problem), but it would be useful to have a better picture of how the residual energy evolves with the concentrations of the different variant types. The problem here is that I do not think there is a unique answer to this question. But then, I guess this is always what will happen when you represent a continuous set of internal variables (the positions of all domain walls) with a finite discrete set (volume concentrations of variant types). My point is that there are certain components of the macro models that can be informed by this type of information.

The short answer to your question would have been that there haven't been many links made between the macroscopic models and phase field models.

Chad

Thanks Ryan

Chad Landis — Mon, 21 Apr 2008 11:34:08 -0400

You are being kind to say I had the right idea. My "ideas" have been refined by your very informative post.

Chad

Chad has the right idea.

ellio167 — Mon, 21 Apr 2008 11:02:36 -0400

Chad has the right idea. Soft phonons are deformation modes in a
crystal that have a low stiffness associated with them. Thus,
they are likely to lead to instabilities. However, these modes
generally involve all atoms in the crystal. They can be described
as plane wave deformations or motions of the atoms. The
distinction between these and more classical continuum deformations,
that most mechanicers are used to, is that the wavelength of the motion
is of the same order of magnitude as the lattice spacing. For
deformations of on this scale, the discrete nature of the crystal leads
to different behavior (wave speed, stiffness, etc.) associated with
different wavelength phonons (eg. dispersion). This is in
contrast to the continuum case where no dispersion exists (there is no
length scale).

The distinction between optical and acoustic phonons is most clear in
the "long wavelength limit". On one hand, long wavelength
acoustic phonons can be directly correlated with "rank-one" uniform
deformations which mechanicers are well acquainted with. Long
wavelength optic phonons, on the other hand, are atomic motions that
correspond to localized excitation of the crystal where the atoms in
each unit cell move about in lock step (Cauchy-Born type
shifts/shuffles) but there is no macroscopic (uniform) deformation.

The phase transformation in Ferroelectrics is often driven by a soft
optic phonon. That is the optic phonon deformation mode is the
"primary order-parameter" (the mode that first goes unstable).
The transformation strains that result are "secondary order-parameters"
that are directly coupled to the primary order-parameter and occur
because of the optic phonon mode deformation. In Ferroelectrics
the primary oder-parameter corresponds to the shifting of the center
atom off of its central location (this happens in each unit cell, thus
illustrating the optic phonon nature of the deformation). As a
result of the atom moving from its central location the lattice spacing
changes which leads to the overall transformation strain.

This type of behavior is not limited to soft optic phonons. In
most Shape Memory Alloys it is a soft "zone-boundary phonon mode"
(neither optic or acoustic, really) that is responsible for the phase
transformation. This type of instability also involves a
"period-extension" where the size of the unit cell changes as a result
of the transformation. So that the martensite crystal structure
has a primitive unit cell with twice as many atoms, say, as the
austenite's primitive unit cell.

Dear Chad, As you've worked

Kaushik Dayal — Mon, 21 Apr 2008 05:07:45 -0400

Dear Chad,

As you've worked with both the phase-field and the macroscopic models, can you provide your perspective on what inputs phase-field calculations can provide to macroscopic models? I'm not aware of work that uses phase-field results in this manner directly, do you know of any references?

Thanks,
Kaushik

Reply to Arash

Chad Landis — Sun, 20 Apr 2008 11:55:22 -0400

Arash,  

In my opinion the best paper on linear piezoelectric fracture mechanics is by Suo and co-workers.  

Suo Z, Kuo CM, Barnett DM, Willis JR. 1992. Fracture mechanics for piezoelectric ceramics. J. Mech. Phys. Solids 40:739–65.    

Ortiz and co-workers have done some work on a cohesive zone model. In my opinion this is a good direction to go in, but more physical understanding is required to determine the coupling between the crack opening displacement and the electric potential drop across the crack faces. I would also be interested to see what such a model predicts for the simple static semi-infinite crack.

Arias I, Serebrinsky S,Ortiz M, 2006. A phenomenological cohesive model of ferroelectric fatigue. Acta Materialia 54:975-984.      

For a more selfish response to your comment, I would refer you to my work on the electromechanical boundary conditions and on small scale domain switching,

http://www.ae.utexas.edu/~landis/Landis/Research_files/IJSS2004b.pdf

http://www.ae.utexas.edu/~landis/Landis/Research_files/ActaMat2004.pdf

http://www.ae.utexas.edu/~landis/Landis/Research_files/JOMMS2006.pdf

If you have a specific topic you are interested in, linear, non-linear, boundary conditions, etc.,then let me know and I can dig up some more specific references.

Chad

Soft optical phonons

Chad Landis — Sun, 20 Apr 2008 11:52:16 -0400

OK, I'll have a go at this, and my understanding comes from some work by Nick Triantafyllidis, John Shaw and Ryan Elliot.

The idea is that the phase transition occurs due to some instability in the lattice. I think that the easiest type of transition for mechanic-ers to understand is one that happens as a result of long wavelength perturbations, or essentially homogeneous deformations of the lattice. For this type of transition one of the eignevalues of the tangent elasticity matrix goes to zero and a transition can occur. This is not the soft optical phonon. However, this is not the only type of instability that can occur. The lattice must also be stable against short wavelength perturbations. My rough mechanistic understanding of this type of instability is that the perturbation of one atom, or a few atoms, from its equilibrium position causes a cascade of the remainder of the crystal to a new lower energy state. I believe that this type of instability is the of the "soft optical phonon" type.

I am really stretching the limits of my knowledge on this, but I think this explanation has some elements of truth to it.

I welcome any additional responses that directly or indirectly point out my ignorance here.

Chad

Reply to Ravi

Chad Landis — Sun, 20 Apr 2008 10:48:09 -0400

Ravi,

Let me be a little more careful. Aside from a recent attempt by Gerold Schneider and his group, the energy release rate has never been unambiguously measured in a piezoelectric/ferroelectric fracture experiment. So when I talk about negative energy release rates I am referring to predictions from linear piezoelectric fracture mechanics (LPFM) with impermeable or semi-permeable crack face boundary conditions. So one point is that the negative energy release rate contributions predicted by LPFM may not actually exist in the experiments. Some possible mechanisms for a reduction in the retarding effects of electric fields include electrical discharge in the crack gap, grain bridges in the crack gap, and large scale domain switching. I have done some theoretical work on the effects of electrical discharge.

http://www.ae.utexas.edu/~landis/Landis/Research_files/IJSS2004b.pdf

In calculations the dielectric permittivity of the crack is usually taken to be that of free space (or air). This assumption is not a requirement in the models, but it is what is most commonly used. Recently Scheinder et al have measured the opening displacement and electric potential drop across a crack and found that the permittivity is higher than that of free space. Again, some possible mechanisms for this observation include crack bridges (although they did not observe any) and electrical discharge.

Schneider GA, Felten F, McMeeking RM. The electrical potential difference across cracks in PZT measured by Kelvin probe microscopy and the implications for fracture. Acta Materialia 51:2235-2241, 2003.

I am not an expert on FRAMs, but my understanding is that defects, generally charge defects in the form of oxygen vacancies or free electrons, are bad. Charge defects act to do two things, 1) they pin the motion of domain walls which leads to an increase in the coercive field and 2) they shield the polarization of domains leading to a reduced current signal during switching. The migration of charge defects leads to the phenomenon of electrical fatigue (very different from mechanical fatigue) in these materials.

Scott JF, Dawber M. Oxygen-vacancy ordering as a fatigue mechanism in perovskite ferroelectrics. Applied Physics Letters 76:3801-3803, 2000.

I will get back to you on the optical phonon question. I do believe that there is a "mechanics" interpretation, but I need to make sure I get my terminology right.

Chad

fracture mechanics papers

Arash_Yavari — Fri, 18 Apr 2008 19:30:26 -0400

Dear Chad:

Thanks for the very interesting discussion. Could you mention a few key
papers that discuss the theoretical issues in fracture mechanics of
ferroelectrics?

Thanks,
Arash

Hi Dr Landis This is a

ravi maranganti — Fri, 18 Apr 2008 19:14:53 -0400

Hi Dr Landis

This is a pretty interesting topic...i was intrigued by the possibility of using electric fields to retard crack growth in ferroelectrics. I do have a couple of questions...

* Is my understanding correct that the "negative" energy release rate should be higher for high dielectric materials like certain perovskites? Also, i noticed in your article that experimental observations have not confirmed the idea of a negative energy release rate due to applied electric fields working against the mechanical positive energy release rate to retard crack growth...while this is a little disheartening, has any body put forward some theoretical explanation for why the effects of the electric field are not being observed?

* In calculations, what is the dielectric permittivity of the crack taken to be? Talking about FRAMs...how defect free should they be?

* On a side note, the phase transition of a ferroelectric from its cubic centrosymmetric phase to a lower symmetry phase is associated with the soft optical phonon becoming unstable at the curie temperature. Is there a "mechanics " interpretation behind this?

Thanks!