Crystallography & Transformations
Crystallography & Transformations
'
I
II
I
I
|
‘ll
_
5%
_“_H_'HW_
J)‘
VWQHM’
M
WP
’_'
w,
M1
I N
My WW
_'>_“
_
___I1
~.&M%M‘W_%.
§;_m)¥M
‘
ll“
in
“_'M_
W _
Il
H M
V_
I
‘W
W
"M
“_Pm
_
>_(_
WM“
~v
‘§V_____’
w____
H _
3M
I_
7Y
‘ V
if
g‘
% _ A
WW
v_
g
“_ A
“__
x‘_ ‘“
an
7;
~‘Q
‘Wm
_v
_
M,
Kh
_
xW
'
_INI_\v’)AMW
§‘X
)_z¥§§§_
_
‘ xq
>V
§§A§g§
aY
_
Y
H
WW"
_
M”
W
*’A“
__
> _
WWW
I
v>‘\
\
MW
WW‘ N
__
A_
“W
iU
k’ w
}_ W
%M‘
NW v
__ __
v‘F
Hp’
‘in _(
)\_%$
' >
H\V
H"
IIL
‘\1AM\
1H“
h
NM
V
AWx
MM
A
I *‘ 6
um
U
__
2
_2%
< k
‘W'
VV
‘ y _
M
M“
“V
W
U“
A"
L_
“HA
WA
I‘ H
_
ax~
N
~W_
w
__>___
WW"
_ _
_W_
NW
_ _ >
d
M
lI
}<
~
WWW Y
*My
W _
WY
W
mu
W_
H
(v
x~‘_
0WW
V
;i‘\
_;
Ni‘
_v
_>
>_\
_\
l“_'_l__ V
%mm0H“§§%%§§_‘§\s
\3\!_\_
M
MW
~_‘
v_
‘_
_v§VvvYALV§§kg!"x§€v}_h\_‘_z_‘$S;;_:§23$
__
V
__
Harshad K. D. H. Bhadeshia
Geometry of crystals,
polycrystals, and phase
transformations
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts
have been made to publish reliable data and information, but the author and publisher cannot assume
responsibility for the validity of all materials or the consequences of their use. The authors and
publishers have attempted to trace the copyright holders of all material reproduced in this publication
and apologize to copyright holders if permission to publish in this form has not been obtained. If any
copyright material has not been acknowledged please write and let us know so we may rectify in any
future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information storage
or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access
www.copyright.com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC),
222 Rosewood Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that
provides licenses and registration for a variety of users. For organizations that have been granted a
photocopy license by the CCC, a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are
used only for identification and explanation without intent to infringe.
Author’ xv
Acronyms xvii
I Basic Crystallography 1
1 Introduction and Point Groups 3
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 The lattice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Primitive representation of Cubic-F . . . . . . . . . . . . . . . . 7
1.3 Bravais lattices . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4 Directions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Number of equivalent indices . . . . . . . . . . . . . . . . . . . . 15
1.5 Planes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.6 Weiss zone law . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.7 Symmetry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.8 Symmetry operations . . . . . . . . . . . . . . . . . . . . . . . . . 17
Five-fold rotation . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.9 Crystal structure . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Structure of graphene . . . . . . . . . . . . . . . . . . . . . . . . 22
1.10 Point group symmetry . . . . . . . . . . . . . . . . . . . . . . . . 23
Point symmetry of chess pieces . . . . . . . . . . . . . . . . . . . 24
Octahedral interstices in iron . . . . . . . . . . . . . . . . . . . 26
1.11 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
References 29
2 Stereographic Projections 31
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Projection of small circle . . . . . . . . . . . . . . . . . . . . . . 31
2.2 Utility of stereographic projections . . . . . . . . . . . . . . . . . 34
2.3 Stereographic projection: construction and characteristics . . . . . 35
Radius of trace of great circle on Wulff net . . . . . . . . . . . . . 39
Traces of plates . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.4 Stereographic representation of point groups . . . . . . . . . . . . 42
vii
viii Contents
References 49
References 59
4 Space Groups 61
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
4.2 Screw axes and glide planes . . . . . . . . . . . . . . . . . . . . . 62
4.3 Cuprite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.4 Location of atoms in cuprite cell . . . . . . . . . . . . . . . . . . . 64
Space group of Fe-Si-U compound . . . . . . . . . . . . . . . . . 66
Cementite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
Diamond and zinc sulphide . . . . . . . . . . . . . . . . . . . . . 69
4.5 Shape of precipitates . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.6 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
References 73
References 89
6.2 Texture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
6.3 Orientation distribution functions . . . . . . . . . . . . . . . . . . 99
Euler angles relating two frames . . . . . . . . . . . . . . . . . . 100
6.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
References 103
References 115
References 125
References 149
x Contents
References 169
References 201
12 Martensite 203
12.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
12.2 Shape deformation . . . . . . . . . . . . . . . . . . . . . . . . . . 203
12.3 Interfacial structure of martensite . . . . . . . . . . . . . . . . . . 206
12.4 Phenomenological theory of martensite crystallography . . . . . . 208
12.5 Stage 1: Calculation of lattice transformation strain . . . . . . . . . 211
Determination of lattice transformation strain . . . . . . . . . . . 212
12.6 Stage 2: Determination of the orientation relationship . . . . . . . . 215
Martensite-austenite orientation relationship . . . . . . . . . . . 215
12.7 Stage 3: Nature of the shape deformation . . . . . . . . . . . . . . 217
Habit plane and the shape deformation . . . . . . . . . . . . . . 218
12.8 Stage 4: Nature of the lattice-invariant shear . . . . . . . . . . . . . 220
Lattice–invariant shear . . . . . . . . . . . . . . . . . . . . . . . 220
12.9 Texture due to displacive transformations . . . . . . . . . . . . . . 222
12.10 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
Contents xi
References 227
13 Interfaces 231
13.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 231
13.2 Misfit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 232
Symmetrical tilt boundary . . . . . . . . . . . . . . . . . . . . . 234
Interface between alpha and beta brass . . . . . . . . . . . . . . 237
13.3 Coincidence site lattices . . . . . . . . . . . . . . . . . . . . . . . 239
Coincidence site lattices . . . . . . . . . . . . . . . . . . . . . . 239
Symmetry and the axis-angle representations of CSL’s . . . . . . 242
13.4 The O-lattice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 246
alpha/beta brass interface using O-lattice theory . . . . . . . . . 248
13.5 Secondary dislocations . . . . . . . . . . . . . . . . . . . . . . . . 249
Intrinsic secondary dislocations . . . . . . . . . . . . . . . . . . 249
13.6 The DSC lattice . . . . . . . . . . . . . . . . . . . . . . . . . . . . 251
13.7 Some difficulties associated with interface theory . . . . . . . . . . 253
13.8 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
References 257
Appendices 259
To state the obvious, crystals contain order. Even when that is disturbed locally, the
perturbations themselves may sometimes form regular patterns. This is evident in
the structure of interfaces where two or more crystals meet in an apparently hap-
hazard manner. Disorder can sometimes be ignored without compromising some of
the consequences of long-range order. Crystalline solid solutions in which atoms
are dispersed at random would fail the strict definition of long-range order, but they
nevertheless show the characteristics of crystals when probed by X-rays. We shall at-
tempt in this book to understand not only the elegance of individual crystals, but also
of clusters of space-filling crystals and transformations between crystalline phases.
Crystallography is hardly a new subject so there are numerous books available,
many of which are beacons of scholarship. So why another text on this much mooted
topic? First, as Buckminster Fuller pointed out, there has been a massive expansion in
human knowledge, with the process continuing at an unabating pace. New subjects
spring up with notorious regularity and some of these have come to be regarded
as essential first steps in the higher education curriculum. Modern undergraduate
students are therefore faced with a much greater palette of distinct courses than has
been the case in the past.
This book is partitioned into two, the first part of which is meant to be self-
contained and deals with what I feel is essential learning for any student in the mate-
rial sciences, physics, chemistry, earth sciences, and the natural sciences in general.
This part, intentionally concise, is based on a set of nine lectures that I give annually
to undergraduate students of the Natural Sciences Tripos in Cambridge University. It
is of generic value and has just sufficient material to deliver concepts. It covers crys-
tals, polycrystals, interfaces and transformations that occur by the disciplined motion
of atoms. I feel that most books on crystallography are too detailed to accommodate
within the schedule of a contemporary undergraduate.
The second part has depth which would be appreciated most in the context of
research. It is a development from a book that I taught and published in 1987 on the
Geometry of Crystals. The treatment is limited to phenomena dominated by crystal-
lography.
The book contains worked examples throughout because crystallography is a sub-
ject that thrives on practice. Video lectures and other electronic materials that can
enhance the content of this book are available freely on
http://www.phase-trans.msm.cam.ac.uk/teaching.html
Crystallography has rules, established by convention, but many of these have ex-
ceptions and there sometimes are multiple conventions. Donnay in his 1943 paper
xiii
xiv Preface and Acknowledgments
H. K. D. H. Bhadeshia
Cambridge
Author’
xv
Acronyms
xvii
Part I
Basic Crystallography
1
Introduction and Point Groups
Abstract
1.1 Introduction
Amorphous solids are homogeneous and isotropic because there is no long range
order or periodicity in their internal atomic arrangement. In contrast, the crystalline
state is characterised by a regular arrangement of atoms over large distances. Crystals
are therefore anisotropic - their properties vary with direction. For example, the inter-
atomic spacing varies with orientation within the crystal, as does the elastic response
to an applied stress.
Crystals can be two or three dimensional; graphene is a two-dimensional crystal if
we neglect the fact that it is one atom thick (Fig. 1.1a), whereas polonium (Fig. 1.1b)
is obviously three-dimensional. There is talk about one-dimensional crystals where
a string of equally spaced atoms is in a straight line, made by confining the atoms
into a carbon nanotube; however, there are interactions between these atoms and the
nanotube which violate the definition.1 Four dimensional crystals which are periodic
in both time and space have been postulated [2]. For example, a cylindrical crystal
with a periodic arrangement of atoms that repeats as time is stepped. Such a structure
would require perpetual motion, i.e. must not radiate rotational energy so there are
difficulties in explaining how such a crystal could exist in a non-stationary ground
state, a subject of discussion in the associated literature. Some experimental evidence
1 To quote Ziman [1], “the hypothetical ordered linear chain is physically unrealisable as a genuine
one-dimensional system”.
3
4 1 Introduction and Point Groups
now exists for discrete time crystals [3, 4]; the question remains as to whether the
‘time crystals’ are everlasting [5].
(a) (b)
FIGURE 1.1
(a) The arrangement of carbon atoms in the two-dimensional structure of graphene.
(b) The arrangement of atoms in the three-dimensional structure of polonium.
The crystalline form usually conjures visions of beautiful facetted forms, but it
is worth noting at the outset that they can be solid or liquid,2 and can have arbitrary
shapes. The familiar elegant shapes are a consequence of the minimisation of surface
energy but in engineering applications the crystals are produced in functional shapes
(Fig. 1.2).
The beautiful shapes of many crystals that are found in nature represent an at-
tempt at the minimisation of surface energy per unit volume of material when at
equilibrium, or because the rate of growth varies with direction. For isotropic mate-
rials, a sphere would represent the equilibrium shape but crystals are not isotropic;
some planes of atoms have a lower surface energy than others so there will be a ten-
dency to maximise the area of those planes, in which case the equilibrium shape is
no longer spherical [7].
Engineering materials usually are space-filling aggregates of many crystals of
varying sizes and shapes; these polycrystalline materials have properties which de-
pend on the nature of the individual crystals, but also on aggregate properties such
as the size and shape distributions of the crystals, and the orientation relationships
between the individual crystals. The degree of randomness in the orientations of the
crystals relative to a fixed frame is a measure of texture, which has to be controlled
in the manufacture of transformer steels, uranium fuel rods and beverage cans.
Connecting a pair of crystals together to form a bi-crystal requires a process in-
volving at least five degrees of freedom; two for the choice of the plane along which
the pair are connected and three for the relative orientations of the two crystals. The
connecting plane is referred to as a ‘boundary’ or an ‘interface’ in general. The crys-
tallography of interfaces connecting adjacent crystals can determine the deformation
2 Liquid crystals contain order that is associated with anisotropic molecules that can flow past each
other and yet maintain orientational order [6]. Normal liquids are isotropic.
1.1 Introduction 5
(a) (b)
FIGURE 1.2
(a) Clusters of naturally occurring pyrite crystals (FeS2 ) from Peru. The structure
belongs to the cubic crystal system. We shall see later that the octahedral shapes with
triangular facets are consistent with the defining symmetry of a cube. (b) A nickel
alloy in the shape of a turbine blade; the part above the spiral is the single crystal. Af-
ter removal of the appendages associated with the manufacturing process, the shape
is consistent with the aerodynamic form required for a turbine engine. Blades of this
kind serve in harsh environments and are subjected to large stresses. This can lead to
slow, stress-driven lengthening of the blade driven by the movement of atoms, even-
tually making it unfit for service. In polycrystalline blades, the boundaries between
the crystals are easy paths for atom movement and hence have shorter service lives.
6 1 Introduction and Point Groups
(a) (b)
FIGURE 1.3
(a) Definition of a unit cell; the lattice points are drawn as blue spheres for clarity
and the red vector is a lattice vector joining two lattice points. (b) The stacking of
unit cells to produce a periodic pattern; the stacking can be continued indefinitely to
fill the space available without leaving unfilled regions.
vectors. For a particular vector, this set of components defines the Miller indices of
that vector.3
Our initial choice of the basis vectors was arbitrary since there is an infinite num-
ber of lattice vectors which could have been used in defining the unit cell. The pre-
ferred choice includes small basis vectors which are as equal as possible, provided
the shape of the cell reflects the essential symmetry of the lattice.
001F
010 F
100 F
FIGURE 1.4
The conventional face-centred cubic unit cell is in black and its primitive version
in blue.
which any crystal within that system must possess as a minimum requirement (Ta-
ble 1.1). The Bravais lattices are illustrated in Fig. 1.5.
TABLE 1.1
The crystal systems. The symmetry axes monad, diad, triad, tetrad and hexad involve
rotations of 360◦ , 180◦ , 120◦ , 90◦ and 60◦ respectively. These particular rotations
leave the lattice in a state that cannot be distinguished from its starting configuration.
ai are the magnitudes of the basis vectors ai (which form a right-handed set). The
defining symmetry cannot be achieved unless the conditions highlighted in red are
satisfied.
of effects such as diffraction and electrical properties, the periodicity of the particles
deposited on substrates is influential.
The structure of a three-dimensional crystal is in general assumed to persist to
the very surfaces of the crystals, but the forces there are different from those in the
bulk of the material. As a consequence, surfaces are known to reconstruct into one of
the five lattice types illustrated in Fig. 1.7. This of course is important in processes
that rely on surface structure, such as catalysis [12].
1.3 Bravais lattices 11
FIGURE 1.5
The fourteen three-dimensional Bravais lattices.
12 1 Introduction and Point Groups
FIGURE 1.6
Projection of the three-dimensional shapes of the fourteen Bravais lattices on to the
a1 -a2 plane. The numbers indicate the coordinates of lattice points relative to the a3
axis; the unlabelled lattice points are by implication located at coordinates 0 and 1
with respect to the a3 axis.
1.3 Bravais lattices 13
FIGURE 1.7
The five two-dimensional lattices.
14 1 Introduction and Point Groups
1.4 Directions
A vector u can be represented as a linear combination of the basis vectors ai of the
unit cell (i = 1, 2, 3):
u = u1 a1 + u2 a2 + u3 a3 (1.1)
and the scalar quantities u1 , u2 and u3 are the components of the vector u with
respect to the basis vectors a1 , a2 and a3 . Once the unit cell is defined, any direction
u within the lattice can be identified uniquely by its components [u1 u2 u3 ], and the
components are called the Miller indices of that direction. They are by convention
enclosed in square brackets (Fig. 1.8).
(a) (b)
(c)
FIGURE 1.8
(a,b) Miller indices for directions and planes. (c) Notice that in non-cubic systems, a
direction with indices identical to those of a plane is not necessarily normal to that
plane.
It is sometimes the case that the properties along two or more different directions
are identical. These directions are said to be equivalent and the crystal is said to
possess symmetry. For example, the [1 0 0] direction for a cubic lattice is equivalent
to the [0 1 0], [0 0 1], [0 1 0], [0 0 1] and [1 0 0] directions; the bar on top of the
number implies that the index is negative.
The indices of directions of the same form are conventionally enclosed in special
1.6 Planes 15
brackets, e.g. ⟨1 0 0⟩. The number of equivalent directions within the form is called
the multiplicity of that direction, which in this case is 6.
Example 1.2: Number of equivalent indices
Given a direction with indices [u v w] such that u ̸ =v ̸=w ̸ = 0, calculate the
number of equivalent directions possible, for the following lattice types: (a) cubic,
(b) tetragonal.
(a) In the cubic system where a1 = a2 = a3 , there are six choices for the first
index, i.e. u, −u, v, −v, w, −w. Suppose that v is selected as the first in-
dex, then the choice for the second index is limited to one of four options
u, −u, w, −w. If the second index is −u, then the choice for the final in-
dex becomes w or −w. The number of equivalent directions of the form
⟨u v w⟩, 6 × 4 × 2 = 48. However, if antiparallel options such as [u v w]
and [−u − v − w] are not counted then the number of equivalent directions
reduces to 24.
Note that if u = v then the number of equivalent directions reduces to 6 ×
1 × 2 = 12 because the choice of the first index leaves only one choice for the
second. An index that is zero also reduces the options by a factor of two since
there is then no negative or positive aspect.Therefore, in the cubic system,
there are 6 variations of ⟨1 1 0⟩: [1 1 0], [1 1 0], [0 1 1], [0 1 1], [1 0 1] and
[1 0 1] if antiparallel directions are not counted.
(b) For the tetragonal system, a1 = a2 ̸ = a3 . Therefore, the choice of the first
index is limited to four options, u, −u, v, −v,; if u is selected as the first
index then the choice of the second is just two (v, −v) and of the third w
or −w. Therefore, the number of equivalent directions of the form ⟨u v w⟩
becomes 6 × 2 × 2 = 24, reducing to 12 if antiparallel options are removed.
1.5 Planes
If a plane intersects the a1 , a2 and a3 axes at distances x1 , x2 and x3 respectively,
relative to the origin, then the Miller indices of that plane are given by (h1 h2 h3 )
where:
h1 = φa1 /x1 , h2 = φa2 /x2 , h3 = φa3 /x3 (1.2)
φ is a scalar which clears the numbers hi off fractions or common factors. Note that
xi are negative when measured in the −ai directions. The intercept of the plane with
an axis may occur at infinity (∞), in which case the plane is parallel to that axis and
the corresponding Miller index will be zero (Fig. 1.8).
Miller indices for planes are by convention written using round brackets:
(h1 h2 h3 ) with braces being used to indicate planes of the same form: {h1 h2 h3 }.
16 1 Introduction and Point Groups
u 1 h1 + u 2 h2 + u 3 h3 = 0 (1.3)
and the law applies to any crystal system. This will be proven when we deal with the
reciprocal lattice (Chapter 5).
FIGURE 1.9
A zone refers to a set of planes which share a common
direction, which in turn is known as a zone axis. The
direction illustrated would satisfy the Weiss law for
all the planes shown. For example, the planes (1 1 1),
(1 1 2) and (1 1 0) all belong to the zone [1 1 0].
1.7 Symmetry
Although the properties of a crystal can be anisotropic, there may be different direc-
tions along which they are identical. These directions are said to be equivalent and
the crystal is said to possess symmetry.
That a particular edge of a cube cannot be distinguished from any other is a
measure of its symmetry; an orthorhombic parallelepiped has lower symmetry, since
its edges can be distinguished by length.
Some symmetry operations are illustrated in Fig. 1.10; in essence, they transform
a spatial arrangement into another that is indistinguishable from the original. The
rotation of a cubic lattice through 90◦ about an axis along the edge of the unit cell is
an example of a symmetry operation, since the lattice points of the final and original
lattice coincide in space and cannot consequently be distinguished.
We have encountered translational symmetry when defining the lattice; since the
environment of each lattice point is identical, translation between lattice points has
the effect of shifting the origin. Dislocations whose Burgers vectors are lattice vectors
therefore accomplish slip without changing the crystal structure. Deformations like
these are known as lattice-invariant deformations because they accomplish a shape
1.8 Symmetry operations 17
change without altering the nature of the lattice. In contrast, partial dislocations have
Burgers vectors that are not lattice vectors – they therefore change the local stacking
fault sequence of the planes on which they glide.
FIGURE 1.10
Some symmetry operations. The mirror
6 6 6 6 6 plane involves reflection symmetry,
translation whereas the glide plane is a combination
of reflection and translation parallel to the
6 mirror. Translational symmetry is intrinsic
6 to all lattices whereas the other elements
may or may not exist for all lattices. For
mirror example, the triclinic lattice has no mirror
plane.
rotation
6 Note that the translations involved in
the screw axis and glide plane are rational
6 fractions of the repeat distance along
glide
the translation direction. The screw axis
6 illustrated involves a rotation of 180◦ (≡
screw diad diad) combined with a translation that is
half of the repeat distance parallel to the
axis.
tation through 120◦ combined with an inversion, the axis being labelled 3. Some of
the symbols representing rotation and rotoinversion axis are as follows, where 2 is
omitted because it is, as will be shown in Chapter 2, equivalent to a mirror operation:
Cl
centre of
CH3 symmetry
H H
CH 3
(a) Cl (b)
FIGURE 1.11
Illustration of centre of symmetry. (a) A molecule with the centre of symmetry iden-
tified by the dot. (b) The common point where the arrows intersect is the centre of
symmetry.
the length of the diagonal concerned, except for diamond glide, where it is a quarter
of the diagonal length.
Example 1.3: Five-fold rotation
Show that a five-fold rotation axis is inconsistent with the translational sym-
metry of a lattice.
In Fig. 1.12, the a is the repeat distance between adjacent lattice points. The
operation of a five-fold rotation axis by θ = 72◦ would lead to lattice points sep-
arated by a distance x. In order for this rotation to be a symmetry operation, x
must equal na, where n is an integer. Since x = a − 2a cos θ, the governing equa-
tion that defines possible values of θ that are consistent with symmetry operations
become:
x = na = a − 2a cos θ
1−n
cos θ = with − 1 ≤ cos θ ≤ 1 (1.4)
2
It follows that for cos θ = −1, − 12 , 0, 12 and 1, θ takes the values 180◦,
120◦ , 90◦ , 60◦ corresponding to two-fold, three-fold, four-fold and six-fold rota-
tion axes. Equation 1.4 shows, therefore, that a five-, or indeed seven-fold rotations
are not symmetry axes. Another way of approaching this problem is on page 112.
FIGURE 1.12
For θ to represent a rotation consistent with the translational symmetry of the lat-
tice, x must equal an integer multiplied by the translation a. The proposed rotation
axis is normal to the plane of the figure.
We note in passing that ordered structures exist that are not periodic but can
yield diffraction patterns that have five-fold or ten-fold symmetries. These struc-
tures are known as quasicrystals; they lack translational symmetry [13].
FIGURE 1.13
Building up the crystal structure of β-brass by placing a motif consisting of an ap-
propriate pair of Cu and Zn atoms at each lattice point of a primitive cubic lattice. (a)
represents the primitive lattice with the motif consisting of a pair of distinct atoms
that have yet to be placed. (b-e) The motif is placed at each lattice point. (f) The
three-dimensional representation of the final structure.
The location of an atom of copper at each lattice point of a cubic-F lattice, gen-
erates the crystal structure of copper, with four copper atoms per unit cell (one per
lattice point), representing the actual arrangement of copper atoms in space.
Consider now a motif consisting of a pair of carbon atoms, with coordinates
[0 0 0] and [ 14 14 41 ] relative to a lattice point. Placing this motif at each lattice point
of the cubic-F lattice generates the diamond crystal structure (Fig. 1.14), with each
unit cell containing 8 carbon atoms (2 carbon atoms per lattice point). The figure also
shows the tetrahedral bonding which makes diamond a giant molecule with strong
covalent bonds in three dimensions. Diamonds do contain dislocations [14] but they
have complex cores, are dissociated and incredibly difficult to move, making it very
hard. The structure of silicon is obtained simply by replacing the carbon atoms by sil-
icon (cubic-F, motif of a pair of Si atoms at 0,0,0 and 41 , 14 , 14 ); the strong directional
bonding also results in a large Peierls barrier to dislocation mobility [15], making
silicon single-crystals rather brittle.
Fig. 1.15 shows how a Cubic-F cell changes into a Cubic-P cell when the nickel
and aluminium atoms order in the classical γ/γ ′ superalloy system.
1.9 Crystal structure 21
(a) (b)
(c)
FIGURE 1.14
(a) Projection of the cubic-F lattice. (b) Projection of cubic-F lattice with a motif of a
pair of carbon atoms at 0, 0, 0 and 14 , 14 , 14 placed at each lattice point. (c) Perspective
of the same structure illustrating the tetrahedral bonding of the carbon atoms.
(a) (b)
FIGURE 1.15
Nickel based superalloys. (a) The face-centred cubic crystal structure of disordered
γ. (b) The primitive cubic crystal structure of γ ′ .
22 1 Introduction and Point Groups
0.142 nm
FIGURE 1.16
Atomic structure of graphene sheet without
defects. Entropy effects mean that optically
visible samples cannot ever be made pristine
[16].
The unit cell, identified blue in Fig. 1.17 is primitive hexagonal with a motif
of two carbon atoms per lattice point. The environment about the black atom is
not identical to that around the red atom. The lattice parameter is 0.142 × 2 ×
cos 30◦ = 0.245 nm.
Copper is cubic close-packed. The arrangement of atoms on the {111} close-
packed plane is hexagonal, primitive. The two dimensional cell has the edges
equal to a2 ⟨110⟩, so the lattice parameter is 0.255 nm. This closely matches that of
graphene.
Given that the unit cell edges of the two-dimensional structures of graphene
and copper match to 4% (100[0.255 − 0.245]/0.255), it is logical that their unit
cell edges are exactly aligned, [10]graphene ∥ [10]Cu and [01]graphene ∥ [01]Cu .
1.10 Point group symmetry 23
FIGURE 1.17
The blue lines represent the unit cell, with
a motif of a pair of carbon atoms coloured
black and red, per lattice point.
(a) (b)
FIGURE 1.18
(a) A molecule with point group symmetry 2m (diad and a mirror plane parallel to the
diad). (b) Vibration modes (symmetrical stretch, asymmetrical stretch and bending)
of the molecule when appropriately stimulated. There is an additional mirror plane
passing through the centres of all three atoms if account is taken of the electron
clouds associated with the molecule.
ror planes, centre of symmetry and inversion axes are permitted. There are 32 point
groups in three dimensions, classified within the seven crystal classes (Table 1.2).
The triclinic, monoclinic and orthorhombic groups do not contain triads, tetrads or
hexads. For those systems, when the point group symbol contains three elements
(e.g. 2mm), then the symbols are presented in the order of the symmetry elements
parallel to the x, y and z axes respectively, as illustrated in Fig. 1.19a for an object
in the orthorhombic class.
TABLE 1.2
Point group symmetries associated with the seven crystal classes.
For crystal systems with higher order axes, the z direction is assigned to that
higher order axis, with the second symbol corresponding to equivalent secondary
axes that are normal to z, and the third also normal to z to equivalent tertiary direc-
tions passing between the secondary ones. Fig. 1.19b shows an example for an object
in the tetragonal class, with the tetrad placed along the z axis, the two mirror planes
with normals along the x and y axes and the additional mirror planes generated by
this symmetry also illustrated.
In the case of the cubic system the triad is always the second symbol given that
the defining symmetry is four triads.
Some further details on the notation are as follows:
• m ≡ 2 represents a mirror plane.
• 2, 3, 4, 6 are 2-fold, 3-fold, 4-fold and 6-fold rotation axes, respectively.
• 1, 2, 3 etc. are inversion axes; the one fold inversion axis is equivalent to a centre
of symmetry; 2 signifies a rotation of 360◦/2 combined with an inversion through
the centre point.
4
• m mm is a 4-fold axis with a mirror normal to it and four mirror planes containing
it. The symbol is usually written without the space as 4/mmm.
• 432 refers to a point group with 3-fold axes which are not parallel to the z axis,
where a unit cell has x, y and z axes. This is limited to the cubic system.
1.10 Point group symmetry 25
zxy
xyz additional
mirrors 4mm
2mm
z
x
y
x
y
z
mirrors
(a) mirrors (b)
FIGURE 1.19
Convention for point group notation. (a) Groups without high order axes. (b) Groups
with high order axes.
The bishop has a three-fold axis passing through the cut corner of the cube
and a mirror parallel to that axis, so the point group symmetry is 3m which fits
the trigonal crystal class. The castle clearly has a single tetrad with two mirrors
parallel and one normal to the tetrad so in this case, the point group is 4/mmm
which is tetragonal class (Table 1.2).
(a) (b)
FIGURE 1.20
Avant garde chess pieces. (a) Bishop. (b) Castle.
Some of the point group notation comes from the classification of the macro-
scopic shapes of crystals. If a crystal exhibits well-developed faces, its point group
symmetry can be derived from its external form. Fig. 1.21 illustrates this for two
cases. In the case of the gypsum there is a diad normal to (010) and a mirror plane
normal to that diad so that the shape has the point group 2/m belonging to mono-
26 1 Introduction and Point Groups
clinic class. With the epsomite there is a diad along the vertical axis of the diagram
and two diads perpendicular to the vertical axis passing through the vertical edges
between the planes (110) and (110), i.e. it belongs to the orthorhombic crystal class.
With appropriate analysis of facetted crystals it therefore becomes possible to
deduce the crystal class simply from symmetry. Fig. 1.22 shows two such examples.
The quartz crystal clearly has a hexad which is the defining symmetry of the hexag-
onal crystal class. The celestite has three mutually perpendicular diads, placing it in
the orthorhombic system.
FIGURE 1.21
(a) Shape of gypsum
(CaSO4 .2H2 O) with point
group symmetry 2/m. (b) Shape
of epsomite (MgSO4 .7H2 O)
with point group symmetry 222.
Note that the angle between the
(110) and (110) faces is 89.37◦C,
because the lattice parameters a,
(a) (b) b, and c are 1.1866, 1.1998 and
0.6855 nm respectively.
(a) (b)
FIGURE 1.22
Each figure shows an idealised model on the left and the corresponding crystal on the
right. (a) Quartz. (b) Celestite. Photographs taken at the Colorado School of Mines.
FIGURE 1.23
Octahedral interstices in austenite (cubic-F) and ferrite (cubic-I). The red objects
are located at a height half, with filled circles representing iron atoms and open
circles the positions of the octahedral interstices. (a) A regular-octahedral inter-
stice in austenite; (b) projection of the austenite unit cell to show the positions of
all such interstices. (a) An irregular-octahedral interstice in ferrite; (b) projection
of the ferrite unit cell to show the positions of all such interstices.
28 1 Introduction and Point Groups
1.11 Summary
It is extraordinary that all crystals without exception can be described in terms of
just fourteen imaginary lattices, the points of which can then be loaded identically
with motifs that may contain one or more atoms that are not necessarily identical.
On the same logic, the apparently great variety of wallpapers available for decora-
tive purposes have just five basic lattices, with the variations generated by changing
the motif. A sixth is of course possible if random two-dimensional patterns can be
printed. One advantage of a random-printed wallpaper is that there would be no effort
required in matching designs at the edges.
The idea that there is an identical configuration of atoms around each lattice point
in a crystal does not in fact sit well with the vast majority of real crystals. Materials
are not pure; random mixtures of atoms in solid solutions cannot result in an identical
environment at each lattice point. Translational periodicity does not then exist, rather,
there is a probability of finding a particular species at a particular location. The con-
sequence is that the crystal must be considered in terms of its averaged properties
and this suffices for most purposes, for example in X-ray diffraction, where an aver-
aged atomic scattering factor is used to calculate intensities. The strain fields around
solute atoms in a random solid solution may then lead to X-ray peak broadening.
The anisotropy of crystals, moderated by symmetry, has major consequences on
how they can be exploited for engineering purposes. The crystallographic direction
along which single-crystal turbine blades are grown is selected in order to minimise
the possibility of vibrations during the operation of gas turbines.
References
29
30 1 References
Abstract
2.1 Introduction
We have seen already that the projection of a three-dimensional crystal structure into
two dimensions, along the z-coordinate, can simplify the perception and representa-
tion of atomic arrangements. There is no loss of information given that the fractional
z-coordinate is identified clearly on the projection. This can be seen in Fig. 1.14
where the three dimensional representation of the structure of diamond lacks clarity
whereas the projected cell is readily visualised. The ease of interpretation using the
projected structure is illustrated for the more complex structure of ε-carbide, which
has a chemical formula between Fe2 C and Fe3 C, in Fig. 2.1. There are six iron atoms
in the cell (3 at z = 21 and the six others that are shared at faces with other cells at
z = 0, 1 and therefore contributing a further 3 to the cell). The illustrated cell con-
tains a full complement of carbon atoms but to achieve the composition Fe2.4 , the
sites coloured blue would only be occupied partially [1].
The projection described in Fig. 2.1 is a straightforward linear operation and there
are no distortions of that linearity. In contrast, Fig. 2.2 shows non-linear projections
of circular arcs, in one case a simple extrapolation onto a horizontal line, and in
the other case via a pole located below the horizontal line where the projections are
recorded. The purpose here is to record angles rather than spatial coordinates.
Imagine now that instead of quadrants of a circle, the diagrams in Fig. 2.2 rep-
resent the surfaces of spheres. A circle drawn on the surface of the sphere would in
general project as an ellipse on the horizontal plane of Fig. 2.2a, whereas it would
project as a true circle in the case of the projection method used in Fig. 2.2b.
31
32 2 Stereographic Projections
(a) (b)
FIGURE 2.1
The structure of hexagonal ε-carbide with lattice parameters a = 0.4767 nm and
c = 0.4353 nm. The large atoms are iron. The small atoms are carbon but not all
the sites designated blue are occupied. (a) The three-dimensional cell with the z-axis
vertical. (b) Projection of the cell along the z-axis, with the fractional z-coordinates
listed.
2.1 Introduction 33
90 80
70
60
50
90 80
70 40
60
50 30
40 20
10
20 0
10
0
90 70 50 30 0
(a) 10 (b)
FIGURE 2.2
Projections of angles made by lines intersecting a circular quadrant onto a horizontal
line.
(a) (b)
FIGURE 2.3
(a) Circle on the surface of a sphere projected on to the equatorial plane via the
south pole. (b) Construction showing that the projection is also a circle.
product phase will form such that the atomic arrangements match as much as
possible, along the parent/product interface. This often leads to a reproducible ori-
entation relationship between the two crystals, one which can be represented on a
stereographic projection in order to understand the mechanism of transformation
or the deformation behaviour of the two-phase mixture. Stereographic projections
can be used to decide whether the orientation relationship is reproducible or oc-
curs by chance [4].
• Diffraction data can be presented on stereographic projections. It is now routine to
to examine both structure and crystallographic information from polycrystalline
samples on a single image [5], with the crystallographic data presented in the
form of corresponding colours on the microstructural image and stereographic
projection.
small
great circle
circle
FIGURE 2.4
Intersections of planes with sphere. (a) Great circle inclined with respect to the north
and south poles of the sphere, (b) small circle; (c) latitude and longitude lines (image
courtesy of vectortemplates.com).
When an intersecting plane does not pass through the centre of the sphere it
results in a small circle at the intersection, as illustrated in Fig. 2.4b. Lines of latitude
are in general small circles, the exception being the equator.
The shortest distance between two points on a sphere is along the great circle
which passes through both of the points.
Fig. 2.5a shows a crystal with a cubic lattice, placed at the centre of a sphere
with the [100] ∥ x-axis, [010] ∥ y-axis and [001] ∥ z-axis of the sphere; this is the
36 2 Stereographic Projections
010 001
100
010
100
(a) (b)
(north) (north)
010 001 010 001
100 100
trace
010 x 010
x
100 100
FIGURE 2.5
Plotting plane normals on a stereographic projection. The insets, top-left in each
case, show a plan view of the equatorial plane. (a) Crystal placed at the centre of the
sphere. (b) Projection of the pole of (011). (c) Projection of (011). (d) Projection of
the trace of the great circle which is (011) on to the equatorial plane. The trace of the
plane in the southern hemisphere is marked as a dashed curve. The crosses identify
the intersections of the projection lines with the equatorial plane.
38 2 Stereographic Projections
trace
90°
pole
(a) (b)
(c) (d)
FIGURE 2.6
The Wulff net. (a) Relation between pole and trace of a plane. (b) The net is rotated
until the poles a and b lie on the same great circle to measure the angle in between.
(c) The geometrical centre c of the small circle is different from its angular centre b.
The angles ab and bd as measured on the great circle are identical. The distances ac
and cd are identical when measured on a ruler. (d) The poles identified with arrows
are both located at angles φ1 to p1 and φ2 to p2 .
2.3 Stereographic projection: construction and characteristics 39
r2 − ro2 = (r − x)2
so it follows that r = (x2 + ro2 )/2x with ro ≤ r ≤ ∞.
Fig. 2.7a,b show the stereographic projection for the cubic system, with the ro-
tational symmetry elements related to the crystallographic directions about which
the rotations occur. Fig. 2.7c shows the stereographic projection being used to illus-
trate how the elastic modulus of a single crystal of ferritic iron varies between 131
and 284 GPa as a function of the crystallographic orientation, in a manner consistent
with its crystalline symmetry. Even aggregates of crystals often do not show collec-
tive macroscopic-properties that are isotropic because the crystals are not distributed
at random. Polycrystalline ferritic steel sheet in which the crystals are aligned with
the ⟨001⟩ direction and with the {110} planes aligned to the rolling plane, exhibit
modulus variation in the range 140 to 210 GPa [6]. The level of anisotropy in aggre-
gates can be controlled, as we shall see in later chapters, to suit specific purposes such
as the formability of metal sheets or the magnetic properties of transformer steels.
40 2 Stereographic Projections
(a) (b)
(c)
FIGURE 2.7
Stereographic projections for the cubic system. (a) Showing the standard way of
drawing the projection with the locations of [100] and [010] on the perimeter and
[001] at the centre. (b) The symmetry elements. The filled squares are tetrads, the
ellipses are diads and triangles are triads. The tetrads are at ⟨001⟩, triads at ⟨111⟩ and
diads at ⟨011⟩. Anything with cubic symmetry must have four triads. (c) Variation
in the elastic modulus of a single crystal of bcc iron as a function of orientation.
Diagram courtesy of Shaumik Lenka, constructed using data from Dieter [7].
2.3 Stereographic projection: construction and characteristics 41
(a) (b)
FIGURE 2.8
(a) A meteorite at the Smithsonian Museum in Washington. (b) The Wid-
manstätten pattern observed in the meteorite, of ferrite plates in an austenitic ma-
trix.
d (112) a
f
-111
- -111
c
e
001
010
b
-
111 111
a d
b f
a 100 a
e
c a
(111) (001)
d (112)
FIGURE 2.9
Stereographic projection referred to the crystallographic axes of the austenite. The
blue lines are traces of {111}F planes, and the dashed lines the planes on which the
macrostructure is observed. The insets at the bottom show what the macrostruc-
tures should look like if all the plates are of identical dimension and shape, as a
function of the plane of observation.
(a)
FIGURE 2.10
(a-c) An arbitrary pole is rotated by 180◦ about the diad normal to the plane of the
diagram, to generate the second pole. This second pole is then inverted through
the centre into the southern hemisphere to to result in inversion diad 2 that is
equivalent to a mirror normal to that axis. (b) Triad, (c) inversion triad and (d)
inversion hexad. Mirror planes in (a) and (d) are identified using bold circles.
(a) (b)
(c) (d)
FIGURE 2.11
Sketch stereograms showing what happens to a general pole on the operation of
the following point groups (bold lines represent mirror planes): (a) trigonal 3m.
(b) The structure of GaN with point group symmetry 6mm, viewed with the 6-fold
axis along z. The gallium atoms are labelled green. (c) Cubic m3. (d) Projection
of the unit cell of pyrite (FeS2 ) with the fractional z coordinates indicated, glide
planes marked with dashed lines and screw diads. Not all the symmetry elements
are illustrated for simplicity. For example, there will be glide planes parallel to the
plane of the diagram at heights 14 z and 34 z.
The point group m3 belongs to the cubic crystal class with the first symbol m
parallel to the z axis (page 24) and the mandatory symbol ‘3’ to follow since the
defining symmetry of a cubic system is four triads, in this case inversion triads.
It is evident from the plotting of a general pole and its symmetry related poles
shown in Fig. 2.11c that there is a centre of symmetry. Pyrite belongs to the cubic
class with point group m3; the actual crystals are illustrated in Fig. 1.2 where
the four triads associated with the regular octahedral shape of the crystals are
obvious. Fig. 2.11d shows a projection of the crystal structure of pyrite, which
2.4 Stereographic representation of point groups 45
is primitive cubic. The structure can be compared with the stereogram showing
the point group symmetry in Fig. 2.11c. Bearing in mind that point groups do
not include translations, the glide planes and screw axes indicated in Fig. 2.11d
become mirrors and diads respectively in Fig. 2.11c.
Crystals which lack a centre of symmetry can be piezoelectric, i.e. they develop
a dipole on being deformed and will change shape on the application of an elec-
trical field. Under normal conditions the crystal has positive and negative electrical
charges which are symmetrically distributed, leaving it in a neutral state when homo-
geneously deformed. The application of a stress causes charge asymmetry, and the
development of a voltage across the crystal. PbTiO3 is such a substance (Fig. 2.12)
where for example, the titanium atom is not located at exactly 21 , 12 , 12 . The structure
is primitive tetragonal, with point group 4mm.1
The stereogram illustrated in Fig. 2.12b shows that the structure does not have a
centre of symmetry. Secondly, if an atom is placed at a general position (marked a)
then there will be seven other identical atoms in the unit cell, inconsistent with the
structure illustrated in Fig. 2.12a. An atom placed on a tetrad (marked c) will generate
just one atom per unit cell; thus, there is just one Ti at 21 , 12 , 0.572. The oxygens at
1 1 1 1
2 , 2 , 0 and 2 , 2 , 1 also lie on a tetrad, and since they are shared between two cells,
there is only one such oxygen in the unit cell. An atom, such as the oxygen at 0, 21 , 12 ,
placed on the mirror plane as illustrated by the poles marked b, would generate two
such atoms in the cell (the four on the vertical faces are each shared by two cells).
This example shows how the placing of an atom within a unit cell must be con-
sistent with the symmetry elements of the cell because there must be identical atoms
at all symmetry related positions.
There are many other consequences of whether or not a crystal has a centre of
symmetry. The free energy of a small volume element containing a one-dimensional
composition gradient varies as follows [9]:
% &2
dc d2 c dc
g = g{c} + κ1 + κ2 2 + κ3 (2.1)
dz dz dz
where c is the concentration, g{c} is the free energy of a homogeneous solution with
composition c, z is the distance and κi are specific functions of the dependence of
free energy on the gradients of concentration as described elsewhere [9]. In this, κ1
is zero for a centrosymmetric crystal since the free energy must then be invariant to
a change in the sign of the coordinate z.
1 That is, by the convention quoted on page 24 the z direction is assigned to the four-fold rotation axis
followed by m for the mirrors with normals parallel to the x and y axes, and the final symbol for the
tertiary mirrors with normals between the x and y axes.
46 2 Stereographic Projections
es indicate
mirror planes
(a) (b)
FIGURE 2.12
(a) Lead titanate (PbTiO3 ) is tetragonal with a = 0.3904 nm and c = 0.4152 nm.
The Pb is located at coordinates (0,0,0.116), the Ti at ( 12 , 21 , 0.572), an oxygen atom
at (0, 12 , 12 ) and another oxygen at ( 12 , 21 , 0). The structure is primitive tetragonal
with point group 4mm. Note that each of the lead atoms is on a cell edge (and hence
shared with three other adjacent cells), not a cell corner so there is only one lead
atom per cell. (b) Stereogram showing the symmetry operations associated with the
point group 4mm.
2.5 Summary
A circle on a sphere projects as a circle on the stereographic projection; angular
relationships are therefore preserved. The method is therefore well suited to look
at angles between planes or directions in a crystal. The symmetry of the crystal in
its three dimensions becomes easy to visualise by plotting the symmetry elements
in their correct angular orientations on a stereographic projection. The placing of
atoms inside the unit cell during attempts to solve structure must clearly be consistent
with the symmetry determined, for example, from the external form of well-formed
crystals.
The stereographic projection can also be used quantitatively, by superimposing
a Wulff net to directly measure angles between a pair of poles. The net consists of
projections of great and small circles, rather like the longitudes and latitudes on a
globe representing the earth, The shortest distance between two points on the surface
of a sphere is when the two points are located on an arc of a circle which has the same
diameter as the sphere. When using a Wulff net to measure the angle, it therefore is
necessary to locate the two poles on a great circle.
We have seen that some quite elementary considerations of symmetry can indi-
2.5 Summary 47
49
3
Stereograms for Low Symmetry Systems
Abstract
Most of the metals in the periodic table have a cubic structure under ambient condi-
tions, but there are many that are in the hexagonal class. The stereographic projection
for non-cubic systems depends on the lattice parameters and the angles between the
basis vectors. And it can no longer be assumed that directions and planes with the
same indices are parallel. We shall see that in the case of the hexagonal system, it
can be useful to use four basis vectors to define directions and planes.
3.1 Introduction
In the cubic system, a plane normal (pole) with indices hkl is parallel to a direction
uvw when h = u, k = v, l = w; thus the normal to (123) is parallel to [123]. This
is generally not the case for less symmetrical systems. Fig. 3.1 illustrates this for a
two–dimensional unit cell in which b > a.
normal to (11)
x
FIGURE 3.1
The direction [11] is not parallel to the normal to the
tra plane (11). On the other hand, the directions parallel
ce
] of ( to the cell edges, i.e. [01] and [10] are parallel to the
[11 11)
y
normals of (01) and (10) respectively.
The principles for the plotting of poles on stereographic projections are never-
theless identical for any lattice. However, the cubic system exhibits extraordinary
symmetry so angular relationships on the stereographic projection are not dependent
on the lattice parameter. This is not the case for other lattices with lesser symmetry,
the stereograms of which will not only appear different from a distance, but for the
same lattice type, the angular positions become dependent on the lattice parameters.
Fig. 3.2a,b show how the projections differ in the case of the cubic and orthorhombic
lattices where the absence of tetrads in the latter is obvious. Fig. 3.2c shows how the
projection changes when one of the lattice parameters of the orthorhombic lattice is
changed relative to that illustrated in Fig. 3.2b.
51
52 3 Stereograms for Low Symmetry Systems
The fact that directions are not in general parallel to plane normals with the same
indices is apparent from a comparison of Fig. 3.2c,d.
(a) (b)
(c) (d)
FIGURE 3.2
Stereographic projections centred on the 001 pole. (a) Cubic, lattice parameter 5
units. (b) Orthorhombic, lattice parameters 2, 3, 8 units. (c) Orthorhombic, lattice
parameters 2, 3, 5 units. (d) The orthorhombic cell as in (c) but with directions plotted
instead of poles.
3.2 Hexagonal system 53
(a) (b)
(c) (d)
FIGURE 3.3
(a) Showing the hexagonal cell with the normals to the (100) and (010) planes iden-
tified with blue-coloured arrows. (b) Stereographic projection of hexagonal structure
(lattice parameters a = 2, b = 2 and c = 3 units with γ = 120◦ ) centred on the 001
pole. (c) Showing the hexagonal cell, now in four index notation. (d) Corresponding
stereogram in four index notation.
It is evident from Fig. 3.3a that the (100), (010) and (110) are planes belong-
ing to the same form {100}, i.e., they are crystallographically equivalent. And yet,
54 3 Stereograms for Low Symmetry Systems
they have different indices. To eliminate this discrepancy, a four index system is
frequently used for hexagonal lattices by introducing a fourth axis (labelled t in
Fig. 3.3c), which gives plane indices as hkil, where from geometry, i = −(h + k).
This has the effect of giving equivalent planes the same combination of digits in
their indices. Thus, (100), (010) and (110) become (1010), (0110) and (1100) re-
spectively. The four components are referred to as the Miller-Bravais indices and the
system is reserved for hexagonal lattices [1].
A similar conversion for directions from uvw to U V JW is a little more difficult
because it is designed in such a way that J = −(U + V ) and at the same time, the
Weiss zone law in which a plane hkl containing a direction uvw is satisfied
hu + kv + lw = 0 (3.1)
so that
hU + kV + iJ + lW should equal zero.
On substituting for i = −(h + k) we get,
U −J =u u = 2U + V U = 13 (2u − v)
V − J = v → v = U + 2V and V = 13 (2v − u) (3.3)
W =w w=W J = − 13 (u + v)
TABLE 3.1
Directions of the same form now have similar combinations of digits in the four
figure notation.
The normal method of taking scalar products between vectors in order to de-
termine angles is difficult in non-cubic systems without a deeper understanding of
real and reciprocal vectors, as explained in Chapter 9. Noting that angles between
planes are also the angles between their normals, the simple geometry illustrated
in Fig. 3.4a,b shows that the required angle is arctan(c/a cos 30) = 61.36◦. Note
that the [0001] direction is normal to (0001) and parallel to the (1010) plane, so
the normals to these planes make angles of 0◦ and 90◦ respectively to [0001].
The slip directions and slip planes are plotted on Fig. 3.4c. The slip direction
lying in a particular slip plane can be identified by applying the Weiss zone rule.
Thus, [2110] lies in each of the planes (0001), (0111) and (0110).
(1011)
[0001]
(1011)
c normal to
(1011)
61.36°
p o p
(a) (b) a cos 30° (c)
FIGURE 3.4
(a) Hexagonal unit cell showing the (1011) plane. (b) The section poq which
contains the normals to the (0001) and (1011) planes. (c) Stereographic projec-
tion showing the poles of the slip planes (filled circles), and the slip directions
(filled squares). The plot is based on the lattice parameters a = 0.2951 nm and
c = 0.4679 nm.
where ϵij with i = j are normal strain and the rest shear strains such that ϵij = ϵji .
Assuming that there is no volume change, ϵ11 + ϵ22 + ϵ33 = 0 so only two of the
normal strains are independent and there are three shear strains. Therefore there
are five independent strains needed in order to produce an arbitrary deformation.
Consider now a polycrystalline material. The application of stress will cause the
individual crystals to respond differently because they have differing orientations
relative to the stress. Their changes in shape are therefore arbitrary as far as neigh-
bouring crystals are concerned, so these differences must be accommodated, with-
out which fracture occurs. Hence the need for five independent slip systems for a
polycrystalline material to maintain structural integrity.
Iron, by far the most used metallic material in the world, under ambient condi-
tions is cubic-I with one iron atom located at each lattice point. The common slip
system is {011}⟨111⟩, each containing two independent slip directions of the form
⟨111⟩. There are therefore, twelve independent slip systems, making the iron ductile
in polycrystalline form. The same can be demonstrated for the cubic-F form of iron.
This makes iron very ductile.
The cubic-I form of iron is stabilised under ambient conditions by its ferromag-
netic properties; ruthenium and osmium, which in terms of their outer-electron struc-
ture are iron analogues in the periodic table, do not exhibit ferromagnetism and
hence have hexagonal crystal structures. In the absence of its ferromagnetic prop-
erties, the stable form of iron would also be hexagonal with a motif of a pair of
identical Fe-atoms, one at 0,0,0 and another at 31 , 23 , 12 . The easy-slip system is then
{0001}⟨1120⟩. Since there is only one (0001) plane in which there are three direc-
tions of the form ⟨1120⟩, polycrystalline hexagonal iron is relatively brittle since the
number of slip systems in each crystal is less than the required five. So in the absence
of ferromagnetism, we would not have civilisation in the form that we know today.
There are caveats to this story. Deformation can occur by mechanisms other than
slip, for example, mechanical twinning, so that the five independent slip system re-
quirement is not sufficient to determine the ductility of a polycrystalline material.
And iron below its Curie temperature is not strictly cubic, but tetragonal because
the magnetic spins in any given crystal tend to be aligned along a ⟨100⟩ axis [2].
However, the tetragonality is quite small and can be neglected in most experiments,
although it does manifest on a macroscopic scale via the magnetostrictive effect,
which is the variation in the length of a sample subjected to a magnetic field [3].
3.2 Hexagonal system 57
Vector scalar products can then be implemented using elementary methods as fol-
lows to calculate the required angles (φ in the case illustrated):
[0 0 c]Fe3C ! [0 b c]orthonormal c2
cos φ = - =- = 37◦
c2 × (b2 + c2 ) c2 × (b2 + c2 )
a c
n
[001] [013]
[010]
[001]
FIGURE 3.5
Discovery of the growth direction of cementite laths in austenitic iron using trace
analysis. The stereogram shows directions in cementite that has the lattice param-
eters a = 0.45165 nm, b = 0.50837 nm and c = 0.67475 nm.
3.3 Summary
Complications arise once the defining symmetry of the cube, i.e., four triads, is lost.
Directions and plane normals with the same indices are then not necessarily parallel
and hence have to be distinguished on a stereographic projection. Furthermore, the
extent of the discrepancy relative to the cubic system depends on the anisotropy of the
lattice parameters and the angles between the basis vectors of the system concerned.
The hexagonal system poses a particular problem of communication because of
the nature of the special symmetry of the hexad. For example, the crystallographi-
cally equivalent planes (100), (010) and (110) have indices that are not permutations
of the same set of digits. The Miller-Bravais four-index system eliminates this am-
biguity but requires additional work to ensure that the Weiss zone rule is satisfied.
Whether the intellectual cost of the four-index system over the usual three-index no-
tation is justified by its elegance is a matter of taste.
References
3. J. P. Joule: ‘On the effects of magnetism, &c. upon the dimensions of iron and
steel bars.’, The London, Edinburgh, and Dublin Philosophical Magazine and
Journal of Science, 1847, 30, 225–241.
4. M. V. Kral, and R. W. Fonda: ‘The primary growth direction of Widmanstätten
cementite laths’, Scripta Materialia, 2000, 43, 193–198.
59
4
Space Groups
Abstract
Some of the symmetry elements of a crystal involve translations that are simple frac-
tions of a repeat distance. These translations have consequences, for example on the
number of symmetry related positions in the unit cell, and hence on the number of
atoms in the unit cell. This chapter deals with space groups, which, unlike point
groups, include screw axes and glide planes.
4.1 Introduction
A point group is a collection of symmetry elements passing through a point, and
therefore, necessarily does not include translations. Space groups, in contrast, in-
clude translations that are fractions of a repeat unit, for example, a 21 axis which
involves a rotation of 180◦ followed by a translation of 21 of the repeat distance along
the axis. Glide planes involve reflections followed by fractional translations. These
translations are small, and hence do not manifest when the point group is determined
from the external shape of well-formed crystals. However, they have consequences
in structure determination and on other properties of crystals. There are 230 space
groups which are made from combinations of the 32 crystallographic point groups
with the 14 Bravais lattices.
The relationships between the placing of atoms in the unit cell and the symmetry
of the structure of lead titanate was described on page 45. The structure of PbTiO3
is primitive tetragonal with point group 4mm. Since there are no screw axes or glide
planes, the space group symbol for the titanate is P 4mm, where the first symbol is
to identify the lattice type as primitive. In this chapter, we will consider a case where
translational symmetry elements1 are present.
1 other than the normal periodicity of a lattice
61
62 4 Space Groups
FIGURE 4.1
An illustration of the operation of a diad
where the structure is restored by a rota-
tion of 180◦ , and a screw diad which in-
volves the same rotation but then a trans-
lation along the axis by 0.5 (hence the no-
tation 21 ) of the repeat distance t.
Glide planes are like mirror planes but a reflection is followed by a translation
in order to recover the symmetry.The terminology for glide planes is explained in
Table 4.1; the translations following reflection are parallel to the glide plane, the
exception begin the diamond glide in a body-centred cubic structure.
TABLE 4.1
Symbols for glide planes. ‘fcc’ and ‘bcc’ are abbreviations for face–centred cubic
and body–centred cubic respectively. The structure of pyrite, which is primitive cubic
and has a point group symmetry m3 is illustrated on page 44, containing axial glide
planes; its space group therefore is P a3 where the symbol a replaces the m in the
point group.
4.3 Cuprite
Consider now the structure of cuprite (Cu2 O) which is primitive cubic, with a motif
of four copper atoms at 000, 21 0 12 , 0 12 21 and 12 21 0, and two oxygen atoms at 14 14 43
and 34 34 41 located at each lattice point.
FIGURE 4.2
Projection along the [001] direction of four unit cells of cuprite (Cu2 O). The frac-
tional height along [001] are indicated, with unlabelled atoms at 0 and 1.
Fig. 4.2 shows a projection of four such unit cells along the [001] direction, with
fractional coordinates along that direction indicated. The environment of each of the
six atoms within the unit cell is clearly different, consistent with the cubic-P lattice
type. There are mirror planes parallel to {110} – the trace of only one of these is
illustrated in Fig. 4.2 for the sake of clarity. Note also that there are 42 axes located
at 14 41 z and 34 43 z and their equivalent positions. A 42 axis involves a rotation of 90◦
followed by a translation by 24 of the repeat distance along that axis. The diagram
also shows that the structure has a centre of symmetry at 12 12 12 .
64 4 Space Groups
The glide plane parallel to {100} (Fig. 4.2) involves a reflection followed by two
orthogonal translations parallel to the plane by 12 the lattice parameter. It is, therefore,
an n-glide plane according to the notation in Table 4.1.
Since cuprite is cubic, the defining symmetry is four triads. However, Fig. 4.3
shows that the triads are in fact inversion triads (3), meaning that a rotation of 120◦
is followed by an inversion through the centre to recover the structure. This is evident
from the projection along ⟨111⟩ that there is no triad since the atoms marked with
z coordinates are in two sets of three, with each set at a different height along z. A
simple three-fold operation would not therefore work. On the other hand, the rotation
combined with the inversion reproduces the structure.
(a) (b)
FIGURE 4.3
(a) The unit cell of cuprite. (b) A projection along ⟨111⟩, with fractional coordinates
along [001] indicated for some of the atoms.
(a) (b)
(c) (d)
FIGURE 4.4
Stereograms consistent with the point group m3m, showing the number of equiva-
lent positions when (a) an atom is placed at general location x, y, z and all equivalent
positions; (b) an atom is placed at general location x, x, z and all equivalent posi-
tions. (c) Projection of the cuprite unit cell with the traces of glide planes marked
as dashed lines and those of mirror planes as heavy lines. Initial placing of atoms
copper in black and oxygen in red. (d) New set of atoms generated by the operation
of the glide planes. The relation between the original and new atoms is marked by
blue lines.
66 4 Space Groups
the symmetry elements 3m passing through it. The operation of the n-glide plane
marked by the horizontal line generates the atom at the opposite face-centre at a
height 12 and the other symmetry elements then generate atoms at all face-centres
and corners of the cube. There are therefore, four equivalent copper atoms in the unit
cell.
Suppose that an oxygen atom is placed at the location 41 , 41 , 34 which has the sym-
metry elements 43m passing through it. The operation of the same n-glide plane on
that oxygen atom leads to an oxygen atom at 34 , 34 , 14 . All other symmetry operations
have the same effect so there are only two oxygen atoms per cell. The locations of the
copper and oxygen atoms are therefore consistent with the chemical formula Cu2 O
and the density of the compound can be measured to verify this. Placing a copper
atom at the general position x, y, z or x, x, z would not lead to the correct density
or stoichiometry. The results are summarise in Table 4.2. The International Union
of Crystallography publishes much more comprehensive space group tables which
list the multiplicity of coordinates as a function of symmetry at that coordinate, for
all possible scenarios and for all space groups. Table 4.2 can be regarded as a much
simplified version of this kind of information.
TABLE 4.2
Space group for cuprite, P n3m. The origin of the cell is set on a copper atom.
1 48 x, y, z x, y, z x, y, z ...
m 24 x, x, z x, x, z x, x, z ...
..
.
..
.
1 1 1 1
3m 4 0,0,0 2, 2, 0 2 , 0, 2 0, 12 , 12
1 1 3 3 3 1
43m 2 4, 4, 4 4, 4, 4
Space groups!P m
FIGURE 4.5
Projection of the unit cell of Fe2 USi normal to the z-axis.
68 4 Space Groups
m
b
a
x
FIGURE 4.6
The crystal structure of cementite, consisting of twelve iron atoms (large) and four
carbon atoms (small). The fractional z coordinates of the atoms are marked.
The lattice type is primitive (P). There are n-glide planes normal to the x-axis,
at 41 x and 34 x involving translations of 2b + 2c . There are mirror planes normal to the
y-axis and a-glide planes normal to the z-axis, at heights 14 z and 34 z with fractional
translations of a2 parallel to the x-axis. The space group symbol is therefore P nma
[1].
The structure of cementite clearly is very anisotropic so its properties, such as
the elastic moduli, also vary strongly with the direction within the crystal [2]. The
shear modulus c44 is exceptionally small, some two times smaller than the cor-
responding term for aluminium. Nevertheless, the cementite has an exceptionally
large ideal shear strength because elastic deformation reduces its symmetry from
orthorhombic to monoclinic (space group P 21 /c), with an accompanying increase
in three dimensional covalent bonding that stiffens the material [3].
4.4 Location of atoms in cuprite cell 69
a
In the context of plastic deformation, the slip systems on which dislocation
glide has been observed include (001)[100], (100)[010], (100)[001], (010)[001]
and (010)[100] [4] but the complex arrangement of atoms means that slip is ac-
tually quite difficult to achieve, making the cementite very hard. There may be
other slip systems that operate when the cementite is forced to deform in a phase
mixture such as pearlite, but the Burgers vectors of slip dislocations would never-
theless be large when compared with dislocations in the allotropic forms of iron
[5].
a It is important to note the labelling of the orthogonal axes of the unit cell, i.e., a = 0.50837 nm,
b = 0.67475 nm and c = 0.45165 nm. since different conventions are used in the published literature.
The indices referred to below have been adjusted from the original sources, for the unit cell as defined
here.
(a) (b)
FIGURE 4.7
Structure projections. (a) Diamond (carbon). The dashed line shows a d-glide
plane involving two orthogonal translations by a quarter of the repeat distance,
parallel to the glide plane. (c) Cubic form of zinc sulphide.
[001]Ω ∥ [111]Al
[100]Ω ∥ [101]Al
Fig. 4.8 shows the stereograms of the symmetry elements for the two crystals in the
4.6 Summary 71
correct relative orientation. It is evident the common point group symmetry is 2/m,
which is consistent with the observed plate shape of the Ω both with respect to the
diad and the mirror plane.
FIGURE 4.8
(a) Stereographic representation of the symmetry elements in the point group m3m
for aluminium. (b) Stereographic representation of the symmetry elements in the
point group mmm for Ω. The symmetry elements shared by Al and Ω are highlighted
in red. (c) The precipitate shape observed, corresponding to the shared symmetry
group 2/m.
4.6 Summary
We have touched on the notion of space groups which represent all the symmetry
elements of the crystal. The appreciation of space groups is essential in the solution
of structures, transformations which lead to changes in symmetry, and the properties
of crystals. A great deal of the work in identifying the number of symmetry related
positions as a function of the space group has already been documented, with the
information routinely available either from the tables published by the International
Union of Crystallography or via the world wide web.
References
73
5
The Reciprocal Lattice and Diffraction
Abstract
The Miller indices for directions and planes are defined differently. For directions
the indices are simply the components of a vector corresponding to that direction
with respect to the basis vectors, but those of a plane are defined rather strangely at
first sight. The intercepts of the plane with the basis vectors are determined, and their
reciprocals form the indices of the plane. The mystery behind this scheme is revealed
in this chapter, whereby a plane is also defined by a vector that is normal to it and
has a magnitude that is the reciprocal of the plane spacing. Diffraction too becomes
easier to visualise when expressed in terms of reciprocal space.
u = u1 a1 + u2 a2 + u3 a3 ,
The term (a1 .a2 ∧ a3 ) represents the volume of the unit cell formed by ai , while the
magnitude of the vector (a2 ∧ a3 ) represents the area of the (1 0 0)A plane (Fig. 5.1).
Since (a2 ∧ a3 ) points along the normal to the (1 0 0)A plane, it follows that a∗1 also
points along the normal to (1 0 0)A and that its magnitude |a∗1 | is the reciprocal of
the spacing of the (1 0 0)A planes (Fig. 5.1).
75
76 5 The Reciprocal Lattice and Diffraction
FIGURE 5.1
The relationship between a∗1 and ai . The vec-
tor a∗1 lies normal to the basal plane wxyz and
the volume of the parallelepiped formed by the
basis vectors ai is given by a1 .a2 ∧ a3 . The
area of the basal plane is |a2 ∧ a3 |. The height
is the spacing between planes parallel to the
basal plane, given by 1/|a∗1 |.
The components of any vector referred to the reciprocal basis represent the
Miller indices of a plane whose normal is along that vector, with the spacing of the
plane given by the inverse of the magnitude of that vector. For example, the vector
(u; A∗ ) = (1 2 3) is normal to planes with Miller indices (1 2 3) and interplanar
spacing 1/|u|. We see that
FIGURE 5.2
Sections of the reciprocal lattice of a cubic-P lattice. (a) Section normal to [001].
(b) Section normal to [011]. (c) Section normal to [111]. Notice how the distance
varies inversely with the spacing of the planes. The origin of the reciprocal lattice
is not labelled, but if necessary, can be designated 000. Each of the spots in these
patterns should strictly be points, but they are weighted here with intensity to show
the relationship with electron diffraction patterns discussed later.
a∧b
c. (5.4)
|a ∧ b|
Every vector in a reciprocal lattice defines a normal to a plane of the same indices
as the components of the reciprocal lattice vector, and the vector has a magnitude
which is the reciprocal of the spacing of those planes. Therefore, for each of the
basis vectors,
a∧b
c∗ = etc. (5.6)
c.a ∧ b
With this equation it is evident that ai .a∗j = 0 when i ̸
= j and ai .a∗ j = 1 when
∗ ◦
i = j. The angle between a and b is clearly 90 . Therefore, if a direction [u v w]
lies in a plane (h k l),
Note that this is a generic proof for the Weiss zone rule since there are no restric-
tions placed on the basis vectors.
A corollary to this example is a general method for taking a scalar product be-
tween two vectors. The magnitude of a vector u is given by
where (u; A) is a row matrix of the components ui of u referred to the basis vectors
ai represented by the basis symbol ‘A’, and [A∗ ; u] is a column matrix containing
the components u∗i of u referred to the reciprocal basis vectors a∗i represented by the
basis symbol ‘A∗ ’. This is an important results that gives a different interpretation
to the scalar (or ‘dot’) product between any two vectors referred to any basis. To
evaluate the dot product, one of the vectors must be referred to the real basis and the
other to the corresponding reciprocal basis.
78 5 The Reciprocal Lattice and Diffraction
FIGURE 5.3
Electromagnetic waves inci-
dent on a set of parallel crys-
tal planes with an interplanar
spacing d. The angle of emer-
gence of the scattered waves
is the same as that of inci-
dence.
k′ − k = g (5.10)
Here g is a reciprocal lattice vector beginning at the origin (0,0,0) and representing
the crystal which is being illuminated by the radiation k. From the geometry of the
triangle, it is evident that sin θ = 12 g/k so that sin θ = (0.5/d)/(1/λ) which on
rearrangement gives the Bragg law.
5.3 Intensities 79
(a) (b)
FIGURE 5.4
(a) The Ewald sphere construction in reciprocal space. k represents the direction of
the incident beam, k′ that of the diffracted beam and g the normal to the diffracting
planes. Both k and k′ have a magnitude λ−1 and g has a magnitude d−1 . The beam is
incident on a crystal which is represented by its reciprocal lattice with origin at 0,0,0.
(b) Another illustration of the Ewald sphere superimposed on the reciprocal lattice.
In this case only the 101 and 101 planes are in Bragg condition. The reciprocal lattice
vectors which do not touch the sphere are not in Bragg orientation.
5.3 Intensities
The electrons from an individual atom can coherently scatter an X-ray beam; the
atomic scattering factor f is the ratio of the amplitude scattered by an atom to that
scattered by one electron. We are interested in coherent scattering by all the atoms
in a unit cell, in which case the resultant amplitude for reflections from hkl planes is
given by a summation known as the structure factor:
n atoms
.
Fhkl = fn exp{2πi(hun + kvn + lwn )} (5.11)
1
n atoms
.
≡ fn [cos 2π(hun + kvn + lwn ) + i sin 2π(hun + kvn + lwn )]
1
where fn is the scattering factor of atom n and the sum is over all atoms in the
unit cell.† The magnitude |F | is now the ratio of the amplitude scattered by a unit
cell to that scattered by one electron. |F |2 is proportional to the scattered intensity.
Notice that for a centrosymmetric system, the sine term can be neglected since that
function is not symmetric about zero. As an example, for the structure of Cu (cubic
† eπi = e3πi = −1, e2πi = e4πi = +1, eπi/2 = i, e3πi/2 = −i, (1 + i)(1 − i) = 2
80 5 The Reciprocal Lattice and Diffraction
11 1 1
close-packed, with four atoms in the cell at 000, 2 2 0, 2 0 2 , and 0 21 12 ),
so that the {100} and {110} reflections would have zero structure factor.
Fig. 5.5 shows how the electron diffraction patterns of austenite, which has a
cubic-F lattice, have zero intensities for {100} and {110} planes, as expected by
substitution of hkl into equation 5.12.
FIGURE 5.5
Electron diffraction patterns from austenite (cubic-F). (a) Zone axis parallel to [001].
(b) Zone axis parallel to [011]. (c) Zone axis parallel to [111]. Notice that in any
pattern, it is necessary to index just two reciprocal lattice vectors in order for all
the others to be determined. For the cases illustrated, the symmetry of the patterns
is consistent with that of the zone axes. Contrast these patterns with the reciprocal
lattice sections presented in Fig. 5.2.
The pattern is rectangular so one possibility is that the axis normal to the
pattern corresponds to a ⟨110⟩, a diad. The g-vectors normal to this direction could
include the normals to {001}, {110}, {112} and {11n} where n is an integer. The
d-spacing of planes in a cubic system is given by
1 h2 + k 2 + l 2
=
d2 a2
where a is the lattice parameter. So the d-spacing decreases in the order {001},
{110}, {112} and {11n}. Fig. 5.6b shows that the intensity diffracted from the
{001} planes would be zero given that there is an atom at the body-centre which
would scatter out of phase with the atoms on the cube faces. However, {002}
planes would scatter in phase. If it is assumed that the g-vector labelled ‘a’ is
a {110} reflection (largest d-spacing and hence closest to the incident beam in
reciprocal space) and that at the reflection at ‘b’ is due to the {002} planes then
5.3 Intensities 81
ratios of the lengths of the ‘a’ and ‘b’ and the 90◦ angle between gives consistent
indexing. Any other choice would not.
Once two of the reflections are labelled, all others are linear combinations of
those two. For example ‘c’ is given by 002 + 110 = 112.
In this particular example, it was easy by examining the structure projection
to show that {001} reflections would be systematically absent. However, if the
indices are of higher order, e.g. {532}, then it is easier to use the structure factor
equation to decide whether these planes would result in diffracted intensity.
(a) (b)
FIGURE 5.6
(a) An electron diffraction pattern from a body-centred cubic crystal. (b) Structure
projection of bcc, with traces of {110}and {200} planes indicated.
The pattern has three-fold symmetry because there are three different g-
vectors if their opposites are not counted. This indicates that there is a triad ⟨111⟩
normal to the pattern. All the g-vectors will be normal to the ⟨111⟩ axis and of
equal length, so the following solutions can be deduced:
• For cubic-P, ‘a’= {011} and ‘b’= {110}. The choice of the first solution as
{110} is arbitrarily made from the three possible, but the second will then be
determined as being at 60◦ from the first. Their right-handed cross product
gives [111] so the beam direction is [111].
• For cubic-I, the solution would be the same as for cubic-P since {011} type
reflections are permitted.
• In the case of cubic-F, {011} is a systematic absence but {022} is not so all
the indices should be doubled.
82 5 The Reciprocal Lattice and Diffraction
(a) (b)
FIGURE 5.7
(a) An electron diffraction pattern from a crystal with a cubic structure. (b) Struc-
ture projection for cubic-F showing that {011} would result in zero diffracted
intensity because of the atoms in-between.
Notice that it is not possible, from the shape of the diffraction pattern alone,
to determine the lattice type. However, if the lattice parameters of the different
forms are known, then the measured d-spacings on the pattern can be compared
with interplanar spacings calculated from the parameters to fix the lattice type.
The relationship between the distance measured on the diffraction pattern and the
camera length of the electron microscope is illustrated in Fig. 5.8. It is evident that
sin θ = R/2L so with the Bragg equation 5.9, it is easy to show that Rd = Lλ,
with the term Lλ referred to as the ‘camera constant’.
FIGURE 5.8
Geometry describing the relationship between the camera length on the transmis-
sion electron microscope and distance measured on an electron diffraction pattern.
(a) (b)
FIGURE 5.9
(a) Al-Fe solid solution (γ). (b) δ-phase. The relative orientations of the two pat-
terns are preserved.
The γ is said to be cubic-F; since the pattern has four-fold symmetry, the direc-
tion normal to the pattern must be parallel to ⟨001⟩ so the spots marked ‘a’ and ‘b’
can be labelled {020} and {200} respectively, with all other reflections obtained
by linear combinations of these two. Therefore, ‘c’ becomes {220}. The shortest
lattice vector in cubic-F is a2 ⟨110⟩. Note that in the cubic-F lattice, reflections of
the type {100} are systematic absences (equation 5.12).
The pattern for the δ-phase also has four-fold symmetry and hence could come
from a cubic or tetragonal lattice. However, it is stated that the lattice has four
triads so it must be cubic-P, cubic-F or cubic-I. The composition and number of
atoms in the unit cell is consistent with cubic-P with a motif of four atoms per
lattice point, iron at 000 and three aluminium atoms at 12 21 0, 0 12 12 and 21 0 12 . As a
result, the spots labelled ‘a’, ‘b’, ‘c’ and ‘d’ in Fig. 5.9b are identified as {010},
{100},{110},{020} respectively.
The structure projections of the two phases are shown in Fig. 5.10, together
84 5 The Reciprocal Lattice and Diffraction
with the shortest lattice vectors along the ⟨110⟩ direction. The significance of this
is that a dislocation with Burgers vector a2 ⟨110⟩ in γ would create faults in the δ
because it would not be a lattice vector in the ordered phase. However, a pair of
γ dislocations moving together as a superdislocation in the δ would achieve slip
without leaving a fault in their wake. The γ and δ are in a ‘cube-cube’ orientation,
i.e.
[100]γ ∥ [100]δ and [010]γ ∥ [010]δ
and have similar lattice parameters, which would facilitate slip that initiates in the
γ to traverse into δ with the caveat that lattice dislocations of γ would need to pair
up when moving into δ.
(a) (b)
FIGURE 5.10
(a) Structure projection of γ. (b) Structure projection of δ.
where f is an average atomic scattering factor for the random solid solution γ.
(a) (b)
FIGURE 5.11
(a) Ewald sphere construction for diffraction from a thin foil. The origin of the re-
ciprocal lattice is at ‘O’ and the layer designated ‘P’ is the next layer normal to the
incident beam, a distance 1/t∗ from the first layer. The arrows mark reciprocal lattice
points which contribute to diffraction even though they are not at the exact Bragg ori-
entation. (b) ⟨111⟩ axis ReO3 showing the higher order Laue zone (the second layer
of the type sketched in (a).) - pattern courtesy of A. Eggeman, T. Chang and P. Midg-
ley.
86 5 The Reciprocal Lattice and Diffraction
FIGURE 5.12
Fe with cubic-F structure. (a) Electron diffraction from a foil which is 200 nm thick,
beam direction [111]. (b) Diffraction from the same foil after thinning it to just 50 nm
thickness, beam direction [111]. Because of this minute thickness, it becomes pos-
sible to pick up reflections from planes of the form {224}. (c) as (a) but with the
beam tilted towards [1 1 0.9]. The Ewald sphere is now tilted away from some of the
reciprocal lattice points on the left hand side.
FIGURE 5.13
The effect of crystal size (in all directions) on an X-ray diffraction pattern. The ma-
terial concerned is halite. The peaks broaden as the crystal gets smaller.
because they now have different neutron scattering factors (for example, hydrogen
and deuterium).
X-rays typically penetrate a few micrometers of the surface of a sample, whereas
neutrons can penetrate several centimetres. Neutrons can therefore be used to deter-
mine the state of residual stress that exists in a macroscopic sample that is otherwise
at equilibrium [2, 3]. Fig 5.14 shows residual stress contours that are determined by
measuring lattice spacings, comparing them with a stress-free sample to calculate
strains, and then converting strains into stresses using the elastic properties of the
material.
FIGURE 5.14
Stress distribution in a steel
weld determined using neu-
tron diffraction [4]. This
would not be possible using
X-rays. The numbers on the
contours represent MPa.
5.6 Summary
It takes a leap of imagination to understand why it is necessary to deal with a space
in which dimensions are inverted. But the fact is that in crystallography, it can be
regarded simply as a lattice in which all reciprocal lattice vectors represent normals
to planes in real space, with the magnitudes of these vectors representing the inverse
of the spacing of those planes. Whereas a plane in real space is defined by a pair of
vectors that lie in that plane, in reciprocal space the plane is represented by a vector.
A plane in reciprocal space represents all real space planes that share a common
direction (a zone axis).
Imagine now that we wish to find a direction that is parallel to a plane normal.
This is trivial for a cubic system where plane normals and directions with the same
indices are parallel. So a [123] direction is parallel to the normal of the (123) plane.
This is not generally true for non-cubic systems, but there exists a metric tensor
which makes it easy to find vectors in reciprocal space that are parallel to those in
real space (and vice versa), a subject reserved for Chapter 9.
References
89
6
Deformation and Texture
Abstract
Single crystals are now used routinely in engineering turbine blades for service in
hottest part of an aeroengine; it is astonishing that the design of the blade permits it
to operate in an environment where the temperature is greater than its melting point.
The deformation of a single crystal on the application of a stress or a system of
stresses is interesting in its own right, but the understanding that emerges from single
crystal studies can useful in explaining the rotation of crystals in a polycrystalline
material during plastic deformation. Real polycrystalline materials rarely consist of
crystals that are randomly oriented relative to the sample frame of reference. Instead,
the long-range order in a polycrystalline material is somewhere between that of a
hypothetical aggregate consisting of randomly oriented crystals, and a single crystal.
This is because the crystallographic axes of the different grains in the polycrystal tend
to align in specific ways. We shall learn in this chapter the methods for representing
such texture.
91
92 6 Deformation and Texture
factors that are smaller. The slip system with the highest resolved shear stress when
the tensile lies in the shaded triangle will therefore be [011](111). The construction
illustrated is known as Diehl’s rule.
(a) (b)
FIGURE 6.1
(a) Deformation of a single crystal. The force is resolved onto the slip plane and
towards the slip direction. Dividing by the area of the slip plane gives the shear stress
along the slip direction. (b) There are 24 stereographic triangles corresponding to the
24 slip systems of the type {111}⟨110⟩. The orientation of a tensile axis is plotted in
the shaded triangle.
FIGURE 6.2
Slip of an unconstrained single-
crystal. Notice how the axis along
which the force is applied rotates as a
consequence of the slip. The amount
of shear implemented is written αs. If
the force axis is constrained to be ver-
tical then the sample will rotate rela-
tive to the testing machine by the ex-
tent illustrated.
94 6 Deformation and Texture
s
d
FIGURE 6.3
Schematic diagram showing the shear of a single crystal. n is a unit normal vector
to the plane on which the shear occurs, d is vector describing the displacement
of the far end of the crystal due to the shear, x is its initial length above the slip
plane and y the corresponding length following the deformation. s is a unit vector
parallel to the shear direction.
A more general method for calculating the elongation during single-crystal deforma-
tion is described on page 175.
6.2 Slip in a single-crystal 95
FIGURE 6.4
Stereographic projection of austenite showing the planes and shear directions cor-
responding to the deformation that leads to the formation of ε martensite.
The shear strain is given by the displacement divided by the height. The height
2a
is the d-spacing of two close packed planes of austenite, i.e., √ 3
. The displacement
1
is the magnitude of a6 ⟨112⟩, i.e. √a6 . Therefore, the shear strain is 2√ 2
.
The habit plane is {111} (four variants) and there are three directions of the
form ⟨112⟩ in each of these planes. Therefore, twelve variants in all.
The strain is a shear deformation on a close-packed plane, but with the shear
direction along ⟨112⟩. Therefore, either using Diehl’s rule, or by calculation, the
plane stimulated would be (111). In that plane, there are three shear directions,
[211], [112] and [211]. Of these, the [211] direction gives the largest Schmid factor.
96 6 Deformation and Texture
6.2 Texture
Imagine that there is a polycrystalline substance in which the individual crystals are
randomly oriented relative to an external frame of reference‡ . Suppose that this poly-
crystalline substance is subjected to rolling deformation. The individual crystals will
then tend to rotate such that the slip planes comply with the external forces in a
manner akin to that illustrated in Fig. 6.2b. This means that the original random dis-
tribution of crystals in the undeformed sample becomes non-random. The deformed
material is said to become crystallographically textured [4].
Most polycrystalline materials show texture due to processing. Texture can arise
during solidification when those crystals which have a fast growth direction parallel
to that of the heat flow will dominate the final structure. It can arise during recrys-
tallisation and phase transformation when selective nucleation leads to the formation
of a biased distribution of crystals.
A convenient method for communicating texture caused by deformation, for ex-
ample, in rolled sheet, is by stating the set {h k l}⟨u v w⟩, the planes which lie
roughly parallel to the rolling plane, and direction in the rolling plane that tends to
be parallel to the rolling direction. The overall texture can be represented as the sum
of components: .
texture = λi {h k l}i ⟨u v w⟩i (6.2)
i
where λ represents the weighting given to a particular type of texture. The ma-
jor components of the deformation texture of austenite are {1 1 0}⟨1 1 2⟩ and
{1 1 2}⟨1 1 1⟩, the so-called brass and copper textures respectively. Recrystalli-
sation in a ferritic steel often leads to a cube texture {1 0 0}⟨0 0 1⟩ cube component
and special thermomechanical processing of the type used to produce magnetically
soft metal is associated with the Goss texture {1 1 0}⟨0 0 1⟩.
Texture can be plotted on a stereogram. A pole figure in this context consists of a
stereogram with its axes defined relative to the external frame of reference, and with
particular hkl poles from each of the crystals in the polycrystalline aggregate plotted
on to it. If the distribution of crystals is random, the pole figure would appear as in
Fig. 6.5a, and if the distribution is non-random then this would be apparent in the
pole figure, as shown in Fig. 6.5b.
An inverse pole figure is one in which the sample frame is plotted relative to the
crystal axes of a reference crystal. Fig. 6.6 illustrates an example where the rolling,
normal and transverse directions (RD, ND and TD respectively) are plotted relative
to the crystal axes as defined by the stereographic triangle. In a cubic system, the 24
stereographic triangles with ⟨100⟩, ⟨101⟩, and ⟨111⟩ at the corner of each triangle, are
crystallographically equivalent. Therefore, the orientation of the rolling direction, for
example, from each of the 24 triangles can be plotted on just one of the stereographic
‡ This frame may, for example, consist of a set of vectors parallel to the rolling direction, transverse
direction and thickness of a rolled plate; alternatively, the principal axes of an applied system of stresses.
6.2 Texture 97
(a) (b)
FIGURE 6.5
(a) Figure showing a random distribution of 100 poles from a polycrystalline sample
containing 500 crystals, plotted relative to the rolling direction (RD) and transverse
direction (T D). Notice that the distribution of poles does not appear random but this
is an artifact due to the angular distortion on a stereographic projection. (b) Corre-
sponding pole figure after texture is introduced by phase transformation. Here there
is a strong tendency for the 100 poles to cluster around the sample axes. The distri-
bution of poles can be plotted as contours rather than the dots illustrated. Diagrams
courtesy of Saurabh Kundu.
triangles in the form of poles, or contour plots representing the design of poles in
angular space.
Crystallographic texture can be measured using a variety of techniques such as
X-ray diffraction or electron back-scattered diffraction (EBSD) [5]. In the former
case, the polycrystalline sample is exposed in a system set to detect X-rays from a
particular reflection (Bragg angle) and the sample is systematically tilted and rotated
in order to capture intensity at a variety of orientations. In EBSD, an electron beam is
rocked about a fixed point on the surface of a grain using a scanning electron micro-
scope. At particular angles, the beam is Bragg diffracted, resulting in a reduction in
intensity picked up by the detector, leading to the formation of channelling patterns
of the type illustrated in Fig. 6.7. These patterns can be interpreted to determine the
orientation of the crystal relative to the sample axes.
98 6 Deformation and Texture
1 1
1 1
FIGURE 6.6
Inverse pole figures for a cubic crystal. The diagrams at the top have poles corre-
sponding to the rolling, normal and transverse directions plotted relative to the crys-
tal axes, whereas those at the bottom have the same data represented as contours
representing the density of poles.
FIGURE 6.7
A simulated electron channelling pattern from silicon carbide. Reproduced from
Winkelmann et al. [6], with permission of Elsevier.
6.3 Orientation distribution functions 99
ND
001
010
normal to ND
(c) (d)
FIGURE 6.8
Operations that define the three (Bunge) Euler angles. (a) Initial relationship between
crystal and sample coordinates. (b) Rotation by ϕ1 about ND to generate RD’ in a
position normal to the plane containing ND and 001. (c) Rotate by Φ about RD’ to
generate ND’ parallel to 001. (d) Rotate by ϕ2 about 001 to bring bring crystal and
sample axes into coincidence. The insets show the rotations relative to an aircraft.
The set of measured Euler angles made by the RD, T D and N D directions for a
100 6 Deformation and Texture
(a) (b)
FIGURE 6.9
(a) An orientation distribution plotted on a cube with axes defined by the Euler an-
gles. Each dot represents a crystal orientation relative to sample axes. (b) Sections of
the cube to show how the density of orientations varies.
particular crystal is plotted as a point in a cube whose axes scale with ϕ1 , Φ and ϕ2 .
Each point inside this cube is a single crystal.
There are very many industrial processes which control and exploit texture on
a grand scale. Beverage cans are incredibly thin even though their bodies are made
without fracture using extreme deformation including deep drawing. The texture of
the material prior to drawing is controlled such that plastic instability is avoided in
the thickness direction. For similar reasons, the steels used in the manufacture of
formed automobile bodies are texture-controlled. Iron alloys used to make electric-
ity transformers have to be magnetically soft. This is achieved by controlling the
crystallographic texture to minimise magnetic losses.
010
010
100
FIGURE 6.10
(a) Note that [001]A ∥ [001]B . (b) Construction showing that the projection is also
a circle.
The Euler angles define the operations needed to bring the two frames of refer-
ence into coincidence. Referring to Fig. 6.8b, since [001]A ∥ [001]B , ϕ1 = 0◦ . For
the same reason, Φ = 0◦ ( Fig. 6.8c) and the remaining angle ϕ2 is clearly 45◦ .
Fig. 6.8b shows the ODF section φ1 = 0, referred to the basis-B as the reference
frame.
6.4 Summary
Crystallographic texture can be understood by a thought experiment in which a large
single crystal is broken up into many smaller ones which are then joined together in
somewhat different orientations that are not random. The material then behaves as if
it has some of the characteristics of a single crystal (such as anisotropy) while other
properties are more a feature of random polycrystals (such as the multiple slip sys-
tems needed per crystal in order to achieve plasticity without the creation of voids).
We have discussed how texture can arise from plastic deformation, but many
kinds of processing can lead to its development. The heat treatment of deformed
materials to induce recrystallisation can lead to recrystallisation textures because not
all ‘nuclei’ grow at the same rate [7]. The nuclei have to be sufficiently different in
orientation from their surroundings to accomplish rapid growth. Similar selection can
occur in phase transformations where the product phase forms at a grain boundary
because certain boundaries make it possible for the product to adopt a favourable
orientation relationship with the parent phase on either side of the interface [8, 9].
102 6 Deformation and Texture
Both magnetic [10] and electrical fields [11] can also lead to orientational order in
appropriate materials.
So how does one produce a polycrystalline material with a random arrangement
of grains? The short answer is, with enormous difficulty. Powder techniques can
partially reach this goal; metallic or ceramic powders that are sintered together tend
to undergo very limited deformation or rotation during fabrication, and hence may
approach a random orientational distribution.
References
1. E. Schmid, and W. Boas: Plasticity of Crystals (translated from the 1935 edition
of Kristalplastizitaet): London, U.K.: F. A. Hughes and Co., 1950.
2. J. W. Christian: ‘Deformation by moving interfaces’, Metallurgical Transac-
tions A, 1982, 13, 509–538.
3. A. Kelly, and K. M. Knowles: Crystallography and Crystal Defects: 2nd ed.,
New York, USA: John Wiley & Sons, Inc., 2012.
4. D. N. Lee: Texture and related phenomena: Seoul, Republic of Korea: Korean
Institute of Metals and Materials, 2006.
103
7
Interfaces, Orientation Relationships
Abstract
7.1 Introduction
Much of the science and technology of polycrystalline materials depends on the na-
ture of interfaces between crystals and the relative orientations of adjacent crystals.
The corrosion resistance of a boundary depends on how coherent it is [1], the conti-
nuity of deformation varies with texture [2] . . . the list is seemingly endless!
Atoms in the boundary between crystals must in general be displaced from po-
sitions they would occupy in the undisturbed crystal, but it is now well established
that many interfaces have a periodic structure. In such cases, the misfit between the
crystals connected by the boundary is not distributed uniformly over every element
of the interface; it is localised periodically into discontinuities that separate patches
of the boundary where the fit between the two crystals is good or perfect. When these
discontinuities are well separated, they may individually be recognised as interface
dislocations which separate coherent patches in the boundary, which is macroscop-
ically said to be semi-coherent. The simple example of a symmetrical tilt boundary
illustrates this.
105
106 7 Interfaces, Orientation Relationships
(a) (b)
FIGURE 7.1
The creation of an interface. (a) Single crystal cut in two halves which are then tilted
through an angle θ about an axis normal to the diagram. (b) The gap within the tilt is
filled with edge dislocations, each with a Burgers vector b and spacing d to create an
interface.
The structure of the interface in terms of dislocations provides a mechanism for
calculating the energy per unit area of the interface as a function of the misorientation
θ. The energy per unit length Ws of a dislocation is given by [3]
3
1 r∞ Gb2 dr
Ws = Wc +
2 r0 2π(1 − ν) r
4 5
Gb2 r∞
= Wc + ln (7.2)
4π(1 − ν) r0
where Wc is the core energy per unit length, covering the region of the dislocation
where elasticity theory fails; G and ν are the shear modulus and Poisson’s ratio,
respectively; r0 is the radius of the dislocation core and r∞ is a cut-off radius beyond
which the dislocation strain field can be neglected or is limited by the size of the
7.3 Symmetrical tilt boundary 107
sample. In arrays of dislocations of the type found in interfaces, the strain fields of
individual dislocations are partly compensated by their neighbours so it is a good
approximation that r∞ ≈ d in which case r∞ ∝ 1/θ. It follows that equation 7.2 can
be written Ws = A(B − ln θ). Therefore, the interfacial energy per unit area σi is
related to the number of dislocations per unit area (1/d):
1
σi = × Ws
d
2 θ
= tan × A(B − ln θ). (7.3)
b 2
This gives a model for the energy of grain boundaries [4], based on the density of
dislocations in the boundary (∝ d−1 ). The variation in energy as a function of mis-
orientation in the symmetrical tilt boundary described above is illustrated in Fig. 7.2.
The model is not valid for large misorientations where the dislocation spacing be-
comes comparable with the magnitude of the Burgers vector, since the dislocation
cores then begin to overlap.
FIGURE 7.2
The interfacial energy as a function of the misorientation for a tilt boundary. The
dashed curve represents cusps when the orientation forms a coincidence site lattice.
(a) (b)
FIGURE 7.3
(a) A rotation of 36.9◦ about ⟨100⟩ defines the orientation relationship between two
primitive cubic crystals, each of lattice parameter 1, leading to a Σ5 coincidence site
lattice. The origin is labelled ‘o’. (b) Simplified illustration of the fact that there are
only one in five of the lattice points of each of the cubic crystals in coincidence – the
four in the middle of the cell are not.
(a) (b)
FIGURE 7.4
(a) Cubic coordinate system with orthogonal basis vectors of equal magnitude. (b)
Tetragonal coordinate system (red) in which the basis vectors on the basal plane end
at the face-centres of the cubic system. The relationship to the cubic cells (black) is
also illustrated.
and transforms the components of vectors referred to the A basis to those referred to
the B basis. The first column of (B J A) represents the components of the basis vector
a1 , with respect to the basis B, and so on.
The components of a vector u can now be transformed between bases using the
matrix (B J A) as follows:
010
010
100
100
FIGURE 7.5
Diagram illustrating the relation between the bases A and B. Note that [001]A ∥
[001]B .
Referring to Fig. 7.5, and recalling that the matrix (B J A) consists of three
columns, each column being the components of one of the basis vectors of A,
with respect to B, we have
u1 u3 (1 − m) + u2 n u2 u3 (1 − m) − u1 n
u3 u3 (1 − m) + m
(7.9)
where m = cos θ and n = sin θ. The right-handed angle of rotation can be obtained
from the fact that
J11 + J22 + J33 = 1 + 2 cos θ (7.10)
and the components of the vector u along the axis of rotation are given by
or
2π 2π 2π 2π
θ = 2π, , , , , (7.14)
2 3 4 6
another way of showing that a five-fold rotation axis is not permitted when transla-
tional symmetry is to be maintained.
7.5 Mathematical method for determining Σ 113
x does not necessarily have integral components in the A basis (i.e. it need not be
a lattice vector of A). CSL vectors, on the other hand, identify lattice points which
are common to both A and B, and therefore are lattice vectors of both crystals. It
follows that CSL vectors have integral indices when referred to either crystal. Hence,
x is only a CSL vector if it has integral components in the basis A. We note that x
always has integral components in B, because a lattice vector of A (such as u) always
deforms into a lattice vector of B.
The meaning of Σ is that 1/Σ of the lattice sites of A or B are common to both
A and B. It follows that any primitive lattice vector of A or B, when multiplied
by Σ, must give a CSL vector. Σx must therefore always be a CSL vector and if
equation 7.15 is multiplied by Σ, then we obtain an equation in which the vector u
always transforms into a CSL vector:
i.e.. given that u is a lattice vector of A, whose components have no common factor,
Σx is a CSL vector with integral components in either basis. This can only be true
if the matrix Σ(A S A) has elements which are all integral since it is only then that
Σ[A; x] has elements which are all integral.
It follows that if an integer H can be found such that all the elements of the matrix
H(A S A) are integers (without a common factor), then H is the Σ value relating A
and B.
Applying this to the problem considered in section 7.3, the rotation matrix corre-
sponding to the rotation 180◦ about [1 1 2]A is given by (equation 7.9)
⎛ ⎞
2 1 2
1⎜ ⎟
(A J A) = ⎜ 1 2 2
⎟ (7.17)
3 ⎝ ⎠
2 2 1
and since 3 is the integer which when multiplied with (A J A) gives a matrix of
integral elements (without a common factor), the Σ value for this orientation is given
by Σ = 3.
114 7 Interfaces, Orientation Relationships
Coincidence site lattices have in recent years become popular because of the
advent of orientation imaging in the scanning electron microscope. The microscope
is usually associated with software that enables the estimation of the Σ value at every
junction between two grains. The accuracy of the technique is limited so such Σ maps
should be taken with a pinch of salt, and large Σ values are not very meaningful as
representations of low-energy boundaries.
7.6 Summary
Interfaces clearly have a structure that can be periodic. Simple boundaries between
misoriented crystals that have identical structures can be described in terms of arrays
of dislocations which become visible during transmission electron microscopy. The
spacing of these boundary dislocations can be related to the degree of misorientation
whereas their line vectors, in the case of a tilt boundary, lie along the tilt axis. The
dislocation structure can be used also to estimate the interfacial energy per unit area
as long as the dislocation description remains valid.
It is possible that at certain values of the misorientations, a significant fraction
of the lattice points between the two crystals that form a boundary, are coincident.
Furthermore, the coincident sites form a lattice that is periodic, known as a coinci-
dence site lattice. It is important to realise that this is a reflection of the orientation
relationships between the crystals rather than the boundary plane. The properties of a
boundary change significantly when the fraction of coincident points, 1/Σ, becomes
large. For example, the interfacial energy per unit area declines dramatically as might
the boundary diffusion coefficient.
References
115
8
Crystallography of Martensitic
Transformations
Abstract
8.1 Introduction
Martensitic transformations, and displacive transformations in general, remain the
only mechanism for the manufacture of bulk nanostructured materials for engineer-
ing applications [1, 2]. The transformation is diffusionless and hence can occur at
incredibly low temperatures [3], or can propagate at the speed of sound in the ma-
terial (a few thousand m s−1 ) [4]. There is no composition change when it grows to
consume the parent phase.
There are, however, a number of peculiarities associated with martensitic trans-
formations, which are important in the design of alloys. The interface plane between
austenite and martensite is known as the habit plane and measurements show that the
crystallographic indices of that plane are irrational and strange (Table 8.1). Another
difficulty is that the martensite-austenite (α/γ) interface has to be glissile, that is, it
must be able to move without diffusion. This requires that the interface must contain
no more than one set of dislocations. Any more than one array of line vectors can lead
to interference between the dislocations which renders the interface sessile. In crys-
tallographic terms, this means that there must be at least one line fully coherent in
the α/γ interface. The transformation strain relating the two lattices must therefore,
as a minimum, be an invariant-line strain.
117
118 8 Crystallography of Martensitic Transformations
TABLE 8.1
Habit plane indices for martensite. The quoted indices are approximate because the
habit planes are in general irrational. There is an interesting study on the regularity
of the habit plane indices as a function of, for example, strain in the austenite [5].
FIGURE 8.1
The measured shape deformation due to martensitic transformation is an invariant-
plane strain with a large shear component (≃ 0.22) and a small dilatational strain
(≃ 0.03) directed normal to the habit plane in steels.
[0 0 1]fcc | |[0 0 1]bcc [1 1 0]fcc | |[1 0 0]bcc [1 1 0]fcc | |[0 1 0]bcc (8.1)
and although the orientations actually observed are consistent with good matching
between the two lattices (Fig. 8.3) but they are inconsistent with the Bain relation-
ship. Typical observed orientation relationships are irrational but are represented ap-
proximately by:
{1 1 1}γ ∥ {0 1 1}α ,
⟨1 0 1⟩γ ∥ ⟨1 1 1⟩α ,
120 8 Crystallography of Martensitic Transformations
(a)
(b) (c)
FIGURE 8.2
The Bain strain. The iron atoms at the centres of the faces of the austenite cell paral-
lel to the plane of the diagram, are omitted for clarity. (a) The lattice correspondence
for formation of martensite from austenite, showing a single carbon atom in an octa-
hedral interstice on the [001]γ axis. (b) Tetragonal unit cell outlined in austenite. (c)
Lattice deformation (compression along c-axis) to form martensite with an appropri-
ate c/a ratio. Note that the martensite in iron alloys can be body-centred tetragonal
(bct) or body-centred cubic (bcc); none of the descriptions in this chapter are sensi-
tive to this.
{1 1 1}γ ∥ {0 1 1}α ,
FIGURE 8.3
Electron diffraction pattern from bcc martensite and fcc austenite lattices in steels.
FIGURE 8.4
(a) and (b) show the effect of the Bain strain on austenite, which when undeformed
is represented as a sphere of diameter wx = yz in three-dimensions. The strain trans-
forms it into an ellipsoid of revolution. (c) Shows the ILS obtained by combining
the Bain strain with a rigid body rotation through an angle θ. a1 , a2 and a3 refer to
[100]γ , [010]γ and [001]γ axes respectively.
122 8 Crystallography of Martensitic Transformations
As can be seen from Fig. 8.4c, there is no rotation which can make B into an
invariant-plane strain since this would require two non-parallel invariant-lines. Thus,
for the fcc→bcc transformation, austenite cannot be transformed into martensite by
a homogeneous strain which is an invariant-plane strain. And yet, the observed shape
deformation leaves the habit plane undistorted and unrotated, i.e., it is an invariant-
plane strain.
The phenomenological theory of martensite crystallography solves this remain-
ing problem (Fig. 8.5) [7–10]. The Bain strain converts the structure of the parent
phase into that of the product phase. When combined with an appropriate rigid body
rotation, the net homogeneous lattice deformation RB is an invariant-line strain (step
a to c in Fig. 8.5). However, the observed shape deformation is an invariant-plane
strain P1 (step a to b in Fig. 8.5), but this gives the wrong crystal structure. If a
second homogeneous shear P2 is combined with P1 (step b to c), then the correct
structure is obtained but the wrong shape since
P1 P2 = RB.
These discrepancies are all resolved if the shape changing effect of P2 is cancelled
macroscopically by an inhomogeneous lattice-invariant deformation, which may be
slip or twinning as illustrated in Fig. 8.5.
The theory explains all the observed features of the martensite crystallography.
The orientation relationship is predicted by deducing the rotation needed to change
the Bain strain into an invariant-line strain. The habit plane does not have rational
indices because the amount of lattice-invariant deformation needed to recover the
correct macroscopic shape is not usually rational. The theory predicts a substruc-
ture in plates of martensite (either twins or slip steps) as is observed experimentally.
The transformation goes to all the trouble of ensuring that the shape deformation
is macroscopically an invariant-plane strain because this reduces the strain energy
when compared with the case where the shape deformation might be an invariant-
line strain.
When austenite transforms martensitically into ε (hcp), the entire change is
achieved by a deformation that is an invariant-plane strain as illustrated in Fig. 8.6.
As a consequence, the habit plane of the martensite is rational, i.e., precisely the
close-packed plane, the orientation relationship is rational:
The strain not only transforms the lattice but the resulting shape of the final transfor-
mation product is as observed, so there is no need for any lattice-invariant deforma-
tion (slip or twinning).
8.4 Summary 123
FIGURE 8.5
The phenomenological theory of martensite crystallography.
FIGURE 8.6
aγ
The fcc→hcp transformation involving a shear of 6 [112] on alternate (111)γ planes
to generate the hcp structure.
8.4 Summary
It is hard to exaggerate the success of the crystallographic theory of martensitic trans-
formations, both in terms of the conceptual leap in explaining away the apparent dif-
124 8 Crystallography of Martensitic Transformations
ficulties accumulated via careful experiments, and in the quantitative closure between
experiment and theory. This has not, in many respects, been achieved for any other
solid-state transformations, particularly those that involve diffusion and composition
change. The theory has generic value because it is not restricted to changes between
particular crystal systems.
References
125
Part II
Abstract
Macroscopic materials consist of aggregates of crystals. The order that typifies sin-
gle crystals is disrupted at boundaries so the properties of polycrystalline samples
can be sensitive to the relative orientations of neighbouring crystals. Or even to the
averaged disposition of the myriads of crystals relative to some macroscopic frame
of reference. Orientation relationships between crystals are useful in understanding
the structure of the interfaces between crystals and the evolution of the shapes of pre-
cipitates as they grow in a solid matrix. An entire subject now exists on grain bound-
ary engineering, where, for example, the corrosion resistance of a stainless steel can
be improved dramatically by ensuring a preponderance of good-fit boundaries. And
there are remarkable experimental techniques that facilitate the crystallographic char-
acterisation of millions of crystals in a long blink of an eye.
9.1 Introduction
A substantial part of research on polycrystalline materials is concerned with the accu-
rate determination, assessment and theoretical prediction of orientation relationships
between adjacent crystals. There are obvious practical applications, such as in the
study of anisotropy and texture, or in defining the role of interfaces in controlling
the mechanical properties. The existence of a reproducible orientation relation be-
tween the parent and product phases might determine the ultimate morphology of
any precipitates, by allowing the development of low interfacial-energy facets. It is
possible that vital clues about nucleation in the solid state might emerge from a de-
tailed examination of orientation relationships, even though these can usually only
be measured when the crystals concerned are well into the growth stage. Needless to
say, the properties of interfaces depend critically on the relative dispositions of the
crystals that they connect.
Perhaps the most interesting experimental observation is that the orientation rela-
tionships that are found to develop during phase transformations (or during deforma-
tion or recrystallisation experiments) usually are not random [1–3]. The frequency of
occurrence of any particular orientation relation usually far exceeds the probability
of obtaining it simply by taking the two separate crystals and joining them up in an
arbitrary way.
129
130 9 Orientation relations
This indicates that there are favoured orientation relations, perhaps because it is
these which allow the best fit at the interface between the two crystals. This would in
turn reduce the interfacial energy, and hence the activation energy for nucleation. Nu-
clei which, by chance, happen to be orientated in this manner would find it relatively
easy to grow, giving rise to the non-random distribution mentioned earlier.
On the other hand, it could be the case that nuclei actually form by the homo-
geneous deformation of a small region of the parent lattice. The deformation, which
transforms the parent structure to that of the product (e.g. Bain strain), would have
to be of the kind which minimises strain energy. Of all the possible ways of accom-
plishing the lattice change by homogeneous deformation, only a few might satisfy the
minimum strain energy criterion - this again would lead to the observed non-random
distribution of orientation relationships. It is a problem to understand which of these
factors really determine the existence of rational orientation relations. In this chap-
ter we deal with the methods for adequately describing the relationships between
crystals.
1 The structure of cementite is described on page 68 but we will assume in what follows that the lattice
The information given indicates parallelisms between vectors in the two lat-
tices. In order to find (α J θ), it is necessary to ensure that the magnitudes of
the parallel vectors are also equal, since the magnitude of a vector must remain
invariant to a coordinate transformation. If the constants k, g and m are defined as
|[1 0 0]θ | a√
k= |[0 1 1]α |
= aα 2
= 1.116120
|[0 1 0]θ | b√
g= |[1 1 1]α |
= aα 3
= 1.024957 (9.2)
|[0 0 1]θ | c√
m= |[2 1 1]α | = aα 6
= 0.960242
then multiplying [0 1 1]α by k makes its magnitude equal to that of [1 0 0]θ ; the
constants g and m can similarly be used for the other two α vectors.
Recalling now the definition of a coordinate transformation matrix, each col-
umn of (α J θ) represents the components of a basis vector of θ in the α basis.
For example, the first column of (α J θ) consists of the components of [1 0 0]θ in
the α basis, [0 k k]α . It follows that (α J θ) can be derived simply by inspection
of the relations 9.1, 9.2, so that
⎛ ⎞
0.000000 1.024957 1.920485
⎜ ⎟
(α J θ) = ⎜ ⎝−1.116120 −1.024957 0.960242⎠
⎟ (9.3)
1.116120 −1.024957 0.960242
132 9 Orientation relations
The transformation matrix can therefore be deduced by inspection when the
orientation relationship 9.1 is stated in terms of the basis vectors of one of the
crystals concerned (in this case, the basis vectors of θ are specified in relation-
ship 9.1). On the other hand, orientation relationships can, and often are, specified
in terms of vectors other than the basis vectors. Furthermore, electron diffraction
patterns may not include direct information about basis vectors. A more general
method of solving for (α J θ) is now presented, a method independent of the
vectors used in specifying the orientation relationship:
From the relations 9.1, 9.2 we see that
6 7
0 k k α = (α J θ)[1 0 0]θ
6 7
g g g α = (α J θ)[0 1 0]θ
6 7
2m m m α = (α J θ)[0 0 1]θ . (9.4)
where the Jij (i = 1, 2, 3 & j = 1, 2, 3) are the elements of the matrix (α J θ).
From equation 9.5, it follows that
⎛ ⎞ ⎛ ⎞
0 g 2m 0 1.024957 1.920485
⎜ ⎟ ⎜ ⎟
(α J θ) = ⎜
⎝k g m ⎠ = ⎝−1.116120 −1.024957 0.960242⎠ . (9.6)
⎟ ⎜ ⎟
⎛ ⎞
0 −3gm 3gm
1 ⎜ ⎟
(θ J α) = ⎜
2mk −2mk −2mk
⎟
6gmk ⎝ ⎠
2gk gk gk
⎛ ⎞
0 −0.447981 0.447981
⎜ ⎟
= ⎜
⎝0.325217 −0.325217 −0.325217⎠
⎟
It should be noted that the determinant of (α J θ) gives the ratio of the volume
of the θ unit cell to that of the ferrite α unit cell. If the coordinate transformation
simply involves a rotation of bases, e.g. when it describes the relation between two
grains of identical structure, then the matrix is orthogonal and its determinant has a
value of unity for all proper rotations (i.e. not involving inversion operations). Such
coordinate transformation matrices are called rotation matrices.
A stereographic representation of the Bagaryatski orientation is presented in
Fig. 9.1, to visualise the angular relationships between poles (plane normals) of crys-
tal planes and give some indication of symmetry; the picture is of course distorted
because distance on the stereogram does not scale directly with angle. Angular mea-
surements on stereograms are in general made using Wulff nets and may typically
be in error by a few degrees, depending on the size of the stereogram. Space and
aesthetic considerations usually limit the number of poles plotted on stereograms,
and those plotted usually have rational indices. Separate plots are needed for cases
where directions and plane normals of the same indices have a different orientation
in space. A coordinate transformation matrix is, by comparison, a precise way of rep-
resenting orientation relationships; angles between any plane normals or directions
can be calculated to any desired degree of accuracy and information on both plane
normals and directions can be derived from just one matrix. With a little nurturing,
it is possible to picture the meaning of the elements of a coordinate transformation
matrix: each column of (α J θ) represents the components of a basis vector of θ in
the basis α, and the determinant of this matrix gives the ratio of volumes of the two
unit cells.
134 9 Orientation relations
(a)
(b) (c)
FIGURE 9.1
(a) Stereographic representation of the Bagaryatski orientation relationship between
ferrite and cementite in steels, where [1 0 0]θ ∥ [0 1 1]α , [0 1 0]θ ∥ [1 1 1]α , and
[0 0 1]θ ∥ [2 1 1]α . All the poles plotted are in the upper hemisphere. (b) Typical
electron diffraction pattern consistent with the Bagaryatski orientation relation with
[1 0 0]θ ∥ [0 1 1]α normal to the plane of the pattern. (c) An actual electron diffraction
pattern corresponding to (b).
9.3 Relations between fcc and bcc crystals 135
⎛ ⎞⎛ ⎞
0 g 2m 2/6 2/6 2/6
⎜ ⎟⎜ ⎟
(X J γ) = ⎜
⎝k g
⎟⎜
m ⎠ ⎝3/6 0 3/6⎠
⎟
In an similar manner,
⎛ ⎞
0.741582 −0.666667 0.074915
aγ ⎜ ⎟
(Y J γ) = ⎜
−0.983163⎠ .
⎟ (9.8)
aα ⎝0.166667 0.074915
0.649830 0.741582 0.166667
it follows that
(X J γ)[γ; u] = (X J Y)[Y; u]
substituting for [Y;u], we get
so that
(X J Y) = (X J γ)(γ J Y)
9.3 Relations between fcc and bcc crystals 137
carrying out this operation, the required X-Y orientation relation is obtained as
follows:
⎛ ⎞
0.988776 0.147308 −0.024972
⎜ ⎟
(X J Y) = ⎜ ⎝−0.024972 0.327722 0.944445 ⎠ .
⎟
(a)
(b)
FIGURE 9.2
(a) Stereographic representation of the Kurdjumov-Sachs orientation relationship.
Note that the positions of the base vectors of the γ lattice are consistent with the ma-
trix (X J γ) derived in example 9.3, equation 9.8, where each column of (X J γ) rep-
resents the components of a basis vector of γ in the basis X, so that [1 0 0]γ , [0 1 0]γ
and [0 0 1]γ are approximately parallel to [1 1 0]α , [1 1 0]α and [0 0 1]α respec-
tively, as seen in the stereographic representation. (b) Stereographic representation of
the Nishiyama-Wasserman orientation relationship which can be generated from the
Kurdjumov-Sachs orientation by a rotation of 5.26◦ about [0 1 1]α . The necessary ro-
tation makes [1 1 2]γ exactly parallel to [0 1 1]α . The Nishiyama-Wasserman relation
is midway between the two variants of the Kurdjumov-Sachs relation which share the
same parallel close-packed plane. The stereograms also show that the Kurdjumov-
Sachs and Nishiyama-Wasserman orientation relationships do not differ much from
the γ/α orientation relationship implied by the Bain strain illustrated in Fig. 8.2.
9.4 Relationships between grains of identical structure 139
where the basis vectors of X and Y define the respective bcc unit cells. Show that
the crystal Y can be generated from X by a right-handed rotation of 60◦ about an
axis parallel to the [1 1 1]X direction.
A rigid body rotation leaves the magnitudes and relative positions of all vec-
tors in that body unchanged. For example, an initial vector u with components
[u1 u2 u3 ]X relative to the X basis, due to rotation becomes a new vector x, with
the same components [u1 u2 u3 ]Y , but with respect to the rotated basis Y. The fact
that x represents a different direction than u (due to the rotation operation) means
that its components in the X basis, [w1 w2 w3 ]X must differ from [u1 u2 u3 ]X .
The components of x in either basis are obviously related by
[Y; x] = (Y J X)[X; x]
in other words,
[u1 u2 u3 ] = (Y J X)[w1 w2 w3 ]. (9.9)
However, if u happens to lie along the axis of rotation relating X and Y, then
not only will [X;u] = [Y;x] as before, but its direction also remains invariant to
the rotation operation, so that [X;x] = [Y;x] and
(Y J X)[X; x] = [Y; x]
so that
(Y J X)[X; u] = [X; u]
and hence
{(Y J X) − I}[X; u] = 0 (9.10)
140 9 Orientation relations
where I is a 3 × 3 identity matrix. Any rotation axis must satisfy an equation of
this form; expanding equation 9.10, we get
1 2 1
− u1 + u2 − u3 =0
3 3 3
1 1 2
− u1 − u2 + u3 =0
3 3 3
2 1 1
u1 − u2 − u3 = 0.
3 3 3
Solving these simultaneously gives u1 = u2 = u3 , proving that the required
rotation axis lies along the [1 1 1]X direction, which also is the [1 1 1]Y direction.
The angle, and sense of rotation can be determined by examining a vector
v which lies at 90◦ to u. If, say, v = [1 0 1]X , then as a result of the rotation
operation it becomes z = [1 0 1]Y = [0 1 1]X , making an angle of 60◦ with v,
giving the required angle of rotation. Since v ∧z gives [1 1 1]X , it is also a rotation
in the right-handed sense.
9.4.1 Comments
(i) The problem illustrates the fact that the orientation relation between two grains
can be represented by a matrix such as (Y J X), or by an axis-angle pair such as
[1 1 1]X and 60◦ . Obviously, the often used practice of stating a misorientation
between grains in terms of just an angle of rotation is inadequate and incorrect.
(ii) The Euler theorem for rigid rotations [17] states that “regardless of the way a
coordinate system is rotated from its original orientation, it always is possible to
find a fixed axis in space about which a single rotation of the initial coordinates
ends at the final orientation” (the original theorem is in Latin - the English quote
is from [18]). If the axis of rotation is constrained to be a unit vector, then only
three independent quantities are needed to define a misorientation between grains:
two components of the axis of rotation, and an angle of rotation. It follows that
a rotation matrix must also have only three independent terms. The relationship
between the axis-angle pair and rotation matrix is given in equation 7.9.
From the definition of a coordinate transformation matrix, each column of (Y J X)
gives the components of a basis vector of X in the basis Y. It follows that
[1 0 0]X ∥ [2 1 2]Y [0 1 0]X ∥ [2 2 1]Y [0 0 1]X ∥ [1 2 2 ]Y .
Suppose now that there exists another ferrite crystal (basis symbol Z), such that
[0 1 0]Z ∥ [2 1 2]Y [1 0 0]Z ∥ [2 2 1]Y [0 0 1]Z ∥ [1 2 2]Y
⎛ ⎞
2 2 1
1⎜ ⎟
(Y J Z) = ⎝2 1
⎜ ⎟
2⎠
3
1 2 2
9.4 Relationships between grains of identical structure 141
with the crystal Y being generated from Z by a right-handed rotation of 70.52◦ about
[1 0 1]Z direction. It can easily be demonstrated that
⎛ ⎞
0 1 0
⎜ ⎟
(Z J X) = ⎜ ⎝1 0 0⎠
⎟ from (Z J X) = (Z J Y)(Y J X)
0 0 1
so that crystal X can be generated from Z by a rotation of 90◦ about [0 0 1]X axis.
However, this clearly is a symmetry operation of a cubic crystal, and it follows that
crystal X can never be distinguished experimentally from crystal Z, so that the matri-
ces (Y J X) and (Y J Z) are crystallographically equivalent, as are the corresponding
axis-angle pair descriptions. In other words, Y can be generated from X either by a
rotation of 60◦ about [1 1 1]X , or by a rotation of 70.52◦ about [1 0 1]X . The two axis-
angle pair representations are equivalent. There are actually 24 matrices like (Z J X)
which represent symmetry rotations in cubic systems. It follows that a cubic bicrys-
tal can be represented in 24 equivalent ways, with 24 axis-angle pairs. Any rotation
matrix like (Y J X) when multiplied by rotation matrices representing symmetry op-
erations (e.g. (Z J X)) will lead to the 24 axis-angle pair representations. The degen-
eracy of other structures is as follows [19, 20]: cubic (24), hexagonal (12), hexagonal
close-packed (6), tetragonal (8), trigonal (6), orthorhombic (4), monoclinic (2) and
triclinic (1). In general, the number N of axis-angle pairs is given by
N = 1 + N2 + 2N3 + 3N4 + 5N6
where N2 , N3 , N4 and N6 refer to the number of diads, triads, tetrads and hexads in
the symmetry elements of the structure concerned.
Fig. 9.3 is an electron diffraction pattern taken from an internally twinned marten-
site plate in a Fe-4Ni-0.4C wt% steel. It contains two ⟨0 1 1⟩ bcc zones, one from
the parent plate (m) and the other from the twin (t). The pattern clearly illustrates
how symmetry makes it possible to represent the same bi-crystal in terms of more
than one axis-angle pair. This particular pattern shows that the twin crystal can be
generated from the parent in at least three different ways: (i) rotation of 70.52◦ about
the ⟨0 1 1⟩ zone axes of the patterns, (ii) rotation of 180◦ about the {1 1 1} plane
normal, and (iii) rotation of 180◦ about the {2 1 1} plane normal. It is apparent that
these three operations would lead to the same orientation relation between the twin
and the parent lattices.
Example 9.4: ‘Double’ twinning
Plates of bcc martensite in Fe-30.4Ni wt% contain {1 1 2} transformation
twins, the two twin orientations X and Y being related by a rotation of 60◦ about
a ⟨1 1 1⟩ axis. Deformation of this martensite at low temperatures leads to the
formation of twins on {5 8 11} planes, the direction of the twinning shear being
⟨5 1 3⟩. This is a very rare mode of twinning deformation in bcc systems; show
that its occurrence may be related to the fact that such twins can propagate without
any deviation, across the already existing transformation twins in the martensite
[21].
142 9 Orientation relations
(a) (b)
FIGURE 9.3
(a) Electron diffraction pattern from a martensite plate (m) and its twin (t). Spots not
connected by lines (e.g. “dd”) arise from double diffraction. (b) Interpretation of the
diffraction pattern.
[5 1 3]X ∥ ⟨5 3 1⟩Y
(5 8 11)X ∥ {5 11 8}Y .
It follows that {5 8 11} deformation twins can propagate without deviation across
the transformation twins, since the above planes and directions, respectively, are
crystallographically equivalent and indeed parallel. This may explain the occur-
rence of such an unusual deformation twinning mode.
a vector referred to the basis A, to those referred to its reciprocal basis A∗ . (The
symbol G is used, rather than J, to follow convention). This matrix, called the metric,
with components Gij can be determined in exactly the same manner as any other
coordinate transformation matrix. Each column of (A∗ G A) thus consists of the
components of one of the basis vectors of A, when referred to the basis A∗ . For
example,
a1 = G11 a∗1 + G21 a∗2 + G31 a∗3
Taking successive scalar dot products with a1 , a2 and a3 respectively on both sides
yields:
G11 = a1 .a1 , G21 = a1 .a2 G31 = a1 .a3
since ai .a∗j = 0 when i ̸= j. The rest of the elements of (A∗ G A) can be determined
in a similar way, so that
⎛ ⎞
a1 .a1 a2 .a1 a3 .a1
⎜ ⎟
(A∗ G A) = ⎜ ⎝a1 .a2 a2 .a2 a3 .a2 ⎠ .
⎟ (9.11)
a1 .a3 a2 .a3 a3 .a3
It is easily demonstrated that the determinant of (A∗ G A) equals the square of the
volume of the cell formed by the basis vectors of A. We note in passing that for
orthonormal coordinates, (Z∗ G Z) =I. A list of metric tensors is given in Table 9.1.
Example 9.5: Plane normals and directions in an orthorhombic structure
A crystal with an orthorhombic structure has lattice parameters a, b and c. If
the edges of the orthorhombic unit cell are taken to define the basis θ, determine
the metric (θ∗ G θ). Hence derive the equation giving the angle φ between a plane
normal (h; θ∗ ) = (h1 h2 h3 ) and any direction [θ; u] = [u1 u2 u3 ].
From the definition of a scalar dot product, h.u/|h||u| = cos φ. Now,
⎛ ⎞ ⎛ ⎞
a2 0 0 a−2 0 0
⎜ ⎟ ⎜ ⎟
(θ∗ G θ) = ⎜⎝ 0 b2 0 ⎠
⎟ (θ G θ∗ ) = ⎜
⎝ 0 b−2 0 ⎠.
⎟
0 0 c2 0 0 c−2
Similarly,
TABLE 9.1
Elements of metric tensor (A∗ G A) and its inverse. The term φ = (−2 cos{β} cos{α} cos{γ} + cos2 {β} + cos2 {α} + cos2 {γ} − 1)−1
(A∗ G A)
Class G11 G12 G13 G21 G22 G23 G31 G32 G33
Cubic a2 0 0 0 a2 0 0 0 a2
Tetragonal a2 0 0 0 a2 0 0 0 c2
2 2
Orthorhombic a 0 0 0 b 0 0 0 c2
Hexagonal a2 −0.5a2 0 −0.5a2 a2 0 0 0 c2
Monoclinic a2 0 ac cos{β} 0 b2 0 ac cos{β} 0 c2
Triclinic a2 ab cos{γ} ac cos{β} ab cos{γ} b2 bc cos{α} ac cos{β} bc cos{α} c2
(A G A ) ∗
Class G11 G12 G13 G21 G22 G23 G31 G32 G33
Cubic a −2
0 0 0 a −2
0 0 0 a−2
Tetragonal a−2 0 0 0 a−2 0 0 0 c−2
Orthorhombic a−2 0 0 0 b−2 0 0 0 c−2
b2 2 2 4
Hexagonal a2 b2 −0.25a4 − a2 −4b2 0 − a2 −4b2 − a2 −4b2 0 0 0 c−2
csc2 {β} csc2 {β}
Monoclinic a2 0 0 b−2 0 0 c2
2 2 2
− cot{β}accsc{β} − cot{β}accsc{β}
cos{γ} cos{γ} cos{γ} cos{γ}
Triclinic 2
−φ sina{α} ab
−φ cos{β} cos{α}−cos{γ} φ cos{β}−cos{α}
ac ab
−φ cos{β} cos{α}−cos{γ} −φ sinb2{β} φ cos{α}−cos{β}
bc φ cos{β}−cos{α}
ac φ cos{α}−cos{β}
bc −φ sinc2{γ}
145
146 9 Orientation relations
u ∧ v = (u1 a + u2 b + u3 c) ∧ (v1 a + v2 b + v3 c)
where [θ; u] = [u1 u2 u3 ] and [θ; v] = [v1 v2 v3 ]. This equation can be expanded to
give:
Since a.b∧c = V , the volume of the orthorhombic unit cell, and since b∧c = V a∗ ,
it follows that
8 9
∗ ∗ ∗
u ∧ v = V (u2 v3 − u3 v2 )a + (u3 v1 − u1 v3 )b + (u1 v2 − u2 v1 )c .
Hence, u ∧ v gives a vector whose components are expressed in the reciprocal basis.
Writing x = u ∧ v, with (x; θ∗ ) = (w1 w2 w3 ), it follows that w1 = V (u2 v3 −
u3 v2 ), w2 = V (u3 v1 − u1 v3 ) and w3 = V (u1 v2 − u2 v1 ). The cross product of two
directions thus gives a normal to the plane containing the two directions. If necessary,
the components of x in the basis θ can easily be obtained using the metric, since
[θ; x] = (θ G θ∗ )[θ∗ ; x]. Similarly, the cross product of two vectors h and k which
are referred to the reciprocal basis θ∗ , such that (h; θ∗ ) = (h1 h2 h3 ) and (k; θ∗ ) =
(k1 k2 k3 ), can be shown to be:
8 9
1
h∧k = (h2 k3 − h3 k2 )a + (h3 k1 − h1 k3 )b − (h2 k1 − h1 k2 )c .
V
Hence, h ∧ k gives a vector whose components are expressed in the real basis. The
vector cross product of two plane normals gives a direction (zone axis) which is
common to the two planes represented by the plane normals.
9.7 Summary
The term ‘orientation relationship’ has more than one meaning. When the relation-
ship is reproducible, for example when precipitates form in a solid matrix away from
heterogeneous nucleation sites, the strong implication is that there is something that
favours that behaviour. The reasons could include:
1. The reduction in interfacial energy between the precipitate and matrix
9.7 Summary 147
149
150 9 References
Abstract
10.1 Introduction
A homogeneous deformation is one in which points that initially are collinear remain
so after the deformation. It follows that a plane will similarly remain a plane in
such circumstances. Each volume element in the solid therefore experiences the same
change in shape. Fig. 10.1 illustrates the difference between a homogeneous and
heterogeneous deformation. In this chapter, we describe deformations that alter the
length and/or direction of vectors, i.e., physical deformations, with the initial and
resultant vectors referred to a fixed reference frame.
Elastic strain energy calculations often are based on the assumption that the de-
formations are homogeneous [1] and this usually is justified given that the elastic
strain fields may extend over large distances whereas atomic perturbations are just
that. Similarly, the stability of crystal lattices has been studied by subjecting them
to homogeneous deformation along the direction of one of its axes, and to another
small deformation along arbitrary direction [2].
151
152 10 Homogeneous deformations
n
atio
initial
ation
FIGURE 10.1
An illustration of the difference between homogeneous deformations and those that
are heterogeneous. Examples of the former include shear, hydrostatic expansion and
uniaxial expansion, whereas slip caused by the motion of a dislocation, or the for-
mation of a crack, or the motion of rigid continental plates, involve heterogeneous
deformation.
η3 = cα′ /aγ
(a) (b)
FIGURE 10.2
(a) Two austenite unit cells, one the conventional cubic-F (basis A) and the other a
body-centred tetragonal cell (basis B). Some of the atoms in face-centering positions
on the vertical faces of the the cubic cell of austenite are not illustrated for the sake of
clarity. (b) How the body centred tetragonal cell can be deformed by the Bain strain
into that of body-centred cubic or body-centred tetragonal cell of martensite.
(A S A) converts the vector [A;u] into a new vector [A;v], with v still referred to the
basis A.
The difference between a co-ordinate transformation (B J A) and a deformation
matrix (A S A) is illustrated in Fig. 10.3, where ai and bi are the basis vectors of the
bases A and B respectively.
FIGURE 10.3
Difference between coordinate transformation (B J γ) and deformation matrix
(γ S γ). The latter is defined with respect to a single reference frame whereas the
former deals with a change in the reference frame.
Vector components before Bain strain Vector components after Bain strain
[1 0 0]A [η1 0 0]A
[0 1 0]A [0 η2 0]A
[0 0 1]A [0 0 η3 ]A
so the matrix (A S A) can be written as follows, with each column consisting of the
components of a vector equal to the basis vector of A, following deformation:
⎛ ⎞
η1 0 0
⎜ ⎟
(A S A) = ⎜ ⎝ 0 η2 0 ⎠ .
⎟ (10.1)
0 0 η3
Each column of the deformation matrix represents the components of the new
vector formed as a result of the deformation of a vector equal to one of the basis
vectors of A.
The strain represented by (A S A) is called a pure strain because it involves
stretching along three orthogonal axes; for this reason, the matrix representation
(A S A) is symmetrical. Stretching means extensions or contractions but no rota-
tions, and the three initially orthogonal directions are referred to as the principal axes
which obviously remain orthogonal and unrotated during the deformation.
The ratio of its final length to initial length is the principal deformation associated
with that principal axis. The directions [1 0 0]B , [0 1 0]B and [0 0 1]B are therefore
the principal axes of the Bain strain, and ηi the respective principal deformations. In
the particular example under consideration, η1 = η2 so that any two perpendicular
lines in the (0 0 1)B plane could form two of the principal axes. Since a3 ∥ b3 and
since a1 and a2 lie in (0 0 1)B , it is clear that the vectors ai must also be the principal
axes of the Bain deformation.
Since the deformation matrix (A S A) is referred to a basis system which co-
incides with the principal axes, the off-diagonal components of this matrix must be
zero. Column 1 of the matrix (A S A) represents the new co-ordinates of the vector
[1 0 0], after the latter has undergone the deformation described by (A S A), and a
similar rationale applies to the other columns. (A S A) deforms the fcc γ lattice into a
bcc or bct α′ lattice, and the ratio of the final to initial volume of the material is sim-
ply η1 η2 η3 ; the ratio is more generally given by the determinant of the deformation
matrix. Finally, it should be noted that any tetragonality in the martensite can readily
be taken into account by writing η3 = cα′ /aγ , where cα′ /aα is the aspect ratio of
the bct martensite unit cell.
Example 10.1: Bain strain: undistorted vectors
Given that the principal distortions of the Bain strain (A S A), referred to the
crystallographic axes of the fcc γ lattice (lattice parameter aγ ), are η1 = η2 =
10.2 Homogeneous deformations 155
1.123883, and η3 = 0.794705, show that the vector
remains undistorted, though not unrotated as a result of the operation of the Bain
strain. Furthermore, show that for x to remain unextended as a result of the Bain
strain, its components xi must satisfy the equation
Now,
|x|2 = (x; A∗ )[A; x] = a2γ (x21 + x22 + x23 ) (10.3)
and,
|y|2 = (y; A∗ )[A; y] = a2γ (y12 + y22 + y32 ). (10.4)
Using these equations, and the numerical values of xi and yi obtained above,
it is easy to show that |x| = |y|. It should be noted that although x remains
unextended, it is rotated by the strain (A S A), since xi ̸
= yi . Using equations 10.3
and 10.4 with yi = ηi xi , yields the required equation 10.2. Since η1 and η2 are
equal and greater than 1, and since η3 is less than unity, equation 10.2 amounts
to the equation of a right-circular cone, the axis of which coincides with [0 0 1]A .
Any vector initially lying on this cone will remain unextended as a result of the
Bain strain.
This process was illustrated in Fig. 8.4 where a a spherical volume of austenite
was deformed into an ellipsoid of revolution. Notice that the principal axes (ai )
remain unrotated by the deformation, and that lines such as xw and yz which
become x′ w′ and y ′ z ′ respectively, remain unextended by the deformation (since
they are all diameters of the original sphere), although rotated through the angle θ.
The lines xw and yz of course lie on the initial cone described by equation 10.2.
Suppose now, that the ellipsoid resulting from the Bain strain is rotated through
a right-handed angle of θ, about the axis a2 , then Fig. 8.4c illustrates that this
rotation will cause the initial and final cones of unextended lines to touch along
yz, bringing yz and y ′ z ′ into coincidence. If the total deformation can therefore
be described as (A S A) combined with the above rigid body rotation, then such
a deformation would leave the line cd both unrotated and unextended; such a
deformation is called an invariant-line strain. Notice that the total deformation,
consisting of (A S A) and a rigid body rotation is no longer a pure strain, since the
vectors parallel to the principal axes of (A S A) are rotated into the new positions
a′i Fig. 8.4c.
must contain an invariant line. Therefore, the Bain strain must be combined with a
suitable rigid body rotation in order to define the total lattice deformation. This ex-
plains why the experimentally observed orientation relationship between martensite
and austenite does not correspond to that implied by Fig. 10.2. The need to have an
invariant line in the martensite-austenite interface means that the Bain Strain does
not in itself constitute the total transformation strain, which can be factorised into
the Bain Strain and a rigid body rotation. It is this total strain which determines the
final orientation relationship although the Bain Strain accomplishes the total fcc to
bcc lattice change.
in λ; the roots of this equation are the three eigenvalues λi . Associated with each of
the eigenvalues is a corresponding eigenvector whose components may be obtained
by substituting each eigenvalue, in turn, into equation 10.7. Of course, since every
vector which lies along the unrotated direction is an eigenvector, if u is a solution of
equation 10.7 then ku must also satisfy equation 10.7, k being a scalar constant. If the
matrix (A S A) is real then there must exist three eigenvalues, at least one of which
is necessarily real. If (A S A) is symmetrical then all three of its eigenvalues are real;
the existence of three real eigenvalues does not however imply that the deformation
matrix is symmetrical. Every real eigenvalue implies the existence of a corresponding
vector which remains unchanged in direction as a result of the operation of (A S A).
Example 10.2: Eigenvectors and eigenvalues
Find the eigenvalues and eigenvectors of
⎛ ⎞
18 −6 −6
⎜ ⎟
(A S A) = ⎜⎝−6 21
⎟
3⎠
−6 3 21
To solve for the eigenvalues, we use equation 10.8 to form the determinant
: :
: :
:18 − λ −6 −6 :
: :
3 :=0
: :
: 6 21 − λ
: :
: −6 3 21 − λ:
which on expansion gives the cubic equation
(12 − λ)(λ − 30)(λ − 18) = 0
with the roots
λ1 = 12, λ2 = 30 and λ3 = 18.
To find the eigenvector u = [A; u] = [u1 u2 u3 ] corresponding to λ1 , we
substitute λ1 into equation 10.7 to obtain
6u1 − 6u2 − 6u3 =0
−6u1 + 9u2 + 3u3 =0
−6u1 + 3u2 + 9u3 = 0.
These equations can be solved simultaneously to show that u1 = 2u2 = 2u3 . The
other two eigenvectors, v and x, corresponding to λ2 and λ3 respectively, can be
determined in a similar way. Hence, it is found that:
1
[A; u] = (6− 2 )[2 1 1]
1
[A; v] = (3− 2 )[1 1 1]
1
[A; x] = (2− 2 )[0 1 1].
All vectors parallel to u, v or x remain unchanged in direction, though not in
magnitude, due to the deformation (A S A).
158 10 Homogeneous deformations
10.3.1 Comments
(i) Since the matrix (A S A) is symmetrical, we find three real eigenvectors, which
form an orthogonal set.
(ii) A negative eigenvalue implies that a vector initially parallel to the corresponding
eigenvector becomes antiparallel (changes sign) on deformation. A deformation
like this is physically impossible.
(iii) If a new orthonormal basis B is defined, consisting of unit basis vectors parallel
to u, v and x respectively, then the deformation (A S A) can be expressed in the
new basis with the help of a similarity transformation. Keeping like basis symbols
adjacent, it follows that
w1 w2 w3 −6 3 21 u3 v3 w3
⎛ ⎞
18 0 0
⎜ ⎟
= ⎝ 0 30 0 ⎠ .
⎜ ⎟
0 0 12
Notice that the off diagonal terms of (B S B) are zero because it is referred to a
basis formed by the principal axes of the deformation - i.e. the three orthogonal
eigenvectors. The matrix representing the Bain strain (page 154) is also diagonal
because it too is referred to the principal axes of the strain. Any real symmetrical
matrix can be diagonalised using the procedure illustrated above. (B S B) is called
the diagonal representation of the deformation and this special representation will
henceforth be identified by placing a bar over the matrix symbol: (B S B).
On the other hand, the shear (Z P2 Z) illustrated in Fig. 10.4a is not a pure de-
formation because it is possible to identify just two mutually perpendicular eigen-
vectors, both of which must lie in the invariant-plane. All other vectors are rotated
by the shearing action. The original object, represented as a sphere, is sheared into
an ellipsoid. The invariant-plane of the deformation contains the z1 and z2 axes. The
deformation can be imagined to occur in two stages, the first one involving simple
extensions and contractions along the y1 and y3 directions respectively (Fig. 10.4b)
and the second involving a rigid body rotation of the ellipsoid, about the axis z2 ∥ y2 ,
through a right-handed angle φ.
(a) (b)
FIGURE 10.4
Factorisation of a simple shear (Z P2 Z) into a pure deformation (Z Q Z) and a right
handed rigid body rotation of φ about z2 . In (a), ac is the trace of the invariant plane.
(Z Q Z) leaves ac undistorted but rotated to a′ c′ and rigid body rotation brings a′ c′
into coincidence with ac. The axes y1 , y2 and y3 are the principal axes of the pure
deformation. The undeformed shape is represented as a sphere in three dimensions.
It was arbitrarily chosen that the pure strain would occur first and be followed by the
rigid body rotation, but the reverse order is equally acceptable,
(Z P2 Z) = (Z Q2 Z)(Z J2 Z)
where in general,
(Z Q2 Z) ̸
= (Z Q Z) and (Z J2 Z) ̸
= (Z J Z).
Any real deformation can in general be factorised into a pure strain and a rigid body
160 10 Homogeneous deformations
0 0 λ3
It is worth repeating that in equation 10.12, λi are the eigenvalues of the matrix
(Z S′ Z)(Z S Z) and ui , vi and wi are the components, in the basis Z of the eigen-
vectors of (Z S′ Z)(Z S Z). The rotation part of the strain (Z S Z) is simply
10.5 Interfaces
A vector parallel to a principal axis of a pure deformation may become extended
but is not changed in direction by the deformation. The ratio η of its final to initial
length is called a principal deformation associated with that principal axis and the
corresponding quantity (η − 1) is called a principal strain. We have seen that when
two of the principal strains of the pure deformation differ in sign from the third,
all three being non-zero, it is possible to obtain a total strain which leaves one line
invariant (page 121). It intuitively seems advantageous to have the invariant-line in
the interface connecting the two crystals, since their lattices would then match exactly
along that line.
A completely undistorted interface would have to contain two non-parallel direc-
tions which are invariant to the total transformation strain. The following example il-
lustrates the characteristics of such a transformation strain, called an invariant-plane
strain, which allows the existence of a plane which remains unrotated and undistorted
during the deformation.
Example 10.3: Deformations and interfaces
A pure strain (A Q A), referred to an orthonormal basis ‘A’ whose basis vectors
are parallel to the principal axes of the deformation, has the principal deformations
η1 = 1.192281, η2 = 1 and η3 = 0.838728. Show that (A Q A) combined
with a rigid body rotation gives a total strain which leaves a plane unrotated and
undistorted.
Because (A Q A) is a pure strain referred to a basis composed of unit vectors
parallel to its principal axes, it consists of simple extensions or contractions along
the basis vectors a1 , a2 and a3 . Hence, Fig. 10.5 can be constructed by analogy for
the Bain strain illustrated on page 121. Since η2 = 1, ef ∥ a2 remains unextended
and unrotated by (A Q A), and if a rigid body rotation (about fe as the axis of
rotation) is added to bring yz into coincidence with y ′ z ′ , then the two vectors
ef and ab remain invariant to the total deformation. Any combination of ef and
ab will also remain invariant, and hence all lines in the plane containing ef and
y ′ z ′ are invariant, giving an invariant plane. Thus, a pure strain when combined
with a rigid body rotation can only generate an invariant-plane strain if two of its
principal strains have opposite signs, the third being zero. Since it is the pure strain
which actually accomplishes the lattice change (the rigid body rotation causes
no further lattice change), any two lattices related by a pure strain with these
characteristics may be joined by a fully coherent interface.
(A Q A) actually represents the pure strain part of the total transformation
strain required to change an fcc lattice to an hcp (hexagonal close-packed) lattice,
assuming that there is no volume change, by shearing on the {1 1 1}γ plane, in
the ⟨1 1 2⟩γ direction, the magnitude of the shear being equal to half the twinning
shear (see Chapter 11). Consistent with the proof given above, a fully coherent in-
terface is observed experimentally when hcp martensite is formed in this manner.
162 10 Homogeneous deformations
It is worth noting that the total deformation, i.e. the combination of the pure
and rotational strains, completely defines the set of crystallographic features of the
hcp martensite [4–6]. The set consists of the habit plane, orientation relationship
and shape deformation associated with the transformation.
(a) (b)
(c)
FIGURE 10.5
(a,b) Illustration of the strain (A Q A), the undeformed crystal represented initially
as a sphere of diameter ef which also is an invariant line since η2 = 1. The strain
(c) illustrates that a combination of (A Q A) with a rigid body rotation about the
axis ef generates another invariant line z ′ y ′ , so that the net deformation becomes an
invariant-plane strain.
The equiaxed grain structure typical of most undeformed materials can be repre-
sented using Kelvin tetrakaidecahedra which accomplish the ‘division of space with a
minimum of partitional area’ [7]. A tetrakaidecahedron has 8 hexagonal and 6 square
faces, Fig. 10.6, with 36 edges, each of length a. All of the edges can be described in
terms of just six vectors, as listed in Table 10.1. It then becomes possible to operate
on these vectors by a deformation matrix in order to calculate consequential changes
in grain parameters.
2 1
FIGURE 10.6
The Kelvin tetrakaidecahedron.
TABLE 10.1
Vectors defining the edges of a tetrakaidecahedron.
Vector Components
1 [a 0 0]
2 [0 a 0]
3 [− a2 − a2 √a ]
2
4 [ a2 − a2 √a ]
2
5 [ a2 a
2
√a ]
2
6 [− a2 a
2
√a ]
2
164 10 Homogeneous deformations
The results are illustrated in Fig. 10.7. For comparison purposes, the results are
10.6 Topology of grain deformation 165
TABLE 10.2
Components of the six vectors listed in Table 10.1, following plane strain or axisym-
metric deformation.
1 [aS11 0 0]
2 [0 aS22 0]
3 [− aS211 − aS222 aS
√33 ]
2
4 [ aS211 − aS222 aS
√33 ]
2
5 [ aS211 aS22
2
aS
√33 ]
2
6 [− aS211 aS22
2
aS
√33 ]
2
(a) (b)
FIGURE 10.7
Calculations for plane strain deformation. The curve represents the data for the
tetrakaidecahedron oriented as illustrated in Fig. 10.6. The small-points are 99 other
cases where the tetrakaidecahedron is randomly oriented relative to S.
Writing the first rolling pass as S and the cross-rolling pass as T, the net deformation
U is given by T × R × S:
⎛ ⎞⎛ ⎞⎛ ⎞ ⎛ ⎞
T11 0 0 0 1 0 S11 0 0 0 T11 0
⎜ ⎟⎜ ⎟⎜ ⎟ ⎜ ⎟
U=⎜ ⎝ 0 1
⎟⎜
0 ⎠ ⎝−1 0 0⎠ ⎝ 0 1
⎟⎜
0 ⎠ = ⎝−S11
⎟ ⎜ 0 0 ⎟ ⎠.
1
0 0 1/T11 0 0 1 0 0 1/S11 0 0 T11 S11
(10.20)
(a) (b)
FIGURE 10.8
Comparison of cross-rolling (identical strains in the two rolling directions), axisym-
metric compression and single-direction rolling.
using equation 10.20. In this graph, the rolling strains are the same in both directions
(T11 = S11 ) and data for axisymmetric compression are also included. To allow a
comparison between these deformation modes, the data are plotted as a function of
the compressive strain. Notice that the axisymmetric compression and cross-rolling
with equal strains in both directions give exactly identical results, illustrating that
the results depend on the final strain components and not on the strain path to reach
them.
The area ratio for single-direction rolling is only slightly larger than for the cross-
rolling when plotted against the compressive strain. This is not surprising given
that for the same rolling reduction, the length along the rolling direction for single-
direction rolling will be much larger than obtained by cross-rolling. In contrast, the
edge ratio becomes much larger for the single-direction scenario, reflecting its greater
microstructural anisotropy.
10.7 Summary
Solid state phase transformations often are limited by the fact that the diffusion of
atoms becomes difficult at low homologous temperatures. Nevertheless, there may
be a tendency for transformation to occur to a lower free energy state, but by a mech-
anism in which the structure of the parent phase is changed into that of the product
by a deformation which is driven not by external forces, but by the favourable change
in free energy. The deformation is real in the sense that the shape of the crystal will
change in order to reflect the difference between the atomic patterns of the parent
and product crystals.
Transformations of this kind are described generically as ‘displacive’ because the
structural change is accomplished by a coordinated shift of atoms [13]. The deforma-
168 10 Homogeneous deformations
tions described in this chapter are homogeneous; in the case illustrated in Fig. 10.9a,
the atoms themselves are placed in the correct final positions by the homogeneous
deformation of the lattice. However, the shuffle transformation in Fig. 10.9 cannot
be described using a homogeneous deformation matrix. The term shuffle refers here
to small displacements within the unit cell. There are cases, as we shall see in Chap-
ter 11 where the transformation is a mixture of lattice distortion and shuffles, so the
homogeneous deformation places some of the atoms in the correct position whereas
others have to move small distances so that overall structure of the product phase is
correct.
(a)
(b)
FIGURE 10.9
An illustration of a lattice distortive and shuffle transformation [13].
References
169
11
Invariant-plane strains
Abstract
11.1 Introduction
The deformation of crystals by the conservative glide of dislocations on a single set
of crystallographic planes causes shear in the direction of the resultant Burgers vector
of the dislocations concerned, a direction which lies in the slip plane. The slip plane
and slip direction constitute a slip system. The material in the slip plane remains
crystalline during slip and since there is no reconstruction of this material during slip
(e.g. by localised melting followed by resolidification), there can be no change in
the relative positions of atoms in the slip plane. The atomic arrangement on the slip
plane is thus completely unaffected by the deformation.
Another mode of conservative plastic deformation is mechanical twinning, in
which the parent lattice is homogeneously sheared into the twin orientation. The
plane on which the twinning shear occurs is again unaffected by the deformation and
forms a coherent boundary between the parent and its twin.
If a material which has a Poisson’s ratio equal to zero is uniaxially stressed below
its elastic limit, then the plane that is normal to the stress axis is unaffected by the
deformation since the only non-zero strain is that parallel to the stress axis. Polycrys-
talline beryllium has a Poisson’s ratio close to zero [1] so that stretching would lead
to little transverse strain.
171
172 11 Invariant-plane strains
FIGURE 11.1
Three kinds of invariant-plane strains. The cubes indicate the shape before deforma-
tion. δ, s and m represent the magnitudes of the dilatational strain, shear strain and
general displacement respectively. p is a unit vector, the shear strain s is parallel to
z1 , whereas δ is parallel to z3 .
planes for the hcp lattice. It has not been possible until recently to prove this, because
such a small region of hcp material gives very diffuse and unconvincing hcp reflec-
tions in electron diffraction experiments. However, the δ component of the fcc-hcp
martensite transformation strain has now been detected to be present for single stack-
ing faults, proving the hcp model of such faults [6]. In Fig. 11.2, the residual contrast
that is visible in part (b) in spite of the fact that the diffraction vector is normal to the
displacement in the plane of the fault, is indicative of a change in interplanar spacing,
i.e. a volume change normal to the fault plane, consistent with the presence of the
hcp phase. The contrast from a fault should otherwise be invisible in the absence of
a change in interplanar spacing.
Turning now to the description of the strains illustrated in Fig. 11.1, we follow
the procedure of Chapter 8, to find the matrices (Z P Z); the symbol ‘P’ in the matrix
representation is used to identify specifically, an invariant-plane strain, the symbol
‘S’ being the representation of any general deformation. Each column of such a ma-
trix represents the components of a new vector generated by deformation of a vector
equal to one of the basis vectors of Z. It follows that the three matrices representing
the deformations of Fig. 11.1a-c are, respectively,
⎛ ⎞ ⎛ ⎞
1 0 0 1 0 s
⎜ ⎟ ⎜ ⎟
(Z P1 Z) = ⎜ ⎝0 1 0 ⎠,
⎟ (Z P2 Z) = ⎜ ⎝0 1 0⎠ ,
⎟
0 0 1+δ 0 0 1
⎛ ⎞
1 0 s
⎜ ⎟
(Z P3 Z) = ⎝0
⎜
1 0 ⎠.
⎟
0 0 1+δ
174 11 Invariant-plane strains
(a) (b)
FIGURE 11.2
Stacking faults in austenitic stainless steel. (a) Bright field image of faults. (b) Image
with g.R equal to an integer for both of the faults illustrated, where g is the diffraction
vector and R is the displacement vector. After Brooks et al. [6], reproduced with
permission of Elsevier.
These have a particularly simple form, because the basis Z has been chosen care-
fully, such that p ∥ z3 and the direction of the shear is parallel to z1 . However, it
is often necessary to represent invariant-plane strains in a crystallographic basis, or
in some other basis X. This can be achieved with the similarity transformation law,
equation 11.2. If (X J Z) represents the coordinate transformation from the basis Z to
X, we have
(X P X) = (X J Z)(Z P Z)(Z J X).
Expansion of this equation gives [7]
⎛ ⎞
1 + md1 p1 md1 p2 md1 p3
⎜ ⎟
(X P X) = ⎜ ⎝ md2 p1 1 + md2 p2
⎟
md2 p3 ⎠ (11.3)
md3 p1 md3 p2 1 + md3 p3
where di are the components of d in the X basis, such that (d; X∗ )[X; d] = 1. The
vector d points in the direction of the displacements involved; a vector which is
parallel to d remains parallel following deformation, although the ratio of its final to
initial length may be changed. The quantities pi are the components of the invariant-
plane normal p, referred to the X∗ basis, normalised to satisfy (p; X∗ )[X; p] = 1.
Equation 11.3 may be simplified as follows:
(referred to the X basis) deform. In order to examine the way in which vectors which
are plane normals (i.e. referred to the reciprocal basis X∗ ) deform, we proceed in the
following manner.
The property of a homogeneous deformation is that points which originally are
collinear remain collinear after the deformation [8]. Lines that initially are coplanar
remain so following deformation. It follows that an initial direction u which lies in
a plane whose normal is initially h, becomes a new vector v within a plane whose
normal is k, where v and k result from the deformation of u and h, respectively. Now,
h.u = k.v = 0, so that:
i.e.
−1
(k; X∗ ) = (h; X∗ )(X P X) (11.5)
Equation 11.5 describes the way in which plane normals are affected by the defor-
mation (X P X). From equation 11.4, it can be shown that
−1
(X P X) = I − mq[X; d](p; X∗ ) (11.6)
011
112 p
011
100
u
441
d
(a) (b) (c) 111
176 11 Invariant-plane strains
FIGURE 11.3
Longitudinal section of the tensile specimen illustrating the (1 1 0) plane. All
directions refer to the crystal basis. The tensile axis rotates towards d, in the plane
containing the original direction of the tensile axis (u) and d.
Fig. 11.3a illustrates the deformation involved; the parent crystal basis α con-
sists of basis vectors which form the conventional bcc unit cell of α-iron. The
effect of the mechanical twinning is to alter the original shape abcd to a′ b′ c′ d′ .
ef is a trace of (1 1 2)α on which the twinning shear occurs in the [1 1 1]α di-
rection. However, as in most ordinary tensile tests, the ends of the specimen are
constrained to be vertically aligned at all times; a′ d′ must therefore be vertical and
the deforming crystal must rotate to comply with this requirement. The configura-
tion illustrated in Fig. 11.3c is thus obtained, with ad and a′ d′ parallel, the tensile
strain being (a′ d′ − ad)/(ad).
As discussed earlier, mechanical twinning is an invariant-plane strain; it in-
volves a homogeneous simple shear on the twinning plane, a plane which is
not affected by the deformation and which is common to both the parent and
twin regions. Equation 11.3 can be used to find the matrix (α P α) describ-
ing the mechanical twinning, given that the normal to the invariant-plane is
1
− 12
(p; α∗ ) = aα 6− 2 (1 1 2), the displacement direction is [α; d] = a−1 α 3 [1 1 1]
1
and m = 2− 2 . It should be noted that p and d respectively satisfy the conditions
(p; α∗ )[α; p] = 1 and (d; α∗ )[α; d] = 1, as required for equation 11.3. Hence
⎛ ⎞
7 1 2
1⎜ ⎟
(α P α) = ⎜ 1 7 2
⎟.
6 ⎝ ⎠
1 1 4
Using this, we see that an initial vector [α; u] = [4 4 1] becomes a new vector
[α; v] = (α P α)[α; u] = 61 [34 34 4] due to the deformation. The need to maintain
the specimen ends in alignment means that v is rotated to be parallel to the tensile
axis. Now, |u| = 5.745aα where aα is the lattice parameter of the ferrite, and
|v| = 8.042aα, giving the required tensile strain as (8.042−5.745)/5.745 = 0.40.
11.1.1 Comments
(i) From Fig. 11.3 it is evident that the end faces of the specimen will also undergo
deformation (ab to a′ b′ ) and if the specimen gripping mechanism imposes con-
straints on these ends, then the rod will tend to bend into the form of an ‘S’. For
thin specimens this effect may be small.
(ii) The tensile axis at the beginning of the experiment was [4 4 1], but at the end is
1
6 [34 34 4]. The tensile direction clearly has rotated during the course of the exper-
iment. The direction in which it has moved is 16 [34 34 4] − [4 4 1] = 16 [10 10 10],
parallel to [1 1 1], the shear direction d. In fact, any initial vector u will be dis-
11.1 Introduction 177
placed towards d to give a new vector v as a consequence of the IPS. Using equa-
tion 11.4, we see that
--
6 5 --- 101 --
b4 d4 111 111
a6 c5
4 ---
b2 c a d1 --- 213
a2 c1 112
d6 2 d4 b4 1 b5 a2 c1
010 -- d4 b4
c5 c3 a3 a6 011 011
b2 d1 c3 a3
a2 d
3
b
b2 d1
b5 d6 c3
5 6
--
111 111
(a) 100 (b) 101
FIGURE 11.4
Stereographic analysis of slip in fcc single-crystals. The dashed curve in (b) is the
trace of (111).
4m 4m
[γ; v] = [γ; u] + √ [0 1 1] = [2 1 3] + √ [0 1 1]
6 6
When duplex slip occurs, v must lie along the intersection of the (1 1 1) and (1 1 0)
planes, the former being the plane on which v is confined to move and the latter
being the boundary between triangles a2 and d4. It follows that v ∥ [1 1 2] and
must be of the form v = [v v 2v]. Substituting this into the earlier equation gives
4m
[v v 2v] = [2 1 3] + √ [0 1 1]
6
11.2 Deformation twins 179
and on comparing coefficients from both sides of this equation, we obtain
[γ; v] = [2 2 4]
(a) (b)
(c) (d)
FIGURE 11.5
Twinning in the fcc austenite crystal structure. The diagrams represent sections in
a
the (1 1 0) plane. The vector v is 6γ [112]. In (c), K2′ and η2′ are the final positions
of the undistorted plane K2 and the undistorted direction η2 , respectively. (d) A
180 11 Invariant-plane strains
close-packed {111} plane showing that a shear along v is not equivalent to one
parallel to −v, since the magnitude of the reverse displacement has to be twice as
large.
s2 = (u/d)2 − 4 (11.7)
From Fig. 11.5a,d we see that the twin lattice could also have been obtained by
displacing the atoms
√ √ in the +1 plane through
√ a distance 2v along [1 1 2] had u been
chosen to equal [ 2 2 0], giving s = 2. This larger shear is of course inconsis-
tent with the hypothesis that the favoured twinning mode involves the smallest shear,
and indeed, this mode of twinning is not observed. To obtain the smallest shear, the
magnitude of the vector v must also be minimised; in the example under considera-
tion, the correct choice of v will connect a lattice site of plane +1 with the projection
of its nearest neighbour lattice site on plane −1. The twinning direction is therefore
expected to be along [1 1 2]. It follows that the operative twin mode for the fcc lattice
1
should involve a shear of magnitude s = 2− 2 on {1 1 1}⟨1 1 2⟩.
The matrix-twin orientation relationship (M J T) can be deduced from the fact
that the twin was generated by a shear which brought atoms in the twin into positions
which are related to the parent lattice points by reflection across the twinning plane
(the basis vectors of M and T define the fcc unit cells of the matrix and twin crystals
respectively). From Fig. 11.5 we note that:
It follows that
⎛ ⎞ ⎛ ⎞⎛ ⎞
1 1 1 J11 J12 J13 1 1 1
⎜ ⎟ ⎜ ⎟⎜ ⎟
1 1⎠ = ⎝J21
⎜ ⎟ ⎜ ⎟⎜ ⎟
⎝1 J22 J23 ⎠ ⎝1 1 1⎠
2 0 1 J31 J32 J33 2 0 1
11.2.1 Comments
(i) Equations like equation 11.7 can be used to predict the likely ways in which dif-
ferent lattices might twin, especially when the determining factor is the magnitude
of the twinning shear [13].
(ii) There actually are four different ways of generating the twin lattice from the
parent crystal: (a) by reflection about the K1 plane on which the twinning shear
occurs, (b) by a rotation of π about η1 , the direction of the twinning shear, (c) by
reflection about the plane normal to η1 and (d) by a rotation of π about the normal
to the K1 plane.
Since most metals are centrosymmetric, operations (a) and (d) produce crystallo-
graphically equivalent results, as do (b) and (c). In the case of the fcc twin discussed
182 11 Invariant-plane strains
above, the high symmetry of the cubic lattice means that all four operations are crys-
tallographically equivalent. Twins which can be produced by the operations (a) and
(d) are called type I twins; type II twins result form the other two twinning opera-
tions. The twin discussed in the above example is called a compound twin, since type
I and type II twins cannot be crystallographically distinguished.
Fig. 11.5b illustrates some additional features of twinning. The K2 plane is the
plane which (like K1 ) is undistorted by the twinning shear, but unlike K1 , is rotated
by the shear. The ‘plane of shear’ is the plane containing η1 and the perpendicular
to K1 ; its intersection with K2 defines the undistorted but rotated direction η2 . In
general, η2 and K1 are rational for type I twins, and η1 and K2 are rational for
type II twins. The set of four twinning elements K1 , K2 , η1 and η2 are all rational
for compound twins. From Fig. 11.5b, η2 makes an angle of arctan(s/2) with the
normal to K1 and simple geometry shows that η2 = [1 1 2] for the fcc twin of
Example 11.2. The corresponding K2 plane which contains η2 and η1 ∧η2 is therefore
(1 1 1), giving the rational set of twinning elements
1
K1 = (1 1 1) η2 = [1 1 2] s = 2− 2 η1 = [1 1 2] K2 = (1 1 1)
In fact, it only is necessary to specify either K1 and η2 or K2 and η1 to completely
describe the twin mode concerned.
The deformation matrix (M P M) describing the twinning shear can be deduced
using equation 11.3 and the information that [M;d] ∥ [1 1 2], (p; M∗ ) ∥ (1 1 1) and
1
s = 2− 2 to give
⎛ ⎞ ⎛ ⎞
7 1 1 5 1 1
1⎜ ⎟ 1⎜ ⎟
(M P M) = ⎜ 1 7 1
⎟ and (M P M)−1 = ⎜ 1 5 1
⎟ (11.8)
6 ⎝ ⎠ 6 ⎝ ⎠
2 2 4 2 2 8
and if a vector u is deformed into a new vector v by the twinning shear, then
(M P M)[M;u] = [M;v] (11.9)
and if h is a plane normal which after deformation becomes k, then
−1
(h; M∗ )(M P M) = (k; M∗ ) (11.10)
These laws can be used to verify that p and d are unaffected by the twinning shear,
and that the magnitude of a vector originally along η2 is not changed by the defor-
mation; similarly, the spacing of the planes initially parallel to K2 remains the same
after deformation, although the planes are rotated.
which are initially coplanar remain coplanar after the strain. Using the data from
equation 11.8, we see that the deformation (M P M) alters the vector [M;u] = [0 0 1]
to a new vector [M;v] = 16 [1 1 4] i.e.
1
(M P M)[0 0 1]M = [1 1 4]M (11.11a)
6
The indices of this new vector v relative to the twin basis T can be obtained using the
coordinate transformation matrix (T J M), so that
1 1
(T J M) [1 1 4]M = [T; v] = [1 1 0]T (11.11b)
6 2
Hence, the effect of the shear stress is to deform a vector [0 0 1]M of the parent lattice
into a vector 12 [1 1 0]T of the twin. Equations 11.11a,b could have been combined to
obtain this result, as follows:
or
(T C M)[M;u] = [T; v]
where
(T J M)(M P M) = (T C M) (11.12)
The matrix (T C M) is called the correspondence matrix; the initial vector u
with indices [M;u] in the parent basis, due to deformation becomes a corresponding
vector v with indices [T;v] in the twin basis. The correspondence matrix tells us that a
certain vector in the twin is formed by deforming a particular corresponding vector of
the parent. In the example considered above, the vector u has rational components in
M (i.e. the components are small integers or fractions) and v has rational components
in T. It follows that the elements of the correspondence matrix (T C M) must also be
rational numbers or fractions. The correspondence matrix can usually be written from
inspection since its columns are rational lattice vectors referred to the second basis
produced by the deformation from the basis vectors of the first basis.
We can similarly relate planes in the twin to planes in the parent, the correspon-
dence matrix being given by
−1
(M C T) = (M P M) (M J T)
where
(h; M∗ )(M C T) = (k; T∗ )
so that the plane (k; T∗ ) of the twin was formed by the deformation of the plane
(h; M∗ ) of the parent.
184 11 Invariant-plane strains
60° [110]
[101]
[010]
[100]
FIGURE 11.6
An alternative correspondence between the austenite, conventionally represented by
the cubic unit cells illustrated, can also be represented by the body-centred triclinic
cell. The main diagram is adapted from [19]. Not all the iron atoms in the cubic cells
are illustrated for clarity. The directions oa, ob and oc are, respectively, parallel to
[101]γ , [110]γ and [013]γ .
11.4 An alternative to the Bain strain 185
a
The basal plane of the triclinic cell consists of vectors oa = 2γ [101] and ob =
aγ
2 [110], making an angle of 60◦ with each other. The normal to the basal plane is
a
therefore parallel to the [111] direction but the vector oc = 2γ [013]. A shear (γ S1 γ)
of s1 = tan 43◦ parallel to [415] on (111) would make oc normal to the basal plane
of the triclinic cell without affecting any vectors in the basal plane defined by the
triclinic cell edges oa and ob. Defining d = √142 [415] as the unit vector parallel
to the shear direction, and p = √13 [111] as the unit normal to the shear plane, and
defining φ = s1 √13 √142 , it follows that
(γ S1 γ) = I + s1 [γ; d](p; γ ∗ )
⎛ ⎞
1 + 4φ −4φ 4φ
⎜ ⎟
= ⎜
⎝ −φ 1+φ
⎟
−φ ⎠ (11.13)
−5φ 5φ 1 − 5φ
This deformation makes the original vector [013] rotated and distorted to
[1.33 1.33 1.33] ∥ [111]. It would need to be compressed by a ratio η3 =
aα /2.3076aγ so that its magnitude becomes the lattice parameter aα of martensite.
To change the angle 60◦ between oa and ob to 90◦ without affecting any vectors
normal to the basal plane requires a second shear (γ S2 γ) of magnitude s2 = tan 30◦
with d = √12 [101] and p = √16 [121]. Taking ψ = s2 √12 √16 :
(γ S2 γ) = I + s2 [γ; d](p; γ ∗ )
⎛ ⎞
1 − ψ −2ψ −ψ
⎜ ⎟
= ⎜
⎝ 0 1 0
⎟.
⎠ (11.14)
ψ 2ψ 1+ψ
This shear leaves [101] unaffected, whereas a vector originally parallel to [110] be-
comes [ 21 1 12 ] which is at 90◦ to [101]. To match the lattice parameter aα , the vector
1
√
2 [101] needs to expanded by η1 = 2aα /aγ whereas 12 [ 12 1 12 ] would need to be ex-
panded by the larger quantity η2 = 1.633aα/aγ because the shear (γ S2 γ) reduces
the length of 21 [110].
Therefore, the total deformation required to convert the triclinic austenite cell
into the body-centred cubic cell of martensite is
⎛ ⎞⎛ ⎞⎛ ⎞
1 + 4φ −4φ 4φ 1 − ψ −2ψ −ψ η1 0 0
⎜ ⎟⎜ ⎟⎜ ⎟
(γ S γ) = ⎜ ⎝ −φ 1 + φ −φ
⎟⎜
⎠ ⎝ 0 1 0
⎟⎜
⎠ ⎝ 0 η2 0
⎟.
⎠
−5φ 5φ 1 − 5φ ψ 2ψ 1+ψ 0 0 η3
Since shears do not lead to any change in volume, the ratio of the volumes of the
martensite and ferrite unit cells is given by the determinant of the distortion matrix,
i.e. η1 η2 η3 = a3α /a3γ as expected. These distortions, even though they are not mea-
sured along the principal axes, are large when compared with the ones associated
186 11 Invariant-plane strains
with the Bain strain. As pointed out by Wayman, it might be concluded on this basis
and and with the backing of other experimental evidence, that the Bain strain is the
homogeneous deformation that leads to a change from fcc→bcc structure.
bI = hmd
The operations outlined above indicate one way of generating the required stepped
interface. They are simply the virtual operations which allow us to produce the re-
quired defect - similar operations were first used by Volterra [23] in describing the
elastic properties of cut and deformed cylinders, operations which were later recog-
nised to generate the ordinary dislocations that metallurgists are so familiar with.
Having defined bI , we note that an initially planar coherent twin boundary can
acquire a step if a dislocation of Burgers vector bm crosses the interface. The height
of the step is given by [10]:
h = bm .p
so that
bI = m(bm .p)d. (11.15)
From equation 11.3, the invariant plane strain necessary to generate the twin from
the parent lattice is given by (M P M) = I + m[M;d](p; M∗ ) so that equation 11.15
becomes
[M;bI ] = (M P M)[M;bm ] − [M;bm ] (11.16)
11.5 Stepped interfaces 187
(a) (c)
FIGURE 11.7
The virtual operations [24] used in determining bI . (a) Single crystal. (b) Crystal cut
along ‘xy’. (c) Gap opens on shearing. The closure failure bI is the Burgers vector
of the final step in the coherent twin boundary, that results on rewelding the gap as
shown in (d).
188 11 Invariant-plane strains
(T C M) = (T J M)(M P M)
(M C T) = (M P M)−1 (M J T)
so that ⎛ ⎞
1 1 1
1⎜ ⎟
(T C M) = (M C T) = ⎜ 1 1⎠ .
⎟
2 ⎝1
2 2 0
The character of a dislocation will in general be altered on crossing an in-
terface. This is because the crossing process introduces a step in the interface,
rather like the slip steps which arise at the free surfaces of deformed crystals. We
consider the case where a dislocation crosses a coherent twin boundary. The in-
terfacial step has dislocation character so that the original dislocation (Burgers
vector bm ) from the parent crystal is in effect converted into two dislocations,
one being the step (Burgers vector bI ) and the other the dislocation (Burgers vec-
tor bt ) which has penetrated the interface and entered into the twin lattice. If the
total Burgers vector content of the system is to be preserved then it follows that
in general, bt ̸= bm , since bm = b+ bt . Using this equation and equation 11.16,
we see that
[M;bt ] = (M P M)[M;bm ]
or
[T; bt ] = (T J M)(M P M)[M;bm ]
so that
[T; bt ] = (T C M)[M;bm ].
Clearly, dislocation glide across the coherent interface will not be hindered if
bt is a perfect lattice vector of the twin. If this is not the case and bt is a partial
dislocation in the twin, then glide across the interface will be hindered because
the motion of bt in the twin would leave a stacking fault trailing all the way from
the interface to the position of the partial dislocation in the twin.
11.5 Stepped interfaces 189
There is an additional condition to be fulfilled for easy glide across the inter-
face; the corresponding glide planes pm and pt of dislocations bm and bt in
the parent and twin lattices respectively, must meet edge to edge in the interface.
Now,
(pt ; T∗ ) = (pm ; M∗ )(M C T).
If the interface plane normal is pi , then the edge to edge condition is satisfied
if pm ∧ pi ∥ pt ∧ pi .
Dislocations in fcc materials usually glide on close-packed {1 1 1} planes
and have Burgers vectors of type a2 ⟨1 1 0⟩. Using the data of Table 11.1 it can
easily be verified that all the close-packed planes of the parent lattice meet the
corresponding glide planes in the twin edge to edge in the interface, which is
taken to be the coherent (1 1 1)M twinning plane. Furthermore, all the a2 ⟨1 1 0⟩
Burgers vectors of glide dislocations in the parent correspond to perfect lattice
dislocations in the twin. It must be concluded that the coherent twin boundary
for {1 1 1} twins in fcc metals does not offer any geometrical restrictions to the
transfer of slip between the parent and product lattices.
TABLE 11.1
Corresponding glide planes and Burgers vectors
Parent Twin
a
2 [1 1 0] a[0 0 1]
a a
2 [1 0 1] 2 [0 1 1]
a a
2 [0 1 1] 2 [1 0 1]
a a
2 [1 1 0] 2 [1 1 0]
a a
2 [1 0 1] 2 [1 0 1]
a a
2 [0 1 1] 2 [0 1 1]
(1 1 1) (1 1 1)
(1 1 1) (1 1 1)
(1 1 1) (0 2 0)
(1 1 1) (2 0 0)
These data (Table 11.1) show also that all dislocations with Burgers vectors in
the (1 1 1)M plane are unaffected, both in magnitude and direction, as a result of
crossing into the twin. For example, a2 [1 1 0]M becomes a2 [1 1 0]T so that |bm | =
|bt |, and using (T J M) it can be demonstrated that [1 1 0]M ∥ [1 1 0]T . This result
is expected because these particular dislocations cannot generate a step in the
(1 1 1)M interface when they cross into the twin lattice (see equation 11.15). Only
dislocations with Burgers vectors not parallel to the interface cause the formation
of steps.
190 11 Invariant-plane strains
The data further illustrate the fact that when bm lies in the (1 1 1)M plane,
there is no increase in energy due to the reaction bm → bI + bt , which occurs
when a dislocation crosses the interface. This is because bI = 0 and bt = bm . For
all other cases bI is not zero and since |bt | is never less than |bm | , bm → bI + bt
is always energetically unfavourable. In fact, in the example being discussed, there
can never be an energy reduction when an a2 ⟨1 1 0⟩ dislocation penetrates the co-
herent twin boundary. The dislocations cannot therefore spontaneously cross the
boundary. A trivial case where dislocations might spontaneously cross a boundary
is when the latter is a free surface, assuming that the increase in surface area (and
hence surface energy) due to the formation of a step is not prohibitive. Sponta-
neous penetration of the interface might also become favourable if the interface
separates crystals with very different elastic properties.
FIGURE 11.8
The passage of a slip dislocation across a coherent twin boundary in a bcc crystal.
1
The twinning system is {1 1 2}⟨1 1 1⟩, s = 2− 2 . The subscripts m and t refer
11.5 Stepped interfaces 191
to the twin and matrix respectively; the arrows at the dislocation symbols indicate
the sense of the Burgers vectors and the dislocation line vectors are all parallel to
[1 1 0]m,t where m refers to the matrix and t to the twin.
The results obtained show that single dislocations can glide into twins in fcc
crystals without leaving a fault; there are no geometrical restrictions to the passage
of slip dislocations across the coherent twin boundaries concerned. It can similarly
be demonstrated that slip dislocations can comfortably traverse the coherent twin
boundaries of {1 1 2} twins in bcc or bct lattices and this has implications on the
interpretation of the strength of martensite [22]. As will be discussed later, the sub-
structure of martensite plates in steels (and in many non-ferrous alloys) often consists
of very finely spaced {1 1 2} transformation twins. It was at one time believed that
the twins were mainly responsible for the high strength of ferrous martensites, be-
cause the numerous twin boundaries should hinder slip – the analysis above clearly
suggests otherwise. Indeed, twinned martensites which do not contain carbon also
do not exhibit exceptionally high strengths and it is now generally accepted that
the strength of virgin ferrous martensites is largely due to interstitial solid solution
hardening by carbon atoms, or in the case of lightly autotempered martensites due
to carbon atom clustering or fine precipitation. Consistent with this, it is found that
Fe-30Ni wt% twinned martensites are not particularly hard.
Finally, it should be mentioned that even when glide across coherent twin bound-
aries in martensites should be unhindered, the boundaries will cause a small amount
of hardening, partly because the corresponding slip systems in the matrix and twin
will in general be differently stressed [22, 25] (simply because they are not necessar-
ily parallel) and partly due to the work necessary to create the steps in the interfaces.
It is emphasised, however, that these should be relatively small contributions to the
strength of martensite. Fig. 11.8 illustrates the passage of a slip dislocation across a
coherent {1 1 2} twin interface in a bcc material.
This can be factorised into a pure strain (Z Q Z) and a rigid body rotation (Z J Z)
(Fig. 10.5). Writing (Z T Z) = (Z P′ Z)(Z P Z), we obtain:
⎛ ⎞
174 30 −6
1 ⎜ ⎟
(Z T Z) = ⎜
30 174 −6
⎟.
144 ⎝ ⎠
−6 −6 102
Notice that the eigenvectors form an orthogonal set and that consistent with the
fact that v lies in the invariant plane, λ2 has a value of unity. u, v and x are also the
eigenvectors of (Z Q Z). The eigenvalues of (Z Q Z) are given by the square roots
of the eigenvalues of (Z T Z); they are 1.192282, 1.0 and 0.838728. Hence, the
maximum extensions and contractions are less than 20% since each eigenvalue is
the ratio of the final to initial length of a vector parallel to an eigenvector. The
maximum extension occurs along u and the maximum contraction along x.
and ⎛ ⎞
0.920562
−0.079438 0.020515
−1 ⎜ ⎟
(Z Q Z) =⎜
⎝−0.079438 0.920562 0.020515⎠
⎟
where the new vector can be shown to have the same magnitude as [1 1 2] but points
in a different direction. The effect of the pure rotation is
Thus, the pure strain deforms [1 1 2]Z into another vector of identical magnitude
and the pure rotation brings this new vector back into the [1 1 2]Z direction, the
net operation leaving it invariant, as expected, since [1 1 2]Z is the shear direction
which lies in the invariant plane. Referring to Fig. 10.5, the direction f e = [1 1 0]Z ,
yz = [1 1 2]Z and y ′ z ′ = [1.230916 1.230916 − 1.723280]Z. y ′ z ′ is brought into
coincidence with yz by the rigid body rotation (Z J Z) to generate the invariant plane
containing f e and yz.
Physically, the fcc to hcp transformation occurs by the movement of a single set
of Shockley partial dislocations, Burgers vector b = a6 ⟨1 1 2⟩γ on alternate close-
packed {1 1 1}γ planes. To produce a fair thickness of hcp martensite, a mechanism
has to be sought which allows Shockley partials to be generated on every other slip
plane. Some kind of a pole mechanism (see for example, p. 310 of [27]) would allow
this to happen, but there is as yet no experimental evidence confirming this. Motion
of the partials would cause a shearing of the γ lattice, on the system {1 1 1}γ ⟨1 1 2⟩γ ,
the average magnitude s of the shear being s = |b|/2d, where d is the spacing of the
1 1 1
close-packed planes. Hence, s = 6− 2 a/2(3− 2 a) = 8− 2 . This is exactly the shear
system we used in generating the matrix (Z P Z) and the physical effect of the shear
on the shape of an originally flat surface is, in general, to tilt the surface (about a
line given by its intersection with the hcp habit plane) through some angle dependant
on the indices of the free surface. By measuring such tilts it is possible to deduce s,
which has been experimentally confirmed to equal half the twinning shear.
In fcc crystals, the close-packed planes stack in the sequence . . . ABCABC . . .;
the passage of a single Shockley partial causes the sequence to change to . . . ABA . . .
194 11 Invariant-plane strains
creating a three layer thick region of hcp phase since the stacking sequence of close-
packed planes in the hcp lattice has a periodicity of 2. This then is the physical
manner in which the transformation occurs, the martensite having a {1 1 1}γ habit
plane – if the parent product interface deviates slightly from {1 1 1}γ , then it will
consist of stepped sections of close-packed plane, the steps representing the Shockley
partial transformation dislocations. The spacing of the partials along ⟨1 1 1⟩γ would
be 2d. In other words, in the stacking sequence ABC, the motion of a partial on B
1
would leave A and B unaffected though C would be displaced by 6− 2 a⟨1 1 2⟩γ
to a new position A, giving ABA stacking. Partials could thus be located on every
alternate plane of the fcc crystal.
Hence, we see that the matrix (Z P Z) is quite compatible with the microscopic
dislocation based mechanism of transformation. (Z P Z) predicts the correct macro-
scopic surface relief effect and its invariant plane is the habit plane of the martensite.
However, if (Z P Z) is considered to act homogeneously over the entire crystal, then
it would carry half the atoms into the wrong positions. For instance, if the habit
plane is designated A in the sequence ABC of close packed planes, then the effect
of (Z P Z) is to leave A unchanged, shift the atoms on plane C by 2sd and those on
plane B by sd along ⟨1 1 2⟩γ . Of course, this puts the atoms originally in C sites into
A sites, as required for hcp stacking. However, the B atoms are located at positions
a
half way between B and C sites, through a distance 12 ⟨1 1 2⟩γ . Shuffles are thus
necessary to bring these atoms back into the original B positions and to restore the
. . . ABA . . . hcp sequence. These atomic movements in the middle layer are called
shuffles because they occur through very small distances (always less than the inter-
atomic spacing) and do not affect the macroscopic shape change [21]. The shuffle
here is a purely formal concept; consistent with the fact that the Shockley partials
glide over alternate close-packed planes, the deformation (Z P Z) must in fact be
considered homogeneous only on a scale of every two planes. By locking the close-
packed planes together in pairs, we avoid displacing the B-site atoms to the wrong
positions and thus automatically avoid the reverse shuffle displacement.
In the particular example discussed above, the dislocation mechanism is estab-
lished experimentally and physically reasonable shear systems were used in deter-
mining (Z P Z). However, in general it is possible to find an infinite number of defor-
mations [8, 21] which may accomplish the same lattice change and slightly empirical
criteria have to be used in selecting the correct deformation. One such criterion could
involve the selection of deformations which involve the minimum principal strains
and the minimum degree of shuffling, but intuition and experimental evidence is al-
most always necessary to reach a decision. The Bain strain which transforms the fcc
lattice to the bcc lattice is believed to be the correct choice because it seems to in-
volve the least atomic displacements and zero shuffling of atoms [20]. The absence
of shuffles can be deduced from the Bain correspondence matrix (α C γ) which can
be deduced from inspection since its columns are rational lattice vectors referred to
the α basis, produced by the deformation of the basis vectors of the γ basis; since
[1 0 0]γ is deformed to [1 1 0]α , [0 1 0]γ to [1 1 0]α and [0 0 1]γ to [0 0 1]α , by the
11.5 Stepped interfaces 195
If u is a vector defining the position of an atom in the γ unit cell, then it can be
verified that (α C γ)[γ; u] always gives a corresponding vector in the α lattice which
terminates at a lattice point. For example, 12 [1 0 1]γ corresponds to 21 [1 1 1]α ; both
these vectors connect the origins of their respective unit cells to an atomic position.
The Bain correspondence thus defines the position of each and every atom in the α
lattice relative to the γ lattice. It only is possible to obtain a correspondence matrix
like this when the primitive cells of each of the lattices concerned contain just one
atom [8].
The primitive cell of the hcp lattice contains two atoms and any lattice corre-
spondence will only define the final positions of an integral fraction of the atoms, the
remainder having to shuffle into their correct positions in the product lattice. This can
be demonstrated with the correspondence matrix for the example presented above.
It is convenient to represent the conventional hcp lattice (basis H) in an alternative
orthorhombic basis (symbol O), with basis vectors:
1 1
[1 0 0]O = [0 1 1]γ = [1 1 2 0]H
2 3
1
[0 1 0]O = [2 1 1]γ = [1 1 0 0]H
2
2
[0 0 1]O = [1 1 1]γ = [0 0 0 1]H
3
The orthorhombic unit cell thus contains three close-packed layers of atoms par-
allel to its (0 0 1) faces. The middle layer has atoms located at [0 31 21 ]O , [1 13 21 ]O
and [ 12 65 12 ]O . The other two layers have atoms located at each corner of the unit cell
and in the middle of each (0 0 1) face, as illustrated in Fig. 11.9.
From our earlier definition of a correspondence matrix, (O C γ) can be written
directly from the relations (between basis vectors) stated earlier:
⎛ ⎞
0 2 1
1⎜ ⎟
(γ C O) = ⎜ −1 −1 1
⎟
2 ⎝ ⎠
1 −1 2
FIGURE 11.9
Representation of bases O, H, and γ. The directions in the hexagonal cell are ex-
pressed in the Weber notation.
vectors of the orthonormal basis Z are parallel to the corresponding basis vectors of
the orthogonal basis γ. It follows that:
⎛ ⎞⎛ ⎞
0 −1 1 13/12 1/12 1/12
⎜ ⎟⎜ ⎟
(O C γ) = ⎜ ⎝2/3 −1/3 −1/3⎠ ⎝ 1/12
⎟⎜
13/12 1/12 ⎠
⎟
We have found that an ILS has two important characteristics: it leaves a line u
invariant and also leaves a plane normal h invariant. If the ILS is factorised into
two IPS’s, then u lies at the intersection of the invariant-planes of these component
IPS’s, and h defines the plane containing the two displacement vectors of these IPS’s.
These results will be useful in gaining a deeper understanding of martensitic trans-
formations.
11.7 Summary
In mathematics, the term invariance generally refers to a property that is unchanged
by a transformation, which in the present context can be a deformation. The deforma-
tions considered here leave planes or lines invariant and have physical consequences
that can be detected experimentally. Ordinary slip in solids leaves traces at free sur-
faces, that for crystalline materials reflect the anisotropy of atomic arrangements but
there is no change in volume. In contrast, solid-state phase transformations are usu-
ally associated with changes in density, but the displacements associated with the
volume change can be directed normal to the broad interface between the parent and
product crystal, so that an invariant plane is preserved. An important result from
the methods described in this chapter is that there is unique connection between the
orientation relationship, invariant plane (habit plane) and the shape deformation.
References
201
202 11 References
Abstract
Martensitic transformations can be quite simple to understand in systems where the
transformation strain is equal to the observed shape deformation. This is the case
when hcp martensite grows in fcc austenite. The odd observation that requires ex-
planation is that during the fcc→bcc or bct transformation, the shape deformation
is inconsistent with the strain needed to change the change the lattice. The crystallo-
graphic theory described in this Chapter deals with this anomaly using a most elegant
formalism that is quantitative.
12.1 Introduction
The basic principles of the characteristics of martensitic transformation have already
been covered in Chapter 8, albeit qualitatively. The purpose here is to explain the
mathematical framework developed by Bowles, MacKenzie, Wechsler, Lieberman
and Read [1–4], a framework that represents one of the most elegant and complete
theories that helped explain an accumulation of findings that could not be reconciled.
And there were predictions made which were only to be verified when advanced in-
strumentation such as transmission electron microscopy became available. The his-
tory of this wonderful theory is described by Wayman [5]; a classic example of a
subject receptive to a new way of thinking when a discipline is in a state of crisis
[6, 7]. The theory has had many consequences, implicitly or directly, in unrelated
subjects such as in the estimation of interfacial structure [8–11], the calculation of
stress-induced crystallographic texture [12–14] and the general understanding of the
relative orientations of crystals that precipitate in the solid state.
203
204 12 Martensite
must be a visible change in the macroscopic shape of the parent crystal during trans-
formation [15]. The shape deformation and its significance can best be illustrated
by reference to Fig. 12.1, where a comparison is made between reconstructive and
diffusionless transformations. For simplicity, the diagram refers to a case where the
transformation strain is an invariant-plane strain and a fully coherent interface exists
between the parent and product lattices, irrespective of the mechanism of transfor-
mation.
Considering the shear transformation first, we note that since the pattern of
atomic arrangement is changed on transformation, and since the transformation is
diffusionless, the macroscopic shape of the crystal changes. The shape deformation
has the exact characteristics of an IPS. The initially flat surface normal to da becomes
tilted about the line formed by the intersection of the interface plane with the surface
normal to da. The straight line ab is bent into two connected and straight segments
ae and eb. Hence, an observer looking at a scratch that is initially along ab and in
the surface abcd would note that on martensite formation, the scratch becomes ho-
mogeneously deflected about the point e where it intersects the trace of the interface
plane. Furthermore, the scratches ae and eb would be seen to remain connected at
the point e. This amounts to proof that the shape deformation has, on a macroscopic
scale, the characteristics of an IPS and that the interface between the parent and prod-
uct lattices does not contain any distortions (i.e. it is an invariant-plane). Observing
the deflection of scratches is one way of deducing the nature of shape deformations
accompanying transformations.
In Fig. 12.1 it also is implied that martensitic transformation is diffusionless;
labelled rows of atoms in the parent crystal are expected to remain in the correct
sequence in the martensite lattice. It is possible therefore to suggest that a particular
atom in the martensite must have originated from a corresponding particular atom in
the parent crystal. A formal way of expressing this property is to say that there exists
an atomic correspondence between the parent and product lattices.
In the case of the reconstructive transformation illustrated in Fig. 12.1, it is ev-
ident that the product phase can be of a different composition from the parent. In
addition, there has been much mixing up of atoms during transformation and the or-
der of arrangement of atoms in the product lattice is different from that in the parent
lattice - the atomic correspondence has been destroyed. Because the transformation
involves a reconstruction of the parent lattice, atoms are able to diffuse around in
such a way that the IPS shape deformation and its accompanying strain energy, do
not arise. The scratch ab remains straight across the interface and is unaffected by
the transformation.
In summary, martensitic transformations always are accompanied by a change
in the shape of the parent crystal, characterised as an invariant-plane strain when
examined on a macroscopic scale. The occurrence of such a shape deformation is
taken to imply the existence of an atomic correspondence between the parent and
product lattices. It is possible to state that a particular atom in the product occupied
a particular corresponding site in the parent lattice.
These results have some interesting consequences. The formation of martensite
in a constrained environment causes distortions in its surroundings. The strain energy
12.2 Shape deformation 205
FIGURE 12.1
Schematic illustration of the mechanisms of diffusional and shear transformations.
elastic accommodation should be taken to be the upper limit of the stored energy
due to the shape change accompanying martensitic transformation. This is because
the plastic accommodation is driven by the shape deformation [18]. In the event that
plastic accommodation occurs, dislocations and other defects may be generated both
in the parent and product lattices, which would make it impossible to reverse the
motion of the interface, thus eliminating the possibility of a shape memory effect.
Martensitic transformation does not require diffusion, can occur at very low tem-
peratures and the interface can move rapidly. The interface between martensite and
the parent phase must therefore be glissile, so we begin with a discussion of what the
appropriate structure of such an interface might be.
activated climb of intrinsic dislocations, so that the interface can only move in a non-
conservative manner, with relatively restricted or zero mobility at low temperatures.
A martensite interface cannot therefore be epitaxially semi-coherent.
In the second type of semi-coherency, the discontinuities discussed above are
screw dislocations, or dislocations whose Burgers vectors do not lie in the interface
plane. This kind of semi-coherency is of the type associated with glissile martensite
interfaces, whose motion is conservative (i.e. the motion does not lead to the cre-
ation or destruction of lattice sites). Such an interface should have a high mobility
since the migration of atoms is not necessary for its movement. Actually, two further
conditions must be satisfied before even this interface can be said to be glissile:
(i) A glissile interface requires that the glide planes of the intrinsic dislocations
associated with the product lattice must meet the corresponding glide planes of
the parent lattice edge to edge in the interface [24], along the dislocation lines.
(ii) If more than one set of intrinsic dislocations exist, then these should either have
the same line vector in the interface, or their respective Burgers vectors must be
parallel [24]. This condition ensures that the interface can move as an integral unit.
It implies in addition, that the deformation caused by the intrinsic dislocations,
when the interface moves, can always be described as a simple shear (caused by
a resultant intrinsic dislocation which is a combination of all the intrinsic dislo-
cations) on some plane which makes a finite angle with the interface plane, and
intersects the latter along the line vector of the resultant intrinsic dislocation.
martensite
ined
martensite
constrained
FIGURE 12.2
The habit plane of martensite (α′ ) under conditions of unconstrained and constrained
transformation, respectively. In the latter case, the dashed line indicates the trace of
the habit plane.
FIGURE 12.3
Schematic illustration of the phenomenological theory of martensite. (a) represents
the austenite crystal and (c) has a bcc structure. (b) has a structure between fcc and
bcc, p is the habit plane unit normal and q is the unit normal to the plane on which
the lattice-invariant shear occurs. The heavy horizontal lines in (e) are coherent twin
boundaries. Note that the vector e is normal to q but does not lie in the plane of the
diagram.
This problem is illustrated in Figs. 12.3. The shape of the starting austenite crystal
with the fcc structure is shown in Fig. 12.3a. On martensitic transformation its shape
alters to that illustrated in Fig. 12.3b and the shape deformation on going from (a)
to (b) is clearly an IPS on the plane with unit normal p and in the unit displacement
direction d. However, the structure of the crystal in Fig. 12.3b is some intermediate
lattice which is not bcc, since an IPS cannot on its own change the fcc structure to
the bcc structure. An invariant-line strain can however transform fcc to bcc, and since
an ILS can be factorised into two invariant-plane strains, it follows that the further
deformation (F Q F) needed to change the intermediate structure of Fig. 12.3b to the
bcc structure (Fig. 12.3c) is another IPS. If the deformation (F Q F) is of magnitude
210 12 Martensite
(F Q F) = I + n[F;e](q; F∗ )
is not affected by these deformations. Furthermore, the invariant normal of the ILS
must define a plane which contains the displacement directions of the lattice invari-
ant shear and of the shape deformation. This ensures that the spacing of this plane is
not affected by (F Q F) or (F P F).
The next example illustrates how the transformation strain can be determined
once the pure deformation which accomplishes the lattice change is deduced. To
ensure that the invariant-line and invariant-normal of the transformation strain are
compatible with the mode of lattice-invariant shear, we first need to specify the lat-
ter. It is assumed here that the plane and direction of the lattice-invariant shear are
(1 0 1)F and [1 0 1]F respectively. One variant of the Bain strain is illustrated in
Fig. 8.2, where we see that [1 0 0]γ is deformed into [1 1 0]α , [0 1 0]γ to [1 1 0]α and
[0 0 1]γ to [0 0 1]α , so that the variant of the Bain correspondence matrix is given by
equation 11.19. We will use this variant of the Bain correspondence matrix through-
out the text, but we note that there are two other possibilities, where [0 0 1]α can be
derived from either [1 0 0]γ or [0 1 0]γ respectively.
0 0 η3
u1 = −u3 (12.3)
Prior to deformation,
|u|2 = (u; F)[F;u] = 1 (23b)
u, as a result of deformation becomes a new vector x with
u21 + u22 + u23 = η12 u21 + η22 u22 + η32 u23 . (12.4)
Equations 12.3 and 12.4 can be solved simultaneously to give two solutions for
undistorted lines:
To solve for the invariant normal of the ILS, we proceed as follows. Writing
(h; F∗ ) = (h1 h2 h3 ), we note that for h to contain [1 0 1]F , its components must
satisfy the equation
h1 = h3 (12.5)
and furthermore,
(h; F∗ )[F∗ ; h] = 1.
h, on deformation becomes a new plane normal l and if |h| = |l| then
so that
h21 + h22 + h23 = (l1 /η1 )2 + (l2 /η2 )2 + (l3 /η3 )2 (12.6)
Solving equations 12.5 and 12.6 simultaneously, we obtain the two possible solu-
tions for the undistorted-normals as
Hence, one way of converting (F B F) into an ILS is to employ a rigid body ro-
tation which simultaneously rotates l into h and x into u. We have found that there
are two undistorted lines and two undistorted normals which satisfy the conditions
of the original question, so there are four ways of choosing pairs of undistorted
214 12 Martensite
lines and undistorted normals, but in this case, the four solutions are crystallo-
graphically equivalent. There are therefore four solutions (different in general) to
the problem of converting (F B F) to (F S F), subject to the condition that the
invariant-line should be in (1 0 1) and that the invariant normal defines a plane
containing [1 0 1]. We will concentrate on the solution obtained using the pair u
and h:
l = (h; F∗ )(F B F)−1 = (0.474554 0.569558 0.671120)
x = (F B F)[F;u] = [−0.762440 − 0.357809 0.539127]
a = u ∧ h = (−0.604053 0.723638 − 0.264454)
b = x ∧ l = (−0.547197 0.767534 − 0.264454)
The required rigid body rotation should rotate x back to u, l back to h and b to a,
giving the three equations:
[F;u] = (F J F)[F;x]
[F;h] = (F J F)[F;l]
[F;a] = (F J F)[F;b]
which can be expressed as a 3 × 3 matrix equation
⎛ ⎞ ⎛ ⎞⎛ ⎞
u 1 h1 a1 J11 J12 J13 w l1 b1
⎜ ⎟ ⎜ ⎟⎜ 1 ⎟
⎜ ⎟=⎜ ⎟⎜ ⎟
⎝u2 h2 a2 ⎠ ⎝J21 J22 J23 ⎠ ⎝w2 l2 b2 ⎠
u3 h3 a3 J31 J32 J33 w3 l3 b3
it follows that
⎛ ⎞ ⎛ ⎞
−0.671120 0.539127 −0.604053 −0.762440 0.474554 −0.547197
⎜ ⎟ ⎜ ⎟
⎝−0.314952 0.647058 0.723638 ⎠ = (F J F) ⎝−0.357808 0.569558 0.767534 ⎠
⎜ ⎟ ⎜ ⎟
[α; u] = (α J γ)[γ; u]
[γ; u] = (γ J α)[α; u]
(h; α∗ ) = (h; γ ∗ )(γ J α)
(h; γ ∗ ) = (h; α∗ )(α J γ).
(α J γ)(γ S γ) = (α C γ)
and
[1 0 1]γ = [−0.932679 − 1.049889 1.061622]α
which means that (1 1 1)γ is very nearly parallel to (0 1 1)α and [1 0 1]γ is about
3◦ from [1 1 1]α . The orientation relationship is illustrated in Fig. 12.4.
12.7 Stage 3: Nature of the shape deformation 217
FIGURE 12.4
Stereographic representation of the orientation relationship between martensite and
austenite, as deduced in equation 12.8. The lattice-invariant shear plane (q) and direc-
tion (-e), and the habit plane (p) and unit displacement vector (d) also are illustrated.
definite link between the initial and final states without being certain of the path in
between.
Example 12.3: Habit plane and the shape deformation
For the martensite reaction considered on page 212, determine the habit plane
of the martensite plate, assuming that the lattice-invariant shear occurs on the sys-
tem (1 0 1)γ [1 0 1]γ . Comment on the choice of this shear system and determine
the nature of the shape deformation.
The lattice invariant shear is on (1 0 1)[1 0 1] and since its effect is to cancel
the shape change due to (F Q F), the latter must be a shear on (1 0 1)[1 0 1].
To solve for the habit plane (unit normal p) it is necessary to factorise (F S F)
into the two invariant-plane strains (F P F) = I + m[F;d](p; F∗ ) and (F Q F) =
I + n[F;e](q; F∗ ).
The transformation strain (F S F) of equation 12.7 was calculated by phe-
nomenologically combining the Bain strain with a rigid body rotation, with the
latter chosen to make (F S F) an invariant line strain, subject to the condition that
the invariant line u of (F S F) must lie in (1 0 1)γ and that the invariant-normal
h of (F S F) must define a plane containing [1 0 1]γ . This is of course compat-
ible with the lattice-invariant shear system chosen in the present example since
u.q = e.h = 0. From equation 12.2,
= (q; F∗ ) − b(p; F∗ )
(p; F∗ ) of course represents the indices of the habit plane of the martensite plate.
As expected, the habit plane is irrational. To completely determine the shape de-
formation matrix (F P F) we also need to know m and d. Using equation 12.9, we
see that
(F S F)[F;e] = [F;e] + m[F;d](p; F∗ )[F;e]
Writing c as the scalar constant c = (p; F∗ )[F;e], we get
and
m = |md| = 0.223435
The magnitude m of the displacements involved can be factorised into a shear
component s parallel to the habit plane and a dilatational component δ normal to
1
the habit plane. Hence, δ = md.p = 0.0368161 and s = (m2 −δ 2 ) 2 = 0.220381.
These are typical values of the dilatational and shear components of the shape
strain found in ferrous martensites.
220 12 Martensite
Using these data, the shape deformation matrix is given by
⎛ ⎞
0.990134 −0.039875 −0.028579
⎜ ⎟
(F P F) = ⎜ ⎝ 0.032037 1.129478 0.092800 ⎠ .
⎟ (12.13)
−0.028583 −0.115519 0.917205
Of course, K must be an integral number, and the non-integral result must be taken
to mean that there will on average be a dislocation located on every 3.7679th slip
plane; in reality, the dislocations will be non-uniformly placed, either 3 or 4 (1 0 1)
planes apart.
The line vector of the dislocations is the invariant-line u and the spacing of the
intrinsic dislocations, as measured on the habit plane is Kd/(u ∧ p.q) where all
the vectors are unit vectors. Hence, the average spacing would be
1
3.7679(aγ 2− 2 )/0.8395675 = 3.1734aγ
It is important to note that the twin plane in the martensite corresponds to a mirror
plane in the austenite; this is a necessary condition when the lattice-invariant shear
involves twinning. The condition arises because the twinned and untwinned regions
of the martensite must undergo Bain strain along different though crystallographi-
cally equivalent principal axes [1, 28].
The above theory clearly predicts a certain volume fraction of twins in each
martensite plate, when the lattice-invariant shear is twinning as opposed to slip. How-
222 12 Martensite
ever, the factors governing the spacing of the twins are less quantitatively established;
the finer the spacing of the twins, the lower will be the strain energy associated with
the matching of each twin variant with the parent lattice at the interface. On the other
hand, the amount of coherent twin boundary within the martensite increases as the
spacing of the twins decreases.
A factor to bear in mind is that the lattice-invariant shear is an integral part of
the transformation; it does not happen as a separate event after the lattice change has
occurred. The transformation and the lattice-invariant shear all occur simultaneously
at the interface, as the latter migrates. It is well known that in ordinary plastic defor-
mation, twinning rather than slip tends to be the favoured deformation mode at low
temperatures or when high strain rates are involved. It is therefore often suggested
that martensite with low MS temperatures will tend to be twinned rather than slipped,
but this cannot be formally justified because the lattice-invariant shear is an integral
part of the transformation and not a physical deformation mode on its own. Indeed, it
is possible to find lattice-invariant deformation modes in martensite which do not oc-
cur in ordinary plastic deformation experiments. The reasons why some martensites
are internally twinned and others slipped are not clearly understood [29]. When the
spacing of the transformation twins is roughly comparable to that of the dislocations
in slipped martensite, the interface energies are roughly equal. The interface energy
increases with twin thickness and at the observed thicknesses is very large compared
with the corresponding interface in slipped martensite. The combination of the rel-
atively large interface energy and the twin boundaries left in the martensite plate
means that internally twinned martensite is never thermodynamically favoured rela-
tive to slipped martensite. It is possible that kinetic factors such as interface mobility
actually determine the type of martensite that occurs.
sample of material. However, circumstances may force variant selection, i.e. where
certain crystals which are more compliant to an external influence, such as an exter-
nally applied stress.
It is the interaction of the applied stress with the shape deformation of individual
martensite plates P (page 217 which determines variant selection. The interaction
energy which provides the mechanical driving force for transformation is given by
[32]:
U = σN ζ + τ s (12.14)
where σN is the stress component normal to the habit plane, τ is the shear stress
resolved on the habit plane in the direction of shear and ζ and s are the respective
normal and shear strains associated with transformation. The energy U can be used as
a rigorous variant selection criterion when the role of any plastic strain is unimportant
[12].
With displacive transformations it is possible to calculate the macroscopic
plastic-strains as a function of texture; they can also be used to characterise the
texture because they are in general anisotropic [12, 33, 34]. An arbitrary vector u
traversing a grain of austenite prior to transformation (Fig. 12.5a) makes an intercept
∆u with a domain of austenite that eventually transforms, after which it becomes a
new vector v given by:
v = P∆u + (u − ∆u) (12.15)
When many plates form in many austenite grains, u traverses a polycrystalline sam-
ple of austenite so this equation is generalised as follows [12]:
. 24
n . % . 24
n . &
v= Pkj ∆ukj + u − ∆ukj (12.16)
k=1 j=1 k=1 j=1
where (S R γk ) is the rotation matrix relating the basis vectors of the kth austenite
grain to the sample axes, and (γk R S) is the inverse of that rotation matrix. In this
way, the calculation described in equation 12.15 can be conducted in the sample
frame of reference.
Calculations illustrating the anisotropy of strains as a function of the number of
crystallographic variants of martensite allowed are illustrated in Fig 12.6a for uni-
axial tension when transformation occurs from a randomly oriented set of austenite
224 12 Martensite
(a) (b)
FIGURE 12.5
The deformation of an initial vector u by the formation of bainite. (a) An austenite
grain prior to transformation, with the ultimate location of a plate of bainite marked.
(b) Following displacive transformation.
grains. That displacive transformations produce highly anisotropic strains when vari-
ant selection is significant has been demonstrated experimentally [33–35].
Fig 12.6b shows that transformation texture is absent when 24 variants form in
each austenite grain, for typical intensities of texture in the austenite; the only strain
visible in these circumstances is the an averaged isotropic volume expansion. The
strength of the transformation texture increases as the number of variants per austen-
ite grain decreases.
It is emphasised here that the calculation of texture as presented here neglects the
fact that favoured variants will in fact have greater volume fractions in the microstruc-
ture. The calculations refer simply to the orientations of the martensite plates, not to
the intensities as might be observed in X-ray experiments where dominant variants
will contribute most. To deal with this problem, it is necessary to couple thermody-
namics and crystallography as described elsewhere [36, 37].
12.10 Summary 225
(a) (b)
FIGURE 12.6
These diagrams show the plastic strain that develops in a sample of 500 austenite
grains placed under tension along the longitudinal axis. The differences in the orthog-
onal strains correlate with the intensity of transformation texture, assuming that each
variant that forms contributes equally to the fraction of transformation. (a) Strains
developed due to transformation along the [1 0 0]S direction (labelled longitudinal,
along the stress axis), and the transverse directions [0 1 0]S and [0 0 1]S . (b) Tensile
stress, but transformation beginning from a variety of starting austenite textures and
illustrating only the longitudinal stress.
12.10 Summary
The crystallographic theory described is one of the most complete methods of pre-
dicting observable features of martensitic transformations, including aspects of the
lattice-invariant deformation (slip or twinning), the orientation relationship, the shape
deformation and even the lenticular shape of the martensite plates. The irrational na-
ture of some of the characteristics of martensite are also predicted to a high precision.
The inputs required are the lattice parameters of the parent and produce unit cells and
possible modes of lattice-invariant deformation.
It is worth emphasising that the theory is based on an elastically accommodated
plate that is not interacting with any other martensite plates. In particular, plastic
deformation driven by the shape deformation, and the impingement of the strain
fields of other plates can alter the morphology and crystallography [38, 39].
References
227
228 12 References
36. H. N. Han, C. G. Lee, C.-S. Oh, T.-O. Lee, and S.-J. Kim: ‘A model for de-
formation behaviour and mechanically induced martensitic transformation of
metastable austenitic steel’, Acta Materialia, 2004, 52, 5203–5214.
37. H. K. D. H. Bhadeshia, A. Chintha, and S. Kundu: ‘Model for multiple stress
affected martensitic transformations, microstructural entropy and consequences
on scatter in properties’, Materials Science and Technology, 2014, 30, 160–165.
38. G. Miyamoto, A. Shibata, T. Maki, and T. Furuhara: ‘Precise measurement
of strain accommodation in austenite matrix surrounding martensite in ferrous
alloys by electron backscatter diffraction analysis’, Acta Materialia, 2009, 57,
1120–1131.
39. J. H. Yang, and C. M. Wayman: ‘Self-accommodation and shape memory
mechanism of epsilon martensite - I. experimental observations’, Materials
Characterization, 1992, 28, 23–35.
13
Interfaces
Abstract
That interfaces have structure, some of which may be periodic, is now established
and observable using techniques that have sufficient spatial resolution. Small misori-
entations between like-crystals that are connected at a boundary, can be related to
arrays of dislocations, which in turn lead to an estimation of the boundary energy
per unit area. In this chapter, we describe a method for estimating the total Burgers
vector content of an arbitrary interface that in general connects crystals with differ-
ent structures, as a function of the degrees of freedom that specify the interface. The
concept of the coincidence site lattice is then generalised as Bollmann’s O-lattice.
The ‘secondary’ dislocations that describe how a boundary may move and yet pre-
serve the level of fit add to the impressive array of tools available in understanding
interfaces. There remain some difficulties, which are highlighted towards the end of
the discussion.
13.1 Introduction
Atoms located at the boundary between crystals must in general be displaced from
positions they would occupy in the undisturbed crystal. Nevertheless, it is now well
established that many interfaces have a periodic structure. In such cases, the misfit be-
tween the crystals connected by the boundary is not distributed uniformly over every
element of the interface; it is localised periodically into discontinuities which sepa-
rate patches of the boundary where the fit between the two crystals is good or perfect.
When these discontinuities are well separated, they may individually be recognised
as interface dislocations which separate coherent patches in the boundary, which is
macroscopically said to be semi-coherent.
Stress-free coherent interfaces can exist only between crystals which can be re-
lated by a transformation strain which is an invariant-plane strain. This transforma-
tion strain may be real or notional as far as the calculation of the interface structure
is concerned, but a real strain implies the existence of an atomic correspondence and
an associated macroscopic shape change of the transformed region between the two
crystals, which a notional strain does not.
Incoherency presumably sets in when the misfit between adjacent crystals is so
231
232 13 Interfaces
13.2 Misfit
The misfit across an interface can formally be described in terms of the net Burgers
vector bt crossing a vector p in the interface [1–3]. If this misfit is sufficiently small
then the boundary structure may relax into a set of discrete interfacial dislocations
where the misfit is concentrated)( which are separated by patches of good fit.
In any case, bt may be deduced by constructing a Burgers circuit across the inter-
face, and examining the closure failure when a corresponding circuit is constructed
in a perfect reference lattice. The procedure is illustrated in Fig. 13.1, where crystal
A is taken to be the reference lattice. An initial right-handed Burgers circuit OAPBO
is constructed such that it straddles the interface across any vector p=OP in the in-
terface; the corresponding circuit in the perfect reference lattice is constructed by
deforming crystal B (of the bi-crystal A-B) in such a way that it is converted into
the lattice of A, thereby eliminating the interface. If the deformation (A S A) con-
−1
verts the reference lattice into the B lattice, then the inverse deformation (A S A)
converts the bi-crystal into a single A crystal, and the Burgers circuit in the perfect
reference lattice becomes OAPP′ B, with a closure failure PP′ , which is identified as
bt . Inspection of the vectors forming the triangle OPP′ of Fig. 13.1b shows that:
Hence, the net Burgers vector content bt crossing an arbitrary vector p in the
interface is given formally by equation 13.1. The misfit in any interface can in general
be accommodated with three arrays of interfacial dislocations, whose Burgers vectors
bi (i = 1, 2, 3) form a non-coplanar set. Hence, bt can in general be factorised into
three arrays of interfacial dislocations, each array with Burgers vector bi , unit line
vector li and array spacing di , the latter being measured in the interface plane. If the
unit interface normal is n, then a vector mi may be defined as (the treatment that
follows is due to Knowles [4] and Read [5]):
mi = n ∧ li /di (13.2)
We note the |mi | = 1/di , and that any vector p in the interface crosses (mi .p)
dislocations of type I (see Fig. 13.1c). Hence, for the three kinds of dislocations we
have
bt = (m1 .p)b1 + (m2 .p)b2 + (m3 .p)b3
13.2 Misfit 233
(a) (b)
(c) (d)
FIGURE 13.1
(a,b) Burgers circuit used to define the formal dislocation content of an interface [3],
(c) the vector p in the interface, (d) relationship between li , mi , ci and n.
We note that the Burgers vectors bi of interface dislocations are generally lattice
translation vectors of the reference lattice in which they are defined. This makes them
perfect in the sense that the displacement of one of the crystals through bi relative
to the other does not change the structure of the boundary. On the basis of elastic
strain energy arguments, bi should be as small as possible. On substituting this into
equation 13.1, we get:
and
(l1 ; A) (A T′ A)[A∗ ; b∗1 ] = 0 (13.6)
/ 01 2
[A∗ ;c1 ]
and
(m1 ; A)[A∗ ; c1 ] = |m1 |2 i.e. m1 .c1 = |m1 |2
These equations indicate that the projection of c1 along m1 is equal to the magnitude
of m1 . Armed with this and the earlier result that c1 is normal to l1 , the vector
diagram illustrated in Fig. 13.1d can be constructed to illustrate the relations between
m1 , l1 , n and c1 . From this diagram, it is evident that
m1 = c1 − (c1 .n)n
so that |m1 ∧ n| = |c1 ∧ n| = 1/d1 (13.7)
l1 ∥ c1 ∧ n
alternatively,
(1/d1 )l1 = c1 ∧ n (13.8)
Relations of this type are, after appropriately changing indices, applicable also to
the other two arrays of interfacial dislocations, so that a method has been achieved
of deducing the line vectors and array spacings to be found in an interface (of unit
normal n) connecting two arbitrary crystals A and B, related by the deformation
(A S A) which transforms the crystal A to B.
Example 13.1: Symmetrical tilt boundary
Given that the Burgers vectors of interface dislocations in low-angle bound-
aries are of the form ⟨1 0 0⟩, calculate the dislocation structure of a symmetrical
tilt boundary formed between two grains, A and B, related by a rotation of 2θ
about the [1 0 0] axis. The crystal structure of the grains is simple cubic with a
unit lattice parameter.
The definition of a tilt boundary is that the unit boundary normal n is per-
pendicular to the axis of rotation which generates one crystal from the other; a
symmetrical tilt boundary has the additional property that the lattice of one crys-
tal can be generated from the other by reflection across the boundary plane.
If we choose the orthonormal bases A and B (with basis vectors parallel to the
cubic unit cell edges) to represent crystals A and B respectively, and also arbitrar-
ily choose B to be the reference crystal, then the deformation which generates the
13.2 Misfit 235
A crystal from B is a rigid body rotation (B J B) consisting of a rotation of 2θ
about [1 0 0]B . Hence, (B J B) is given by (see equation 7.9):
⎛ ⎞
1 0 0
⎜ ⎟
(B J B) = ⎜ ⎝0 cos 2θ − sin 2θ⎠
⎟
0 sin 2θ cos 2θ
−1
and (B T B) = I − (B J B) is given by
⎛ ⎞
0 0 0
⎜ ⎟
(B T B) = ⎜
⎝0 1 − cos 2θ
⎟
sin 2θ ⎠
0 − sin 2θ 1 − cos 2θ
so that ⎛ ⎞
0 0 0
⎜ ⎟
(B T′ B) = ⎜
⎝0 1 − cos 2θ − sin 2θ ⎠ .
⎟
0 sin 2θ 1 − cos 2θ
Taking [B; b1 ] = [1 0 0], [B; b2 ] = [0 1 0] and [B; b3 ] = [0 0 1], from
equation 13.4 we note that [B∗ ; b∗1 ] = [1 0 0], [B∗ ; b∗2 ] = [0 1 0] and [B∗ ; b∗3 ] =
[0 0 1]. Equation 13.6 can now be used to obtain the vectors ci :
6 7
[B∗ ; c1 ] = (B T′ B)[B∗ ; b∗1 ] = 0 0 0
6 7
[B∗ ; c2 ] = (B T′ B)[B∗ ; b∗2 ] = 0 1 − cos 2θ sin 2θ
6 7
[B∗ ; c3 ] = (B T′ B)[B∗ ; b∗3 ] = 0 − sin 2θ 1 − cos 2θ
Since c1 is a null vector, dislocations with Burgers vector b1 do not exist in the
interface; furthermore, since c1 is always a null vector, irrespective of the bound-
ary orientation n, this conclusion remains valid for any n. This situation arises
because b1 happens to be parallel to the rotation axis, and because (B J B) is an
invariant-line strain, the invariant line being the rotation axis; since any two crys-
tals of identical structure can always be related by a transformation which is a rigid
body rotation, it follows that all grain boundaries (as opposed to phase boundaries)
need only contain two sets of interface dislocations. If b1 had not been parallel to
the rotation axis, then c1 would be finite and three sets of dislocations would be
necessary to accommodate the misfit in the boundary.
To calculate the array spacings di it is necessary to express n in the B∗ basis.
A symmetrical tilt boundary always contains the axis of rotation and has the same
indices in both bases. It follows that
so that
[B; l2 ] = [1 0 0]
and
d2 = 1/2 sin θ
similarly,
(1/d3 )l3 = c3 ∧ n = [0 0 0]
so that dislocations with Burgers vector b3 have an infinite spacing in the inter-
face, which therefore consists of just one set of interface dislocations with Burgers
vector b2 .
These results are of course identical to those obtained from a simple geometri-
cal construction of the symmetrical tilt boundary [3]. We note that the boundary is
glissile (i.e. its motion does not require the creation or destruction of lattice sites)
because b2 lies outside the boundary plane, so that the dislocations can glide con-
servatively as the interface moves. In the absence of diffusion, the movement of the
boundary leads to a change in shape of the ‘transformed’ region, a shape change
described by (B J B) when the boundary motion is towards the crystal A.
If on the other hand, the interface departs from its symmetrical orientation (with-
out changing the orientation relationship between the two grains), then the boundary
ceases to be glissile, since the dislocations with Burgers vector b3 acquire finite spac-
ings in the interface. Such a boundary is called an asymmetrical tilt boundary. For
example, if n is taken to be (n; B∗ ) = (0 1 0), then
The edge dislocations with Burgers vector b3 lie in the interface plane and there-
fore have to climb as the interface moves. This renders the interface sessile.
Finally, we consider the structure of a twist boundary, a boundary where the axis
of rotation is parallel to n. Taking (n; B∗ ) = (1 0 0), we find:
actually have a small edge component. This is because the dislocations, where they
mutually intersect, introduce jogs into each other so that the line vector in the region
between the points of intersection does correspond to a pure screw orientation, but
the jogs make the macroscopic line vector deviate from this screw orientation.
Example 13.2: Interface between alpha and beta brass
The lattice parameter of the fcc alpha phase of a 60 wt% Cu-Zn alloy is 3.6925
Å, and that of the bcc beta phase is 2.944 Å [6]. Two adjacent grains A and B
(alpha phase and beta phase respectively) are orientated in such a way that
[1 1 1]A ∥ [1 1 0]B
[1 1 2]A ∥ [1 1 0]B
[1 1 0]A ∥ [0 0 1]B
and the grains are joined by a boundary which is parallel to (1 1 0)A . Assuming
that the misfit in this interface can be fully accommodated by interface disloca-
tions which have Burgers vectors [A; b1 ] = 21 [1 0 1], [A; b2 ] = 12 [0 1 1] and
[A; b3 ] = 21 [1 1 0], calculate the misfit dislocation structure of the interface. Also
assume that the smallest pure deformation which relates fcc and bcc crystals is the
Bain strain.
The orientation relations provided are first used to calculate the coordinate
transformation matrix (B J A), using the procedure given in Example 9.3. (B J A)
is thus found to be:
⎛ ⎞
0.149974 0.149974 1.236185
⎜ ⎟
(B J A) = ⎜ ⎝0.874109 0.874109 −0.212095⎠ .
⎟
The smallest pure deformation relating the two lattices is stated to be the Bain
strain, but the total transformation (A S A), which carries the A lattice to that of
B, may include an additional rigid body rotation. Determination of the interface
structure requires a knowledge of (A S A), which may be calculated from the
equation (A S A)−1 = (A C B)(B J A), where (A C B) is the correspondence
matrix. Since (B J A) is known, the problem reduces to the determination of the
correspondence matrix; if a vector equal to a basis vector of B, due to transforma-
tion becomes a new vector u, then the components of u in the basis A form one
column of the correspondence matrix.
For the present example, the correspondence matrix must be a variant of the
Bain correspondence. The vector that [1 0 0]B must become, as a result of trans-
formation, either a vector of the form ⟨1 0 0⟩A , or a vector of the form 12 ⟨1 1 0⟩A ,
although we do not know its final specific indices in the austenite lattice. How-
ever, from the matrix (B J A) we note that [1 0 0]B is close to [0 0 1]A , making it is
reasonable to assume that [1 0 0]B corresponds to [0 0 1]A , so that the first column
of (A C B) is [0 0 1]. Similarly, using (B J A) we find that [0 1 0]B is close to
238 13 Interfaces
[1 1 0]A and [0 0 1]B is close to [1 1 0]A , so that the other two columns of (A C B)
are [ 12 21 0] and [− 21 21 0]. Hence, (A C B) is found to be:
⎛ ⎞
0 1 1
1⎜ ⎟
(A C B) = ⎜ 0 1 1
⎟
2 ⎝ ⎠
2 0 0
−1
so that (A S A) = (A C B)(B J A) is given by:
⎛ ⎞
0.880500 −0.006386 −0.106048
−1 ⎜ ⎟
(A S A) = ⎜ ⎝ −0.006386
⎟
0.880500 −0.106048⎠
0.149974 0.149974 1.236183
and since
−1
(A T A) = I − (A S A) ,
′
so that (A T A) becomes:
⎛ ⎞
0.119500 0.006386 −0.149974
⎜ ⎟
(A T′ A) = ⎜
⎝0.006386 0.119500 −0.149974⎠ .
⎟
and from equation 13.8 and the fact that (n; A∗ ) = 3.6925(−0.707107 0.707107 0),
Grain A
! "# $
ABCDEF ABC DEF ABCDEF ABC DEF ABC DEF
| | | | | | [1 1 2]A →
# B C D E$!F A B C D" C
A # B A F E D C B A F$!E D C B A F E D C B"
Grain A Grain B
The first sequence represents a stack of (1 1 2)A planes, as do the first 9 layers of the
second sequence. The remainder of the second sequence represents a stack of (1 1 2)B
planes (note that B lattice is obtained by reflection of A about the 9th layer of the second
sequence). Since the two sequences have the same origin, comparison of the first sequence
with the B part of the second sequence amounts to allowing the two crystals to interpene-
trate in order to identify coincidences. Clearly, every 3rd layer of grain B coincides exactly
with a layer from the A lattice (dashed vertical lines), giving Σ = 3.
A boundary parallel to (1 1 2)A will be fully coherent; at least 1 in 3 of the sites in any
other boundary, such as (0 2 1)A will be coincidence sites.
x does not necessarily have integral components in the A basis (i.e. it need not be
a lattice vector of A). CSL vectors, on the other hand, identify lattice points which
are common to both A and B, and therefore are lattice vectors of both crystals. It
follows that CSL vectors have integral indices when referred to either crystal. Hence,
x is only a CSL vector if it has integral components in the basis A. We note that x
always has integral components in B, because a lattice vector of A (such as u) always
deforms into a lattice vector of B.
13.3 Coincidence site lattices 241
The meaning of Σ is that 1/Σ of the lattice sites of A or B are common to both
A and B. It follows that any primitive lattice vector of A or B, when multiplied
by Σ, must give a CSL vector. Σx must therefore always be a CSL vector and if
equation 13.9 is multiplied by Σ, then we obtain an equation in which the vector u
always transforms into a CSL vector:
i.e. given that u is a lattice vector of A, whose components have no common factor,
Σx is a CSL vector with integral components in either basis. This can only be true
if the matrix Σ(A S A) has elements which are all integral since it is only then that
Σ[A; x] has elements which are all integral.
It follows that if an integer H can be found such that all the elements of the matrix
H(A S A) are integers (without a common factor), then H is the Σ value relating A
and B.
The rotation matrix corresponding to the rotation 180◦ about [1 1 2]A is given by
(equation 7.9) ⎛ ⎞
2 1 2
1⎜ ⎟
(A J A) = ⎜ ⎟
3 ⎝1 2 2⎠
2 2 1
and since 3 is the integer which when multiplied with (A J A) gives a matrix of in-
tegral elements (without a common factor), the Σ value for this orientation is given
by Σ = 3 . For reasons of symmetry (see Chapter 9), the above rotation is crystal-
lographically equivalent to a rotation of 60◦ about [1 1 1]A with the rotation matrix
given by ⎛ ⎞
2 2 1
1⎜ ⎟
(A J A) = ⎜ ⎟
3 ⎝1 2 2⎠
2 1 2
◦
so it is not surprising that a rotation of 60 about [1 1 1]A also corresponds to a Σ = 3
value.
Finally, we see from equation 13.10 that if the integer H (defined such that
H(A S A) has integral elements with no common factor) turns out to be even, then
the Σ value is obtained by successively dividing H by 2 until the result H’ is an
odd integer. H’ then represents the true Σ value. This is because (equation 13.10) if
H[A; x] is a CSL vector and if H is even, then Hx has integral even components in
A, but H ′ x would also have integral components in A and would therefore represent
a smaller CSL vector. From page 191, the transformation strain relating fcc-austenite
and hcp-martensite is given by:
⎛ ⎞ ⎛ ⎞
1.083333 0.083333 0.083333 13 1 1
⎜ ⎟ 1 ⎜ ⎟
(Z P Z) = ⎜ ⎝ 0.083333 1.083333 0.083333⎠ = 12 ⎝ 1
⎟ ⎜
13 1⎠
⎟
fcc
1 2/ 0
ABCABCABCABCABCABCABCAB
| | | | | | [1 1 1]γ ∥ [0 0 0 1]hcp →
A
/ BCAB
01 C A B C2 A
/ CACAC A01C A C A C A C2
fcc hcp
so that Σ = 11 (A is the basis symbol representing one of the cubic crystals). The
24 symmetry operations of the cubic lattice are given by:
13.3 Coincidence site lattices 243
This procedure can be used to derive all the axis-angle pair representations of
any Σ value, and the table below gives some of the CSL relations for cubic crys-
244 13 Interfaces
tals, quoting the axis-angle pair representations which have the minimum angle of
rotation, and also those corresponding to twin axes.
TABLE 13.1
Some CSL relations for Cubic crystals [3]
3 60.0 ⟨1 1 1⟩ ⟨1 1 1⟩, ⟨1 1 2⟩
5 36.9 ⟨1 0 0⟩ ⟨0 1 2⟩, ⟨0 1 3⟩
7 38.2 ⟨1 1 1⟩ ⟨1 2 3⟩
9 38.9 ⟨1 1 0⟩ ⟨1 2 2⟩, ⟨1 1 4⟩
11 50.5 ⟨1 1 0⟩ ⟨1 1 3⟩, ⟨2 3 3⟩
13a 22.6 ⟨1 0 0⟩ ⟨0 2 3⟩, ⟨0 1 5⟩
13b 27.8 ⟨1 1 1⟩ ⟨1 3 4⟩
15 48.2 ⟨2 1 0⟩ ⟨1 2 5⟩
17a 28.1 ⟨1 0 0⟩ ⟨0 1 4⟩, ⟨0 3 5⟩
17b 61.9 ⟨2 2 1⟩ ⟨2 2 3⟩, ⟨3 3 4⟩
19a 26.5 ⟨1 1 0⟩ ⟨1 3 3⟩, ⟨1 1 6⟩
19b 46.8 ⟨1 1 1⟩ ⟨2 3 5⟩
21a 21.8 ⟨1 1 1⟩ ⟨2 3 5⟩, ⟨1 4 5⟩
21b 44.4 ⟨2 1 1⟩ ⟨1 2 4⟩
repulsive interaction at short range), coincidence site lattices may not exist. For
example, in the case of a Σ3 twin in a ‘hard’ bcc material, with a {1 1 2} coherent
a
twin plane, it is found that a small rigid translation (by a vector 12 ⟨1 1 1⟩) of the
twin lattice lowers the energy of the interface (Fig. 13.2) [3, 14]. Because of this
relaxation, the lattices no longer have a common origin and so the coincidences
vanish. Nevertheless, boundaries which contain high densities of coincidence sites
before relaxation may be expected to represent better fit between the lattices, and
thus have low energies relative to other boundaries. This is because the periodic
nature and the actual repeat period of the structure of the interface implied by the
CSL concept is not destroyed by the small translation. The rigid body translations
mentioned above have been experimentally established in the case of aluminium;
although not conclusively established, the experiments suggest that the translation
may have a component outside the interface plane, but the atomistic calculations
cannot predict this since they always seem to be carried out at constant volume
[14].
(iv) The physical significance of CSL’s must diminish as the Σ value increases,
because only a very small fraction of atoms in an interface can then be common
to both the adjacent crystals.
(v) The Σ = 3 value is independent [3] of symmetry considerations. This is conve-
nient since such matrices do not uniquely relate the grains; it is usually necessary
to impose criteria to allow physically reasonable choices to be made.
(a) (b)
FIGURE 13.2
{112} coherent twin boundary in a bcc material, initially with an exact Σ3 CSL. The
figure on the right illustrates the structure after a rigid body translation is included
[3, 14]
246 13 Interfaces
[A; y] = (A S A)[A; x]
If the point y is crystallographically equivalent to the point x, in the sense that it has
the same internal coordinates as x, then y is also a point of the O-lattice, designated
O. This means that y = x + u, where u is a lattice vector of A. Since y is only an
O-point when y = x + u, we may write that y is an O-lattice vector O if [15, 16]
or in other words,
[A; u] = (A T A)[A; O] (13.11)
−1
where (A T A) = I − (A S A) . It follows that
−1
[A; O] = (A T A) [A; u]. (13.12)
By substituting for u the three basis vectors of A in turn [3], we see that the columns
of (A T A)−1 define the corresponding base vectors of the O-lattice.
Since O-points are points of perfect fit, mismatch must be at a maximum in be-
tween neighbouring O-points. When O is a primitive O-lattice vector, equation 13.12
states that the amount of misfit in between the two O-points connected by O is given
by u. A dislocation with Burgers vector u would thus accommodate this misfit and
localise it at a position between the O-points, and these ideas allow us to consider a
dislocation model of the interface in terms of the O-lattice theory.
Three sets of dislocations with Burgers vectors bi (which form a non-coplanar
set) are in general required to accommodate the misfit in any interface. If the Burgers
vectors b1 , b2 and b3 are chosen to serve this purpose, and each in turn substituted
into equation 13.12, then the corresponding O-lattice vectors O1 , O2 and O3 are
obtained. These vectors Oi thus define the basis vectors of an O-lattice unit cell ap-
propriate to the choice of bi . If the O-points of this O-lattice are separated by ‘cell
walls’ which bisect the lines connecting neighbouring O-points, then the accumulat-
ing misfit in any direction can be considered to be concentrated at these cell walls
[15, 16]. When a real interface (unit normal n) is introduced into the O-lattice, its
line intersections with the cell walls become the interface dislocations with Burgers
vectors bi and unit line vectors li parallel to the line of intersection of the interface
with the cell walls.
The three O-lattice cell walls (with normals O∗i ) have normals parallel to O∗1 =
O2 ∧ O3 , O∗2 = O3 ∧ O1 and O∗3 = O1 ∧ O2 , so that the line vectors of the
dislocations are given by l1 ∥ O∗1 ∧ n, l2 ∥ O∗2 ∧ n and l3 ∥ O∗3 ∧ n. Similarly,
1/di = |O∗i ∧ n| . These results are exactly equivalent to the theory developed at
248 13 Interfaces
the beginning of this chapter (equation 13.1-13.8) but the O-lattice theory perhaps
gives a better physical picture of the interface, and follows naturally from the CSL
approach [3]. The equivalence of the two approaches arises because equation 13.11
is identical to equation 13.1 since u and O are equivalent to bt and p respectively.
Example 13.5: alpha/beta brass interface using O-lattice theory
Derive the structure of the alpha/beta brass interface described on page 237
using O-lattice theory.
From page 237, the matrix (A T A) is given by:
⎛ ⎞
0.119500 0.006386 0.016048
⎜ ⎟
(A T A) = ⎜ ⎝ 0.006386 0.119500 0.016048 ⎠
⎟
so that ⎛ ⎞
−52.448 −61.289 −51.069
−1 ⎜ ⎟
(A T A) =⎜
⎝−61.289 −52.448 −51.064⎠
⎟
The matrix (A J3 A) describing the deviation from the exact CSL is given by [16]
−1
(A J3 A) = (A J2 A) (A J A), so that
⎛ ⎞
0.996829 0.071797 −0.034313
⎜ ⎟
(A J3 A) = ⎜ ⎝−0.070528 0.996829 0.036850 ⎠ .
⎟
The secondary dislocation structure can now be calculated using the procedures
described on page 233:
[A∗ ; b∗1 ] = [1 1 1]
[A∗ ; b∗2 ] = [1 1 1]
[A∗ ; b∗3 ] = [1 1 1]
so that
001 001
001 001
(a) 010 010 (b) 010 010
FIGURE 13.3
Σ5 coincidence system for fcc crystals [22]. Filled symbols are lattice A, unfilled
ones lattice B and coincidence sites are a mixture of the two. lattice sites in the
plane of the diagram are represented as circles whereas those displaced by 12 [100]
are represented as triangles. The [100] axis is normal to the plane of the diagram.
1
Fig. 13.3b is obtained by displacing lattice B by [A; b] = 10 [0 3 1] relative to
lattice A and it is obvious that the basic pattern of lattice sites and CSL sites remains
unaffected by this translation, despite the fact that b is not a lattice vector of A or B.
It is thus possible for secondary dislocations to have Burgers vectors which are not
lattice translation vectors, but are vectors of the DSC lattice. The DSC lattice, or the
Displacement Shift Complete lattice, is the coarsest lattice which contains the lattice
points of both A and B, and any DSC lattice vector is a possible Burgers vector for
252 13 Interfaces
a perfect secondary dislocation. We note that the displacement b causes the original
coincidences (Fig. 13.3a) to disappear and be replaced by an equivalent set of new
coincidences (Fig. 13.3b), and this always happens when b is not a lattice translation
vector. This shift of the origin of the CSL has an important consequence on the
topography [22] of any boundary containing secondary dislocations with non-lattice
Burgers vectors.
Considering again Fig. 13.3, suppose that we introduce a boundary into the CSL,
with unit normal [A; n] ∥ [0 1 2], so that its trace is given by XX on Fig. 13.3a. The
effect of the displacement b of crystal B relative to A, due to the presence of a sec-
ondary dislocation, is to shift the origin of the CSL; if the boundary originally at XX
is to have the same structure after the displacement then it has to shift to the position
YY in Fig. 13.3b. Because a dislocation separates slipped from unslipped regions,
the shift of the boundary occurs at the position of the secondary dislocation so that
the boundary is stepped at the core of this dislocation. One such step is illustrated in
Fig. 13.4.
b
001
010
001
010
FIGURE 13.4
The presence of a step [23] in a Σ5, (3 1 0)A boundary of an fcc crystal containing
a secondary interface dislocation with [A; b] = (a/10)[1 3 0]. The symbolism is
identical to that of Fig. 13.3.
The following further points about the DSC lattice and its consequences should
be noted:
(i) The DSC lattice can be constructed graphically simply by inspection, bearing in
mind that it is the coarsest lattice containing lattice sites from both the crystals
orientated at an exact CSL orientation. Rather detailed analytical methods for
computing the basis vectors of the DSC lattice have been presented elsewhere
13.7 Some difficulties associated with interface theory 253
13.8 Summary
The Burgers vector content of an interface between two crystals that do not neces-
sarily have the same structure, can be obtained by the closure failure of a Burgers
circuit constructed across the interface, with the comparison made against the same
circuit constructed after transforming the bi-crystal into a single reference-crystal.
The Burgers vector content can then be de-convoluted into arrays of interfacial dis-
locations, with each array defined by a Burgers and line vector. This process requires
a prior choice of individual dislocation Burgers vectors but once this is done, the
calculations permit a comparison against experimental observations.
There exist special orientations, the coincidence site lattices, where there is a
fraction of lattice points that is common to both crystals when their lattices are al-
lowed to interpenetrate and fill all space, assuming they have a common origin. The
single crystal may be regarded at Σ = 1 CSL, and when divided to create a bi-crystal
with a small misorientation, leads to the classical dislocation description of the tilt
or twist boundary. Other CSL orientations such as Σ = 3 etc. can also be used as
reference lattices, with deviations from the exact CSL described in terms of dislo-
cation arrays. Furthermore, a boundary located on a CSL plane can be displaced so
that it maintains the coincidences using secondary dislocations as described in the
‘complete pattern shift’ or DSC lattice.
The idea of coincidences of lattice points can be generalised to coincidences of
points that have the same internal coordinates but are not necessarily lattice sites,
as in Bollmann O-lattice theory. This again is a useful concept that permits the in-
terface structure to be calculated after making assumptions about the nature of the
dislocations in the interface.
All the theory described here refers to pure interfaces that have no segregants and
which are not relaxed in the sense of rigid body translations parallel to the interface
plane. Furthermore, when the crystals are solid solutions, there will be a chemical
component of interfacial energy that would need to be taken into account.
References
257
258 13 References
27. A. P. Sutton, and V. Vitek: ‘On the structure of tilt grain boundaries in cubic
metals I. symmetrical tilt boundaries’, Philosophical Transactions of the Royal
Society A, 1983, 309, 1–36.
Appendices
259
A
Matrix methods
A.1 Vectors
Quantities (such as force, displacement) which are characterised by both magnitude
and direction are called vectors; scalar quantities (such as time) only have magnitude.
A vector is represented by an arrow pointing in a particular direction, and can be
identified by underlining the lower–case vector symbol (e.g. u). The magnitude of
u (or |u|) is given by its length, a scalar quantity. Vectors u and v are only equal if
they both point in the same direction, and if |u| = |v|. The parallelism of u and v is
indicated by writing u ∥ v. If x = −u, then x points in the opposite direction to u,
although |x| = |u|.
Vectors can be added or removed to give new vectors, and the order in which
these operations are carried out is not important. Vectors u and x can be added by
placing the initial point of x in contact with the final point of u; the initial point of
the resultant vector u + x is then the initial point of u and its final point corresponds
to the final point of x. The vector mu points in the direction of u, but |mu|/|u| = m,
m being a scalar quantity. A unit vector has a magnitude of unity; dividing a vector
u by its own magnitude u gives a unit vector parallel to u.
It is useful to refer vectors to a fixed frame of reference; an arbitrary reference
frame would consist of three non-coplanar basis vectors a1 , a2 and a3 . The vector u
could then be described by means of its components u1 , u2 and u3 along these basis
vectors, respectively, such that
If the basis vectors ai form an orthogonal set (i.e. they are mutually perpendicu-
lar) then the magnitude of u is:
|u|2 = (u1 |a1 |)2 + (u2 |a2 |)2 + (u3 |a3 |)2
A dot or scalar product between two vectors u and x (order of multiplication not
261
262 A Matrix methods
important) is given by u.x = |u| × |x| cos θ, θ being the angle between u and x. If x
is a unit vector then u.x gives the projection of u in the direction x.
The cross or vector product is written u ∧ x = |u| × |x| sin θy, where y is a unit
vector perpendicular to both u and x, with u,x and y forming a right–handed set. A
right–handed set u, x, y implies that a right–handed screw rotated through an angle
less than 180◦ from u to x advances in the direction y. The magnitude of u ∧ x gives
the area enclosed by a parallelogram whose sides are the vectors u and x; the vector
y is normal to this parallelogram. Clearly, u ∧ x ̸= x ∧ u.
If u, x and z form a right–handed set of three non-coplanar vectors then u ∧ x.z
gives the volume of the parallelepiped formed by u, x and z. It follows that u ∧ x.z =
u.x ∧ z = z ∧ u.x.
The following relations should be noted:
u ∧ x = −x ∧ u
u.(x ∧ y) = x.(y ∧ u) = y.(u ∧ x)
u ∧ (x ∧ y) ̸
= (u ∧ x) ∧ y
u ∧ (x ∧ y) = (u.y)x − (u.x)y (A.1)
A.2 Matrices
A.2.1 Definition, addition, scalar multiplication
A matrix is a rectangular array of numbers, having m rows and n columns, and is
said to have an order m by n. A square matrix J of order 3 by 3 may be written as
⎛ ⎞ ⎛ ⎞
J11 J12 J13 J11 J21 J31
⎜ ⎟ ⎜ ⎟
J=⎜ ⎝J21 J22 J23 ⎠
⎟ and its transpose J′ = ⎜⎝J12 J22 J32 ⎠
⎟
where the vertical lines imply a determinant. The cofactor jij of the element Jij is
then given by multiplying the minor of Jij by (−1)i+j . The determinant (∆) of J is
thus
.n
det J = J1j j1j withj = 1, 2, 3 (A.4)
j=1
264 A Matrix methods
= > = >
2 3 ′ 2 1
A= A =
1 4 3 4
= >
−1 0.8 −0.6
det A = 5 A =
−0.2 0.4
B
General rotation matrix
The simple matrix representing a right-handed rotation of θ = 45◦ about [001] was
written on page 111 as: ⎛ ⎞
m n 0
⎜ ⎟
R=⎜
⎝−n m
⎟
0⎠
0 0 1
where m = cos θ and n = sin θ. The matrix was derived by inspection of Fig. 7.5.
This was possible because a rotation about one of the basis vectors is easy to visu-
alise. At the same time, it was stated without derivation that a general rotation about
an arbitrary unit vector [u1 u2 u3 ] can be described by the matrix g R = JRJ′ , where
the rotation matrix J represents the new coordinate system:
⎛ ⎞
u1 u1 (1 − m) + m u1 u2 (1 − m) + u3 n u1 u3 (1 − m) − u2 n
g
⎜ ⎟
R=⎜ ⎝ u1 u2 (1 − m) − u3 n u2 u2 (1 − m) + m u2 u3 (1 − m) + u1 n ⎠
⎟
u1 u3 (1 − m) + u2 n u2 u3 (1 − m) − u1 n u3 u3 (1 − m) + m
(B.1)
The right-handed angle of rotation can be obtained from the fact that
and the components of the vector u along the axis of rotation are given by
267
268 13 General rotation matrix
where for clarity, only the first two elements of g R are listed.
• The components of the rotation axis in the basis R are [001], which in the co-
ordinates of J become [J13 J23 J33 ] ≡ [u1 u2 u3 ]. And since each row and each
2 2 2
column of a rotation matrix is a unit vector, J11 + J12 = 1 − J13 . If follows that
g
the element R11 becomes u1 u1 (1 − m) + m as in equation B.1.
• Since JJ′ = I, the identity matrix, if follows that J11 J21 + J12 J22 + J13 J23 = 0
and since J13 = u1 and J23 = u2 , the term m(J11 J21 + J12 J22 ) = −u1 u2 m.
• The determinant of J is
J33 (J11 J22 − J12 J21 ) + J23 (J31 J12 − J11 J32 ) + J13 (J21 J32 − J22 J31 ) = 1
u3 (J11 J22 − J12 J21 ) + u2 (J31 J12 − J11 J32 ) + u1 (J21 J32 − J22 J31 ) = 1
Since u21 + u22 + u23 = 1, it follows that J11 J22 − J12 J21 = u3 etc. so considering
the last two listed-items together, g R12 = u1 u2 (1 − m) + u3n as in equation B.1.
The remaining terms can be proven similarly.
Finally, the Euler angles described on page 99 can be converted into a rotation
matrix. An individual rotation matrix is created for each of the three operations,
followed by taking their product in the correct sequence, to give an overall rotation
matrix (and therefore an axis-angle pair) for the equivalent orientation relation-
ship.
13.0 INDEX 269
Index
Symmetry, 16
axis-angle pairs, 242
Tetrakaidecahedron, 162
Texture, 96
displacive transformations, 222
Euler angles, 99
inverse pole figures, 96
orientation distribution, 99
pole figures, 96
representation, 96
Topology, 162
Transformation twinning, 221
Translations, 23
Turbine blade, 4
Twinning
deformation, 179
double twinning, 141
orientation relationship, 141
slip across twins, 186
transformation twins, 221
Unit cell
primitive, 6
Zinc sulphide, 69
Organized into a two-part structure aimed at readers of
differing experience levels, Geometiy of Crystals, Poly
crystals, and Phase Transformations is accessible to both
newcomers and advanced researchers within the field of
crystallography. The rst part of the text covers what any
reader in the materials sciences, physics, chemistry, earth
sciences "and natural sciences in general“ should know
about crystallography. lt is intentionally concise and covers
sufficient material to form a firm foundation. The second
part is aimed at researchers and discusses phase transfor- ‘X
.,.
concepts T
YKX
Supported by slides, video lectures, and other
,»?<
..,¢
electronic ancillaries freely accessible at: 5
4??
http.//www.phase trans.msm.cam.ac.uk/teaching.html
§¢
,<k
32> 1;’:
9?.
at
gi
we Q
gs
0/