Metal and Alloy Bonding: An Experimental Analysis
Metal and Alloy Bonding: An Experimental Analysis
An Experimental Analysis
R. Saravanan M. Prema Rani
•
123
Dr. R. Saravanan M. Prema Rani
Research Centre and PG Department Research Centre and PG Department
of Physics of Physics
The Madura College The Madura College
Madurai 625 011 Madurai 625 011
Tamil Nadu Tamil Nadu
India India
e-mail: [email protected]; e-mail: [email protected]
[email protected]
This book has been written based on the experimental results obtained through
several experimental techniques, especially the powder X-ray diffraction method,
on various metals and alloys we encounter frequently. An analysis of the
interactions of electrons in different atoms has been discussed.
Metals have useful properties including strength, ductility, high-melting points,
thermal and electrical conductivity and toughness. The key feature that distin-
guishes metals from non-metals is their bonding. The existence of free electrons in
metals has a number of profound consequences for the properties of metallic
materials. There are a large number of possible combinations of different metals
and each has its own specific set of properties. The physical properties of an alloy,
such as density, reactivity, Young’s modulus and electrical and thermal conduc-
tivity, may not differ greatly from those of its elements, but engineering properties,
such as tensile strength and shear strength, may be substantially different from
those of the constituent materials. Metals and their alloys make today’s manu-
facturing industry, agriculture, construction and communication systems, trans-
portation, defense equipments, etc. possible. Some of the major reasons for the
continuing advancements in alloys are the availability of materials, new manu-
facturing techniques and the ability to test alloys before they are produced. Most
modern alloys are, in fact, preplanned using sophisticated computer simulations,
which help to determine what properties they will display. Semiconductors have
been studied extensively due to their importance in applications. These materials
receive much attention because physical properties such as the band gap, mobility
and lattice parameter can be continuously controlled. Having such continuous
control is of importance in applications such as electronic and optical devices.
Metals and alloys have high-melting temperatures because of the heavy
bonding between the atoms. There are a variety of applications for metals and
alloys. Due to the importance of these materials, a study of their bonding
interactions has been carried out in this monograph using experimentally observed
X-ray diffraction data.
Today’s technological evolution results in developing new and sophisticated
materials of immense use in domestic, technical and industrial applications.
v
vi Preface
Usually, the synthesis of new materials, especially metals and alloys, results in
single-phase materials, but often not in single crystalline form. Hence, a complete
analysis of the structure, local distribution of atoms and electron distribution in
core, valence and bonding region is necessary using powder diffraction methods,
in addition to single crystal diffraction results, since most of the recent materials
will be initially obtained in powder form. Since one can make efforts to grow
single-crystals from powders, a prior analysis is required using powders to proceed
for single crystal growth.
In this context, we have taken some simple metals (Al, Cu, Fe, Mg, Na, Ni, Te,
Ti, Sn, V, Zn) and alloys (AlFe, CoAl, FeNi, NiAl) and collected powder XRD data
sets or used single crystal XRD data sets from the literature, to study the structure in
terms of the local and average structural properties using pair distribution function
(hereafter PDF), electron density distribution between atoms using Maximum
Entropy Method (hereafter MEM) and bonding of core and valence electron dis-
tribution using multipole technique. Particularly, the PDF analysis requires data
sets of very high values of Q (=4pSinh/k) which is achievable only through syn-
chrotron studies, but not accessible for common crystallographer/material scien-
tists. The present work gives reasonable results obtained through single crystal
work or through high Q data sets, using only powder samples. Also, a study on the
electronic structure of metals using the most versatile currently available tech-
niques like MEM and multipole method is worthwhile. If the tools available for
analysis yield highly precise information, then it is appropriate to apply it to precise
data sets available as in this work, and thus the methodology can also be tested. In
order to elucidate the distribution of valence electrons and the contraction/expan-
sion of atomic shells, multipole analysis of the electron densities was also carried
out. Recently, multipole analysis of the charge densities and bonding has been
widely used to study the electronic structure of materials.
Bonding studies in crystalline materials are very important, especially in metals,
because of their extensive use. These studies can reveal the qualitative nature of
bonding as well as the numerical values of mid-bond densities which indicate the
strength of the material under study. With the advent of versatile methods like
MEM and multipole method, bonding studies gained impetus because of the
accuracy of these methods and the fact that the experimental data can be used with
these methods to accurately determine the actual bonding between atoms.
The precise study of bonding in materials is always useful and interesting, yet
no study can reveal the real picture as no two sets of experimental data are
identical. This problem is enhanced when the model used for the evaluation of
electron densities is not entirly suitable. Fourier synthesis of electron densities can
be of use in picturing bonding between two atoms, but it suffers from the major
disadvantages of series termination error and negative electron densities which
prevent the clear understanding of bonding between atoms; the factor intended to
be analysed. The advent of MEM solves many of these problems. MEM electron
densities are always positive and even with limited number of data, one can
determine reliable electron densities resembling true densities. Currently, the
multipole analysis of charge densities has been widely used to study crystalline
Preface vii
materials. This synthesizes the electron density of an atom into core and valence
parts and yields an accurate picture of bonding in a crystalline system.
In this research monograph on metals and alloys, a complete analysis of
bonding has been made on 11 important metals and four alloys. Powder X-ray
diffraction data as well as single crystal data sets have been used for the purpose.
Charge density analysis of materials provides a firm basis for the evaluation of
the properties of the materials. Designing and engineering of new combinations of
metals requires firm knowledge of the intermolecular features. Recent advances in
technology and high-speed computation has put the crystal X-ray diffraction
technique on a firm pedestal as a unique tool for the determination of charge
density distribution in molecular crystals. Methods have been developed to make
experimental probes to unravel the features of charge densities in the intra and
intermolecular regions in crystal structures. In this report the structural details have
been elucidated from the X-ray diffraction technique through Rietveld technique.
The charge density analysis has been carried out with MEM and multipole method,
and the local and average structure analysis by atomic PDF.
This research work reveals the local and average structural properties of some
technologically important materials, which are not studied along these lines. New
understandings of the existing materials have been gained in terms of the local and
average structures of the materials. The electron density, bonding, and charge
transfer studies analysed in this work will give fruitful information to researchers in
the fields of physics, chemistry, materials science, metallurgy, etc. These properties
can be properly utilized for the proper engineering of these technologically
important materials.
Chapter 1 introduces the significance and applications of metals, alloys and
semiconductors studied in this research work. The objectives of this book are
presented. The essential mechanism of ball milling which has evolved to be a
simple and useful method for the formation of nano crystalline materials is
discussed. The current state of art of non-destructive characterisation techniques
such as X-ray diffraction and scanning electron microscope are discussed.
Chapter 2 provides a survey of the current applications of X-ray diffraction
techniques in crystal structure analysis, with focus on the recent advances made in
the scope and potential for carrying out crystal structure determination directly
from diffraction data. The basic concepts of crystal structure analysis, Rietveld
refinement and the concepts used for the estimation and analysis of charge density
in a crystal are discussed. The more reliable models for charge density estimation
like multipole formalism and MEM are discussed in detail. The local structural
analysis technique and atomic PDF is also discussed.
Chapter 3 presents the results and discussions of this research work. A detailed
account of the results of the materials analysed are presented in the subsections.
Section 3.1 (Sodium and Vanadium Metals) describes about the nature of
bonding and the charge distribution in sodium and vanadium metals are analysed
using the reported X-ray data of these metals. MEM and multipole analysis used
for bonding in these metals are elucidated and analysed. The mid-bond densities in
sodium and vanadium are found to be 0.014 and 0.723 e/Å3 respectively, giving an
viii Preface
indication of the strength of the bonds in these materials. From multipole analysis,
the sodium atom is found to contract more than the vanadium atom.
Section 3.2 (Aluminium, Nickel and Copper) describes the average and local
structures of simple metals Al, Ni and Cu are elucidated for the first time using
MEM, multipole and PDF. The bonding between constituent atoms in all the above
systems is found to be well pronounced and clearly seen from the electron density
maps. The MEM maps of all the three systems show the spherical core nature of
atoms. The mid-bond electron density profiles of Al, Ni and Cu reveal the metallic
bonding nature. The local structure using PDF profile of Ni has been compared
with that of the reported results. The R value in this work using low Q XRD data
for the PDF analysis of Ni is close to the value reported using high Q synchrotron
data. The cell parameters and displacement parameters were also studied and
compared with the reported values.
Section 3.3 (Magnesium, Titanium, Iron, Zinc, Tin and Tellurium) describes the
average and local structures of magnesium, titanium, iron, zinc, tin and tellurium
are analysed using the MEM, and PDF. The structural parameters of the metals
were refined with the well-known Rietveld powder profile fitting methodology.
One-, two- and three-dimensional electron density distributions of Mg, Ti, Fe, Zn,
Sn and Te have been mapped using the MEM electron density values obtained
through refinements. The mid-bond density in Ti is the largest value along [110]
direction among the six metal systems. From PDF analysis the first neighbour
distance is observed to decrease as the atomic number increases for all the metals.
Section 3.4 (Cobalt Aluminium and Nickel Aluminium Metal Alloys) describes
the precise electron density distribution and bonding in metal alloys CoAl and NiAl
is characterized using MEM and multipole method. Reported X-ray single-crystal
data used for this purpose. Clear evidence of the metal bonding between the con-
stituent atoms in these two systems is obtained. The mid-bond electron densities in
these systems are found to be 0.358 and 0.251 e/Å3 respectively, for CoAl and NiAl
in the MEM analysis. The two-dimensional maps and one-dimensional electron
density profiles have been constructed and analysed. The thermal vibration of the
individual atoms Co, Ni and Al has also been studied and reported. The contraction
of atoms in CoAl and expansion of Ni and contraction of Al atom in NiAl is found
from multipole analysis, in line with the MEM electron density distribution.
Section 3.5 (Nickel Chromium (Ni80Cr20)) describes the alloy Ni80Cr20 was
annealed and ball milled to study the effect of thermal and mechanical treatments
on the local structure and the electron density distribution. The electron density
between the atoms was studied by MEM and the local structure using PDF. The
electron density is found to be high for ball-milled sample along the bonding
direction. The particle sizes of the differently treated samples were realized by
SEM and through XRD. Clear evidence of the effect of ball milling is observed on
the local structure and electron densities.
Section 3.6 (Silver doped in NaCl (Na1-xAgxCl)) describes the alkali halide
Na1-xAgxCl, with two different compositions (x = 0.03 and 0.10) is studied with
regard to the Ag impurities in terms of bonding and electron density distribution.
X-ray single crystal data sets have been used for this purpose. The analysis focuses
Preface ix
on the electron density distribution and hence the interaction between the atoms is
clearly revealed by MEM and multipole analysis. The bonding in these systems is
studied using two-dimensional MEM electron density maps on the (100) and (110)
planes and one-dimensional electron density profiles along the [100], [110] and
[111] directions. The mid-bond electron densities between atoms in these systems
are found to be 0.175 and 0.183 e/Å3, respectively, for Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl. Multipole analysis of the structure is performed for these two
systems, with respect to the expansion/contraction of the ion involved.
Section 3.7 (Aluminium Doped with Dilute Amounts of Iron Impurities (0.215
and 0.304 wt% Fe)) describes the electronic structure of pure and doped alu-
minium with dilute amounts of iron impurities (0.215 and 0.304 wt % Fe) has been
analysed using reported X-ray data sets and the MEM. Qualitative as well as
quantitative assessment of the electron density distribution in these samples is
made. The mid-bond characterization leads to a conclusion about the nature of
doping of impurities. An expansion of the size of the host aluminium atom was
observed with Fe impurities.
Chapter 4 presents the conclusion of the results of the reported work.
A complete analysis on the electron density of important metals and alloys is
presented in this book. This book will be highly useful for scientists and
researchers working in the areas of metallurgy, materials science, crystallography,
chemistry and physics.
Acknowledgments
The author Dr. R. Saravanan, acknowledges his family for their kind support, help
and for making the atmosphere conducive during the course of the compilation of
this book.
The author Ms. M. Prema Rani, wishes to thank her family, husband and espe-
cially her children for their support and for motivating her in writing this book.
The authors thank the various finding agencies in India, the University Grants
Commission (UGC), Council of Scientific and Industrial Research (CSIR) and
Department of Science and Technology (DST), though they did not fund the
compilation of this book directly. But, the authors believe that the various research
tasks accomplished during the course of the work for the book may involve usage
of the resources arising out of the funds by the above agencies and hence these
agencies are gratefully acknowledged.
The authors wish to render their cordial thank to the authorities of the Madura
College, Madurai, 625 011, India for their generous support in the various research
efforts by the authors which led to the successful compilation of this book.
Research of high quality needs good support from various people including the
authorities in the concerned institutions from where the research efforts originate.
In that respect, the authors thank the principal and the board of management of the
Madura College, Madurai, 625 011, India, particularly the secretary, Mr. M.S.
Meenakshi Sundaram, The Madura College Board, Madurai, 625 011, India for his
support and encouragement in the academic and research efforts of the authors.
Editing a book on a special topic like the present one involves help, support,
and constant motivation by a large number of clause of people, right from clerical
level and up to intellectual level. The authors wish to acknowledge all those people
who could not find a place in this page of this book but who rendered their cordial
help for successfully editing this book.
The authors dedicate this book for real hard working people with real positive
qualities.
Dr. R. Saravanan
M. Prema Rani
xi
Contents
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1 Introduction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Significance of the Present Work . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Objectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4 Metals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4.1 Sodium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4.2 Vanadium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4.3 Magnesium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4.4 Aluminium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.5 Titanium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.6 Iron . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4.7 Nickel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4.8 Copper . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4.9 Zinc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4.10 Tin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.4.11 Tellurium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.5 Significance of Alloys . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.5.1 Alloys in Nuclear Reactors . . . . . . . . . . . . . . . . . . . . 14
1.5.2 Alloy Wheels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.6 Significance of the Alloys Dealt With in this Research Work. . . 14
1.6.1 Cobalt Aluminium . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.6.2 Nickel Aluminium . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.6.3 Nickel Chromium. . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.6.4 Iron–Nickel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.6.5 Sodium Chloride Doped with Silver . . . . . . . . . . . . . . 17
1.6.6 Aluminium Doped with Iron . . . . . . . . . . . . . . . . . . . 18
xiii
xiv Contents
4 Conclusion. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
4.1 Sodium and Vanadium Metals . . . . . . . . . . . . . . . . . . . . . . . . 148
4.2 Aluminium, Nickel and Copper . . . . . . . . . . . . . . . . . . . . . . . 148
4.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium . . . . . . . . 149
4.4 Cobalt Aluminium and Nickel Aluminium Metal Alloys. . . . . . 149
4.5 Nickel Chromium (Ni80Cr20) . . . . . . . . . . . . . . . . . . . . . . . . . 150
4.6 Silver Doped in NaCl (Na1-xAgxCl) . . . . . . . . . . . . . . . . . . . 150
4.7 Aluminium Doped with Dilute Amounts of Iron
Impurities (0.215 and 0.304 wt% Fe) . . . . . . . . . . . . . . ..... 150
Chapter 1
Introduction
Abstract The properties of a material are a direct result of its internal structure.
The ability to control structures through processing, and to develop new struc-
tures through various techniques, requires qualitative and quantitative analysis of
the atomic and electronic structure. The average and local structure of some
significant metals and important alloys have been analyzed and reported in this
book. This introduction chapter deals with the significance and applications of
metals, alloys and semiconductors. The essential mechanism of ball milling
which has evolved to be a simple and useful method for the formation of nano
crystalline materials is discussed. The current state-of-the-art of non-destructive
characterisation techniques such as X-ray diffraction and scanning electron
microscope are discussed.
1.1 Introduction
R. Saravanan and M. Prema Rani, Metal and Alloy Bonding: An Experimental Analysis, 1
DOI: 10.1007/978-1-4471-2204-3_1, Springer-Verlag London Limited 2012
2 1 Introduction
1.3 Objectives
Though the materials studied and reported in this research book are all metals and
alloys, the work has been divided into several parts for the sake of convenience.
They have been given as below.
1. The average and electronic structure of the following elemental metals using
Rietveld (Rietveld 1969), multipole (Hansen and Coppens 1978) and MEM
(Collins 1982) by single-crystal XRD data.
• Sodium (Na)
• Vanadium (V)
2. a. The local, average and electronic structure of the following elemental metals
using Rietveld (1969), and MEM (Collins 1982) by powder XRD data.
• Magnesium (Mg)
• Aluminium (Al)
• Titanium (Ti)
• Iron (Fe)
• Nickel (Ni)
• Copper (Cu)
• Zinc (Zn)
• Tin (Sn)
• Tellurium (Te)
b. The local structural information by analyzing the atomic pair distribution
function (PDF) (Proffen and Billinge 1999).
3. The average and electronic structure of the following metal alloys using
multipole (Hansen and Coppens 1978) and MEM (Collins 1982) by single-
crystal XRD data.
• cobalt aluminium (CoAl)
• nickel aluminium (NiAl)
• Iron nickel (FeNi)
4. a. The annealing and ball milling of the alloy nickel chromium (Ni80Cr20)
b. The study of the local, average and electronic structure of the annealed and
ball milled alloy Ni80Cr20 using Rietveld (Rietveld 1969) and MEM (Collins
1982) by powder XRD data.
c. To study the local structure using PDF (Proffen and Billinge 1999).
4 1 Introduction
1.4 Metals
Metals account for about two-thirds of all the elements and about 24% of the mass
of the planet. Metals have useful properties including strength, ductility, high-
melting points, thermal and electrical conductivity and toughness. The key feature
that distinguishes metals from nonmetals is their bonding. (Gallagher and Ingram
2001) Metallic materials have free electrons that are free to move easily from one
atom to the next. The existence of these free electrons has a number of profound
consequences in the properties of metallic materials (Kittel 2007).
The local and average structures of some technologically important metals such
as sodium, magnesium, aluminium, titanium, vanadium, iron, nickel, copper, zinc,
tin and tellurium are analyzed in this work and their typical properties and uses are
presented below.
1.4.1 Sodium
Sodium is a soft, silvery-white, highly reactive metal having only one stable isotope;
23
Na. Sodium ion is soluble in water in nearly all of its compounds. Sodium metal is
so soft that it can be cut with a knife at room temperature (Zumdahl 2007)
Sodium compounds are important for the chemical, glass, metal, paper,
petroleum, soap, and textile industries. A sodium–sulphur battery is a type of
molten metal battery constructed from sodium and sulphur. This type of battery
has a high-energy density, high efficiency of charge/discharge (89–92%) and long
cycle life, and is fabricated from inexpensive materials (Oshima et al. 2004).
NaS batteries are a possible energy storage technology to support renewable
energy generation, specifically in wind farms and solar generation plants. In the
case of a wind farm, the battery would store energy during times of high wind but
low-power demand. This stored energy could then be discharged from the batteries
during peak load periods. In addition to this power shifting, it is likely that sodium
sulfur batteries could be used throughout the day to assist in stabilising the power
output of the wind farm during wind fluctuations (Walawalkar et al. 2007). Due to
its high-energy density, the NaS battery has been proposed for space applications
(Auxer 1986).
1.4 Metals 5
1.4.2 Vanadium
Pure vanadium is a bright white metal, and is soft and ductile. It has good
corrosion resistance to alkalis, sulphuric and hydrochloric acid and salt water. The
metal has good structural strength and a low fission neutron cross section, making
it useful in nuclear applications. (Lynch 1974).
Vanadium is used in producing rust resistant, spring and high-speed tool steels.
It is an important carbide stabiliser in making steels. Vanadium is also used in
producing superconductive magnets with a field of 175,000 gauss (Lide 1999).
The role of vanadium complexes in catalytically conducted redox reactions
(Crans et al. 2004) and potential medicinal applications, such as in the treatment of
diabetes type I and type II (Crans 2000), has stimulated interest in the stereo-
chemistry and reactivity of its coordination compounds (Monfared et al. 2010).
Vanadium oxides and vanadium oxide-related compounds have a wide range of
practical applications such as catalysts, gas sensors and cathode materials for
reversible lithium batteries, electrochemical and optical devices, due to their
structural, novel electronic and optical properties (Zhang et al. 2010). Mixed metal
oxides find applications in a variety of fields due to the wide variation in their
dielectric and electrical properties. The vanadium-based oxide ceramics have high-
dielectric constant, low-dissipation factor and high-quality factor, which favour the
use of these ceramics in many fields (Nithya and Kalaiselvan 2011).
Secure and reliable power is essential in areas such as telecommunications and
information technology to safeguard the vast computer networks that have been
established. Uninterruptible power systems have incorporated battery technology
to allow smooth power feeding switch-over in the case of a power failure. In such
systems lead-acid batteries are commonly being used until generators come online
or for safe computer shutdown. The vanadium redox battery provides many
advantages over conventional batteries for emergency back-up applications. This
system stores all energy in the form of liquid electrolytes which are re-circulated
around the battery system. The electrolytes can be recharged for indefinite number
of times, or the system can be instantly recharged by mechanically exchanging the
discharged solution with recharged solution (Kazacos and Menictas 1997).
Figure 1.1 shows a vanadium redox battery. A vanadium redox battery consists of a
power cell in which two electrolytes are kept separated by an ion exchange membrane.
Both the electrolytes are vanadium based. Vanadium redox batteries are based on the
ability of vanadium to exist in four different oxidation states (V2, V3, V4 and V5), each
of which holds a different electrical charge. The electrolyte in the negative half-cell
has V3+ and V2+ ions, while the electrolyte in the positive half-cell contains V3+ and
V2+ ions. During charging, reduction in the negative half-cell converts the V3+ ions
into V2+ ions. During discharge, the process is reversed, oxidation in the negative half-
cell converts V2+ ions back to V3+ ions. The typical open-circuit voltage created
during discharge is 1.30 V at 25C (Skyllas-Kazacos 2003).
6 1 Introduction
Other useful properties of Vanadium flow batteries are their quick response to
changing loads and their extremely large overload capacities. Their extremely
rapid response times also make them perfectly well suited for UPS-type applica-
tions, where they can be used to replace lead-acid batteries and even diesel
generators.
1.4.3 Magnesium
and rear axes of cars, propeller shaft tunnel, pedals, dashboards, elements of seats,
steering wheels, elements of timer-distributors, air filters, wheel bands, oil sumps,
elements and housings of the gearbox, framing of doors and sunroofs and others
(Dobrzánski et al. 2007).
In recent times, the increased environmental concerns and the rising costs of oil
have again made magnesium and its alloys a material of interest for the automotive
industry. Considering the characteristics of low density of magnesium, its exten-
sive use in structural body parts of vehicles will offer major reductions of weight
and hence reduction in fuel consumption. Such weight reduction provides a sig-
nificant contribution to reducing the carbondioxide emission. It is estimated that an
average new car produces 156 g CO2/km travelled. This could be reduced to
around 70 g CO2/km through the application of magnesium technology (Mehta
et al. 2004).
The advantages of magnesium and magnesium alloys are, lowest density of all
metallic constructional materials, high-specific strength, good castability, suitable
for high-pressure die-casting, can be turned or milled at high speed, good wel-
dability under controlled atmosphere, much improved corrosion resistance, readily
available, better mechanical properties, resistant to ageing, better electrical and
thermal conductivity and recyclability (Mordike and Ebert 2001a, b).
Magnesium alloys have attracted increasing interest in the past few years due to
their potential as implant materials. Magnesium and its alloys are degradable
during their time of service in the human body. Magnesium alloys offer a property
profile that is very close or even similar to that of human bone (Hort et al. 2010).
1.4.4 Aluminium
Aluminium has been the dominant material in the aircraft industry for more than a
half century due to its attractive combination of light weight, strength, ductility,
corrosion resistance, ease of assembly and low cost (Dorward and Pritchett 1988).
Aluminium foam sandwiches (AFS) due to their flexible process ability and
potential of cost reduction find application in space components. Currently, these
light-weight materials find some first applications in particular fields of mechanical
engineering such as race cars, and small series of other land-based vehicles
(Schwingel et al. 2007). The use of high-strength aluminium alloys in automotive
and aircraft industries allows reducing significantly the weight of the engineering
constructions. In these fields, very often the main requirements for the components
include high fatigue and wear-resistance (Lonyuk et al. 2007). Aluminium solar
mirrors are an alternative for solar concentrators. The aluminium reflectors often
offer an initial reflectance of 85–91% for solar irradiance. They have good
mechanical properties and are easy to recycle (Almanza et al. 2009). The high
strength-to-weight advantage of aluminium alloys has made it the material of
choice for building airplanes and sometimes for the construction of land-based
structures. For marine applications, the use of high-strength, weldable and
8 1 Introduction
1.4.5 Titanium
Titanium has many desirable physical properties. The pure metal is relatively soft
and weak, but it becomes much stronger when mixed with other metals to form
alloys. The high-melting point of titanium (1,668C) shows that it is an ideal
material for the construction of high-speed aircraft and space vehicles.
Due to the exceptional strength-to-weight ratios, toughness, high stiffness and
excellent biocompatibility, titanium and its alloys are used extensively in aero-
space, chemical and biomedical applications (Kartal et al. 2010). Titanium alloys
are widely used in the aerospace industry due to their excellent fatigue/crack
propagation behaviour, and corrosion resistance (Markovsky and Semiatin 2010).
The alloys of titanium represent significant advantages over most other
engineering materials used for a variety of industrial applications due to their
resistance to corrosion, oxidation and erosion. Titanium and its alloys have high-
chemical durability as well as high strength. Their use is significant in nuclear
industry, since the mechanical strength, high-heat proof and radiation proof are
desired by many components, such as the steam condenser tubes, the irradiation
targets for transmuting radioactive wastes and the overpacks for geological
disposal of high level radioactive wastes (Setoyama et al. 2004).
Titanium is the preferred choice for surgical instrumentation due to its lighter
weight, bacterial resistance and durability. High strength-to-weight ratio, corrosion
resistance, non-toxic state and non-ferromagnetic property has made titanium ‘‘the
metal of choice’’ within the field of medicine. It is also durable and long-lasting.
When titanium cages, rods, plates and pins are inserted into the body, they can
last for more than 15 years. And dental titanium, such as titanium posts and
implants, can last even longer. Osseo integration is a unique phenomenon where
the body’s natural bone and tissue actually bonds to the artificial implant. This
firmly anchors the titanium dental or medical implant into place. Titanium is the
only metal that allows this integration. Titanium and its alloys are widely used to
replace failed hard tissues, such as artificial hip joints and dental implants
(Li et al. 2008). Mechanical properties such as high strength, ductility and fatigue
resistance, as well as a low modulus make titanium and its alloy suitable for
applications in jet propulsion systems and human body implant (Heinrich et al.
1996). Titanium has long been used as an implant material in different medical
applications, showing excellent performance in forming a close contact to the
surrounding tissues (Petersson et al. 2009).
1.4 Metals 9
1.4.6 Iron
Pure iron is silvery metal with a shining surface. It is a good conductor of heat and
electricity. Iron is used to make bridges, automobiles and support for buildings,
machines and tools. It is mixed with other elements to make alloys, the most
important of which is steel (Sparrow 1999).
Iron-based glassy alloys seem to be one of the most interesting materials due to
their soft magnetic properties including high-saturation magnetisation. They are
suitable materials for many electrical devices such as electronic measuring and
surveillance systems, magnetic wires, sensors, band-pass filters, magnetic
shielding, energy-saving electric power transformers (Nowosielski et al. 2008).
1.4.7 Nickel
Nickel is a silvery-white lustrous metal with a slight golden tinge. It is one of the
four ferromagnetic elements that exist around room temperature, the other
three being iron, cobalt, and gadolinium. Its Curie temperature is 355C. Nickel is
non-magnetic above this temperature (Kittel 1996). Nickel belongs to the transi-
tion metals and is hard and ductile. The isotopes of nickel range from 48Ni to 78Ni.
The isotope of nickel with 28 protons and 20 neutrons 48Ni is ‘‘double magic’’ and
therefore unusually stable (Audi 2003).
The metal is corrosion-resistant, finding many uses in alloys, as plating, in the
manufacture of coins, magnets, common household utensils, rechargeable batter-
ies, electric guitar strings, as a catalyst for hydrogenation, and in a variety of other
applications. Enzymes of certain life-forms contain nickel as an active centre,
which makes the metal an essential nutrient for those life-forms. It is also used for
plating and as a green tint in glass. In the laboratory, nickel is frequently used as a
catalyst for hydrogenation. Nickel is often used in coins, or occasionally as a
substitute for decorative silver.
Rechargeable nickel batteries are one type of alkaline storage cylindrical bat-
tery and classified as secondary batteries.
Nickel battery has a positive electrode made of active material-nickelous
hydroxide. Because of the perfectly, sealed construction and the efficient charge/
discharge characteristics, nickel batteries provide superior features and practical
values in long service life, high-rate discharge and stable performance. As a result,
they are widely used in many fields such as communication and telephone equip-
ment, office equipment, tools, toys and emergency devices and consumer
applications.
A typical jet engine today contains about 1.8 tonnes of nickel alloys and includes
a long list of tailor-made nickel-based-alloys to meet specific needs (Nickel
Magazine 2007). Pure nickel is a strong candidate for protective coating in bio-
diesel storage applications, due to its high resistance to the corrosive nature of
10 1 Introduction
biodiesel and its vapors and minimal catalytic effects on the oxidation of biodiesel
(Boonyongmaneerat et al. 2011).
1.4.8 Copper
Copper is a ductile metal with very high thermal and electrical conductivity. Pure
copper is rather soft and malleable, and a freshly exposed surface has a pinkish or
peachy color. It is used as a thermal conductor, an electrical conductor, a building
material, and a constituent of various metal alloys.
Copper is the most widely used metal because of high conductivity. Copper and
copper-based alloys are unique in their physical and mechanical properties. They
have excellent corrosion resistance, high resistance to fatigue and relative ease of
joining by soldering available in wide variety of forms (Pillai 2007).
Copper is easily worked, being both ductile and malleable. The ease with which
it can be drawn into wires makes it useful for electrical work in addition to its
excellent electrical properties. Copper can be machined, although it is usually
necessary to use an alloy for intricate parts, such as threaded components, to get
really good machinability characteristics. Good thermal conduction makes it
useful for heat sinks and in heat exchangers. Copper has good corrosion resistance;
it has excellent brazing and soldering properties and can also be welded, although
best results are obtained with gas metal arc welding (Sambamurthy 2007). Copper
as both metal and pigmented salt has a significant presence in decorative art.
1.4.9 Zinc
In operation, a mass of zinc particles form a porous anode, which is saturated with
an electrolyte. Oxygen from the air reacts at the cathode and forms hydroxyl ions
which migrate into the zinc paste and form zincate, releasing electrons to travel to the
cathode. The zincate decays into zinc oxide and water returns to the electrolyte. The
water and hydroxyls from the anode are recycled at the cathode, so the water is not
consumed. The reactions produce a theoretical 1.65 V, but this is reduced to 1.4–
1.35 V in available cells. Zinc-air batteries have some properties of fuel cells as well
as batteries, with zinc as the fuel, the reaction rate can be controlled by varying the air
flow, and oxidised zinc/electrolyte paste can be replaced with fresh paste. Metallic
zinc could be used as an alternative fuel for vehicles, in a zinc-air battery (Noring
et al. 1993). Zinc-air batteries are considerably more safer in combating situations
and more environmental friendly than lithium batteries (http://www.defense-update.
com/products/z/zinc-air-battery-new.htm).
A Switzerland-based company, ReVolt uses zinc-air battery technology for
hearing aids. ReVolt’s battery claims to store three times more energy than lith-
ium–ion by volume, and could incur just half the costs (http://www.goodcleantech.
com/2009/11/a_powerful_rechargeable_zinc-a.php).
1.4.10 Tin
1.4.11 Tellurium
Due to the remarkable physical properties of tellurium such as low band gap and
transparency in the infrared region, Te is used extensively in various technological
areas. Te thin films find use in microelectronic devices such as gas sensor (Shashwati
et al. 2004; Tsiulyanu et al. 2004) and optical information storage (Josef et al. 2004).
Tellurium-based thin films, suitable for applications in environmental monitoring
with considerably short-response time and high sensitivity to nitrogen dioxide at
room temperature have been reported (Tsiulyanu et al. 2001).
Tellurium, with a low band gap of 0.32 eV, is one of the most promising
materials for a shield in a passive radiative cooling (Engelhard et al. 2000).
Radiative cooling is the one among today’s challenges in materials science
research. It occurs when a body gets cold by loosing energy through radiative
processes. The phenomenon of radiative cooling uses the fact that the thermal
energy emitted by a clear sky in the ‘‘window region’’ (8–13 mm) is much less
than the thermal energy emitted by a blackbody at ground air temperature in this
wavelength range. Hence, a surface on the earth facing the sky experiences an
imbalance of outgoing and incoming thermal radiation and cools to below the
ambient air temperature. While this concept can work well at night, assuming a
1.4 Metals 13
relatively dry atmosphere, the solar energy input during the day, which is normally
much greater than that radiated out, causes heating of the system. To prevent this,
a shield is required to cover the radiating surface in order to block solar radiation
during the day as well as to prevent convective mixing in the cooled space. An
ideal radiation shield should completely reflect solar radiation, but allow complete
transmission in the ‘‘atmospheric-window’’ region. Solar radiation should prefer-
ably be reflected, as any absorbed radiation will be converted to heat somewhere in
the system (Dobson et al. 2003). The different approaches for the design of shield
are the introduction of optical scattering materials into shield substrates, and
coating of shield substrates with high-solar reflector films. Te thin films with
thickness of 111–133 nm having high-IR transmission across the full 8–13 lm
band region suitable for solar radiation shield devices have been prepared by
chemical vapor deposition method (Tian et al. 2006). Tellurium is used in
photocopiers to enhance picture quality.
Technology is reshaping our day-to-day life. Metals and their alloys make today’s
manufacturing industry, agriculture, construction and communication systems,
transportation, defense equipments, etc., possible. Our manufacturing industries
are using different metals and alloys as raw materials for their finished goods.
There are a large number of possible combinations of different metals and
each has its own specific set of properties. Alloying one metal with other metal or
non-metal often enhances its properties. For example, steel is stronger than iron, its
primary element. The physical properties, such as density, reactivity, Young’s
modulus and electrical and thermal conductivity, of an alloy may not differ greatly
from those of its elements, but engineering properties, such as tensile strength
(Francis 2008) and shear strength may be substantially different from those of the
constituent materials. This is sometimes due to the size of the atoms in the alloy,
since larger atoms exert a compressive force on neighboring atoms, and smaller
atoms exert a tensile force on their neighbors, helping the alloy resist deformation.
Sometimes alloys may exhibit marked differences in Behaviour even when small
amounts of one element occur (Hogan 1969; Zhang and Suhl 1985). Some of the
major reasons for the continuing advances in alloys are the availability of
materials, new manufacturing techniques, and the ability to test alloys before they
are ever produced. Most modern alloys are, in fact, preplanned using sophisticated
computer simulations, which help determine what properties the alloy will display.
14 1 Introduction
Structural materials employed in reactor systems must possess suitable nuclear and
physical properties and must be compatible with the reactor coolant under the con-
ditions of operation. The most common structural materials employed in reactor
systems are stainless steel and zirconium alloys. Zirconium alloys have favourable
nuclear and physical properties (Kutty et al. 1999), whereas stainless steel has
favourable physical properties. Aluminium is widely used in low-temperature test
and research reactors; zirconium and stainless steel are used in high-temperature
power reactors. Zirconium is relatively expensive, and its use is therefore confined to
applications in the reactor core where neutron absorption is important.
Alloy wheels are automobile (car, motorcycle and truck) wheels which are made
from an alloy of aluminium or magnesium (or sometimes a mixture of both). They
are typically lighter for the same strength and provide better heat conduction and
improved cosmetic appearance. Lighter wheels can improve handling by reducing
vehicle mass which helps to reduce fuel consumption. Better heat conduction can
help dissipate heat from the brakes, which improves braking performance in more
demanding driving conditions and reduces the chance of brake failure due to
overheating (Nunney 2006). Alloy wheels are also purchased for cosmetic pur-
poses as the alloys used are largely corrosion-resistant. This permits the use of
attractive bare-metal finishes, with no need for paint or wheel covers, and the
manufacturing processes allow intricate, bold designs. Magnesium alloy wheels, or
‘‘mag wheels’’ are sometimes used on racing cars, in place of heavier steel or
aluminium wheels, for better performance.
Some of the commercially and industrially important metals and alloys have been
analyzed in this work and their significance are discussed briefly in the following
lines.
Intermetallic alloys, including CoAl and NiAl, are of great importance since these
materials not only have good strength-to-weight ratio but also has excellent
corrosion and oxidation resistance, which make them good candidates for high-
1.6 Significance of the Alloys Dealt in this Research Work 15
temperature and soft magnetic applications (Wan et al. 2010). Among the
intermetallics, Cobalt aluminides are of a considerable technological interest for
high-temperature applications (Lee et al. 2004). Cobalt forms a stable B2 (CsCl)
structure with Al in a wide concentration range (Botton et al. 1996). This structure
can be considered as two interpenetrating primitive cubic sublattices, where each
Co atom has eight Al atoms as nearest neighbors and vice versa. Depending on the
concentration, CoAl alloys exhibit different mechanical and magnetic properties
(Kudryavtsev et al. 1998). CoAl is especially attractive for epitaxial growth due to
its low lattice mismatch with respect to GaAs (Wan et al. 2010). In combination
with semiconductor hetero structures, they are able to be integrated in a number of
devices, e.g., spin injectors and mirrors.
Cobalt-based super alloys consume most of the produced cobalt. The
temperature stability of these alloys makes them suitable for use in turbine
blades for gas turbines and jet aircraft engines. Cobalt-based alloys are also
corrosion and wear-resistant. Cobalt-based alloys are used in total joint
replacement applications where high strength and corrosion resistance are
necessary. Special cobalt–chromium–molybdenum alloys are used for prosthetic
parts such as hip and knee replacements (Michel et al. 1991). Cobalt alloys are
also used for dental prosthetics, where they are useful to avoid allergies to
nickel. Some high-speed steels also use cobalt to increase heat and wear-
resistance. The special alloys of aluminium, nickel, cobalt and iron, known as
Alnico, are used in permanent magnets (Luborsky et al. 1957).
Nichrome wire, an alloy typically made of 80% nickel and 20% chrome, has
specific characteristics that make it ideal for certain uses. It is silvery-grey in color,
is corrosion resistant and has a high-melting point of about 1,400C. Its ability to
conduct electricity, withstand heat, its strength, flexibility and melting point
determine its uses and applications.
Nichrome wire is used to make coils for heating elements. Due to its relatively
high resistivity and resistance to oxidation at high temperatures, it is widely used
in heating elements, such as in hair dryers, electric ovens, toasters and irons.
Industrial uses include a wide variety of heating elements for ovens of all sizes,
and devices that must apply heat directly to surfaces, such as for sealing plastic
packages. Medical laboratories use nichrome wire loops somewhat like small
spoons for handling specimens because they withstand frequent and repeated
sterilisation, then cool rapidly. In industrial use, fine nichrome wire meshes serve
as filters for liquids at high temperatures. Since the wire heats up rapidly with even
a small amount of electricity, it makes safe igniters (sometimes called electric
matches) to light the fuses of fireworks and model rockets. A piece of nichrome
wire about an inch-and-a-half long connected to a length of lead wire allows the
operator to safely detonate the device from a distance with only a battery and a
switch. Fireworks professionals use nichrome igniters with timers at large displays.
NiCr thin films are widely used in several applications in microelectronics such as
thin film resistors, filaments and humidity sensors because of their relatively large
resistivity, more resistant to oxidation and a low-temperature coefficient of resis-
tance (Kwon et al. 2005). Ni80Cr20 matrix offers excellent high-temperature oxi-
dation/corrosion resistance and essential mechanical strength (Ding et al. 2007).
For dental clinical applications, the NiCr-based casting alloys are developed as an
alternative to gold-based alloys due to their corrosion resistance in oral environ-
ment (Huang 2003). Nichrome Ni80Cr20 wt% films are used for strain gauge
applications because of their high resistivity, low-temperature coefficient of
resistance, commercial availability and low-temperature dependence of gauge
factor (Chen et al. 2008).
1.6 Significance of the Alloys Dealt in this Research Work 17
1.6.4 Iron–Nickel
Iron group metals and binary alloys have a number of important industrial
applications (Fabio and Mascaro 2006). Fe–Ni alloys have attracted much atten-
tion due to their interesting mechanical and magnetic properties (Karayannis and
Moutsatsou 2006). Permalloy is a Fe–Ni alloy used in soft magnetic read/write
heads (Cacciamani et al. 2010). Fe–Ni-based alloy powders are interesting in their
applications as soft magnetic materials with low coercivity and high permeability
(Gheisari et al. 2009).
One of the important challenges of electrical motors is to increase the efficiency
to obtain more power with the same electrical energy consumed. The use of iron–
nickel alloys increases the efficiency by reducing the magnetic losses significantly
(Frederic et al. 2007). Fe–Ni thin films and multilayered structures are
materials for high-frequency devices, such as inductors or magneto impedance
effect (MI)-based magnetic field sensors (Kurlyandskaya et al. 2011).
Typical applications that are based on the low coefficient of thermal expansion
of Fe–Ni alloys include thermostatic bimetals, glass sealing, integrated circuit
packaging, cathode ray tube shadow masks, composite molds/tooling and mem-
branes for liquid gas tankers (Cacciamani et al. 2010). Applications based on the
soft magnetic properties include read-write heads for magnetic storage, magnetic
actuators, magnetic shielding and high-performance transformer cores. Due
to their unique low coefficient of thermal expansion and soft magnetic properties,
Fe–Ni alloys are used in several industrial applications (McCrea et al. 2003).
Silver alkali halides provide interesting model systems for the study of decom-
position processes in ionic solids. This is not only due to the absence of structural
phase transitions but also due to the invariance of the anion sublattice, which is not
involved in the demixing process. The phase separation is entirely confined to the
cationic system and the anions exhibit an almost rigid frame. Along with the strong
polarizability of silver ions, this feature guarantees that even single crystals are not
destroyed during demixing (Elter et al. 2005). The dynamics of Na1-xAgxCl show
that in the homogeneous phase, the doping of NaCl with silver ions leads to a
considerable softening of the lattice (Caspary et al. 2007).
Semiconductors in confined surroundings experience a very significant devel-
opment considering their importance in optoelectronic technologies. The semi-
conductors concerned are those having a direct and wide band gap conferring them
with a radiative character and consequently a rather considerable output of
photoluminescence, contrary to indirect gap semiconductors with a weaker output.
Photosensitivity is observed in ionic semiconductors such as sodium chloride
doped with silver (Madani et al. 2004).
18 1 Introduction
A ball mill is a type of grinder and a cylindrical device used in grinding (or
mixing) materials like ores, chemicals and ceramic raw materials. Ball mills rotate
around a horizontal axis, partially filled with the material to be ground and the
grinding medium (Fig. 1.2). Different materials are used as media, including
ceramic balls, alumina balls and stainless steel balls. An internal cascading effect
reduces the material to a fine powder. Industrial ball mills can operate continu-
ously fed at one end and discharged at the other end. Large-to-medium sized ball
mills are mechanically rotated on their axis, but small ones normally consist of a
cylindrical capped container that sits on two drive shafts (pulleys and belts are
used to transmit rotary motion) (Rajagopal 2009).
Ball mills are also used in pyrotechnics and the manufacture of black powder,
but cannot be used in the preparation of some pyrotechnic mixtures such as flash
powder because of their sensitivity to impact. High-quality ball mills are poten-
tially expensive and can grind mixture particles to as small as 5 nm, enormously
increasing surface area and reaction rates. The grinding works on principle of
critical speed. The critical speed can be understood as the speed after which the
balls (which are responsible for the grinding of particles) start rotating along the
direction of the cylindrical device; thus causing no further grinding (Malik and
Singh 2010). Ball mills are used extensively in the Mechanical alloying process in
which they are not only used for grinding but also for cold welding as well, with
the purpose of producing alloys from powders.
The mechanism for formation of nano crystalline materials by the ball milling
technique has been summarised as the phenomenology of the grain size reduction
into three stages (Fecht 1995).
First stage: Plastic deformation is produced in the crystal lattices of the ball-
milled powders by slip and twinning. This deformation is localised in shear bands
containing a high-dense network of dislocation. The local shear instability of a
20 1 Introduction
crystal lattice can be triggered by the material heterogeneity and enhance insta-
bilities. These instabilities result from a non-uniform heat transfer during the
mechanically induced deformation of the milled powders. During this stage of
milling, the atomic level strain increases as a result of increasing the dislocation
density.
Second stage: Due to the successive accumulation of the dislocation density,
the crystals are disintegrated into sub-grains that are initially separated by low-
angle grain boundaries. The formation of these sub-grains is attributed to the
decrease of the atomic level strain.
Third stage: Further ball milling time leads to further deformation occurring in
the shear bands located in the unstrained parts of the powders which leads to sub-
grain size reduction so that the orientation of final grains become random in
crystallographic orientations of the numerous grains and hence, the direction of
slip varies from one grain to another.
In ball mills, the useful kinetic energy can be applied to the powder particles of
the reactant materials by
• Collision between the balls and the powders (Fig. 1.2)
• Pressure loading of powders pinned between milling media or between the
milling media and the liner
• Impact of the falling milling media
• Shear and abrasion caused by dragging of particles between moving milling
media
• Shock wave transmitted through crop load by falling milling media.
There are many types of milling media suitable for use in a ball mill, each material
having its own specific properties and advantages (Rajagopal 2009). Common in
some applications are stainless steel balls. While usually very effective due to their
high density and low contamination of the material being processed, stainless steel
balls are unsuitable for some applications, such as black powder and other flam-
mable materials require non-sparking lead, antimony, brass or bronze grinding
media. In some applications ceramic or flint grinding media is used. Ceramic
media are also very resistant to corrosive materials. High-density alumina media
are widely used to grind clay bodies, frits, glazes and other ingredients. It is more
expensive than silica media but is more efficient.
1.7 Ball Milling 21
The laboratory ball mill is a low-energy ball mill which is used in the present
study. It is less expensive and operates with minimum maintenance requirements.
The vial is a cylinder made of stainless steel with a radius of 2.5 cm. Stainless steel
balls of different radii have been used as the milling media. The rotation speed is
about 200 rpm. It produces homogenous and uniform powders. The laboratory ball
mill used in our work is shown in Fig. 1.3.
The modern understanding of metals and alloys, their structures, defects and
various properties would not be possible if their crystal structures had not been
revealed by XRD studies. XRD is a non-destructive technique for analyzing a wide
range of materials, including crystals, metals, minerals, polymers, thin film coating
and ceramics. Crystalline materials are characterised by the orderly, periodic
arrangements of atoms. The unit cell is the basic repeating unit that defines a
crystal. Parallel planes of atoms intersecting the unit cell are used to define
directions and distances in the crystal. These crystallographic planes are identified
by miller indices (Warren 1990). Diffraction from different planes of atoms
produces a diffraction pattern, which contains information about the atomic
arrangement within the crystal.
XRD is based on constructive interference of monochromatic X-rays and a
crystalline sample. These X-rays are generated by a cathode ray tube, filtered to
produce monochromatic radiation, collimated to concentrate and directed towards
the sample. The interaction of the incident rays with the sample produces
constructive interference (and a diffracted ray) when conditions satisfy Bragg’s
Law (nk = 2d sin h) (Bragg 1912). This law relates the wavelength of
22 1 Introduction
Classically, the two main ways of studying metals and alloys were metallography
(the examination of polished and etched surfaces) and cooling curves (looking for
discontinuities that indicated some sort of phase change). Both these methods
involved considerable skill and experience, and the results were not always
unambiguous. The introduction of XRD provided a much clearer, simpler and
more objective way of investigation. XRD is now a common technique for the
study of crystal structures and atomic spacing.
In terms of the specimen handled, two methods can be identified,
1. Powder diffraction
2. Single-crystal diffraction
In the former, the specimen is a collection of crystallites. Since these fragments
are completely randomly oriented the incident X-ray beam meets with every
possible lattice plane, oriented in all directions. Whereas in the single-crystal
method, the whole specimen is a single piece, without any discontinuity in the
lattice arrangements (Tareen and Kutty 2001). All diffraction methods are based
on generation of X-rays in an X-ray tube. These X-rays are directed to the sample,
and the diffracted rays are collected. A key component of all diffraction is the
angle between the incident and diffracted rays.
1.8 X-Ray Diffraction 23
1.8.2 Diffractometers
The powder method essentially has two different ways of registering the diffracted
X-rays (Tareen and Kutty 2001).
1. The whole diffraction pattern is recorded simultaneously on a photographic film
called the powder photographic method.
2. The diffraction pattern is scanned by a counter device, or a solid-stated semi-
conductor detector. The counter or detector registers the diffracted beam in
successive stages, away from the direct beam.
In X-ray diffractometers, X-rays are generated in a cathode ray tube by
heating a filament to produce electrons, accelerating the electrons towards a
target by applying a voltage, and bombarding the target material with electrons.
When electrons have sufficient energy to dislodge inner shell electrons of the
target material, characteristic X-ray spectra are produced. These spectra consist
of several components, the most common being Ka and Kb. Ka consists, in part,
of Ka1 and Ka2. Ka1 has a slightly shorter wavelength and twice the intensity as
Ka2 (Stout and Jensen 1989). The specific wavelengths are characteristic of the
target material (Cu, Fe, Mo, and Cr). The monochromatic radiation required for
the powder method is usually the Ka1a2 doublet, monochromated by crystal
reflection or by the use of a filter whose K absorption wavelength falls between
the Ka and the Kb wavelengths. The Cu Ka1a2 doublet (k = 1.542 Å ´ ) with a
´
Ni filter kk = 1.488 Å is probably used more than any other source (Warren
1990). These X-rays are collimated and directed onto the sample. The sample is
mounted on a goniometer (Shirane et al. 2002) and gradually rotated while
being bombarded with X-rays, producing a diffraction pattern of regularly
spaced spots known as reflections. As the sample and detector are rotated, the
intensity of the reflected X-rays is recorded. When the geometry of the incident
X-rays impinging the sample satisfies the Bragg Equation, constructive inter-
ference occurs and a peak in intensity is seen. A detector records and processes
this X-ray signal and converts the signal to a count rate which is then output to
a device such as a printer or computer monitor.
24 1 Introduction
The foremost essential criteria for single-crystal XRD are to obtain an adequate
crystal of the material under study. The crystal should be sufficiently large (typ-
ically larger than 0.1 mm in all dimensions), pure in composition and regular in
structure, with no significant internal imperfections such as cracks or twinning.
The crystal is placed in an intense beam of X-rays, usually of a single wave-
length (monochromatic X-rays), producing the regular pattern of reflections.
Molybdenum is the most common target material for single-crystal diffraction,
with MoKa radiation = 0.7107Å. These X-rays are collimated and directed onto
the sample. When the geometry of the incident X-rays impinging the sample
satisfies the Bragg Equation, constructive interference occurs. A detector records
and processes this X-ray signal and converts the signal to a count rate which is
then output to a device such as a printer or computer monitor. Modern single-
crystal diffractometers use CCD (charge-coupled device) technology to transform
the X-ray photons into an electrical signal which are then sent to a computer for
processing.
Single-crystal diffractometers use either three- or four-circle goniometers.
These circles refer to the four angles (2h, v, u, and X) (Shirane et al. 2002) that
define the relationship between the crystal lattice, the incident ray and the detector
as shown in Fig. 1.5. Samples are mounted on thin glass fibers using an epoxy or
cement. The thin glass fiber is attached to brass pins and mounted onto goniometer
heads. The goniometer head and sample are then affixed to the diffractometer.
Samples are centered by viewing the sample under an attached microscope and are
placed under the cross-hairs for all crystal orientations. Once the crystal is centred,
a preliminary rotational image is often collected to screen the sample quality and
to select parameters for later steps. As the crystal is gradually rotated, previous
reflections disappear and new ones appear, the intensity of every spot is recorded
at every orientation of the crystal. Multiple data sets may have to be collected, with
each set covering slightly more than half a full rotation of the crystal and typically
containing tens of thousands of reflections (Massa 2004).
An automatic collection routine can then be used to collect a preliminary set
of frames for determination of the unit cell. Reflections from these frames are auto-
indexed to select the reduced primitive cell and calculate the orientation matrix
(which relates the unit cell to the actual crystal position within the beam). These
data are combined computationally with complementary chemical information to
produce and refine a model of the arrangement of atoms within the crystal and
converted to the appropriate crystal system and Bravais lattice (Giacovazzo 2002).
The final, refined model of the atomic arrangement is the essential crystal
structure.
1.8 X-Ray Diffraction 25
After the data have been collected, corrections for instrumental factors,
polarisation effects of X-ray absorption (Glusker et al. 1994) and (potentially)
crystal decomposition must be applied to the entire data set. This integration
process also reduces the raw frame data to a smaller set of individual integrated
intensities. These correction and processing procedures are typically part of the
software package which controls and runs the data collection.
Single-crystal XRD is most commonly used for precise determination of a unit
cell, including cell dimensions and positions of atoms within the lattice. Bond-
lengths and angles are directly related to the atomic positions. A single, robust,
optically clear sample, generally between 50 and 250 microns in size is essential.
Data collection generally requires between 24 and 72 h.
The size of the grains in a polycrystalline material has more effect on the prop-
erties of the material, for example, the hardness of a metal or alloy increases with
decrease in the grain size (Cullity and Stock 2001). This dependence of properties
on grain size makes the measurement of grain size important in the control of most
metal forming operations.
In this book the grain size has been reported for some metals which has been
analyzed from XRD. The grain morphology was examined by SEM. The software
GRAIN written by Dr.R.Saravanan was used to estimate approximate grain sizes
from XRD (Saravanan, GRAIN software). The grain size is analyzed using full
26 1 Introduction
width at half maximum of the powder XRD peaks. The Debye–Scherrer formula
given in Eq. 1.1has been used to calculate the particle size.
Kk
s¼ ð1:1Þ
b cos h
where K is the shape factor, k is the X-ray wavelength, typically 1.54 Å, b is the
line broadening at half the maximum intensity (FWHM) in radians and h is the
Bragg angle (Patterson 1939) s is the mean size of the ordered (crystalline)
domains, which may be smaller or equal to the grain size. The dimensionless shape
factor has a typical value of about 0.9, but varies with the actual shape of the
crystallite. The Scherrer equation is limited to nano-scale particles. It is not
applicable to grains larger than about 0.1 lm, which precludes those observed in
most metallographic and ceramographic microstructures.
It is important to realise that the Scherrer formula provides a lower bound on
the particle size. The reason for this is that a variety of factors can contribute to the
width of a diffraction peak; besides particle size, the most important of these are
usually inhomogeneous strain and instrumental effects. If all of these other
contributions to the peak width were zero, then the peak width would be deter-
mined solely by the particle size and the Scherrer formula would apply. If the other
contributions to the width are non-zero, then the particle size can be larger than
that predicted by the Scherrer formula, with the ‘‘extra’’ peak width coming from
the other factors.
Accelerated electrons in SEM carry significant amounts of kinetic energy, and this
energy is dissipated as a variety of signals are produced by electron-sample
interactions when the incident electrons are decelerated in the solid sample. These
signals include secondary electrons (that produce SEM images), backscattered
electrons (BSE), diffracted backscattered electrons (EBSD that are used to
determine crystal structures and orientations of minerals), photons (characteristic
X-rays that are used for elemental analysis and continuum X-rays), visible light
(cathodoluminescence-CL) and heat. Secondary electrons and backscattered
electrons are commonly used for imaging samples, secondary electrons are most
valuable for showing morphology and topography of samples and backscattered
electrons are most valuable for illustrating contrasts in composition in multiphase
samples (i.e., for rapid phase discrimination). X-ray generation is produced by
inelastic collisions of the incident electrons with electrons in discrete orbitals
(shells) of atoms in the sample. As the excited electrons return to lower energy
states, they yield X-rays that are of a fixed wavelength (that is related to the
difference in energy levels of electrons in different shells for a given element).
Thus, characteristic X-rays are produced for each element in a mineral that is
‘‘excited’’ by the electron beam (Goldstein 2003).
SEM analysis is considered to be ‘‘non-destructive’’, i.e., X-rays generated by
electron interactions do not lead to volume loss of the sample, so it is possible to
analyze the same materials repeatedly.
References
Caspary D, Eckold G, Elter P, Gibhardt H, Güthoff F, Demmel F, Hoser A, Schmidt W, Ding CH,
Yang ZM, Zhang HT, Guo YF, Zhoua JN (2007) Compost A 38:48–352
Chen R, Zuo D, Sun Y, Li D, Lu W, Wang M (eds) (2008) Key engineering materials, vol 375.
Trans Tech Publications, Switzerland, pp 690–694
Choudry MS, Dollar M, Eastman JA (1998) Mater Sci Eng A256:25–33
Collins DM (1982) Nature 298:49–51
Crans DC (2000) J Inorg Biochem 80:123–131
Crans DC, Smee JJ, Gaidamauskas E, Yang L (2004) Chem Rev 104:849–902
Cullity BD, Stock SR (2001) Elements of X-ray diffraction, Chap. 14, 3rd edn. Prentice Hall,
Upper Saddle
Ding CH, Yang ZM, Zhang HT, Guo YF, Zhoua JN (2007) Compos Part A 38:348–352
Dobrzánski LA, Tánski T, Čížek L, Brytan Z (2007) J Mater Process Technol 192:567–574
Dobson KD, Hodes G, Mastai Y (2003) Sol Energ Mater Sol Cells 80:283–296
Dorward RC, Pritchett T (1988) Mater Des 9(2):63–69
El-Eskandarany MS (2001) Mechanical alloying for fabrication of advanced engineering
materials, Chap. 1. William Andrew, New York
Elter P, Eckold G, Gibhardt H, Schmidt W, Hoser A (2005) J Phys Condens Matter 17:6559–
6573
Engelhard T, Jones ED, Viney I, Mastai Y, Hodes G (2000) Thin Solid Films 370:101–105
Eveloy V, Ganesan S, Fukuda Y, Wu J, Pecht MG (2005) IEEE Trans Compon Packag Technol
28:884–894
Fabio RB, Mascaro LM (2006) Surf Coat Technol 201:1752–1756
Fecht HJ (1995) NanoStruc Mater 6:33–42
Francis J (ed) (2008) Philosophy of mathematics. Global publishing house, New York, p 56
Frederic B, Thierry W, Herve F (2007) Electrical insulation conference and electrical
manufacturing expo, pp 394–401
Futsuhara M,Yoshioka K,Takai O (1998) Thin Solid Films 322:274–281
Gallagher R, Ingram P (eds) (2001) New coordination science, chemistry for higher tier ingram,
3rd edn. Oxford University Press, New York, p 56
Gaoa Q, Guob JT, Huaib KW, Zhanga JS (2005) Mater Lett 59:2859–2862
Gheisari KH, Javadpoura S, Oh JT, Ghaffari M (2009) J Alloys Compd 472:416–420
Giacovazzo C (2002) Fundamentals of crystallography issue 7 of international union of
crystallography texts Oxford science publication, Chap. 5, 2nd edn. Oxford University Press,
Oxford
Gleiter H (1981) Deformation of polycrystals. In: Hansen N, Horsewell A, Leffers T, Lilholt H
(eds) Proceedings second riso intrnational symposium metallurgy and materials science,
Roskilde, p 15
Gleiter H (1989) Prog Mat Sci 33:223–315
Glusker JP, Lewis M, Rossi M (eds) (1994) Crystal structure analysis for chemists and biologists,
volume 16 of methods in stereochemical analysis. Wiley, New York, p 261
Goldstein J (2003) Scanning electron microscopy and X-ray microanalysis, Chap. 6. Kluwer
Adacemic/Plenum Pulbishers, New York
Hansen NK, Coppens P (1978) Acta Cryst A34:909–921
Hauptmann H, Walter W (1962) Chem Rev 62:347–404
Heinrich G, Sgler T, Rosiwal SM, Singer RF, Ley L (1996) Diam Relat Mater 5:304–307
Hogan CM (1969) Phys Rev 188:870–874
Hort N, Huang Y, Fechner D, Störmer M, Blawert C, Witte F, Vogt C, Drücker H, Willumeit R,
Kainer KU, Feyerabend F (2010) Acta Biomater 6:1714–1725
http://itp.nyu.edu/physcomp/uploads/battery_nimh.jpg)
http://en.wikipedia.org/wiki/File:Vanadinite2_sur_goethite_ (Maroc).jpg
http://image.made-in-china.com/2f0j00DvhtLbGPOekg/Grinding-Ceramic-Balls.jpg
http://www.shubhmets.in/full-images/663342.jpg
Huang HH (2003) Biomaterials 24:1575–1582
References 29
Hunziker W, Voyt W, Melchior H (1996) Proceedings of the 1996 IEEE 46th electronic
components and technology conference, Orlando, Florida, pp 8–12
Inoue A (1994) Mater Sci Eng A 179:57–61
Ivey DG (1998) Micron 29(4):281–287
Josef P, Dana P, Jaroslav B et al (2004) Anal Appl Pyrolysis 71:739–746
Kainer KU, Sb (2000) International congress magnesium alloys and their application, Mnichov,
pp 534–608
Karayannis VG, Moutsatsou AK (2006) J Mater Process Technol 171:295–300
Kartal G, Timur S, Urgen M, Erdemir A (2010) Surf Coat Technol 204:3935–3939
Kazacos MS, Menictas C (1997) Telecommunications energy conference, 19th international,
pp 463–471
Kear BH, Strutt PR (1995) Nanostruct Mater 6:227–236
Kittel C (1996) Introduction to solid state physics, 7th edn. Wiley, India
Kittel C (2007) Introduction to solid state physics, Chap. 6, 7th edn. Wiley, India
Kudryavtsev YV et al (1998) J Appl Phys 83:1575–1581
Kurlyandskaya GV, Bhagat SM, Svalov AV, Fernandez E, Arribas AG, Barandiaran JM (2011)
Solid State Phenom 168:257–260
Kutty TRG, Ravi K, Ganguly C (1999) J Nucl Mater 265:91–99
Kwon Y, Kim N, Choi G, Lee W, Seo Y, Park J (2005) Microelectron Eng 82:314–320
Lee YS, Lee CH, Kim KW, Shin HJ, Lee YP (2004) J Magn Magn Mater 272:2151–2153
Li ZG, Smith DJ (1989) Appl Phys Lett 55:919–921
Li SJ, Cui TC, Hao YL, Yang R (2008) A Biomater 4:305–317
Lide DR (ed) (1999) In: CRC handbook of chemistry and physics, 80th edn. CRC Press,
New York, pp 4–32
Lonyuk B, Apachitei I, Duszczyk J (2007) Surf Coat Technol 201:8688–8694
Lu K, Wei WD, Wang JT (1995) Scr Metall Mater 24:2319–2323
Luborsky FE, Mendelsohn LI, Paine TO (1957) J Appl Phys 28:344–351
Lynch CT (1974) CRC handbook of materials science. Vol 1,CRC press
Madani S, Hocine K, Madjid B, Abdelhamid C, Mohamed AB, Sofiane M (2004) Physica E
23:217–220
Maeshima T, Nishida M (2004) Mater Trans 45:1096–1100
Malik HK, Singh AK (2010) Engineering physics, Tata McGraw-Hill, New Delhi, p 6
Markovsky PE, Semiatin SL (2010) J Mater Process Technol 210:518–528
Massa W (2004) Crystal structure determination, Chap. 6, 2nd edn. Springer, Berlin
McCrea JL, Palumbo G, Hibbard GD, Erb U (2003) Rev Adv Mater Sci 5:252–258
McHenry ME, Laughlin DE (2000) Acta Mater 48:223–238
Mehta DS, Masood SH, Song WQ (2004) J Mater Process Technol 156:1526–1531
Michel R, Nolte M, Reich M, Löer F (1991) Systemic effects of implanted prostheses made of
cobalt–chromium alloys, Archives of orthopaedic and trauma surgery, 110:61–74
Monfared HH, Alavi S, Bikas R, Mahedpour M, Mayer P (2010) Polyhedron 29:3355–3362
Mordike BL, Ebert T (2001a) J Mater Process Technol 117:381–385
Mordike BL, Ebert T (2001) Mater Sci Eng A 302:37–45
Nazmy M, Noseda C, Staubli M (1995) Switzerland, United States (19), Patent Number 5411702,
May 2
Newbury DE (1986) Advanced scanning electron microscopy and X-ray microanalysis, Chap. 2.
Springer, Berlin
Ng C, Simkin BA, Crimp MA (1997) Mater Sci Eng A239:150–156
Nickel Magazine (2007), 22(2) March issue
Nithya D, Kalaiselvan R (2011) Physica B 406:24–29
Noring J et al (1993) In: Proceedings of the symposium on batteries and fuel cells for stationary
and electric vehicle applications vols 93–98 of proceedings (electrochemical society), The
Electrochemical Society, pp 235–236
Nowosielski R, Babilas R, Ochin P, Stoktosa Z (2008) Arch Mater Sci Eng 30:13–16
30 1 Introduction
Nunney MJ (ed) (2006) Light and heavy vehicle technology, 4th edn. Butterworth-Heinmann,
London, p 454
Oshima T, Kajitani M, Okuno A (2004) Int J Appl Ceram Technol 1269–1276
Paik J, Veen SVD, Duran A, Collette M (2005) Thin-Walled Str 43:1550
Patterson A (1939) Phys Rev 56:978–982
Pawlikowski M (1981) Infrared Phys 21:181–187
Petersson IU, Löberg JEL, Fredriksson AS, Ahlberg EK (2009) Biomaterials 30:4471–4479
Pillai SO (2007) Solid state physics, New age international (p) limited, New Delhi, Reprint, p 288
Poth U (2008) Automotive coatings formulation: chemistry, physics and practices European
coatings tech files, Vincentz, p 152
Proffen TH, Billinge SJL (1999) J Appl Cryst 32:572–575
Qin Y, Wilcox GD, Liu C (2010) Electrochim Acta 56:183–192
Rajagopal K (2009) Textbook of Engineering Physics, Part II, Publisher-PHI learning Pvt. Ltd,
New Delhi, p 168
Reimer L (1998) Scanning electron microscopy, physics of image formation and micro analysis,
Chap. 1. Springer, Berlin
Rietveld HM (1969) J Appl Crystallogr 2:65–71
Ryan VW, Schutz RJ, Warren (1993) United States Patent (19), Patent No: 5243221 September 7
Sambamurthy K (2007) Pharmaceutical engineering, new age international, p 427
Saravanan R http:\\www.saraxraygroup.net. Grain software
Scheppe F, Sahm PR, Hermann W, Paul U, Preuhs J (2002) Mater Sci Eng A329:596–601
Schwingel D, Seeliger HW, Vecchionacci C, Alwes D, Dittrich J (2007) Acta Astronaut 61:326–330
Setoyama D, Matsunaga J, Muta H, Uno M, Yamanaka S (2004) J Alloys Compds 381:215–220
Shashwati S, Muthe KP, Joshi N (2004) Sens Actuators B 98:154–159
Shirane G, Shapiro SP, Tranquada JM (2002) Neutron scattering with a triple-axis spectrometer:
basic techniques, Cambridge University Press, Cambridge, p 88
Skyllas-Kazacos M (2003) J Power Sour 124:299–302
Sparrow G (1999) Iron—the elements-group 1, the elements series. Marshall Cavendish, p 4
Stout GH, Jensen LH (1989) X-ray structure determination, Chap. 1, 2nd edn. Wiley-
Interscience, New York
Suryanarayana C (2004) Mechanical alloying and milling, Chap. 4, vol 180 of Mechanical
engineering. Marcel Dekker, New York
Sutou Y, Yamauchi K, Takagi T, Maeshima T, Nishida M (2006) Mater Sci Eng A 438:64–69
Tareen JAK, Kutty TRN (2001) A basic course in crystallography, chemistry series, Universities
Press, p 149
Tian MY, Qing GZ, Hong XW, Jian H (2006) Trans Nonferrous Met Soc China 16:693–699
Tsiulyanu D, Mariana S, Miron V, Liess HD (2001) Sens Actuators B 73:35–39
Tsiulyanu D, Marian S, Liess HD, Eisele I (2004) Sens Actuators B 100:380–386
Walawalkar R, Apt J, Mancini R (2007) Energ Policy 35(4):2558–2568
Wan Q, Hey R, Trampert A (2010) 16th International conference on microscopy of
semiconducting materials, IOP Publishing J Phys Conference Series 209, 012023
Wang RM, Eliezer A, Gutman E (2002) Mater Sci Eng A344:279–287
Warren BE (1990) X-ray diffraction, 2nd edn. reprint, Courier Dover Publications, New York,
p 16
Wolf MJ, Engelmann G, Dietrich L, Reichl H (2006) Nucl Instrum Method Phys Res A 565:290–
295
Zhang X, Suhl H (1985) Phys Rev A 32:2530–2533
Zhang Y, Liu X, Xie G, Yu L, Yi S, Hua M, Huang C (2010) Mater Sci Eng B 175:164–171
Zhao KY, Li CJ, Tao JM, Ng DHL, Zhu XK (2010) J Alloys Compds 504S:306–310
Zumdahl SS (2007) Chemical principles,6th edn. Cengage Publications, p 33
Chapter 2
Charge Density Analysis from X-Ray
Diffraction
2.1 Introduction
R. Saravanan and M. Prema Rani, Metal and Alloy Bonding: An Experimental Analysis, 31
DOI: 10.1007/978-1-4471-2204-3_2, Springer-Verlag London Limited 2012
32 2 Charge Density Analysis from X-Ray Diffraction
correlation of structural features with physical properties. This can include optical,
magnetic, and conductive effects as well as the structural behavior. The resulting
knowledge increases our understanding of the complex, underlying processes,
ultimately aiding the design of new materials in which the desired chemical or
physical properties are enhanced.
The field of crystallography developed to great heights after the discovery of X-rays
by Wilhelm Conrad Röntgen in 1895 (Novelline and Squire 2004). The first X-ray
experiment was performed by Max von Laue in Munich in 1912 (Purrington 1997).
At the time of von Laue’s work, there was some experimental evidence that X-rays
were high energy particles; other data indicated that X-rays might be waves. Von
Laue surmised that, if X-rays were waves, they would have rather short wavelengths
(on the order of 1 9 10-10 m) and the dimensions of the objects in crystals would be
the appropriate size to produce the phenomenon of diffraction. He exposed a crystal
of copper sulfate to X-rays and recorded the diffraction pattern on a piece of pho-
tographic film, the experimental setup similar to that shown in Fig. 2.1.
This experiment proved the wave nature of X-rays and began the exploration of
molecular structure by X-ray diffraction methods. Modern X-ray data collection is
substantially the same as von Laue’s experiment though various kinds of electronic
detectors are used rather than photographic film.
An X-ray photo contains two important pieces of information about the crystal
structure, the shape and size of the unit cell and the distribution of electron density
throughout the unit cell. A useful explanation of von Laue’s experiment was
formulated by Bragg (1913) in Cambridge, in 1912. Bragg considered the dif-
fractions to arise from ‘reflections’ of separate X-ray waves from parallel planes of
electron density.
X-ray diffraction results from the interaction between X-rays and electrons of
atoms. Depending on the atomic arrangement, interfaces between the scattered
rays are constructive when the path difference between two diffracted rays differs
by an integral number of wavelengths. This selective condition is described by
Bragg’s equation also called ‘‘Bragg’s law’’ (Bragg 1912).
Crystals consist of a periodic arrangement of atoms or molecules that form a
crystal lattice. In such an arrangement of atoms numerous planes run in different
directions through the lattice points (atoms, molecules), horizontally, vertically,
and diagonally which are the lattice planes. All the planes parallel to a lattice plane
are also lattice planes and are set at a distance apart from each other which is
2.2 X-Ray Diffraction 33
called the lattice plane distance ‘d’. Bragg diffraction occurs when electromagnetic
radiation or subatomic particle waves with wavelength comparable to atomic
spacings are incident upon a crystalline sample, scattered in a specular fashion by
the atoms in the system, undergo constructive interference. When the scattered
waves interfere constructively, they remain in phase since the path length of each
wave is equal to an integer multiple of the wavelength. The path difference
between two waves undergoing constructive interference is given by 2d Sin h,
where h is the scattering angle. This leads to Bragg’s law which describes the
condition for constructive interference from successive crystallographic planes
(h,k,l) of the crystalline lattice (Azaroff 1968).
In Fig. 2.2, two monochromatic X-ray beams of wavelength ‘k’ are represented
by wave 1 and wave 2. The spacing between the atomic planes occurs over the
distance, d. Wave 2 reflects off at the upper atomic plane at an angle h equal to its
angle of incidence. Similarly, wave 1 reflects off at the lower atomic plane at the
same angle h. although wave 1 is in the crystal, it travels a distance of 2d farther
than wave 2. Hence the path difference between waves 1 and 2 is 2d. If this
distance 2d is equal to an integral number of wavelengths nk, then waves 1 and 2
will be in phase on their exit from the crystal and constructive interference will
occur. If the distance 2d is not an integral number of wavelengths, then destructive
interference will occur and the waves will not be as strong as when they entered
the crystal. Thus, the condition for constructive interference to occur is nk = 2d.
In figure, BC ¼ d sin h; BD ¼ d sin h
The path difference ¼ BC þ BD
ð2:1Þ
BC þ BD ¼ 2d sin h
Thus, nk ¼ 2d Sin h
This is known as Bragg’s Law (Bragg 1913) for X-ray diffraction.
On the basis of Bragg’s Law, by measuring the angle ‘h’, the wavelength ‘k’,
the chemical elements can be determined, if the lattice plane distance ‘d’ is known,
or, if the wavelength ‘k’ is known, the lattice plane distance ‘d’ and thus the
crystalline structure can be determined. This provides the basis for two measuring
34 2 Charge Density Analysis from X-Ray Diffraction
Fig. 2.2 Schematic representation of X-ray reflection from two parallel planes of a crystal
W $ W2 ¼ q ð2:2Þ
wave function probability
ðprobability amplitudeÞ ðor electron densityÞ
A ‘particle’ electron can also be described by a ‘wave’ function or orbital. The
orbital, whether atomic or molecular, covers a region of space and does not move.
A moving electron looks like a stationary cloud of charges. The physical
2.2 X-Ray Diffraction 35
The crystal structure, which is determined by the lattice parameters and the atomic
positions within the unit cell, is an essential element for the diffraction intensities.
Each atom scatters according to its electronic distribution following dependence on
sin h/k. When there is more than one atom per unit cell, the interference between the
waves scattered by each atom is considered as the diffraction unit. The addition of
the waves scattered by each of the atoms of the unit cell, considering that there is
coherence in time and space is called the structure factor (Marin and Dieguez 1999).
The structure factor Fhkl is the resultant of j waves scattered in the direction of
the reflection hkl by the j atoms in the unit cell (Stout and Jensen 1968).
The structure factor may be expressed as
X
Fhkl ¼ Fhkl expðiahkl Þ ¼ fj exp 2pi hxj þ kyj þ lzj ð2:3Þ
j
X X
Fhkl ¼ fj cos 2p hxj þ kyj þ lzj þ i fj sin 2p hxj þ kyj þ lzj ð2:4Þ
j j
where the sum is over all atoms in the unit cell, xj, yj, zj are the positional
coordinates of the jth atom, fj is the scattering factor of the jth atom, and ahkl is the
phase of the diffracted beam.
The scattering factor of the atom is,
amplitude of the radiation scattered from the atom
fj ¼ ð2:6Þ
amplitude of the radiation scattered from a single electron
The structure factor describes the way in which an incident beam is scattered by
the atoms of a crystal unit cell, taking into account the different scattering power of
the elements through the term fj. Since the atoms are spatially distributed in the
unit cell, there will be a difference in phase when considering the scattered
amplitude from two atoms. This phase shift is taken into account by the complex
exponential term. The atomic form factor or scattering power, of an element
depends on the type of radiation considered as electrons interact with matter
through different processes.
If the summation over discrete atoms in the following structure factor
expression
36 2 Charge Density Analysis from X-Ray Diffraction
X
Fhkl ¼ fj e2piðhxj þkyj þlzj Þ ð2:7Þ
This implies that the electron density is the Fourier Transform of the Structure
Factor. Likewise, the Structure Factor is the Fourier Transform of the electron
density.
Z
qð x:y:zÞ ¼ Fhkl e2piðhxþkyþlzÞ dV ð2:9Þ
V
According to Eq. 2.10 the electron density, q, can be calculated at any point (x,
y, z) by constructing a Fourier Series which has coefficients that are equal to the
Structure Factors (Warren 1990). This is the basic equation of crystallography.
It enables the calculation of a three-dimensional electron density map
throughout the unit cell. Maxima in the electron density map define the locations
of individual atoms.
Once a structure solution has been achieved, there are actually two structure
models, a calculated model based on the approximate co-ordinates obtained from
interpretation of a three-dimensional electron density map,
X
Fhkl ¼ fj e2piðhxj þkyj þlzj Þ ð2:11Þ
and an observed model based on the calculation of structure factors from exper-
imental intensities,
pffiffiffiffiffiffi
jFhkl j ¼ K Ihkl ð2:12Þ
The task now is to adjust the various atomic parameters so that the calculated
structure factors match the observed structure factors as closely as possible. One of
the ways to measure the agreement between the observed and calculated models is
with the Residual Index (more commonly referred to as the R-factor)
P
jFobs Fcalc j
R¼ P ð2:13Þ
jFobs j
In addition to the approximate atomic x, y, z-coordinates, there is another factor
influencing the magnitudes of the structure factors, the thermal motion. Although
the structure determination is done on the solid phase, the atoms still have some
thermal motion (vibration and rotation). This thermal effect is introduced into the
structure factor equation by a factor that serves to attenuate the atomic scattering
factor,
ðsin2 hÞ
fj ¼ fo eB k2 ð2:14Þ
B ¼ 8p2 u2 ð2:15Þ
As B increases, the scattering power of atom j(fj) decreases. To represent
anisotropic thermal motion, a total of six parameters per atom are necessary. Three
of these parameters provide the orientations of the principal axes of the ellipsoid
produced by anisotropic thermal motion; the other three parameters represent the
magnitudes of displacement along the ellipsoid axes.
With thermal effects included, the complete expression for the calculated
structure factor is
38 2 Charge Density Analysis from X-Ray Diffraction
X ðsin2 hÞ
Fhkl ¼ fj e2piðhxj þkyj þlzj Þ eB k2 ð2:16Þ
In order to refine the structure, the coordinates and temperature factors must be
adjusted so that these calculated Fhkl’s match, as closely as possible, the observed
Fhkl’s derived from the experimentally measured intensities.
In terms of the specimen handled, two methods can be identified, viz., the single
crystal method and the powder method.
1. Single-crystal X-ray diffraction
2. Powder X-ray diffraction
The single crystal diffraction technique, using relatively large crystals of the
material, gives a set of separate data from which the structure can be obtained.
However, most materials of technical interest cannot grow large crystals, so one
has to resort to the powder diffraction technique using material in the form of very
small crystallites. In reality, the two techniques are highly complementary, have
their own strengths and weaknesses and domains of applicability, and one will
never supercede the other. However, the domain where powder diffraction is
having an impact is certainly growing and diversifying. In general, the structure
solution process poses more significant challenges, which has prompted much
research in recent years on the development of new strategies and techniques in
this field.
2.4 Methods in X-Ray Crystallography 39
300
200
100
result in R-values lower than 2.5%. wR refers to squared F-values. GooF refers to
the ‘‘goodness of fit’’ of the solution. In addition to the difference in F values, the
GooF also takes into account the number of observed reflections and the param-
eters used. At the end of refinement, the GooF should approach 1.
The possibility of using powder diffraction methods to study materials was rec-
ognized shortly after the discovery of X-ray diffraction by Laue and Von Knipping
in 1910. Although, the powder method was developed as early as 1916 by Debye
and Scherrer, for more than 50 years its use was almost exclusively limited to
qualitative and semi-quantitative phase analysis and macroscopic stress mea-
surements. In 1922, Bain using the powder photograph method of Debye and Hull
carried out some very important pioneering work (Weiss 1925). Within a few
years many others including Bragg and Pauling, had exploited the powder method
to study a wide range of materials, including metals, minerals, and simple organic
solids. Quantitative analysis of the pattern using modern computers and software
yields a wealth of additional information about the sample structure. Figure 2.3
shows a typical powder diffraction profile.
Powder diffraction has played a central role in structural physics, chemistry,
and materials science over the past 20 years. Important advances in structural
studies of materials ranging from high temperature superconductors and fullerenes
to zeolites and high-pressure research have relied heavily on the powder diffraction
technique.
Although most of the structures determined from powder diffraction data have
been solved in the last few years, notable contributions have been from Zacha-
riasen (1948) and Berg and Werner (1977). Zachariasen, in particular, used a
number of ingenious methods to solve crystal structures from powders. His work
on b-plutonium (Zachariasen and Ellinger 1963), for example, utilized differential
2.4 Methods in X-Ray Crystallography 41
thermal expansion to resolve Bragg peak overlap. Many of the early zeolite
structures that were solved from powder diffraction data involved model building
and significant chemical intuition (Breck et al. 1956; Kokatalio et al. 1978), and
these concepts are now implemented in computer algorithms.
Solving a structure from powder diffraction data has developed rapidly over the
last 10 years. Prior to 1990, very few unknown crystal structures had been
determined directly from powder diffraction data. Today, numerous crystal
structures, both organic and inorganic, have been solved from powder data.
Developments in instrumentation, computer technology, and powder diffraction
methodology have all contributed to this increased success rate.
The first choice is the selection of radiation source and instrumentation geometry.
The first step in data analysis is indexing of the diffraction pattern. From the unit
cell and space group the diffraction intensity associated with each reflection should
be determined by applying whole-profile fitting. The intensity extraction can be
performed using either a least-squares method (Pawley 1981) or an iterative
approach (Le Bail et al. 1988).
Rietveld proposed a simple summation approach to the evaluation of an
observed structure factor magnitude for partially and completely overlapped
reflections, if the calculated diffraction pattern is good but the observed and cal-
culated structure factor magnitudes are different from one another.
The peak area is proportional to the square of the structure factor magnitude and
the problem reduces to finding the peak area. For an isolated peak, the observed
peak area is easy to evaluate. For overlapping peaks, the contribution for a given
reflection weighted by the calculated peak contribution for that reflection is
divided by the sum of the calculated peak values for each contributing reflection.
The estimated standard deviations of observed intensities are not normally
included as part of the standard Le Bail approach, they may nevertheless be
evaluated.
Pawley (1981) declared a method for determining Bragg peak intensities from
powder diffraction data in the absence of structural model. The variables associ-
ated with peak positions and widths are the same, the variables associated with
peak areas are peak areas themselves. Although the Pawley method was intro-
duced some 6 years earlier than its counterpart, LeBail method is currently still the
more popular approach.
The final step in the structure determination maze is the completion of the
structure and the refinement of the structural parameters using the Rietveld method
(Rietveld 1969). The structural proposal from the structure determination can be
confirmed only when the refinement has been brought to a successful conclusion.
Chemical information and intuition play important roles in guiding the path in the
maze.
42 2 Charge Density Analysis from X-Ray Diffraction
The drawback of the conventional powder method is that the data grossly overlap,
thereby preventing proper determination of the structure. The major breakthrough
in the value of the powder method as a quantitative tool was the development of
the Rietveld method (Rietveld 1969), a technique for crystal structure refinement
which, for the first time, made use of the entire powder pattern instead of analyzing
individual, non-overlapped, Bragg reflections separately. This approach minimizes
the impact of overlapped and degenerate peaks by calculating the entire powder
pattern of a crystalline model, including various experimental and sample
dependent peak-broadening effects. Parameters in the model such as atomic
positions, lattice parameters, and experimental factors that affect peak-shape and
background are varied, using a least-squares approach, until the agreement
between the calculated and measured diffraction profiles are optimized and the
model is refined by iterative procedure.
The method was first reported for the diffraction of monochromatic neutrons
where the reflection-position is reported in terms of the Bragg angle 2h. The
method was quickly extended from reactor neutron data, with its nice Gaussian
line-profiles and lack of atomic form-factor, to in-house X-ray powder diffraction,
synchrotron powder diffraction, and time-of-flight neutron data from pulsed
spallation sources, and to refinements of incommensurate and magnetic structures.
The application of the Rietveld method to neutron data in the early 1970s was soon
followed by its extension to laboratory X-ray diffractometer data (Malmros and
Thomas 1977; Young et al. 1977). The technique has been applied to a wide range
of solid-state problems and has been reviewed by several authors during the last
25 years (Cheetham and Taylor 1977; Hewat 1986).
The problem of the more complex peak shape was resolved by employing
alternative peak-shape functions, such as the Lorentzian and the pseudo-Voigt.
Other problems that can plague X-ray studies include preferred orientation and
poor powder averaging (graininess), both of which arise from the fact that X-rays
probe a smaller sample volume than neutrons; these were addressed by paying
closer attention to the data collection strategy. The accuracy and precision of a
structure refinement from X-ray data can normally be optimized by collecting
high-resolution data at a synchrotron source (Cox et al. 1983). The resolution of
the powder diffractometers at second and third generation sources is so good that
sample imperfections now play a major role in determining the shape of the Bragg
peaks.
Powder diffraction patterns are collected in a step scan mode. Intensity and two-
theta axes of Bragg diffraction are collected for the crystal structure determination
(Young 1995). The only wavelength and technique independent scale is in
2.5 The Rietveld Method 43
where Hk is the full width at half peak height (full-width half-maximum), 2hk is the
center of the reflex, and Ik is the calculated intensity of the reflex (determined from
the structure factor, the Lorentz factor, and multiplicity of the reflection).
At very low diffraction angles the reflections may acquire an asymmetry due to
the vertical divergence of the beam. Rietveld used a semi-empirical correction
factor, As to account for this asymmetry.
" #
sPð2hi hk Þ2
AS ¼ 1 ð2:21Þ
tan hk
The width of the diffraction peaks are found to broaden at higher Bragg angles.
This angular dependency (Caglioti et al. 1958) was originally represented by
where U, V, and W are the half-width parameters and may be refined during the fit.
44 2 Charge Density Analysis from X-Ray Diffraction
where Iobs is the intensity expected for a random sample, G is the preferred
orientation parameter, and a is the acute angle between the scattering vector and
the normal of the crystallites.
2.5.1.4 Refinement
1 calc 2
obs
M¼ Wi yi yi ð2:24Þ
i
c
where Wi is the statistical weight and c is an overall scale factor such that
ycalc = cyobs.
Rietveld (1969) method calculates the entire powder pattern using a variety of
refinable parameters and improves the selection of these parameters by minimizing
the weighted sum of the squared differences between the observed and the cal-
culated pattern using least square methods. Thus the systemic and accidental peak
overlap is overcome.
The method has been so successful that nowadays the structure of materials, in
the form of powders, is routinely being determined, nearly as accurately as the
results obtained by single-crystal diffraction techniques. An even more widely
used application of the method is in determining the components of chemical
mixtures. This quantitative phase analysis is now routinely used in industries
ranging from cement factories to the oil industry. The success of the method can be
gauged by the publication of more than a thousand scientific papers yearly using it.
The Rietveld refinement was performed using the softwares JANA 2000 (Petříček
et. al. 2006) and an improvised version JANA 2006 (Petříček et. al. 2006). JANA
2.5 The Rietveld Method 45
2000 and JANA 2006 are crystallographic programs for structure analysis of
crystals periodic in three or more dimensions from diffraction data. It is focused to
solution, refinement, and interpretation of difficult, especially modulated struc-
tures. It calculates structures having up to three modulation vectors from powder
as well as single-crystal data measured with X-ray or neutron diffraction. The input
diffraction data can be unlimitedly combined, the combination of powder neutron
data with single-crystal X-ray data being a typical example. The structure solution
can be done using the builtin charge flipping algorithm or by calling an external
direct methods program. Jana can handle multiphase structures (for both powder
and single-crystal data), twins with partial overlap of diffraction spots, commen-
surate, and composite structures. It contains powerful transformation tools for
symmetry (groupsubgroup relations), cell parameters, and commensurate—super
cell relations. A wide scale of constrains and restrains is available including a
powerful rigid body approach and the possibility to define a local symmetry
affecting only part of the structure.
In the refinement using the software, the observed profiles are matched with the
profiles constructed similarly by using pseudo-voigt (Wertheim 1974) profile
shape function of Thompson, Cox and Hastings (1987) modified to some extent
that accommodate various Gaussian FWHM parameters and Scherer coefficient for
Gaussian broadening. The profile asymmetry is introduced by employing multi-
term Simpson rule integration devised by Howard (1982) that incorporates sym-
metric profile shape function with different coefficient for weights and peak shift.
Jana 2006 also employs the correction for preferred orientation which is inde-
pendent of diffraction geometry using March—Dollase function (March 1932;
Dollase 1986). The calculated profiles thus evolved are compared with the
observed ones. Finally, the structure factors evolved from the Rietveld refinements
were further utilized for the estimation of charge density in the unit cell.
The distribution of positive and negative charges in a crystal fully defines physical
properties like the electrostatic potential and its derivatives, the electric field, and
the gradient of the electric field. The electrostatic potential is of importance in the
study of intermolecular interactions, and has received considerable attention
during the past two decades. It plays a key role in the process of molecular
recognition, including drug-receptor interactions, and is an important function in
the evaluation of the lattice energy. The evaluation of the electrostatic potential
and its derivatives may be achieved directly from the structure factors, or indi-
rectly from the experimental electron density by the multipole formalism.
The ability to measure the experimental charge distribution in crystals from the
intensities of the scattered X-rays was realized almost immediately after the dis-
covery of X-ray diffraction. Notwithstanding this early recognition, the technical
developments of the 1960s and beyond, which occurred in difftractometry,
46 2 Charge Density Analysis from X-Ray Diffraction
F ð hÞ ¼
" X X #
atoms X
X symmetry h 4
h l
hp
Pc fcore ðhÞ þ Pv fvalence 00 þ ukl 00 Pklm yklm
k p
k l¼0
k m¼1
h
exp 2pih:rkp Tk ðhÞ
ð2:27Þ
where fcore and fvalence are the Fourier transforms of qcore and qvalence respectively,
and ukl is the Fourier transform of Rkl defined as:
Z1
l
ukl ðhÞ ¼ 4pi Rkl ðr Þjl ð2phr Þ r2 dr ð2:28Þ
0
The temperature factor Tk is the Fourier transform of tk and hp and rlp are
symmetry transformed scattering and position vectors. We note that the model
reverts to the free-atom model when j0 = 1 and the Pklm0 s are zero. It is in this
respect, similar to Hirshfeld’s deformation model (Harel and Hirshfeld 1975), but
different from the single Slater-type models used by a series of authors (Cromer
et.al. 1976; Price and Maslen 1978; Chen et.al. 1977).
Although an aspherical atom model (multipole) refinement based on Eq. 2.25 may
lead to a lowering of agreement indices, the improvement is not necessarily sig-
nificant because of the large number of parameters involved, and the results are not
necessarily physically meaningful as they may be affected by parameter correla-
tion or systematic errors in the measurements. The following tests of the multipole
refinements have therefore been applied.
A successful refinement should give parameters with minimal bias due to atomic
asphericity. The rationale for this requirement is twofold. First, deformation maps
(i.e., total minus ‘sum of spherical atom’ densities) are quite sensitive to errors in
positional and thermal parameters, and second, physical properties to be calculated
from the refinement results, such as dipole and quadrupole moments, depend on
the nuclear positions as well as charge distribution.
One of the justifications for performing the multipole refinement is the experi-
mental determination of derived physical properties (Stewart 1972). Properties
50 2 Charge Density Analysis from X-Ray Diffraction
such as net molecular charges, dipoles, and quadrupole moments may also be
obtained by direct-space integration of the experimental charge distribution.
The raw data is refined for their structural parameters with extinction, absorption,
and TDS corrections. Structure factors are refined to yield the smallest possible
reliability indices. The results are compared for the thermal parameters obtained
using MEM. The refined structure factors are further refined for the population and
k-parameters using the formalism proposed by Hansen and Coppens (1978). In this
refinement, the canonical Hartree–Fock atomic orbitals of the free atoms nor-
malized to one electron are used for the construction of qcore and qvalence charge
densities.
model has significant advantages and its introduction has greatly contributed to the
increasing application of experimental results in charge-density analysis,
1. Experimental noise is generally not fitted by the model functions and therefore
effectively filtered out.
2. Thermal motion is treated separately and de-convoluted from the final result.
3. The resulting static density provides an effective level of comparison with
theoretical results, especially if the latter has been filtered through the model by
refinement of theoretical structure factors.
4. Notwithstanding the development of alternative formulations, including bond-
charge models and orbital-based algorithms, no generally competitive alter-
native has been developed. While the pseudo atom model is widely used in
experimental density analysis, it is important not to lose sight of the implied
assumptions.
5. The results are dependent on the adequacy of the thermal motion formalism
used (Mallinson et al. 1988), which generally is limited to the harmonic
approximation.
The accurate electron density distribution could be obtained if all the structure
factors are known without any ambiguities. It is however impossible to collect
exact values of all the structure factors by X-ray diffraction methods. The number
of observed structure factors by the experiment is always limited and has some
errors. The uncertainties in the results due to the incompleteness of the experi-
mental information must be rectified. Maximum entropy method is one of the
appropriate methods in which the concept of entropy is introduced to handle the
uncertainties properly. The principle of MEM is to obtain electron density dis-
tribution which is consistent with the observed structure factors and to leave the
uncertainties to a minimum. The mathematical formalism of MEM is given in the
following lines.
The probability of a distribution of N identical particles over m boxes, each
populated by ni particles, is given by
N!
P¼ ð2:35Þ
n1 !n2 !n3 !. . .. . .nm !
As in statistical thermodynamics, the entropy is defined as ln P. Since the
numerator is constant, the entropy is, apart from a constant, equal to
X
S¼ ni ln ni ð2:36Þ
i
where both q0 ðr Þ and prior probability s0 ðr Þ are related to the actual electron
density in a unit cell as
qðr Þ sðr Þ
q 0 ðr Þ ¼ P and s 0 ðr Þ ¼ P ð2:40Þ
r qðr Þ r sðr Þ
where qðr Þ and sðr Þ are the electron density and prior electron density at a fixed r
in a unit cell respectively. In the present theory, the actual densities are treated
hereafter instead of normalized densities, and q0 ðr Þ becomes s0 ðr Þ when there is no
information. The q0 ðr Þ and s0 ðr Þ are normalized as
X X
q0 ðr Þ ¼ 1 and s0 ð r Þ ¼ 1 ð2:41Þ
where N is the number of reflections used for MEM analysis, rðH Þ; standard
deviation of Fobs ðH Þ; the observed structure factor and Fcal ðH Þ is the calculated
structure factor given by
X
Fcal ðH Þ ¼ V qðr Þ exp ð2piH:r ÞdV ð2:43Þ
r
ln x ¼ x 1 we get,
( " # )
kF000 X 1
qðri Þ ¼ sðri Þ exp jFobs ðH Þ Fcal ðH Þj exp ð2pjH:r Þ
N rðH Þ2
ð2:46Þ
where F000 ¼ Z; the total number of electrons in a unit cell. Eq. 2.46 cannot be
solved as it is, since Fobs ðH Þ is defined on qðr Þ: In order to solve Eq. 2.46 in a
simple manner, we introduce the following approximation which replaces Fcal ðH Þ
as
X
Fcal ðkÞ ¼ V sðr Þ exp ð2piH:r ÞdV ð2:47Þ
where kf and ki are the momenta of the scattered and incident photons or neutrons,
respectively, and h is the diffraction angle. For X-ray scattering the quasi-elastic
scattering intensity usually includes the inelastic scattering intensity due to phonons,
since unless a special high resolution setup is utilized it cannot be separated from the
elastic scattering, while the single-electron inelastic (Compton) scattering intensity
is excluded. For neutron scattering the phonon is either excluded or included only
approximately. The average, or total, structure factor can be defined as
h
i
I ðQ Þ hf ðQ Þi 2 f ð Q 2 Þ
Sð Q Þ ¼ þ ð2:52Þ
hf ðQÞi2 hf ðQÞi2
where G(r) gives the number of atoms in a spherical shell of unit thickness at a
distance r from a reference atom. It peaks at characteristic distances separating
pairs of atoms and thus reflects the atomic structure. G(r) is the Fourier transform
of the experimentally observable total structure function, S(Q) that is
Z
Q¼max
2
Gðr Þ ¼ Q½SðQÞ 1 sinðQr ÞdQ ð2:56Þ
ð
Q¼0
The structure function is related to the coherent part of the total diffracted
intensity of the
P 2
~ ci fi Q
I coh Q ~
S Q~ ¼1þ P 2 ð2:57Þ
~
ci f i Q
where I coh ðQÞ is the coherent scattering intensity per atom in electron units and ci
and fi are the atomic concentration and X-ray scattering factor, respectively, for the
atomic species of type i. G(r) is simply another representation of the diffraction
data. However exploring the diffraction data in real space is advantageous, espe-
cially in the case of materials with significant structural disorder.
In order to refine an experimental PDF one needs to calculate a PDF from a
structural model. This can be done using the relation
" #
1 X X bi bj
Gcalc ðr Þ ¼ a r rij 4ðr~n0
€ ð2:58Þ
r i j hbi2
where the sum goes over all pairs of atoms i and j within the model separated by
rij : The scattering power of atom i is bi and hbi is the average scattering power of
the sample. In case of neutron scattering, bi is simply the scattering length, in case
of X-rays, it is the atomic form factor evaluated at a user defined value of Q.
-0
G(r)(Å) -2
-1
-2
-3
-4
Difference
-5
-6
4.0 10.5 17.0 23.5 30.0
r(Å)
such a way as to improve the agreement between the calculated and experimental
PDFs. This is done with, or without, observing predefined constraints imposed by
the symmetry of the space group of the crystal structure being tested (Proffen and
Billinge 1999). In this way, local distortions away from the average structure or
lower, unresolved, symmetries can be modeled. The PDF is very sensitive to the
coordination environment of atoms over short (\5 Å) and intermediate (5–20 Å)
ranges. The approach has proved to be quite successful in determining the structure
of various crystalline materials exhibiting different degrees of structural disorder
(Gutmann et al. 2000; Proffen et al. 1999; Petkov et al. 1999a, b).
In Fig. 2.4 the result of PDF refinement is plotted against G(r) (the number of
atoms in a spherical shell of unit thickness) and r (the distance from a reference
atom). The solid line is the calculated PDF. The fitted circles are the observed data.
The difference between the calculated and observed PDF is plotted below the data.
Since PDF contains Bragg and diffuse scattering, the information about local
arrangements is preserved. The PDF can be understood as a bond—length dis-
tribution between all pairs of atoms i and j within the crystal (up to a maximum
distance), however each contribution has a weight corresponding to the scattering
power of the two atoms involved.
The study of a measured PDF ranges from a simple peak width analysis revealing
information about correlated motion (Jeong et al. 1998) to the full profile refine-
ment of the PDF based on a structural model either using the reverse Monte Carlo
technique (Toby and Egami 1992) or least-squares regression (Billinge 1998) as
implemented in the program PDFFIT (Proffen and Billinge 1999).
The experimental PDF was obtained by initially extracting the coherently
scattered intensities from X-ray diffraction pattern by applying appropriated cor-
rection for flux, Compton scattering, and sample absorption. The intensity was
2.8 Pair Distribution Function 61
normalized in absolute electron unit, reduced to atomic PDF. All data procession
was done using the program PDFgetX (Jeong et al. 2001).
PDFgetX (Jeong et al. 2001) is a program used to obtain the observed atomic
pair distribution function (PDF) from measured X-ray powder diffraction data. The
observed and calculated PDF has been obtained from the graphical software
PDFgui (Farrow et al. 2007), which is a graphical environment for PDF fitting.
This allows for powerful usability features such as real-time plotting and remote
execution of the fitting program while visualizing the results locally. Modeling of
PDF was done using the software PDFFIT (Proffen and Billinge 1999) to yield
structural parameters.
The program PDFFIT is designed for the full profile structural refinement of the
atomic pair distribution function (PDF). In contrast to conventional structure
refinement based on Bragg intensities, the PDF probes the local structure of the
studied material. The program allows the refinement of atomic positions, aniso-
tropic thermal parameters, and site occupancies as well as lattice parameters and
experimental factors. By selecting individual atom types one can calculate partial
and differential PDFs in addition to the total PDF. Furthermore, one can refine
multiple data sets and/or multiple structural phases. The program is controlled by a
command language, which includes a FORTRAN style interpreter supporting
loops and conditional statements. This command language is also used to define
the relation between refinement parameters and structural or experimental infor-
mation, allowing virtually any constraint to be implemented in the model. PDFFIT
is written in Fortran-77. The basic concept, command language, and some file
formats of PDFFIT are taken from the diffuse scattering and defect structure
simulation program DISCUS (Proffen and Neder 1997).
Using the PDFFIT (Proffen and Billinge 1999) software the PDF refinement
was executed with the refinement of structural parameters like lattice parameters,
phase scale factor, linear atomic correlation factor, quadratic atomic correlation
factor, low r peak sharpening, peak sharpening cut-off and cut-off, for profile setup
functioning to get the absolute phase. The data configuration parameters are PDFfit
range with step size, data scale factor, upper limit for Fourier transform to obtain
data PDF, resolution peak broadening factor, data collection temperature etc.,
which can be refined to get accurate PDF fitting. In the end the observed and
calculated PDFs are visualized and compared.
• We would like to be able to sit on an atom and look at our neighborhood. The
PDF method allows us to do that.
• The PDF gives different information about different length-scales. We can see
the structure within a domain at low-r and between domains at high-r.
• PDF gives both local and average structure information.
62 2 Charge Density Analysis from X-Ray Diffraction
References
Azaroff LV (1968) Elements of x-ray crystallography. Mc Graw hill book company, New York,
p 79
Bader RFW (1991) Atoms in molecules, a quantum theory. Oxford University Press, Oxford
Berg JE, Werner PE (1977) Z. Kristallogr 145:310–320
Billinge SJL (1998) Local structure from diffraction. In: Billinge SJL, Thorpe MF. Plenum, New
York, p 137
Billinge SJL, DiFrancesco RG, Kwei GH, Neumeier JJ, Thompson JD (1996) Phys Rev Lett
77:715–718
Billinge SJL, Proffen TH, Petkov V (2000) Phys rev B 1 62(Number 2):1203–1211
Bragg WL (1912) Nature 90:410–410
Bragg WL (1913) The diffraction of short electromagnetic waves by a crystal. Proc Cambridge
Philoso Soc 17:43–57
Breck DW, Everslole WG, Milton RM, Reed TB, Thomas TL (1956) J Am Chem Soc 78:5963–
71
Bricogne G (1988) Acta Cryst A44:517–545
Caglioti G, Paeletti A, Ricci FP (1958) Nucl Instrum 3:223–228
Chandler GS, Figgis BN, Reynolds PA, Wolff SK (1994) Chem Phys Lett 225:421–426
Cheetham AK, Taylor JC (1977) J Solid State Chem 21:253–257
Chen R, Trucano P, Stewart RF (1977) Acta Cryst A33:823–828
Collins DM (1982) Nature 298:49–51
Coppens P, Yang YW, Blessing RH, Cooper WF, Larsen FK (1977) J Am Chem Soc 99:760–766
Coppens P, GuruRow TN, Leung P, Stevens, ED, Becker PJ, Yang YW (1979) Acta Cryst A
35:63–72
Cotton FA, Stanley TT (1977) Inorg Chem 16:2668–2671
Cox DE, Hastings JB, Thomlinson W, Prewill CT (1983) Nucl Instrum Methods 208:573–578
Cromer DT, Larsen AC, Stewart RF (1976) J Chem Phys 65:336–349
Cullity BD, Stock SR (2001) Elements of X-ray diffraction, Pearson education. 3rd edn. Prentice
Hall, Upper Saddle River, p 558
Dawson B (1967) Proc Royal Soc London Ser A 298:255–263
De Marco JJ, Weiss RJ (1965) Phys Rev 137 A1:869–871
Dinnibier RE (2008) Powder diffraction: theory and practice. In: Dinnibier RE, Billinge SJL
(eds.) Royal society of chemistry, p 5
Dollase WAJ (1986) Appl Crystallogr 19:267–272
Egami T (1990) Mater Trans JIM 31:163–176
Egami T, Billinge SJL, Thorpe MF (eds) (1998) Local structure from diffraction. Plenum, New
York
Farrow CL, Juhas P, Liu JW, Bryndin D, Bozin ES, Bloch J, Proffen T, Billinge SJL (2007) J
Phys Condens Matter 19:335219
Figgis BN, Iverson BB, Larsen FK, Reynolds PA (1993) Acta Cryst B49:794–806
Fox AG, Tabernor MA (1991) Acta Metall Mater 39:669–678
George W (2006) Powder diffraction—The Rietveld method and two stage method chapter 1.
Springer, Berlin
Goeta AE, Howard JAK (2004) Chem Soc Rev 33:490–500
Gull SF, Daniell GJ (1978) Nature 272:686–690
Guru Row TN (1983) J Chem Sci 92:4–5
Gutmann M, Billinge SJL, Brosha EL, Kwei GH (2000) Phys Rev B 61:11762–11768
Hansen NK, Coppens P (1978) Acta Cryst A34:909–921
Hansen NK, Schneider JR, Yellon WB (1987) Pearson WH Acta Cryst A43:763–769
Harel M, Hirshfeld FL (1975) Acta Cryst B31:162–172
Hewat AW (1986) Chem Scripta 26A:119–130
Hirshfeld FL (1971) Acta Cryst B 27:769–781
References 63
Shmueli U (2008) IUCr series, International tables of crystallography. Volume 2, chapter1, 3rd
edn. Springer, Berlin
Stevens ED, Coppens P (1976) Acta Cryst A32:915–917
Stewart RF (1969) J Chem Phys 51:4569–4577
Stewart RF (1972) J Chem Phys 57:1664–1668
Stewart RF (1976) Acta Cryst A32:545–574
Stout GH, Jensen LH (1968) X-ray structure determination—a practical guide. The Macmillan
Company Collier-Macmillan, London, p 217
Takata M, Umeda B, Nishibori E, Sakata M, Saito Y, Ohno M, Shinohara H (1995) Nature
377:46–49
Takata M, Nishibori E, Shinmura M, Tanaka H, Tanigaki K, Kosaka M, Sakata M (2001) Mater
Sci Eng A312:66–71
Thompson P, Cox DE, Hastings JB (1987) J Appl Crystallogr 20:79–83
Toby BH, Egami T (1992) Acta Cryst A48:336–346
Wagner CNJ (1978) J Non Cryst Solids 31:1–40
Warren BE (1990) X-ray diffraction, chapter 3. Dover publications, New York
Weiss H (1925) Proceedings of the royal society of london series A 108(748):643–654
Wertheim GK, Butler MA, West KW, Buchanan DNE (1974) Rev Sci Instrum 45:1369–1371
Young RA (1995) The rietveld method, Volume 5 of International union of crystallography
momographs on crystallography, Chap 1. Oxford Science publications, Reprint Oxford
University press
Young RA, Mackie PE, Von Dreele RBJ (1977) Appl Crystallogr 10:262–269
Zachiarsen WH (1948) Acta Crstallogr 1:265–8, 277–87
Zachiarsen WH, Ellinger FH (1963) Acta Crystallogr 16:369–75
Chapter 3
Results and Discussion on Metals
and Alloys
3.1.1 Introduction
R. Saravanan and M. Prema Rani, Metal and Alloy Bonding: An Experimental Analysis, 65
DOI: 10.1007/978-1-4471-2204-3_3, Ó Springer-Verlag London Limited 2012
66 3 Results and Discussion on Metals and Alloys
The precise study of bonding in materials is always useful and interesting yet no
study can reveal the real picture, because no two experimental data are identical.
This problem is enhanced when the model used for the evaluation of electron
densities is not entirely suitable. Fourier synthesis of electron densities can be of
use in picturing the bonding between two atoms, but it suffers from the major
disadvantage of series termination error and negative electron densities which
prevent the clear understanding of the bonding, a factor intended for analysis. The
advent of Maximum Entropy Method (MEM) has solved many of these problems.
MEM electron densities are always positive and even with limited data, one can
determine reliable electron densities resembling the true densities. Currently,
multipole analysis of charge densities has been widely used to study crystalline
materials (Coppens and Volkov 2004; Poulsen et al. 2004; Frils et al. 2004; Pillet
et al. 2004; Marabello et al. 2004). The multipole technique of synthesising the
electron density of an atom into core and valence parts yields an accurate picture
of the bonding in a crystalline system. In this work, such an analysis has been
carried out using the formalism of Hansen and Coppens (1978).
Bonding in materials has been analysed by Israel et al. (2002), Saravanan et al.
(2002a, b), Kajitani et al. (2001) using MEM technique and reported X-ray data
and they obtained precise information about bonding in materials. The nature of
bonding in these materials as analysed are found to be ionic, mixed covalent and
ionic, and ‘oxide’ bonding.
The nature of bonding and the charge distribution in sodium and vanadium metals
have been analysed using the reported X-ray data of these metals. MEM and
multipole analysis have been used for elucidation and analysis of bonding in these
metals. The mid-bond densities in sodium and vanadium are found to be 0.014 and
0.723 e/Å3 respectively, giving an indication of the strength of the bonds in these
materials. From the multipole analysis, the sodium atom is found to contract more
than the vanadium atom.
In this study on sodium and vanadium, we have attempted to study the metal
bonding in these metals using MEM and multipole techniques. The X-ray structure
factor data of sodium has been taken from Field and Medlin (1974) and that of
vanadium from Linkoaho (1972). Other relevant details can be found from these
papers. Sodium has been chosen in this present study because of its very large
thermal vibration (Field and Medlin 1974) and because it would be interesting to
analyse the highly reactive nature of the metal in terms of strength of bonding and
3.1 Sodium and Vanadium Metals 67
where tk ð~
uÞ is the Gaussian term and indicates a convolution. The atomic charge
density is then defined as
X4 X l !
! ! 03 0 003 00 r
q r ¼ Pc qcore r þ Pv k qvalence ðk r Þ þ k Rl ðk r Þ Plm Ylm
l¼0 m¼l
r
ð3:2Þ
where Pc, Pv and Plm are population coefficients. Canonical Hartree–Fock atomic
orbitals of the free atoms normalised to one electron can be used for the con-
struction qcore and qvalence but the valence function is allowed to expand and
contract by the adjustment of the variable parameters k0 and k00 : The effect of the
temperature can be distinguished from the convoluted and the deconvoluted form
of thermal contribution to the charge density as dynamic and static multipole
deformation maps. The deformation density in these maps are characterised by
! 1 X ! ! h ! i
Dqmultipole deformation r ¼ F h F h exp 2 h ; !
r
v multipole spherical atom
ð3:3Þ
where, Fmultipole is the Fourier transform of the multipole charge density with or
without the convolution of thermal contribution, where the Fourier components are
terminated at the experiment resolution.
Both the reported data sets have been analysed using the GSAS (Larson and
Von Dreele 2004) programme of the Los Alamos National Laboratory, USA.
The reliability index wR for sodium data is found to be 2.2% after successful
68 3 Results and Discussion on Metals and Alloys
refinement. The refined Debye–Waller factor is found to be BNa = 7.899 (96) Å2.
This value is very large and comparable with that reported [BNa = 7.86 Å2, Field
and Medlin (1974)], the difference being only 0.039 Å2. The reliability index for
the vanadium refinement is wR = 1.1% and the Debye–Waller factor is found to
be 0.404(24) Å2. This value is also comparable to that reported by Linkoaho
(1972), i.e., BV = 0.37 Å2, the difference being 0.034 Å2.
The reliability indices and the matching of thermal parameters with the
corresponding reported values for both the atoms indicate that the quality of the
data for the electron density analysis is suitable enough. The same number of data
for both the crystal systems has been used for the present analysis in order to make
a comparison also of the bonding (only insignificant changes in the density maps
and numerical results has been observed on including all the reflections for
sodium).
The MEM analysis has been carried out by dividing the cell into 643 pixels. In
the case of sodium data, convergence condition C = 1.0 was obtained after 1464
cycles of refinement; the reliability indices are found to be RMEM = 2.157% and
wRMEM = 2.356%. The Lagrange parameter k used was 0.024 for sodium. For
vanadium, convergence was obtained after 236 cycles with R values as
RMEM = 0.392% and wRMEM = 0.498%. A value of 0.024 was used for k. The
prior densities assumed for sodium and vanadium were 0.2785 and 1.6634 e/Å3,
respectively. The resolution of the map was 0.067 and 0.047 Å/pixel for sodium
and vanadium, respectively. The analysis on the MEM results has been given in
the following section. Some important results from the MEM refinements are
given in Table 3.1.
The MEM electron density maps in the density rich regions of sodium and
vanadium on the (100) plane have been given in Fig. 3.1a, b, respectively.
The density distribution of sodium and vanadium on the (110) plane has been
shown in Fig. 3.2a, b, respectively.
3.1 Sodium and Vanadium Metals 69
3.0
2.0
[010]
[010]
2.0
1.0
1.0
Na V
Fig. 3.1 a MEM electron density of sodium on (100) plane. Contour lines are between 0.0 and
1.0 e/Å3. Contour interval is 0.05 e/Å3. b MEM electron density of vanadium on (100) plane.
Contour lines are between 0.2 and 4.5 e/Å3. Contour interval is 0.224 e/Å3
3.0 2.0
[001]
[001]
2.0
1.0
1.0
Na V
0.0
1.0 2.0 3.0 4.0 5.0 6.0 1.0 2.0 3.0 4.0
[110] [110]
Fig. 3.2 a MEM electron density of sodium on (110) plane. Contour lines are between 0.0 and
1.0 e/Å3. Contour interval is 0.05 e/Å3. The plane has been shifted half a unit cell along the z
direction. b MEM electron density of vanadium on (110) plane. Contour lines are between 0.04
and 6.0 e/Å3. Contour interval is 0.298 e/Å3. The plane has been shifted half a unit cell along the
z direction
1.275 1.850 2.425 3.000 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6
Distance (Å) Distance (Å)
Fig. 3.3 a One-dimensional electron density variation of sodium along [100], [110] and [111]
directions. b One-dimensional electron density variation of vanadium along [100], [110] and
[111] directions
Number of reflections
5
of number of reflections
plotted in the error space of 2 4
vanadium
3
1 2
0 0
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
[Fobs-FMEM ]/σ (K) [Fobs- FMEM ]/σ (K)
of the body-centred atom is seen at the centre of the (100) plane of the cell as seen
from Fig. 3.1a. It also represents the very large thermal vibration of the sodium
atom.
The corners show that the core regions are perfectly symmetric and spherical.
The four voids seen on the (100) plane conform to the crystal system. The conical-
shaped contour lines seem to be attracted towards the centre from the corners. This
may be due to the fact that the core of the centre atom attracts the valence electrons
of the eight corner atoms. Since the atomic number of sodium is only 11, this core–
valence attraction is very flexible and leads to large smearing of the electrons.
Moreover, it should be remembered that the thermal vibration of the sodium atom
being very large (BNa = 7.899(96) Å2), makes this electrostatic interaction pro-
cess highly feasible. The large thermal parameter of sodium and also the smaller
atomic number of sodium makes the body-centred sodium atom possible to reveal
itself at the face centre.
3.1 Sodium and Vanadium Metals 71
Figure 3.1b shows the density distribution of vanadium on the (100) plane. The
centre of the outer contour lines are bent towards the core of each atom, showing
strong attraction of valence electrons by the core, which is due to the large atomic
number of vanadium and its smaller thermal vibration parameter. Two humps are
seen at the outer contour lines of each of the corner atoms. This may be due to the
overlapping of charges between the cores in the bonding direction.
Figure 3.2a, b show the densities of sodium and vanadium respectively, on the
(110) plane (shifted by half the cell distance along z direction). Large thermal
smearing of sodium atom is seen (Fig. 3.2a) as seen from the body-centred atom.
The core of one atom attracting the valence of other atom is also seen. In vanadium
also, this process is taking place as seen from Fig. 3.2b.
From Fig. 3.3a, which shows the one-dimensional variation of the electron
densities of sodium along the three directions, [100], [110] and [111], the mini-
mum density is found at a distance of 1.732 Å at which the electron density is
found to be 0.014 e/Å3. The electron density along the bonding region of the [111]
direction is flat with low level of density similar to the density along the [110]
direction. But along [100] direction, the density is slightly higher than that along
the other two directions. The large thermal vibration of sodium atom and the
smallest inter atomic distance along [100] direction leads to higher density value
along [100]. The inter atomic distances are larger along [110] and [111] directions
and the lone valence electron of sodium atom leads to lower densities along these
directions. As far as the [100] direction is concerned, the amount of electrons
contributing to the bond electron density comes from the corner atom along [100]
direction and the body-centred atom too (due to the large thermal smearing), to
decide the electron density along the shortest [100] direction.
Figure 3.3b shows the one-dimensional variation of electron density of vana-
dium along the three crystallographic directions. There is a slight increasing
undulation of electron density along [111] direction. This starts beyond a distance
of about 0.7 Å from the centre. The electron density at 1.309 Å is found to be
0.723 e/Å3. This is an indication of the bond strength in vanadium compared to
that of sodium, which has a density of 0.014 e/Å3. The small ellipses of charge
density in the [110] plane as seen from Fig. 3.8b and the void in Fig. 3.8a confirm
these magnitudes. The reported (Wilkinson) radii for sodium and vanadium are
1.57 and 1.22 Å. In this work, the radii are estimated from the electron densities
along the [110] direction, because, along [100] direction, the atoms are too close to
decide the inter atomic interactions and along [111] direction, the bonding effects
are predominant. Thus, the estimated radii are 1.517 and 1.136 Å respectively, for
sodium and vanadium, very close to those reported (Wilkinson).
Figure 3.4a, b representing the error analysis of the data used for the present
analysis conform to the quality of good data sets and justify their usage for the
analysis.
The results of the multipole analysis have been given in Table 3.2.
The results of multipole refinement as given in Table 3.2 show the thermal
vibration parameters, the multipole parameters and the reliability indices. The
large j0 parameter of the valence part of the sodium atom (one electron) denotes
72 3 Results and Discussion on Metals and Alloys
the contraction needed to restrict the lone valence electron to its nearly ideal value
of one. In this work, the number of core and valence electrons has also been
refined in the multipole formalism, for both sodium and vanadium, which leads to
very close and reasonable values. For sodium, the core and valence electrons are
refined to be 9.94 and 1.13 e- the total being 11.07 e- well within the experi-
mental error limits. Only a large j0 parameter will restrict the large thermal
vibration of sodium atom to regain this amount of electrons in the valence region.
In vanadium, the value of j0 is 1.11, again showing contraction of the atom, but not
as much as in sodium. The number of core and valence electrons in vanadium is
found to be 17.876 and 4.573 e-. There is a deficit of around 0.551 e- in vana-
dium, which actually shows up as a share in the bonding electron and shoots up the
mid-bond density. In sodium there is no deficit in the total number of electrons
leading to meagre share to the bonding and a decrease in the mid-bond density
value. Figure 3.2a, b clearly show these features; clear and density-free regions in
sodium and charge accumulated mid-bond regions in vanadium. Thus MEM and
multipole results are complementary to each other.
The static multipole deformation maps are shown in Fig. 3.5a, b for sodium and
vanadium respectively on the (100) plane. The static multipole deformation maps
are shown in Fig. 3.6a, b for sodium and vanadium respectively on the (00)
plane.
Figure 3.5a shows positive electron contours at the centre of the (100) plane,
and at the edge centres and corners, although the values of these densities are very
small. Yet, right from the core, the positive densities are seen showing enhanced
thermal vibration of sodium atom and the contracted spherical core in line with the
j0 parameter of sodium.
The (00) section (Fig. 3.6a) also shows the sodium atom with large thermal
vibration of the valence region. Figure 3.5b shows the nullified core and residual
valence region (in terms of thermal vibration) of vanadium atom on (100) plane.
The core is not seen to be as much contracted as in the case of sodium as seen from
Fig. 3.6b.
3.1 Sodium and Vanadium Metals 73
Fig. 3.5 a Static multipole deformation map of sodium on the (100) plane. Dotted lines indicate
negative electron densities. b Static multipole deformation map of vanadium on the (100) plane.
Dotted lines indicate negative electron densities
Fig. 3.6 a Static multipole deformation map of sodium on the (00) plane. Dotted lines indicate
negative electron densities. b Static multipole deformation map of vanadium on the (00) plane.
Dotted lines indicate negative electron densities
3.1.6 Conclusion
The nature of bonding and the charge distribution in sodium and vanadium metals
has been analysed using the reported X-ray data of these metals. MEM and the
multipole analysis have been used for this purpose. The bonding in these metals
has been elucidated and analysed. The nature of bonding are clearly revealed by
74 3 Results and Discussion on Metals and Alloys
the two-dimensional MEM maps plotted on (100) and (110) planes and the mid-
bond densities are clearly revealed by the one-dimensional electron density along
the [100], [110] and [111] directions. The mid-bond densities in sodium and
vanadium are found to be 0.014 and 0.723 e/Å3 respectively, giving an indication
of the strength of the bonds in these materials. The sodium atom is found to
contract more than the vanadium atom from the multipole analysis.
3.2.1 Introduction
The average and local structures of simple metals Al, Ni and Cu have been
elucidated for the first time using MEM, multipole and PDF. The bonding between
constituent atoms in all the above systems is found to be well pronounced and is
clearly seen from the electron density maps. The MEM maps of all the three
systems show the spherical core nature of atoms. The mid-bond electron density
profiles of Al, Ni and Cu reveal the metallic bonding nature. The local structure
using PDF profile of Ni has been compared with that of the reported results. The R
value in this work using low Q XRD data for the PDF analysis of Ni is close to the
value reported using high Q synchrotron data. The cell parameters and displace-
ment parameters were also studied and compared with the reported values.
The powder X-ray intensity data was collected at the Regional Research Laboratory
(RRL), Council of Scientific and Industrial Research (CSIR), Thiruvananthapuram,
India, using X-PERT PRO (Philips, Netherlands) X-ray diffractometer with a
monochromatic incident beam, which offers pure Cu Ka1 radiation. The wavelength
used for X-ray intensity data collection is 0.154056 nm with 2h range of data col-
lection from 10 to 1208. The raw intensities were refined using the software pro-
gramme JANA 2000 (Petřiček et al. 2000), considering the fcc unit cell of Al, Ni and
Cu with 4 atoms/unit cell in the space group Fm3m. The fitted profiles of observed
and calculated relative intensities with the Bragg peaks along with the difference are
given in Fig. 3.7a, b, c for Al, Ni and Cu respectively.
The results of the refinements are given in Table 3.3.
The refined structure factors are given in Table 3.4a, b, c respectively, for Al,
Ni and Cu.
The refined structure factors were used for MEM analysis using the methodologies
discussed in earlier works, e.g., Saravanan et al. (2003), Israel et al. (2003, 2004).
In the present calculation, the unit cell was divided into 643 pixels and the initial
electron density at each pixel was fixed uniformly as Z/a30, where Z is the number
of electrons in the unit cell. The electron density is evaluated by carefully selecting
the Lagrange multiplier in each case such that the convergence criterion C
becomes unity after performing the minimum number of iterations. The resulting
parameters of MEM computations have been given in Table 3.5.
The electron density distribution of Al, Ni and Cu has been mapped using the
MEM electron density values obtained through refinements. Figure 3.8a, b, c show
76 3 Results and Discussion on Metals and Alloys
Fig. 3.7 a The Rietveld refined powder profile for Al along with difference between observed
and calculated relative intensities. b The Rietveld refined powder profile for Ni along with
difference between observed and calculated relative intensities. c The Rietveld refined powder
profile for Cu along with difference between observed and calculated relative intensities
the MEM electron density distribution of Al, Ni and Cu on the (100) plane of the
unit cell.
Figure 3.8d, e, f show the MEM electron density distribution of Al, Ni and Cu
on the (110) plane.
Figure 3.9a, b, c show the one-dimensional variation of the electron densities of
Al, Ni and Cu along the three directions [100], [110] and [111] of the unit cell
respectively.
The electron densities at the saddle points along these three directions have
been given in Table 3.6.
3.2 Aluminium, Nickel and Copper 77
Table 3.4 Observed and calculated structure factors for a Aluminium b Nickel c Copper
hkl Fo Fc Fo - Fc r(Fo)
a Aluminium
111 35.07 35.05 0.03 0.22
0 02 32.87 32.84 0.03 0.16
202 26.63 26.82 -0.19 0.36
113 22.87 23.18 -0.32 0.11
222 22.34 22.17 0.17 0.58
004 18.83 18.64 0.19 0.83
313 15.96 16.47 -0.51 0.22
204 16.23 15.82 0.41 0.37
b Nickel
111 68.16 67.93 0.23 0.42
002 61.40 61.96 -0.57 0.59
202 44.38 45.76 -1.39 0.68
113 38.60 37.59 1.00 0.56
222 34.01 35.38 -1.36 1.31
c Copper
111 78.59 78.49 0.10 0.53
002 72.71 72.65 0.06 0.76
202 55.34 55.46 -0.12 0.87
113 47.00 46.73 0.27 0.75
222 44.06 44.33 -0.27 1.67
004 35.95 36.58 -0.64 2.99
Table 3.6 Electron density profiles along the three directions of the unit cells
Direction metals [100] [110] [111]
Position (Å) ED (e/Å3) Position (Å) ED (e/Å3) Position (Å) ED (e/Å3)
Al 2.026 0.1467 1.433 0.2636 3.509 0.1467
Ni 1.763 0.2296 1.247 1.1010 3.054 0.2296
Cu 1.813 0.3318 1.282 0.7649 3.139 0.3317
ED Electron density
78 3 Results and Discussion on Metals and Alloys
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
[001]
[001]
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[010] [010]
(c) (d)
1.0
1.0
0.9
0.9
0.8 0.8
0.7 0.7
[001]
[001]
0.6 0.6
0.5 0.5
0.4 0.4
0.3 0.3
0.2
0.2
0.1
0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[010] [110]
(e) (f)
1.0 1.0
0.9 0.9
0.8 0.8
0.7 0.7
[001]
[001]
0.6 0.6
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[110] [110]
Fig. 3.8 a MEM electron density distribution of Al on the (100) plane. Contour range is
0.0–2.0 e/Å3 and contour interval is 0.1 e/Å3. b MEM electron density distribution of Ni on the
(100) plane. Contour range is 0.0–2.0 e/Å3 and contour interval is 0.1 e/Å3. c MEM electron
density distribution of Cu on the (100) plane. Contour range is 0.0–2.0 e/Å3 and contour interval
is 0.1 e/Å3. d MEM electron density distribution of Al on the (110) plane. Contour range is 0.0 to
3.0 e/Å3 and contour interval is 0.15 e/Å3. e MEM electron density distribution of Ni on the (110)
plane. Contour range is 0.0–3.0 e/Å3 and contour interval is 0.15 e/Å3. f MEM electron density
distribution of Cu on the (110) plane. Contour range is 0.0–3.0 e/Å3 and contour interval is
0.15 e/Å3
3.2 Aluminium, Nickel and Copper 79
Electron density(e/Å )
0.7
3
3
1.4
0.6
1.2
0.5
1.0
0.4
0.8
0.3
0.6
0.2
0.4
0.1
0.0 0.2
1.0 1.5 2.0 2.5 3.0 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 2.3
Distance along [100] direction (Å) Distance along [110] direction (Å)
(c) 1.4
Cu
1.2 Ni
Electron density(e/Å )
3
Al
1.0
0.8
0.6
0.4
0.2
Fig. 3.9 a One-dimensional electron density profiles of Al, Ni and Cu along the [100] direction
of the unit cell. b One-dimensional electron density profiles of Al, Ni and Cu along the [110]
direction of the unit cell. c One-dimensional electron density profiles of Al, Ni and Cu along the
[111] direction of the unit cell
For materials whose structure is not reflected in the long-range order of the crystal,
an alternate structural analysis called the Pair Distribution Function (PDF)
approach is used. This method is sometimes called the real-space structure
determination method, because the PDF is modelled in real-space, rather than in
the reciprocal space. The PDF reflects the short-range ordering in a material or the
probability to find an atom at a distance r from another atom. In other words PDF
gives the bond length distribution of the material under consideration (Proffen and
Neder 1997). It well with the inter atomic distances computed from a crystallo-
graphic model, when there are no short-range deviations from the average struc-
ture (Tobi and Egami 1992; Peterson et al. 2003).
This approach has been widely used for studying the structure of glasses and
liquids (Debay and Menki 1930; Bowran and Finney 2003). More recently, it has
been applied to disordered crystalline and partially crystallised materials.
80 3 Results and Discussion on Metals and Alloys
G(r) (Å -2)
10
0.00
0
-10
-20
-3.25 Difference
Difference
-30
-40
4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20
r(Å) r(Å)
0.3875
G(r) (Å -2)
-0.0750
-0.5375
Difference
4 6 8 10 12 14 16 18 20
r(Å)
Fig. 3.10 a The observed and calculated PDFs of Al with difference profile. b The observed and
calculated PDFs of Ni with difference profile. c The observed and calculated PDFs of Cu with
difference profile
Table 3.7 Distances from System First nearest neighbour distance (Å)
the PDF profiles
Observed Calculateda
Al 2.86 2.863
Ni 2.50 2.492
Cu 2.38 2.564
a
Wyckoff 1963
fitted observed and calculated PDF profiles along with the difference between them
for Al, Ni and Cu respectively, and Table 3.7 gives the bond length distribution of
atoms.
The fitted profiles of Rietveld refinements (Fig. 3.7a, b, c for Al, Ni and Cu
respectively) give a clear picture of the quality of the sample as well as the data.
The accuracy of refinements can be judged from the difference curve between the
observed and calculated intensity, which is almost a straight line. The results of the
82 3 Results and Discussion on Metals and Alloys
Fig. 3.11 a Static Multipole Deformation map of Al on (100) plane. Dotted lines indicate
negative electron densities. b Static Multipole Deformation map of Ni on (100) plane. Dotted
lines indicate negative electron densities. c Static Multipole Deformation map of Cu on (100)
plane. Dotted lines indicate negative electron densities
3.2 Aluminium, Nickel and Copper 83
structural refinement (Petřiček et al. 2000) given in Table 3.3 show reasonable
values of the Debye–Waller factor for Al, Ni and Cu obtainable from a powder
data set, which is comparable to those reported in Peng et al. (1996). The
experimental cell constant of Al, Ni and Cu are highly comparable to the reported
values (Wyckoff 1963) and the differences are 0.0026, 0.0028 and 0.0096 Å for
Al, Ni and Cu respectively, which are very low indicating the accuracy of
refinements and also the precision of the observed data sets. The reliability indices
for all the powder samples are very low indicating the correctness of the refine-
ments (RP = 9.84, 6.21 and 6.77 for Al, Ni and Cu respectively).
The MEM electron density distribution map of Al as given in Fig. 3.8a for
(100) plane reveals the core of the Al atom being spherical and this sphericity
persists even at slightly larger distances away from the centre. The distribution of
charges all over the plane, indicates the distribution of electrons and the charges at
the edge centres on (100) plane show the distribution from the perpendicular face-
centred atoms, indicating the extension of the spread of the charges. The electron
density map of Al on (110) plane shows (Fig. 3.8b) highly concentrated charges at
the core of the atoms, and the distribution of charges on places other than the
atomic positions, due to the valence, free electrons.
The core of Ni and Cu atom as seen from the electron density map on (100)
plane (Fig. 3.8b, c), show the perfect spherical nature of the electronic charge
clouds. The edge centres show the distribution of electron density of the atoms
located at the face centres of the planes perpendicular to the paper. Since the
atomic numbers of nickel and copper are higher than that of aluminium, more
diffuse distribution of charges are seen in the electron density maps of Ni and Cu
on the (100) plane. The electron density map of Ni and Cu on the (110) plane
(Fig. 3.8e, f), show similar trends as observed in (100) plane.
The one-dimensional profiles of electron density constructed along [100], [110]
and [111] directions for Al, Ni and Cu are shown in Fig. 3.9a, b and c respectively.
The positions of minimum electron densities and density values have been given in
Table 3.6. The mid-bond density for Al is found to be 0.263 e/Å3 at a distance of
1.433 Å along the bonding direction [110]. The mid-bond electron density
between the Al atoms along [100] direction is 0.147 e/Å3 at a distance of 2.026 Å.
The first minimum and the mid-bond electron density between the atoms along
[111] direction are 0.119 e/Å3, 0.147 e/Å3 at a distance of 1.755 and 3.509 Å
respectively.
The mid-bond density for Ni is found to be 1.101 e/Å3 at a distance of 1.247 Å
along the bonding direction [110]. Table 3.6 shows that the mid-bond density
along [100] direction of Ni occurs at a distance 1.763 Å with a value 0.229 e/Å3.
The first and second minimum electron densities between the atoms along [111]
directions are 0.841 and 0.229 e/Å3 at a distance of 1.527 and 3.054 Å.
The mid-bond density for Cu is found to be 0.765 e/Å3 at a distance of 1.282 Å
along the bonding direction [110]. The mid-bond density along [100] direction
occurs at a distance 1.813 Å with a value 0.332 e/Å3. The first and second min-
imum electron densities between the atoms along [111] direction are 0.664 and
0.332 e/Å3 at a distance of 1.374 and 3.139 Å respectively. The mid-bond density
84 3 Results and Discussion on Metals and Alloys
in Ni is the largest value (Table 3.6) among the three metal systems we have
studied. But along the [100] and [111] directions (non-bonding), the element with
the higher atomic number has the higher electron density. Hence obviously the
bonding interaction is stronger in Ni as evidenced by the higher electron density
along [110] direction.
Hence, from the qualitative and quantitative analysis of the MEM electron
densities, the bonding in Al, Ni and Cu is predicted to be predominantly metallic
and the interaction of charges along other nonb-onding directions for Ni and Cu
seems to be heavier due to their metallic nature. (Because, if the systems are ionic
one cannot expect higher mid-bond values as in the present cases and if the
systems are covalent, one can expect values higher than the present values due to
overlapping of charges). Also the strength of the electron density is higher in the
case of Ni, which has lower number of electrons in the unit cell compared to Cu.
The electron density of Al in all the three directions [100], [110] and [111] are
relatively low compared to Ni and Cu, which conforms to the loosely packed
structure of Al. (The interaction of atomic charges will be less and hence electron
densities along bonding and directions other than bonding are expected to be
minimum). !
The refined pair distribution function (the Fourier transform of S Q ; the
reduced structure factor) given in Fig. 3.10a, b, c for Al, Ni and Cu show nice
matching of observed and calculated PDFs. The PDF refinement can be considered
as equivalent to matching the observed X-ray powder data with a model with too
little structure parameters. Moreover, high Q data are required for this type of
analysis on short-range order and local structure of materials (Egami 1990). In
spite of these factors, an attempt has been made to use the X-ray powder data for
this analysis, to test how well these powder data sets can be used for purposes like
PDF analysis.
The peaks in the PDF profile correspond to the nearest neighbour distances
given in Table 3.7. In Al, the difference between observed and calculated first
nearest neighbour distances turns out to be 0.003 Å. In Ni, the difference is
0.008 Å and for Cu it turns out to be 0.184 Å. These differences are not too high
considering the fact that we have used only (i) powder data with (ii) limited Q
values. There may also be local undulations in the structure of these three systems
which lead to these differences. The local disorder is expected to affect the repeat
distances, because of the diffuse scattering, which does not have sharp, single
pointed X-ray diffraction phenomenon. Hence, it is reflected in the local structural
analysis using suitable tools like PDF. From these differences (0.003 Å for Al,
0.008 Å for Ni and 0.184 Å for Cu), one can conclude that (i) the local disorder
increases as the atomic number increases due to more concentration of charges and
(ii) the disorder propagates at a longer distance from the origin leading to
enhanced differences in neighbour distances. In the case of Cu, an additional
disorder is possible due to the fact that CuKa radiation has been used for the
powder data collection which gives rise to resonance effects and possible
fluorescence.
3.2 Aluminium, Nickel and Copper 85
The PDF profile of Ni has been compared with that obtained with reported data
(Proffen and Billinge 1999). Very close resemblance is seen between these two
refinement results. The R values are also close to each other [18% in this work
compared to 14% in the work by Proffen and Billinge (1999)].
The multipole refinement shows no expansion/contraction for Al as seen from
the j0 value given in Table 3.8 [j0 = 1.007(0.052)]. In the case of Cu, there is a
slight contraction of the atom, but in Ni there is appreciable contraction compared
to the other two systems. There is a possibility that the contraction is effected to
obey Pauli’s principle and to make the system with the strongest mid-bond
interaction stable (as evidenced by the mid-bond electron density of Ni).
The results of SMD calculations given in Table 3.8 show the positive and
negative (maximum and minimum) electron densities for all the three systems.
Since these calculations are based on difference densities, the lowest electron
densities are seen for all the systems, i.e., the maximum positive difference is
0.05 e/Å3 for Al and the maximum negative density is -0.08 e/Å3 for Ni. These
values indicate the accuracy with which the multipole analysis has been carried
out. Moreover, the positive difference densities are the same for both Ni and Cu
(0.02 e/Å3), thereby justifying the fact that no biassed results are obtained as far as
the comparison of electron densities of Ni and Cu is concerned.
3.2.5 Conclusion
A clear understanding of the structure of the metals A1, Ni and Cu has been gained
from this work, in terms of cell parameters, thermal vibrations and electron density
distribution in the bonding region using MEM and local structural analysis using
PDF. Multipole model analysis reveals the clear qualitative and quantitative pic-
ture of the charge densities both in core and valance regions with the expansion
and contraction of the valence region.
3.3.1 Introduction
The applications of metallic materials are spread into various areas of science and
technology including electronics, aeronautics and space travel, to name a few.
Hence a study on some technologically important metals is essential in terms of
the local and average structures which are completely different. The usual methods
of analysis using structural refinement of X-ray or neutron data will give only the
average structure of the materials under investigation. There is only limited
information available about the investigations of materials in terms of the local
86 3 Results and Discussion on Metals and Alloys
structure. In the present investigation, powder data sets of six elemental metals
magnesium, titanium, iron, zinc, tin and tellurium were collected to study the
structure in terms of the electron density distribution between atoms using MEM
and the local and average structural properties using PDF. The PDF analysis
requires data sets of very high values of Q (= 4k sin h/k) which is achievable only
through synchrotron studies, which is not always accessible. But this work gives
reasonable results, which could be obtained through high Q data sets.
The average and local structures of the metals magnesium, titanium, iron, zinc, tin
and tellurium have been analysed using MEM, and PDF. The bonding between the
constituent atoms in all these systems is found to be well pronounced and is clearly
seen from the electron density maps. The mid-bond electron density profiles reveal
the metallic nature of the bonding. The local structure of Mg, Ti, Fe, Zn, Sn and Te
using the low Q XRD have given results comparable to high Q synchrotron data.
The cell and displacement parameters are also studied.
For the accurate charge density study, the MEM (Sakata and Sato 1990;
Collins 1982) has been used. It is a non-linear calculation that has its roots in
information and probability theory, and was originally designed to reconstruct the
most probable and least biassed probability distribution in an underdetermined
situation. The MEM can yield a high-resolution density distribution from a limited
number of diffraction data without using a structural model. The obtained density
distribution gives detailed structural information in the inter atomic region, e.g.
bonding electron densities. The ability of MEM in terms of a model-free recon-
struction of the charge densities from measured X-ray diffraction data can be
interpreted as ‘imaging of diffraction data’. Combining MEM with the Rietveld
method (Rietveld 1969) has been successful in developing a sophisticated method
of structure refinement in charge density level, that is, the MEM/Rietveld (Takata
et al. 1995). The MEM/Rietveld analysis is an iterative method combining the
MEM and Rietveld analyses. The MEM charge density at certain iteration step
normally succeeds in providing a better structural model for the Rietveld of the
next iteration. In this method, the final MEM electron density distribution derived
is compatible with the structural model used in Rietveld refinement. The method
has been successfully applied to the structural studies of metals and alloys (Robert
et al. 2009; Saravanan and Prema Rani 2007), semiconductors (Saravanan et al.
2008), magnetic materials (Syed Ali et al. 2009]) In this study this method has
been used to elucidate the electron density distribution of elemental metals.
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 87
Gðr Þ ¼ 4pr½qðr Þ q0
Z
2 !h ! i ! !
¼ Q S Q 1 sin Q r d Q ;
p
where G(r) is the atomic pair distribution function and qðr Þ corresponds to the
(atomic) number density at a distance r from the average atom. The atomic pair
distribution function, obtained from powder diffraction data is thus a valuable tool
for the study of the local atomic arrangements in a material, since both Bragg and
diffuse scattering information about local arrangements are preserved in PDF.
High quality powder samples of Mg, Ti, Fe, Zn, Sn and Te were purchased and
used for this work. The general properties of these metals are given in Table 3.9.
88
Fig. 3.12 a The Rietveld refined powder profile for Mg along with difference between observed
and calculated relative intensities. b The Rietveld refined powder profile for Ti along with
difference between observed and calculated relative intensities. c The Rietveld refined powder
profile for Fe along with difference between observed and calculated relative intensities. d The
Rietveld refined powder profile for Zn along with difference between observed and calculated
relative intensities. e The Rietveld refined powder profile for Sn along with difference between
observed and calculated relative intensities. f The Rietveld refined powder profile for Te along
with difference between observed and calculated relative intensities
90 3 Results and Discussion on Metals and Alloys
Table 3.10 Observed and calculated structure factors for a Magnesium b Titanium c Iron d Zinc
e Tin f Tellurium
hkl Fo Fc Fo-Fc r(Fo)
a Magnesium
220 130.19 130.30 -0.11 3.41
101 109.31 109.82 -0.51 0.95
200 158.93 158.44 0.49 0.73
220 130.19 130.30 -0.11 3.41
211 93.60 93.52 0.08 1.34
301 78.01 78.25 -0.24 2.57
112 103.02 103.72 -0.70 2.13
400 95.89 101.90 -6.01 3.17
321 71.97 72.39 -0.42 2.05
420 96.46 93.20 3.26 3.03
411 64.56 62.40 2.17 2.03
312 89.12 85.19 3.93 2.15
501 65.04 54.87 10.17 3.40
431 57.19 48.25 8.94 2.99
103 51.03 47.87 3.16 2.16
332 70.87 66.37 4.50 2.76
440 69.55 64.68 4.87 2.57
521 53.62 47.79 5.83 2.20
213 40.20 43.10 -2.90 2.14
600 56.13 60.17 -4.03 2.99
303 34.93 40.05 -5.12 1.22
512 48.75 55.56 -6.81 1.64
620 47.57 54.20 -6.63 1.60
611 33.15 37.77 4.63 1.12
b Titanium
100 8.85 8.93 -0.09 0.20
002 17.47 17.46 0.01 0.07
101 14.83 14.79 0.04 0.04
102 7.64 7.63 0.01 0.16
2–10 13.60 13.76 -0.16 0.34
103 11.00 11.17 -0.17 0.21
200 5.97 6.14 -0.17 0.61
2–12 12.51 12.10 0.41 0.30
201 10.17 10.31 0.14 0.34
004 11.44 11.56 -0.12 0.55
202 5.53 5.43 0.10 0.10
104 5.28 5.20 0.09 0.24
203 7.76 8.13 -0.36 0.41
3–10 4.52 4.49 0.02 0.85
3–11 7.66 7.57 0.09 0.35
2–14 8.08 8.52 -0.44 0.38
105 6.85 6.99 -0.15 0.28
3–12 3.97 4.03 -0.06 0.16
(continued)
92 3 Results and Discussion on Metals and Alloys
The powder X-ray intensity data were collected at the National Institute of
Interdisciplinary Science and Technology (NIIST), Trivandrum, India, using
X-PERT PRO (Philips, Netherlands) X-ray diffractometer with a monochromatic
incident beam, which offers pure Cu Ka1 radiation. The wavelength used for X-ray
intensity data collection is 0.154056 nm with 2h range of data collection from 10
to 120° and 0.05° step size.
The structural parameters of the metals were refined with the well-known
(Rietveld 1969) powder profile fitting methodology using the software package
JANA 2006 (Petřiček et al. 2006). The initial cell parameters and space group were
fixed as given in Table 3.9. The fitted profiles and the positions of Bragg peaks for
the six metals have been shown in Fig. 3.12a–f. In these figures, the dots represent
the observed powder patterns and the continuous lines represent the calculated
powder profiles of the respective samples. The small vertical lines indicate Bragg
positions. The difference between the observed profile and fitted calculated profile
is shown at the bottom of each figure. The refined profiles show good matching
between the observed and calculated profiles for all the three systems.
The refined structure factors of the samples given in Tables 3.10a–f show very
good agreement of the observed and calculated structure factors with very low
statistical errors.
The refined structural parameters from the software JANA 2006 are given in
Table 3.11. The refinement indices are low which show the accuracy of the
refinement (Rp = 10.03, 8.92, 20.89, 9.08, 15.84, 17.11% respectively, for sam-
ples Mg, Ti, Fe, Zn, Sn and Te).
The refined structure factors were used for MEM analysis. In the present calcu-
lation, the unit cell was divided into 643 pixels and the initial electron density at
each pixel was fixed uniformly as Z/a30, where Z is the number of electrons in the
unit cell. The electron density is evaluated by carefully selecting the Lagrange
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 97
multiplier in each case such that the convergence criterion C becomes unity after
performing minimum number of iterations. The resulting parameters of MEM
computations are given in Table 3.12.
The electron density distributions of Mg, Ti, Fe, Zn, Sn and Te have been
mapped using the MEM electron density values obtained through refinements.
Three-dimensional MEM electron density isosurface of Mg, Ti, Fe, Zn, Sn and Te
are shown in Fig. 3.13a–f respectively.
Figure 3.14a–f show the MEM electron density distributions of Mg, Ti, Fe, Zn,
Sn and Te on the (100) plane of the unit cell.
Figure 3.15a–f show the MEM electron density distributions of Mg, Ti, Fe, Zn,
Sn and Te on the (110) plane.
Figure 3.16a, b, c show the one-dimensional variation of the electron densities
of Mg, Ti, Fe, Zn, Sn along the three directions [100], [110] and [111] of the unit
cell, respectively.
98 3 Results and Discussion on Metals and Alloys
Fig. 3.14 a MEM electron density of magnesium on (100) plane. Contour lines are between 0.02
and 2.0 e/Å3. Contour interval is 0.33 e/Å3. b MEM electron density of titanium on the (100)
plane. Contour lines are between 0.05 and 5.0 e/Å3. Contour interval is 0.33 e/Å3. c MEM
electron density of iron on (100) plane. Contour lines are between 0.03 and 5.0 e/Å3. Contour
interval is 0.33 e/Å3. d MEM electron density of zinc on the (100) plane. Contour lines are
between 0.05 and 5.0 e/Å3. Contour interval is 0.33 e/Å3. e MEM electron density of tin on (100)
plane. Contour lines are between 0.05 and 5.0 e/Å3. Contour interval is 0.33 e/Å3. f MEM
electron density of tellurium on the (100) plane. Contour lines are between 0.05 and 5.0 e/Å3.
Contour interval is 0.001 e/Å3
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 99
Fig. 3.15 a MEM electron density of magnesium on (110) plane. Contour lines are between 0.0
and 2.0 e/Å3. Contour interval is 0.02 e/Å3. b MEM electron density of titanium on the (110)
plane. Contour lines are between 0.05 and 5.0 e/Å3. Contour interval is 0.001 e/Å3. c MEM
electron density of iron on (110) plane. Contour lines are 0.0 and 5.0 e/Å3. Contour interval is
0.03 e/Å3. d MEM electron density of zinc on (110) plane. Contour lines are between 0.00 and
10.0 e/Å3. Contour interval is 0.05 e/Å3. e MEM electron density of tin on (110) plane. Contour
lines are between 0.00 and 5.0 e/Å3. Contour interval is 0.05 e/Å3. f MEM electron density of
tellurium on the (110) plane. Contour lines are between 0.00 and 5.0 e/Å3. Contour interval is
0.05 e/Å3
100 3 Results and Discussion on Metals and Alloys
(a) 0.4 Fe Sn
(b) 0.3 Sn
Sn Ti
Fe Fe
Te Mg
0.3 Mg Te
Zn Zn
Ti 0.2
Fe Ti
Te
0.2 Ti
Te 0.1
Mg
0.1
Mg
Sn
Zn
Zn
0.0 0.0
0.40 1.05 1.70 2.35 3.00 0 1 2 3
(c) 0.6 Fe Sn
Ti
0.5 Fe
Zn Zn
Te
Electron density (e/Å 3 )
Sn Mg
0.4 Te
0.3
0.2 Ti
Mg
0.1
0.0
0.20 0.85 1.50 2.15 2.80
Fig. 3.16 a One-dimensional electron density profiles of Mg, Ti, Fe, Zn, Sn, Te along the [100]
direction of the unit cell. b One-dimensional electron density profiles of Mg, Ti, Fe, Zn, Sn, Te
along the [110] direction of the unit cell. c One-dimensional electron density profiles of Mg, Ti,
Fe, Zn, Sn, Te along the [111] direction of the unit cell
3.3.7.1 MEM
The fitted profiles of Rietveld refinements (Fig. 3.12a–e for Mg, Ti, Fe, Zn, Sn and
Te respectively) give a clear picture of the quality of the sample as well as the data.
The accuracy of refinements can be judged from the difference curve between the
observed and calculated intensity, which is almost a straight line. The results of the
structural refinement (Petřiček et al. 2000) given in Table 3.10 shows reasonable
values of the Debye–Waller factor for Mg (1.46 Å2) (Peng et al. 1996), Fe
(0.350 Å2) (Mohanlal 1979), Zn (1.26 Å2) (Bashir et al. 1988) and Te (0.74 Å2)
(Bashir et al. 1988]) obtainable from a powder data set, which is comparable to
those reported. The reliability indices for all the powder samples are very low
indicating the accuracy of the refinements (RP = 10.03, 8.92, 20.89, 9.08, 15.84,
17.11% for Mg, Ti, Fe, Zn, Sn and Te, respectively).
The MEM electron density distribution map of Mg as given in Fig. 3.14a for
(100) plane reveals the core of the Mg atom to be spherical and shows that the
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 101
Table 3.13 Electron density profiles along the three directions of the unit cells
Direction metal [100] [110] [111]
Position (Å) ED (e/Å3) Position (Å) ED (e/Å3) Position (Å) ED (e/Å3)
Mg 0.6527 0.06 0.6751 0.06 0.5741 0.10
Ti 1.4772 0.19 1.4772 0.19 2.7721 0.19
Fe 1.4346 0.16 1.5850 0.13 1.2424 0.52
Zn 0.4593 0.00 0.5011 0.00 3.6121 0.12
Sn 1.6319 0.33 2.1796 0.02 2.1987 0.02
Te 1.1857 0.13 0.9764 0.05 1.1601 0.17
ED Electron density
valence region of Mg is not spherical but smearing of charges is seen. The MEM
electron density distribution map of Ti as given in Fig. 3.14b for (100) plane
reveals the core of the Ti atom to be spherical but the valence region of Ti shows a
concave nature. The MEM electron density distribution map of Fe as given in
Fig. 3.14c for (100) plane reveals the core of the Fe atom being spherical and this
sphericity persists even at slightly larger distances away from the centre. The
distribution of charges all over the plane indicates the distribution of electrons and
the charges at the edge centres on (100) plane showing the distribution from the
perpendicular face-centred atoms and indicating the extension of the spread of the
charges. The MEM electron density distribution map of Zn as given in Fig. 3.14d
for (100) planes reveals the distortion of the core atom which may be due to the
high Debye–Waller factor as seen in Table 3.11. The distortion extends to the
valence region which is also concave. The MEM electron density distribution map
of Sn as given in Fig. 3.14e for (100) plane reveals the distortion in the core region
which is again substantiated by the high Debye–Waller factor as seen in Table 3.11.
The atomic number of Sn is higher and hence, more diffuse distribution of charges
is seen in the electron density map. Due to high-density of charges, there is
smearing of charges in the valence region. The MEM electron density distribution
map of Te as given in Fig. 3.14d for (100) planes reveals the distorted core region
which may again be due to the high Debye–Waller factor as seen in Table 3.11. The
atomic number of Te is high which can be visualised from the dense distribution of
charges in the electron density map. Since the atomic numbers of Sn and Te are
higher than that of Mg, Ti, Fe and Zn, more diffuse distribution of charges are seen
in the electron density maps of Sn and Te on the (100) plane.
The electron density map of Fe on (110) plane shows (Fig. 3.15c) highly
concentrated charges at the core of the atoms, and the distribution of charges on
places other than the atomic positions, due to the valence, free electrons. The
electron density map of Mg, Ti, Zn, Fe on the (110) plane (Fig. 3.15a, d), show
similar trends as observed in (100) plane. The electron density map of Sn and Te
on the (110) plane (Fig. 3.15e, f), show less dense charge distribution compared to
(100) plane.
The one-dimensional profiles of electron density constructed along [100], [110]
and [111] directions for Mg, Ti, Fe, Zn, Sn and Te are shown in Fig. 3.16a–f,
102 3 Results and Discussion on Metals and Alloys
respectively. The positions of minimum electron densities and density values have
been given in Table 3.14.
The electron densities at the saddle points along these three directions have
been given in Table 3.13.
The mid-bond density in Ti is the largest value (Table 3.13) along [110]
direction among the six metal systems studied. The Debye–Waller factor is less for
Ti when compared with Mg, Fe, Zn, Sn and Te as shown in Table 3.11. Since the
thermal vibration is less the bonding interaction is more significant in Ti. Hence,
obviously, the bonding interaction is stronger in Ti as evidenced by the higher
electron density along [110] direction. From Table 3.13, Zn has minimum mid-
bond density of 0.00 10-3 e/Å3 at a distance of 0.4593 Å, among the metals
studied along [110] and this may be due to the high Debye–Waller factor as
observed in Table 3.11. Sn has less mid bond density 0.02 10-3 e/Å3 at a distance
of 2.1796 Å, along [110] due to the high Debye–Waller factor as observed in
Table 3.11. This substantiates the two-dimensional electron density map of tin
(Fig. 3.15e) where sparse distribution of electron density is observed along [110]
direction. From Tables 3.11 and 3.13, it can be observed that the metals with
higher Debye–Waller factors have low mid-bond density along [110] direction
which confirms the higher thermal vibration.
But along the [100] (nonbonding), Sn has higher electron density of 0.3310-3
e/Å3 at a distance of 1.6319 Å whose atomic number is higher than Mg, Ti, Fe, Zn.
This can also be visualized from Fig. 3.14e which gives the MEM electron density
distribution of Sn on the (100) plane.
Along [111] direction Fe has higher electron density of 0.52 e/Å3 at a distance
of 1.2424 Å which is also visualized in Fig. 3.15c which gives the MEM electron
density distribution of Fe on the (110) plane.
Hence, from the qualitative and quantitative analysis of the MEM electron
densities, the bonding in Mg, Ti, Fe, Zn, Sn and Te is predicted to be predomi-
nantly metallic and the interaction of charges along other non-bonding directions
seems to be heavier due to the metallic nature.
3.3.7.2 PDF
In this study, the observed PDFs of Mg, Ti, Fe, Zn, Sn and Te were obtained from
the raw intensity data using a software programme PDFgetX (Jeong et al. 2001)
after performing corrections such as multiple scattering correction, polarisation
correction, absorption correction, normalisation correction and Compton correc-
tion. A comparison between the observed and calculated PDF is carried out using
the software package PDFFIT (Farrow et al. 2007). Figure 3.17a–f give the
observed and calculated PDF profiles along with the difference between them for
Mg, Ti, Fe, Zn, Sn and Te respectively.
The bond lengths deduced from the PDF analysis have been compared with
those calculated using software programme Grenoble Thermal Ellipsoids Plot
Programme (GRETEP) is a Windows interactive programme that displays the list
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 103
of calculated bonds, bond length between two atoms and the distances and angles
around an atom) (GRETEP). The refined pair distribution function (the Fourier
transform of, the reduced structure factor) given in Fig. 3.17a–f for Mg, Ti, Fe, Zn,
Sn and Te show good matching of observed and calculated PDFs. The difference
between the observed and calculated PDFs are also shown in the figures. The peaks
in the PDF profile correspond to the nearest neighbour distances given in
Table 3.14.
The first neighbour distance as observed from Table 3.15 decreases as the
atomic number increases for all the metals. The difference between observed and
calculated nearest neighbour distances are also given in the table. These differ-
ences are not too high considering the fact that we have used only data with limited
Q values. There may be local undulations in the structure of these six systems that
will lead to slightly higher neighbour differences which correspond to lattice repeat
distances. The local disorder would have arisen due to several factors viz., X-ray
diffraction instrumental error or fine impurity in the sample. This is reflected in the
local structural analysis only like PDF.
104 3 Results and Discussion on Metals and Alloys
G(r)(Å -2 )
G(r)(Å-2 )
-0 -2
-1
-4
-2
-3 -6
-4
-8
-5
-6 -10
5 10 15 20 25 30 2 6 10 14 18 22 26 30
r(Å) r(Å)
G(r)(Å-2 )
G(r)(Å -2 )
-0
0 -2
-4
-6
-10 -8
-10
-12
-20
2 6 10 14 18 22 26 30 4.0 10.5 17.0 23.5 30.0
r(Å) r(Å)
(e) 8
Obs (f) 3
Cal
Cal Obs
6 Diff 2 Diff
4 1
2 -0
G(r)(Å -2 )
G(r)(Å-2 )
-0
-1
-2
-2
-4
-3
-6
-8 -4
-10 -5
-12 -6
2 6 10 14 18 22 26 30 3 5 7 9 11 13 15 17 19 21 23 25
r(Å) r(Å)
Fig. 3.17 a The observed and calculated PDFs of Mg with difference profile. b The observed
and calculated PDF of Ti with difference profile. c The observed and calculated PDFs of Fe with
difference profile. d The observed and calculated PDFs of Zn with difference profile. e The
observed and calculated PDFs of Sn with difference profile. f The observed and calculated PDFs
of Te with difference profile
3.3.8 Conclusion
A clear understanding of the structure of Mg, Ti, Fe, Zn, Sn and Te has been gained
from this work, in terms of the cell parameters, thermal vibrations and electron
density distribution in the bonding region using MEM and local structural analysis
using PDF. The electron density, bonding and charge transfer studies analysed in
this work will give fruitful information to researchers in the fields of physics,
3.3 Magnesium, Titanium, Iron, Zinc, Tin and Tellurium 105
chemistry, materials science, metallurgy, etc. These properties can be properly used
for the proper engineering of these technologically important materials.
3.4.1 Introduction
Metal alloys are highly useful because they have specific properties and charac-
teristics that are more attractive than those of pure, elemental metals. For example,
some alloys possess high strength, some have low melting points, some are
refractory with high melting temperatures, some are especially resistant to cor-
rosion and some have desirable magnetic, thermal or electrical properties. These
characteristics arise from both the internal and electronic structure of the alloys.
Hence, a study on the electronic structure of alloys using the most versatile
techniques like MEM (Collins 1982) and multipole method (Hansen and Coppens
1978) available today is worthwhile. If the tool available for the analysis yields
highly precise information, then it is appropriate to apply it to precise data sets
available as has been done in this work, thereby the methodology can also be
tested. Hence, a study on the electronic structure of the metal alloys CoAl and
NiAl has been carried out using accurate single crystal X-ray data sets reported
earlier by Weiss (1966). Electron density distribution and bonding in different
materials using MEM have been analysed for different materials e.g., GaAs
(Saravanan et al. 2003), InP (Israel et al. 2004), Oxides (Israel et al. 2003), ionic
systems (Israel et al. 2003). Research works have also been reported by Sakata and
Sato (1990), Saka and Kato (1986), Yamamura et al. (1998), Ikeda et al. (1998),
Yamamoto et al. (1996). In order to confirm the results and to elucidate the
distribution of valence electrons, and the contraction/expansion of atomic shells,
multipole analysis of the electron densities was carried out using the software
package JANA 2000 (Petřiček et al. 2000). Recently, the multipole analysis of the
charge densities and bonding has been widely utilised to study the electronic
structure of materials (Coppens and Volkov 2004; Poulsen et al. 2004; Frils et al.
2004; Pillet et al. 2004; Marabello et al. 2004). In this work, the multipole model
proposed by Hansen and Coppens (1978) has been used for elucidating the elec-
tronic structure of CoAl and NiAl.
The precise electron density distribution and bonding in metal alloys CoAl and
NiAl is characterised using MEM and multipole method. Reported X-ray single
crystal data have been used for this purpose. Clear evidence of the metal bonding
between the constituent atoms in these two systems has been obtained. The mid-
106 3 Results and Discussion on Metals and Alloys
bond electron densities in these systems are found to be 0.358 and 0.251 e/Å3
respectively, for CoAl and NiAl in the MEM analysis. The two-dimensional maps
and one-dimensional electron density profiles have been constructed and analysed.
The thermal vibration of the individual atoms Co, Ni and Al has also been studied
and reported. The contraction of atoms in CoAl and expansion of Ni and con-
traction of Al atom in NiAl is found from the multipole analysis, in line with the
MEM electron density distribution.
Weiss (1966) has reported single crystal X-ray structure factor data of CoAl and
NiAl collected using AgKa and CoKa radiations. Among these highly precise data
sets, structure factor sets collected using AgKa were utilised in this work, because,
more number of Bragg reflections is available using this radiation.
The two sets of data were subjected to preliminary data analysis and the
structure of CoAl and NiAl was refined using the software GSAS (Larson and Von
Dreele 2004) (General Structure Analysis System). Further, the same data were
refined using a code written by Dr. R. Saravanan using an algorithm which uses the
standard least-squares refinement. The individual thermal vibration parameters of
atoms were refined both in GSAS (Larson and Von Dreele 2011) and also using
the code by Dr. R. Saravanan (hereafter referred as RS refinements). Exact
matching of the structure factors and structural parameters were realised in both
the refinements. The refined structure factors were used for MEM analysis. Mul-
tipole refinements were carried out using the refinement software package JANA
2000 (Petřiček et al. 2000) to refine the thermal parameters and the j0 parameter
which gives information about the valence electronic status of the atomic shells.
CoAl and NiAl belong to the CsCl structure in which there are 2 atoms/cell and
the fractional atomic coordinates are 0, 0, 0, and , , . The cell constants are
2.8619 and 2.8864 Å respectively, for CoAl and NiAl, (Weiss 1966). Figure 3.18
shows the structure of NiAl.
It has been reported that no order–disorder transition occurs in these alloys and
there is a possibility of squeezing of charge densities and directionality in the outer
electron charge density. A total of 11 Bragg reflections have been reported (Weiss
1966) for both these systems. The data were found to be precise and the reported r
(Fobs) values are very low indicating the precision of the data. Using these data,
first, the refinements were carried out with GSAS and the results of these refine-
ments are given in Table 3.15 both for CoAl and NiAl. Then RS refinements were
carried out, the results of which are also given in Table 3.15. The quantity
3.4 CoAl and NiAl Metal Alloys 107
Table 3.15 The Debye–Waller factors in NiAl and CoAl in GSAS and RS Refinements
Refinement type System
NiAl CoAl
GSAS BNi (Å2) 0.221(0.070) BCo (Å2) 0.368(0.049)
BAl (Å)2 0.421(0.144) BAl (Å2) 0.591(0.088)
R (%) 1.5 R (%) 1.1
wR (%) 1.9 wR (%) 1.2
RS BNi (Å2) 0.279(0.073) BCo (Å2) 0.353(0.082)
BAl (Å)2 0.389(0.120) BAl (Å2) 0.528(0.130)
Boverall (Å)2 0.288(0.035) Boverall (Å)2 0.382(0.042)
R (%) 1.413 R (%) 1.134
B Debye–Waller factor
[001]
[001]
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 Co 0.1 Ni
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[010] [010]
(c) 1.0 (d) 1.0
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
[001]
[001]
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 Co 0.1 Ni
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[010] [010]
(e) 1.0 (f) 1.0
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
[001]
[001]
Al Al
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 Co 0.1 Ni
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[010] [010]
Fig. 3.19 a MEM high electron density distribution of CoAl on the (100) plane. Contour Range
is between 0.0 and 15.0 e/Å3. Contour interval is 0.75 e/Å3. b MEM high electron density
distribution of NiAl on the (100) plane. Contour range is between 0.0 and 15.0 e/Å3. Contour
interval is 0.75 e/Å3. c MEM low electron density distribution of CoAl on the (100) plane.
Contour range is between 0.0 and 2.5 e/Å3. Contour interval is 0.125 e/Å3. d MEM low electron
density distribution of NiAl on the (100) plane. Contour range is between 0.0 and 2.5 e/Å3.
Contour interval is 0.125 e/Å3. e MEM electron density distribution of CoAl on the (100) plane
showing the extent of electron clouds of Al. Contour range is between 0.0 and 0.2 e/Å3. Contour
interval is 0.008 e/Å3. f MEM electron density distribution of NiAl on the (100) plane showing
the extent of electron clouds of Al. Contour range is between 0.0 and 0.2 e/Å3. Contour interval is
0.008 e/Å3
3.4 CoAl and NiAl Metal Alloys 109
1.0 1.0
(a) Al (b) Al
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
[001]
[001]
Co Ni
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[110] [110]
1.0 1.0
(c) Al (d) Al
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
[001]
[001]
Co Ni
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
[110] [110]
Fig. 3.20 a MEM high electron density distribution of CoAl on the (110) Plane. Contour range:
0.0–6.6 e/Å3. Contour interval is 0.33 e/Å3. The plane has been shifted half the unit cell along the
z direction. b MEM high electron density distribution of NiAl on the (110) Plane. Contour range:
0.0–6.6 e/Å3. Contour interval is 0.33 e/Å3. The plane has been shifted half the unit cell along the
z direction. c MEM low electron density distribution of CoAl on the (110) Plane. Contour range:
0.0–1.2 e/Å3. Contour interval is 0.06 e/Å3. The plane has been shifted half the unit cell along the
z direction. d MEM low electron density distribution of NiAl on the (110) Plane. Contour range:
0.0–1.2 e/Å3. Contour interval is 0.06 e/Å3. The plane has been shifted half the unit cell along the
z direction
P
minimised in RS refinements was D ¼ hkl Whkl ðjFobs j jkFcal jÞ2 ; where whkl is
the weight assigned on each observation, Fobs is the observed structure factor and
Fcal is the calculated structure
h P factor. Thei estimated standard deviation in any
2
parameter is given by rp ¼ bii wr¼1 ðDF rÞ
mn ; where bii is the diagonal element of
the inverse matrix, wr is the weight of the rth factor, where DFr is the difference
between Fobs and Fcal, m is the number of observations
P and n is the number of
jjFobs jjFcal jj
parameters. The reliability index given by R ¼ P of these refinements is
jFobs j
also given in Table 3.15. The average thermal vibration of the overall system is
studied and reported in Table 3.15. In RS refinements, the analytical coefficients of
the five-parameter model by Waasmaier and Kirfel (1995) were used for the
calculation of the atomic scattering factors.
MEM studies on the electron density distribution of these two systems were
carried out by dividing both the unit cells into 64 9 64 9 64 pixels. The initial
electron densities were fixed as (F000/a3o), where F000 is the total number of
110 3 Results and Discussion on Metals and Alloys
(a) (b)
2.0
400 High density profiles 1.8 Low density profiles
1.6
Electron Density (e/Å 3)
200 1.0
0.8
100 0.6
0.4
0 0.2
Co Co
0.0
0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 0.5500 0.9875 1.4250 1.8625 2.3000
Distance (Å) Distance (Å)
(c) (d)
500 2.0
High density profiles Low density profiles
1.8
400 111
1.6 111
Electron Density (e/Å3)
110
Electron Density (e/Å 3)
110
100 100
1.4
300
1.2
1.0
200
0.8
0.6
100
0.4
0 Ni 0.2
Ni
0.0
0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 0.5500 0.9875 1.4250 1.8625 2.3000
Distance (Å) Distance (Å)
Fig. 3.21 a One-dimensional high electron density distribution of CoAl along the [100], [110]
and [111] directions of the unit cell. b One-dimensional low electron density distribution of CoAl
along the [100], [110] and [111] directions of the unit cell. c One-dimensional high electron
density distribution of NiAl along the [100], [110] and [111] directions of the unit cell. d One-
dimensional low electron density distribution of NiAl along the [100], [110] and [111] directions
of the unit cell
3.4 CoAl and NiAl Metal Alloys 111
Table 3.17 Profile electron densities along the three directions of the unit cells of CoAl and
NiAl
Direction system [100] [110] [111]
Position (Å) ED (e/Å3) Position (Å) ED (e/Å3) Position (Å) ED (e/Å3)
CoAl 1.431 0.079 2.024 0.083 1.317 0.358
NiAl 1.443 0.050 2.041 0.078 1.328 0.251
ED Electron density
5
[Fobs(K)-FMEM(K)]/ Fobs (K)
(a)
4
3
2
1
0
-1
-2
-3
-4
-5
0.1 0.2 0.3 0.4 0.5 0.6 0.7
sin (Å-1)
5
[Fobs(K)-FMEM(K)]/ Fobs (K)
(b)
4
3
2
1
0
-1
-2
-3
-4
-5
0.1 0.2 0.3 0.4 0.5 0.6 0.7
-1
sin (Å )
Fig. 3.22 a The distribution of errors in the reciprocal space of CoAl. b The distribution of
errors in the reciprocal space of NiAl
electrons in the respective unit cells and ao is the cell dimension. The resolution of
the MEM maps turns out to be 0.045 Å/pixel along the three directions of the unit
cells. The Lagrange parameter k was chosen so that 188 and 223 cycles of
refinement were required for CoAl and NiAl respectively, to achieve the con-
vergence criterion C = 1; C being,
112 3 Results and Discussion on Metals and Alloys
Fig. 3.23 a Static Multipole Deformation map of CoAl on (100) plane. Dotted lines indicate
negative electron densities. b Static Multipole Deformation map of CoAl on (00) plane. Dotted
lines indicate negative electron densities. c Dynamic Multipole Deformation map of CoAl on
(100) plane. Dotted lines indicate negative electron densities. d Dynamic Multipole Deformation
of CoAl on (00) plane. Dotted lines indicate negative electron densities
3.4 CoAl and NiAl Metal Alloys 113
Fig. 3.24 a Static Multipole Deformation map of NiAl on (100) plane. Dotted lines indicate
negative electron densities. b Static Multipole Deformation map of NiAl on (00) plane. Dotted
lines indicate negative electron densities. c Dynamic Multipole Deformation map of NiAl on
(100) plane. Dotted lines indicate negative electron densities. d Dynamic Multipole Deformation
map of NiAl on (00) plane. Dotted lines indicate negative electron densities
1 X
C¼
ð3:5Þ
N
rðFobs ðkÞÞ
Figure 3.19e, f show the extent to which the aluminium atom is distributed as
evidenced by the electron clouds on the (100) plane.
Figure 3.20a, b show the MEM high electron density distribution on the (110)
plane of CoAl and NiAl respectively, wherein the shifted atoms (half the unit cell
along the Z direction) have been shown.
Figure 3.20c, d show the MEM low electron density distribution on the (110)
plane of CoAl and NiAl.
Figure 3.21a, b show the one-dimensional variation of the high and low elec-
tron densities along the three directions ([100], [110] and [111]) of the unit cell of
CoAl and Figs. 3.21c, d represent the one-dimensional variation of high and low
electron densities of NiAl.
The electron densities at the saddle points along these three directions are given
in Table 3.17.
Based on these two-dimensional MEM maps and one-dimensional profiles, a
discussion is formulated and given in the following section.
The errors in the observations have also been studied in this work and repre-
sented in Fig. 3.22a, b, which show the distribution of errors in the reciprocal
space CoAl and NiAl respectively.
The results of the multipole refinements done using JANA 2000 are presented
in Table 3.18.
Hansen and Coppens (1978) proposed a modified electron density model with
the option that allows the refinement of population parameters at various orbital
levels where the atomic density is described as a series expansion in real spherical
harmonic functions through fourth order Ylm. According to this model, the charge
density in a crystal is written as the superposition of harmonically vibrating
aspherical atomic density distribution convolving with the Gaussian thermal dis-
placement distribution as
X
atom
qðr Þ ¼ qk !r !
r k!
u tk !
u ; ð3:6Þ
k
where tk ð~
uÞ is the Gaussian term and indicates a convolution. The atomic charge
density is then defined as
X4 Xl !
r
r ¼ Pc qcore !
q ! r þ Pv k03 qvalence ðk0 r Þ þ k003 Rj ðk00 :Þ P Y
¼l lm lm r
l¼0
ð3:7Þ
where Pc, Pv and Plm are population coefficients. Canonical Hartree–Fock atomic
orbitals of the free atoms normalised to one electron can be used for the con-
struction of qcore and qvalence but the valence function is allowed to expand and
contract by the adjustment of the variable parameters k0 and k00: The effect of the
temperature can be distinguished from the convolution and the deconvolution of
3.4 CoAl and NiAl Metal Alloys 115
the valence region of Co is clearly not spherical and the four corner atoms on (100)
plane show concave outer charges. It appears that the eight Co atoms give way for
the thermal vibration of the Al atom located at the body centre. It is noted from
Table 3.15 that BAl in CoAl is higher than BCo. In the case of NiAl (Fig. 3.19d), no
such concave outer charges are seen, since the thermal vibration of both Ni and Al
atoms are lower compared to Co and Al in CoAl. Figure 3.19e, f show very low
electron densities on (100) plane of CoAl and NiAl. The body-centred Al charges
touching the (100) plane of the unit cell is visible in Fig. 3.19e. The valence charges
of both the atoms (Co and Al) do not repel each other. In fact, charge clouds of both
the atoms touch each other, but with no resemblance of an attraction. This should be
due to the characteristics of metal bonding. Moreover, in multipole analysis, the
expansion/contraction parameter j0 is estimated to be 0.953 and 0.92 for Co and Al
atoms, which means expansion of both the atoms. The expansion is more in Al atom.
Hence, the above said trends in the low-density regions are visible as seen from
Fig. 3.19e. In NiAl (Fig. 3.19f), the charge accumulation is less dense due to
reduced thermal vibration of both the atoms in NiAl. Again from multipole
refinement, the j0 parameter is found to be 0.984 and 1.258 for Ni and Al atoms
respectively. This clearly shows that Ni atom is expanded slightly (not as in Co), but
the Al atom is contracted to a large extent (j0 = 1.258). These facts are clearly
visualized in the MEM maps. Thus, the results from one analysis (MEM) could be
correlated to the other (multipole). Figure 3.20a, b show the electron density dis-
tribution of CoAl and NiAl in the high-density regions of the (110) plane. The core
regions of both the atoms are seen to be spherical both in CoAl and NiAl. Again, no
repulsion of charges is seen in these systems. The outer regions of both Co and Al
atoms are more smeared in CoAl compared to that of Ni and Al atoms in NiAl, a fact
due to the increased thermal vibrations in CoAl. Figure 3.20c, d show the low
electron density distribution of CoAl and NiAl on the (110) plane. Charges are seen
around the periphery of Co and Al atoms in CoAl (Fig. 3.20c) and Ni and Al atoms in
NiAl (Fig. 3.20d). The bonding charges show a metallic nature in both the systems.
Figure 3.21a–d show the one-dimensional electron density distribution along
the [100], [110] and [111] directions of the CoAl and NiAl unit cells respectively.
Both high- and low-density variation of charges is shown in these figures. The
high-density peaks in Fig. 3.21a, c show the peak densities along these three
directions which are larger in NiAl. The low-density maps show that the density
minima do not occur at the same points in both the cases. Moreover, the distri-
bution of charges along the three directions does not overlap. These facts suggest
that the outer regions are not perfectly spherical. The profile along [100] direction
is perfectly symmetrical due to the distribution of charges from like atoms. The
distribution along [111] direction shows the bonding nature between unlike atoms.
The previous discussion clears that there is attraction/repulsion like interaction
between Co and Al and Ni and Al. Yet, the mid-bond charges are considerably
higher as seen from Table 3.17. CoAl has a mid-bond density value of 0.358 e/Å3
3.4 CoAl and NiAl Metal Alloys 117
and NiAl has a value of about 0.251 e/Å3. If the systems are ionic one cannot
expect higher mid-bond values as in the present cases. If the systems are covalent,
one can expect values higher than the present values due to overlapping of charges.
Hence, from the qualitative and quantitative analysis of the MEM electron den-
sities, the bonding in CoAl and NiAl is predicted to be predominantly metallic.
The strength of the electron density is higher in the case of CoAl, which has lower
number of electrons in the unit cell compared to NiAl. This is due to the relatively
larger thermal vibration of aluminium and smearing of charges in CoAl. The
increase in the mid-bond density is about 0.107 e/Å3 consistent with the increased
thermal vibration of atoms, particularly aluminium, in CoAl. Hence, it is important
to study precisely the thermal vibration of atoms in electron density analyses,
particularly in comparing the bonding densities in systems with similar constituent
atoms.
Figure 3.22a, b shows the distribution of errors in the observations in the
reciprocal space of CoAl and NiAl respectively. The error in the form of [Fobs
(K)-FMEM(K)]/r(Fobs(K)), is distributed uniformly around ±2 indicating that
there is no bias in the experimental data or the calculations and justifies the use of
these experimental data for the present analysis.
The B values from the multipole refinement are in close agreement with those
obtained from other refinements. Although all the B values are comparable, those
from the multipole analysis give accurate interpretation of the electron density
distribution in both the systems.
The SMD and DMD maps obtained from multipole refinements show that
(Figs. 3.23, 3.24) more pronounced deformation in density is at the valence
regions than at the core regions for both the systems. Figure 3.23a (SMD—
CoAl—100 plane) shows features similar to those seen in Fig. 3.19e for CoAl. The
expansion of Al atom and the concave nature of charge clouds of Co atom are
visible in this map. Figure 3.24a (NiAl—SMD 100 plane) does not show any
concave nature/repulsive interaction, because of the heavy contraction of Al atom.
The map in Fig. 3.23b (SMD—CoAl— 0 0 plane) shows the expansion of both
the atoms and Fig. 3.24b (SMD—NiAl—00 plane) shows nothing but Al atom
because of its contracted nature. Further, the expansion of Ni atom (j = 0.984) is
not large compared to Co (0.953) atom. Hence, only the distribution of Al atom is
seen in this map. The DMD maps show trends in line with the expansion/con-
traction of atoms in the two systems. Comparing Fig. 3.23b (SMD of CoAl) and
Fig. 3.23d (DMD of CoAl), it is seen that the SMD map shows a number of
positive contours at the edge centres and at the centre of (100) plane. These are the
residual densities due to the thermal vibration of aluminium atom. The core of
aluminium atom in the Fig. 3.23b is full of negative charges. The DMD map
shows (Fig. 3.23d) the uncompensated thermal vibration of the valence regions
that show up as the positive contours starting from the core. This suggests slight
expansion of the aluminium atom.
118 3 Results and Discussion on Metals and Alloys
Figure 3.24b (SMD of NiAl) and 3.24d (DMD of NiAl) can be compared in a
similar manner. Figure 3.24b shows positive electron densities not only at the core
but also at far-valence regions, even after the deconvolution process, suggesting
shrinkage of the aluminium atom in NiAl. The DMD (Fig. 3.24d) shows a well
compensated core region (without negative charge densities) but uncompensated
positive islands of valence region confirming the contraction of aluminium atom in
NiAl.
3.4.6 Conclusion
3.5.1 Introduction
Alloying one metal with other metals or non-metals often enhances its properties.
The physical properties of an alloy may not differ greatly from those of its ele-
ments, but engineering properties may be substantially different from those of the
constituent materials.
Whether in the form of element or alloy with other metals, nickel materials
have made significant contributions to present-day society. The favourable phys-
ical and mechanical properties as well as corrosion resistance characteristics of
nickel and its alloys provide aids to the solution of many industrial problems
(Teeple 1954). Alloying with nickel increases the impact strength and ductility of
sintered materials without adversely affecting their strength (Enríquez and Mathew
2003).
Nickel alloys have good corrosion resistance and heat resistance, the standard
alloy being used for electrical-resistance for heaters and electrical appliances is
Nickel–Chromium (Murray 1997). Chromium–nickel steel alloys have good
mechanical properties (Radomysel’’skii and Kholodnyi 1975). Conductive
3.5 Nickel Chromium (Ni80Cr20) 119
nichrome probe tips have been developed as functionalized nanotools for nano-
scale electrical testing and nano-welding, and are advantageous in competition
with other scanning probe tips because of their unique combination of high
resistance to oxidation, high hardness and relatively high resistivity (Peng et al.
2009). Ni80Cr20 is a commonly used material for application in thin-film strain
gauges (Kazi et al. 2003, 2006]. NiCr coatings on stainless steel substrates are used
as a high velocity oxy-fuel technique for corrosion applications (Ak et al. 2003).
NiCr cermet coatings are commonly utilised in the industry as a means to provide
wear resistance in high temperature (Wirojanupatump et al. 2001). NiCr-coated
super alloys show better oxidation resistance due to the formation of a compact
and adhesive thin Cr2O3 scale on the surface of the coating during oxidation
(Kamal et al. 2009).
For practical use in different applications, alloys need to be shaped at elevated
temperatures according to the requirements. It is therefore, essential to study the
changes in the microstructure after heat treatment of materials. One of the popular
heat treatment processes is annealing (Al-Quran and Al-Itawi 2010). Annealing
induces ductility, softens material, relieves internal stresses, refines the structure
by making it homogeneous and improves cold working properties. The resistance
characteristics of NiCr improve in the double annealing process (Jeng et al. 1990).
Annealing processes appear to be responsible for an increase in passivity of Ni/Cr
alloys toward oxidation (Kim et al. 2005).
Ball milling may be referred to as the breaking down of relatively coarse
materials to the ultimate fineness. At present, it has evolved as an important
method for the preparation of either materials with enhanced physical and
mechanical properties or indeed, new phases or new engineering materials. The
morphology of the powders is modified when they are subjected to collisions.
The initial ball-powder-ball collision causes ductile metal powders to flatten
into flake-like structures (Eskandarany 2001). The high-energy ball milling
technique is an effective method for producing powders for metal matrix
composites (Zwick and Lugscheider 2004). There exist methods to control the
microstructure of pellet and among others, mechanical milling of powder is
simple and easy to control. With this method, the characteristics of powder and
the resultant compaction behaviour and sinterability are properly controlled
(Na et al. 2002).
The alloy Ni80Cr20 was annealed and ball-milled to study the effect of thermal and
mechanical treatments on the local structure and on the electron density distri-
bution. The electron density between the atoms was studied by MEM and the local
structure using PDF. The electron density is found to be high for ball-milled
sample along the bonding direction. The particle sizes of the differently treated
120 3 Results and Discussion on Metals and Alloys
samples were realised by SEM and through XRD. Clear evidence of the effect of
ball milling is observed on the local structure and electron densities.
In this work, the purchased alloy Ni80Cr20 was subjected to different types of
mechanical and thermal treatments viz., ball milling and annealing. The three
different processes are as follows:
i. annealing the sample (hereafter sample A)
ii. ball milling (hereafter sample BM)
iii. ball milling and annealing (hereafter sample BMA).
The charge density and local structure variations of these three samples have
been analysed in this work. This analysis would reveal the properties of Ni80Cr20
suitable for its wide applications. The electronic structure analysis has been carried
out by the most versatile technique, MEM (Gull and Daniel 1978) and local and
average structure properties have been analysed using PDF (Petkov et al. 2002;
Xiangyun et al. 2005).
3.5.3 Experimental
Powder sample of Ni80Cr20 (99.9% purity) was purchased from Alfa Aesar,
Johnson’s Mathew Company. To anneal the sample, 20 mg of the sample was
placed in a small silica crucible and was annealed at 1023 K for 3 h in a pro-
grammable furnace (sample A). To prepare the ball-milled sample, the alloy was
milled at room temperature in a low energy laboratory ball mill for 40 h with a
milling speed of 170 rpm, using stainless steel balls (sample BM). Further, a part
of the ball-milled powder was annealed at 1023 K for 3 h (sample BMA). The low
energy laboratory ball mill used to mill the sample and the furnace in which the
samples were annealed are shown in Figs. 3.25a, b, respectively.
To analyse the bonding charge densities (structural behaviour) of Ni80Cr20
under various experimental conditions, powder X-ray diffraction data were col-
lected individually for sample A, sample BM, sample BMA. X-ray powder
diffraction measurements were performed using X-PERT PRO (Philips, Nether-
lands) at the National Institute of Interdisciplinary Science and Technology (NI-
IST), Trivandrum, India using CuKa1 radiation with a 2h range from 10 to 120°
and 0.05° step size.
Particle size is an important parameter that influences the physical properties
of materials. A particle may be made of several different crystallites. A crys-
tallite size is the minimum part of the material that diffracts coherently. Several
techniques could be used for the investigation of the particle size. X-ray
diffraction (XRD) and Scanning electron microscopy (SEM) techniques are
used in this work to obtain the particle size of the samples. From XRD data the
particle size was obtained as the particle size is related to the diffraction peak
3.5 Nickel Chromium (Ni80Cr20) 121
Fig. 3.26 a SEM picture of sample A. b SEM picture of sample BM. c SEM picture of sample
BMA
125.00
250 (a) (111)
A (b) (002) A
BM BM
BMA BMA
Intensity (counts)
Intensity (counts)
200 93.75
150
62.50
100
50 31.25
0
0.00
43.0 43.5 44.0 44.5 45.0 45.5 50.7 50.9 51.1 51.3 51.5 51.7 51.9 52.1 52.3
2 (degrees) 2 (degrees)
75.0000 65.0000
64.2857
(c) (202) A
54.1667
(d) (113)
A
BM
Intensity (counts)
BM
Intensity (counts)
BMA
BMA
53.5714
43.3333
42.8571
32.5000
32.1429
(222)
21.6667
21.4286
10.7143 10.8333
0.0000 0.0000
74 75 76 77 78 89 90 91 92 93 94 95 96 97 98 99 100
2 (degrees) 2 (degrees)
Fig. 3.27 a Comparison of (111) XRD peaks of the samples A, BM, BMA respectively.
b Comparison of (002) XRD peaks of the samples A, BM, BMA respectively. c Comparison of
(202) XRD peaks of the samples A, BM, BMA respectively. d Comparison of (113, 222) XRD
peaks of the samples A, BM, BMA respectively
BMA is with uniform morphology with particle size of 59.35 nm (Fig. 3.26c).
The obtained particle sizes from XRD and SEM are given in Table 3.19.
The sizes of the samples measured by XRD line broadening are about one-third
of the particle sizes calculated from SEM. The individual particles of the samples
may have been composed of coherently diffracting domains rather than each being
a single coherently diffracting crystal.
From Fig. 3.27a–d, it is observed that the intensity of samples decrease in the
order of A, BMA, BM respectively.
The broadening of the peaks is also observed in the case of sample BM. The
full width at half-maximum (FWHM) as observed in Table 3.20 indicates that
peaks are broader for sample BM compared to samples A and BMA. The
reduction in particle size due to ball milling as observed from SEM pictures is
substantiated.
3.5 Nickel Chromium (Ni80Cr20) 123
The structural parameters of the samples (A, BM, BMA) of Ni80Cr20 which have
FCC structure with four atoms/unit cell in the space group Fm-3 m were refined
with the well-known Rietveld (1969) powder profile fitting methodology using the
software package JANA 2006 (Petřiček et al. 2006). The fitted profiles and the
positions of Bragg peaks for the three samples have been shown in Fig. 3.28a–c.
In these figures, the dots represent the observed powder patterns and the con-
tinuous lines represent the calculated powder profiles of the respective samples.
The small vertical lines indicate Bragg positions. The difference between the
observed profile and fitted calculated profile is shown at the bottom of each figure.
The refined profiles show good matching between the observed and calculated
profiles for all the three systems. The refined structure factors of the samples given
in Table 3.21a–c shows very good agreement of the observed and calculated
structure factors with very low statistical errors. The composition of the sample
was also refined and was determined as Ni80.1Cr19.9 for the annealed sample
(sample A).
The refined structural parameters from the software JANA 2006 are given in
Table 3.22. The refinement indices are low which show the accuracy of
the refinement (Rp = 8.17, 8.35, 10.16% respectively for samples A, BM, BMA).
The cell parameters and thermal vibration parameters are observed to reduce in the
order of sample BM, A, BMA respectively. Ball milling and subsequent annealing
appears to have reduced the cell constant.
The refined structure factors of the three samples were used for the MEM analysis.
The software package PRIMA (Ruben and Fujio 2004) was used for the numerical
MEM computations. The two- and three-dimensional representation of the electron
densities were picturised using the software VESTA (Momma and Izumi 2006).
The MEM refinements were carried out by dividing the unit cell into
64 9 64 9 64 pixels.
The Lagrange parameter is suitably chosen so that the convergence criterion
C = 1 is reached after minimum number of iterations. The MEM parameters are
given in Table 3.23.
The electron density distributions of samples A, BM, BMA are mapped using
the electron density values obtained through refinements. Figure 3.29a–f shows the
three-dimensional electron densities of samples A, BM, BMA in the unit cell in
two different orientations.
Figure 3.30a–c shows the two-dimensional electron densities on (100) plane in
the unit cell for the three samples respectively.
124 3 Results and Discussion on Metals and Alloys
Fig. 3.28 a The fitted XRD profiles and the positions of Bragg peaks for sample A. b The fitted
XRD profiles and the positions of Bragg peaks for sample BM. c The fitted XRD profiles and the
positions of Bragg peaks for sample BMA
3.5 Nickel Chromium (Ni80Cr20) 125
Table 3.21 Observed and calculated structure factors for a Sample A, b Sample BM, c Sample
BMA
hkl Fo Fc Fo - Fc r (Fo)
a Sample A
111 65.62 65.67 -0.05 0.94
002 59.48 59.67 -0.19 0.93
202 43.77 42.81 0.96 0.78
113 34.36 34.60 -0.24 0.79
22 2 31.55 32.40 -0.85 1.31
b Sample BM
111 65.10 65.15 -0.06 0.93
0 02 59.01 59.09 -0.08 1.11
202 41.87 42.07 -0.20 1.03
113 33.74 33.80 -0.06 0.86
22 2 30.75 31.59 -0.84 1.59
c Sample BMA
111 66.85 67.11 -0.26 0.96
00 2 61.16 61.43 -0.27 1.21
202 45.17 45.37 -0.20 1.11
113 37.33 37.49 -0.16 0.92
222 34.99 35.36 -0.37 1.68
Fig. 3.29 a, b Three-dimensional MEM electron density isosurface of sample A in two different
orientations. c, d Three-dimensional MEM electron density isosurface of sample BM in two
different orientations. e, f Three-dimensional MEM electron density isosurface of sample BMA in
two different orientations
density in sample BM. This enhancement reflects the high thermal vibration in this
sample. The two-dimensional MEM electron density distributions on (100) plane
shown in Fig. 3.30a–c exhibit the large thermal vibration in sample BM. This is
observed from the smearing of charges at the corner atoms and valence electrons
in the face-centred atoms. The thermal vibration is comparatively less in sample A.
Expansion of the charge region is observed in sample BMA. Expansion of charge
clouds of the face-centred atom and their extension towards the corner atoms is
observed. This reveals the restriction in the thermal vibration resulting in the low
Debye–Waller factor in this case. The increase in the atomic sizes can also be
visualized. Ball milling appears to have increased the average cell constant, which
is then reduced by subsequent annealing of the ball-milled sample.
The one-dimensional electron density profiles shown in Fig. 3.31a along [100]
direction illustrates that the electron density is higher for sample BM initially but
at the mid-bond region the density for sample A is increased. This can be attributed
to the removal of dislocations and vacancies due to annealing. Along the [110]
direction (Fig. 3.32b), which is the bonding direction, the electron density of
sample BM is observed to be higher. Ball milling has effectively increased the
bond strength in Ni80Cr20 along the bonding direction than in the other directions.
The flattening of the particles as seen from SEM picture (Fig. 3.26b) is in support
of this conclusion. Again, the electron density along [111] direction (Fig. 3.32c),
3.5 Nickel Chromium (Ni80Cr20) 127
Fig. 3.30 a Two-dimensional electron density distributions on the (100) plane for sample A.
Contour lines are drawn between 0.0 and 1.5 e/Å3 with 0.1 e/Å3 interval. b Two-dimensional
electron density distribution on the (100) plane for sample BM. Contour lines are drawn between
0.0 and 1.5 e/Å3 with 0.1 e/Å3 interval. c Two-dimensional electron density distribution on the
(100) plane for sample BM. Contour lines are drawn between 0.0 and 1.5 e/Å3 with 0.1 e/Å3
interval
show that the density is higher for sample A. In general, sample BMA has lower
electron density. Ball milling and subsequent annealing has reduced the charge
density in the alloy Ni80Cr20. The strain produced in the sample BM (as evidenced
from the flat region in electron density in Fig. 3.32b) is relieved along the other
directions, particularly the [111] direction as seen from Fig. 3.32c. This figure
shows undulations in the ball-milled sample which is much higher than the other
two samples.
128 3 Results and Discussion on Metals and Alloys
Fig. 3.31 a Two-dimensional electron density distribution on (110) plane for sample A. Contour
lines are drawn between 0.05 and 1.0 e/Å3 with 0.05 e/Å3 interval. b Two-dimensional electron
density distribution on (110) plane for sample BM. Contour lines are drawn between 0.05 and
1.0 e/Å3 with 0.05 e/Å3 interval. c Two-dimensional electron density distribution on (110) plane
for sample BMA. Contour lines are drawn between 0.05 and 1.0 e/Å3 with 0.05 e/Å3 interval
In this study, the observed PDF was obtained from raw X-ray powder data using
software programme PDFgetX (Jeong et al. 2001). A comparison between
observed and calculated PDFs has been carried out using software package
PDFFIT (Farrow et al. 2007). The results of the PDF studies are presented in
Table 3.25.
The refined pair distribution function given in Fig. 3.33a–c for samples A, BM,
BMA respectively, show good matching of observed and calculated PDFs. The
PDF refinement can be considered as equivalent to matching the observed X-ray
powder data with too little structure parameters. Moreover, high Q data are required
for these types of analyses on short-range order and local structure of materials.
Powder XRD data have been used in this work only for comparison purposes.
The bond lengths deduced from the PDF analysis have been compared with those
calculated using software programme GRETEP is a Windows interactive pro-
gramme that displays the list of calculated bonds, bond length between two atoms
and the distances and angles around an atom) (Laugier and Bochu). On comparing
the observed bond lengths obtained from the PDFs of the three samples, it is found
that maximum deviation in the observed bond lengths is noted in sample BM.
3.5 Nickel Chromium (Ni80Cr20) 129
1.2 A 1.175
1.0
0.8 0.850
0.6
0.4 0.525
0.2
0.0 0.200
0.5 1.0 1.5 2.0 2.5 3.0 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
Distance(Å)along [100] direction Distance(Å) along [110] direction
(c) 1.500
Electron density (e/Å3)
A
1.175 BM
BMA
0.850
0.525
0.200
0.5 1.5 2.5 3.5 4.5 5.5
Distance along [111] direction(Å)
Fig. 3.32 a One-dimensional MEM electron density profiles of samples A, BM,BMA along
[100] direction. b One-dimensional MEM electron density profiles of samples A, BM, BMA
along [110] direction. c One-dimensional MEM electron density profiles of samples A, BM,
BMA along [111] direction
Table 3.24 Electron density profiles along the three directions of the unit cell
Sample [100] [110] (Bonding direction) [111]
Position (Å) ED (e/Å3) Position (Å) ED (e/Å3) Position (Å) ED (e/Å3)
A 1.7731 0.484 1.2537 0.378 3.0711 0.484
BM 1.7767 0.460 1.2561 0.605 3.0774 0.460
BMA 1.7691 0.337 1.2509 0.479 3.0641 0.337
ED Electron density
(a) (b)
3 A cal 3 BM cal
obs obs
2 2
1 1
-0 -0
G(r)(Å)-2
G(r)(Å)-2
-1
-1
-2
-2
-3
-3
-4 Difference
-4
-5 Difference
-5
-6
-6
4.0 10.5 17.0 23.5 30.0 4.0 10.5 17.0 23.5 30.0
r(Å) r(Å)
(c)
2 BMA cal
obs
1
-0
-1
G(r) (Å-2)
-2
-3
-4 Difference
-5
-6
Fig. 3.33 a The observed and calculated pair distribution function of sample A. b The observed
and calculated pair distribution function of sample BM. c The observed and calculated pair
distribution function of sample BMA
3.5.5 Conclusion
The effects of ball milling and annealing on the local structure have been studied in
this work. The ball-milled sample shows enhanced electron density along bonding
direction, which decreases on annealing. The cell parameters and the bond lengths
appear to be larger when the sample is ball-milled, but after annealing these
parameter values are reduced due to relieving of strain in the sample.
3.6.1 Introduction
doping. Many studies have been conducted on alkali halides, for example,
Linkoaho (1969), Banerjee and Mitra (1969), Shahi and Wagner (1981), Valvoda
(1970), Burley (1967). The effect of doping on the physical properties of alkali
halides is a subject of sustained interest. While the changes in electrical and optical
properties accompanying doping of alkali halides crystals have received consid-
erable attention, lattice strains and lattice deformations arising out of the intro-
duction of foreign atoms have not been studied extensively. Metals containing salts
are among the inorganic substances that act as antimicrobial agents. Also, the
metal salts are used in gradient index of optical fibre communication, development
of passive fault current limiter in parallel biassing, mass spectroscopic analysis of
metals etc. The present analysis is carried out on the accurate electron density
distribution using MEM and multipole analysis. This research work is on single
crystal data sets of Na0.97Ag0.03Cl and Na0.90Ag0.10Cl (Mohanlal and Pathinettam
Padiyan 1989). When the system is dilute doped, it is possible to study only the
average electronic structure. But, when suitable tools like MEM and multipole
techniques are available, they can be used effectively in conjunction to elucidate
the fine features of bonding as done in this work. In order to elucidate the valence
electrons distribution, MEM and to analyse the contraction/expansion of atomic
shells, multipole analysis of the electron density was carried out.
The alkali halide Na1-xAgxCl, with two different compositions (x = 0.03 and
0.10), was studied with regard to the Ag impurities in terms of the bonding and
electron density distribution. X-ray single crystal data sets have been used for the
purpose. The present analysis focused on the electron density distribution and
hence the interaction between the atoms is clearly revealed by MEM and multipole
analyses. The bonding in these systems has been studied using two-dimensional
MEM electron density maps on the (100) and (110) planes and one-dimensional
electron density profiles along the [100], [110] and [111] directions. The mid-bond
electron densities between atoms in these systems are found to be 0.175 and
0.183 e/Å3, respectively, for Na0.97Ag0.03Cl and Na0.90Ag0.10Cl. Multipole anal-
ysis of the structure has been performed for these two systems, with respect to the
expansion/contraction of the ion involved.
The two data sets (Na0.97Ag0.03Cl and Na0.90Ag0.10Cl) reported by Mohanlal and
Pathinettam Padiyan (1989) have been analysed using Rietveld (1969) structural
refinement using the software package Jana 2000 (Petřiček et al. 2000). The
refined multipole structure factors were used for the MEM computations of the
132 3 Results and Discussion on Metals and Alloys
electron densities. Table 3.26 represents the results of the Rietveld refinements for
Na0.97Ag0.03Cl and Na0.90Ag0.10Cl respectively.
Fig. 3.35 a MEM electron density distribution of Na0.97Ag0.03Cl on the (100) plane. b MEM
electron density distribution of Na0.90Ag0.10Cl on the (100) plane
134 3 Results and Discussion on Metals and Alloys
Fig. 3.36 a MEM electron density distribution of Na0.97Ag0.03Cl on the (110) plane. b MEM
electron density distribution of Na0.90Ag0.10Cl on the (110) plane
0.6 0.6
(a) (b)
110 111
0.5 0.5 110
100
100
Electron density (e/Å3)
Electron density (e/Å3)
111
0.4 0.4
0.3 0.3
0.2 0.2
Na Cl
0.1 Na Cl 0.1
1.2 1.6 2.0 2.4 2.8 3.2 1.3 1.8 2.3 2.8 3.3
Distance along direction (Å) Distance along direction (Å)
Fig. 3.37 a One-dimensional variation of electron density along [100], [110] and [111]
directions of the Na0.97Ag0.03Cl unit cell. b One-Dimensional variation of electron density along
[100], [110] and [111] directions of the Na0.90Ag0.10Cl unit cell
Na0.90Ag0.10Cl
Electron density (e/Å3)
80 Na0.90Ag0.10Cl
70 0.23
60 0.22
50
0.21
40
30 0.20
20 0.19
10
0.18
0
0.17
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8
Distance along [100] direction (Å) Distance along [100] direction (Å)
Fig. 3.38 a The peak electron density profiles of Na (Ag) along the [100] direction for the two
compositions, (Na 0.97Ag0.03Cl and Na0.90Ag0.10Cl). b The mid-bond electron density of Na (Ag)
along the [100] direction for Na0.97Ag0.03Cl and Na0.90Ag0.10Cl
3.6 Silver Doped in NaCl (Na1-xAgxCl) 135
Fig. 3.39 a Dynamic multipole deformation (DMD) map of Na0.97Ag0.03Cl on the (100) plane.
b Dynamic multipole deformation (DMD) map of Na0.90Ag0.10Cl on the (100) plane
Hansen and Coppens (1978) proposed a modified electron density model with the
option that allows the refinement of population parameters at various orbital levels
where the atomic density is described as a series expansion in real spherical
harmonic functions through fourth order Ylm.
Multipole electron density computation was carried out according to Hansen
and Coppens (1978) in this work. The electron density from multipole analysis is
represented in Fig. 3.39a, b as DMD maps on the (100) plane for Na0.97Ag0.03Cl
and Na0.90Ag0.10Cl respectively.
Figure 3.40a, b represent the SMD maps on (100) plane for Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl respectively.
136 3 Results and Discussion on Metals and Alloys
Fig. 3.40 a Static multipole deformation map of Na0.97Ag0.03Cl on the (100) plane. b Static
multipole deformation map of Na0.90Ag0.10Cl on the (100) plane
Fig. 3.41 a Dynamic multipole deformation map of Na0.97Ag0.03Cl on the (110) plane.
b Dynamic multipole deformation map of Na0.90Ag0.10Cl on the (110) plane
Fig. 3.42 a Static multipole deformation map of Na0.97Ag0.03Cl on the (110) plane. b Static
multipole deformation map of Na0.90Ag0.10Cl on the (110) plane
Figure 3.41a, b represent DMD maps on (110) plane for Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl respectively.
Similarly Fig. 3.42a, b represent the SMD maps on (110) plane for
Na0.97Ag0.03Cl and Na0.90Ag0.10Cl respectively. The results of the multipole
analysis are given in Table 3.29.
The effect of impurity doping has been analysed in this work, in terms of the
electron density distribution using the two versatile techniques MEM and multi-
pole analysis. The experimental structure factor data sets of Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl were subjected to Rietveld (1969) refinements to obtain the indi-
vidual thermal vibration parameters of the atoms Na(Ag) and Cl. The refined
thermal vibration parameters as given in Table 3.26 shows that the Debye–Waller
factor of Na(Ag) increases with the concentration of Ag, whereas the Debye–
Waller factor value of Cl decreases. When the NaCl is doped with Ag impurity, the
local distortion increases with impurity concentration. Hence, the Debye–Waller
factor of Na(Ag) increases with Ag. Since the electrostatic interaction between the
Na(Ag) and Cl ions increases with Ag concentration, the Cl atom has less thermal
vibration as reflected in Table 3.26. The Rietveld refinements give very low
reliability indices as given in Table 3.26. The refined results were used for the
MEM computations. The three-Dimensional electron density of Na0.97Ag0.03Cl
and Na0.90Ag0.10Cl in terms of isosurfaces as shown in Fig. 3.34a, b can be
138 3 Results and Discussion on Metals and Alloys
visually compared for the ionic sizes and the interionic distances. It is seen from
these figures that the electronic charge clouds are more diffuse in Na0.90Ag0.10Cl
and the separation between charges decreases as the concentration of Ag increases.
Since the charge clouds of Cl expand due to more electrostatic attraction in
Na0.90Ag0.10Cl, the thermal vibration of Cl is restricted, resulting in smaller
Debye–Waller factor of Cl.
The increase in the ionic sizes can be visualized in the two-Dimentional electron
density contour maps on (100) and (110) plane (Figs. 3.35a, b and 3.36a, b) for
Na0.97Ag0.03Cl and Na0.90Ag0.10Cl respectively. Particularly, Fig. 3.36b for
Na0.90Ag0.10Cl shows increased ionic sizes compared to those in Fig. 3.36a. Hence,
the effect of impurity addition in the host lattice is substantiated.
From the one-dimensional analysis of electron density it is found that there is an
increase in the mid-bond density in Na0.90Ag0.10Cl compared to that in
Na0.97Ag0.03Cl in all the three directions. The minimum density along [100]
direction occurs at 1.4251 and 1.4227 Å respectively, for Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl indicating an increase in the ionic radius in Na0.90Ag0.10Cl.
Moreover, the electron density at these positions is 0.175 and 0.183 e/Å3 for
Na0.97Ag0.03Cl and Na0.90Ag0.10Cl respectively. This clearly indicates the addition
of Ag ion in the host lattice. A similar trend i.e., increased mid-bond position and
density is observed along [110] direction also, in Na0.90Ag0.10Cl. Along [111]
direction, there is an increase in electron density in Na0.90Ag0.10Cl compared to
Na0.97Ag0.03Cl, but the position of mid-bond density is decreased in
Na0.90Ag0.10Cl. But one can observe horizontal lines along [111] direction and
only the minimum numerical value decides the electron density, but the bonding
characteristics is decided by the electron density distribution along [100] direction
which supports the addition of Ag impurities in the host lattice. The [100] profiles
provide an estimation of the ionic radii of atoms Na(Ag) and Cl.
The DMD maps on (100) plane (Fig. 3.39a, b) for Na0.97Ag0.03Cl and
Na0.90Ag0.10Cl clearly indicate that the residual densities are more concentrated
and contracted in Na0.90Ag0.10Cl than Na0.97Ag0.03Cl. The SMD maps in
Fig. 3.40a, b show positive clusters (contour line) between the atom indicative of
more electrostatic interaction in Na0.90Ag0.10Cl compared to Na0.97Ag0.03Cl.
A similar trend is observed on (110) planes (Fig. 3.41a, b). Figure 3.42b shows
(for Na0.90Ag0.10Cl) charge clouds of larger radius between the atoms compared to
the electron clouds of Na0.90Ag0.10Cl. The multipole analysis gives the multipole
parameters j and j0 which are as follows. The j and j0 values of Na(Ag) in
Na0.97Ag0.03Cl and Na0.90Ag0.10Cl was refined keeping the Cl ion fixed. For
Na0.90Ag0.10Cl the j and j0 parameters are 1.805(0.605) and 1.0703(0.5411)
respectively. These values are 1.776(0.5040) and 1.0009(0.0544) respectively for
Na0.90Ag0.10Cl. Both values for Na0.90Ag0.10Cl reveal an expansion of Na(Ag) ion
in Na0.90Ag0.10Cl. Table 3.29 shows that the residual density is larger in
Na0.90Ag0.10Cl as analysed by DMD method.
3.6 Silver Doped in NaCl (Na1-xAgxCl) 139
3.6.7 Conclusion
Both MEM and multipole analyses give a clear picture of the bonding in the doped
system NaCl. Hence, these two methods can be used in combination to study the
electronic charge distribution in pure and doped systems which will help to
understand the physical and conducting phenomena of these materials.
3.7.1 Introduction
The analysis of bonding and electronic structure of metals is useful because of the
variety of applications of metals, particularly when they are doped with other
elements. There has been tremendous impetus in the study of electronic structure
of solids, both theoretically and experimentally due to necessity. When the system
is dilute doped, it is possible to study only the average electronic structure. But,
when suitable tools like MEM are available, they can be effectively used to elu-
cidate the fine features of bonding as attempted in this work. Fourier synthesis of
electron densities can be of use in picturing bonding between two atoms, but, it
suffers from the major disadvantage of series termination error and negative
electron densities that prevent the clear understanding of the fine nature of bonding
between atoms, the factor intended to be analysed. The advent of MEM solves
many of these problems. MEM electron densities are always positive and even
with limited data, one can determine reliable electron densities resembling true
densities. In this work, such an analysis has also been carried out using the for-
malism of Collins (1982) using the data sets reported by Mohanlal and Pathinettam
Padiyan (1987).
In our series of research studies on bonding in materials, we have analysed LiF,
NaF, GaAs, CdTe and some oxides (Israel et al. 2002; Israel et al. 2003; Saravanan
et al. 2002a, b; Kajitani et al. 2001) using MEM and reported X-ray data and
obtained precise information about bonding in the above materials. The nature of
bonding in these materials as analysed are found to be ionic, mixed covalent and
ionic and ‘oxide’ bonding.
The electronic structure of pure and doped aluminium with dilute amounts of iron
impurities (0.215 and 0.304 wt% Fe) has been analysed using reported X-ray data
sets and MEM. Qualitative as well as quantitative assessment of the electron
140 3 Results and Discussion on Metals and Alloys
density distribution in these samples has been made. The mid-bond character-
ization leads to a conclusion about the nature of doping of impurities. An expan-
sion of size of the host aluminium atom is observed with Fe impurities.
In this study on aluminium with iron impurities, we have attempted to study
metal bonding and bring out the details of doping using MEM. The X-ray structure
factor data sets have been taken from Mohanlal and Pathinettam Padiyan (1987)
who reported precise structure factors collected from a home-built manual X-ray
diffractometer. Data sets for pure aluminium, Al ? 0.189, Al ? 0.215 and
Al ? 0.304 wt% Fe are available from this report Mohanlal and Pathinettam
Padiyan (1987). Other relevant details can be found from this paper. Among these
data sets, the one corresponding to Al ? 0.189 wt% Fe has been excluded from the
present analysis after examining the preliminary results from MEM and multipole
methods (the low doping level does not allow any meaningful conclusions of the
obtained results). The remaining three data sets, i.e., those of pure aluminium,
Al ? 0.215 and Al ? 0.304 wt% Fe were used in this study. It was observed that
the Beff and the expected number of vacant sites in Al ? 0.215 wt% Fe are much
better than those in Al ? 0.189 wt% Fe, although slight variations of these
quantities is observed, compared to pure aluminium and Al ? 0.304 wt% Fe [it
has been reported (Mohanlal and Pathinettam Padiyan (1987) that the number of
vacant sites created in the Al lattice increases with the number of Fe impurities and
in the case of Al ? 0.189 wt% Fe, less vacancies are created than expected].
Hence, a data set for Al ? 0.215 wt% Fe has been included in the present analysis
in order to have a composition of iron between pure aluminium and Al ? 0.304
wt% Fe and to compare the results.
All the reported data sets have been analysed using JANA 2000 (Petřiček et al.
2000) software package. The refined structure factors were used for MEM
refinements of the electron densities. The characteristics of the structural refine-
ments are given in Table 3.30.
The reliability indices (R-factors) and the matching of thermal parameters (B)
with the corresponding reported values for all three crystal structures indicate that
the quality of the experimental data is suitable enough for the electron density
analysis. The same number of data for all three compositions samples has been
utilised for the present analysis in order to make a comparison of the results.
MEM calculations were done by dividing the cell into 1283 pixels. For
numerical MEM computations, the software package PRIMA (Izumi and Dilanian
2002) was used. For two- and three-dimensional representations of the electron
densities, the programme VESTA (Momma and Izumi 2006) package was used. In
all three cases, the convergence condition C = 1.0, was obtained after several
cycles of MEM refinement. The main characteristics of the refinements are listed
in Table 3.31.
3.7 Aluminium Doped with Dilute Amounts of Iron Impurities 141
Table 3.32 Maxima of the MEM electron density revealed in different one-dimensional sections
of a unit cell for different compositions
Composition Direction
[100] [110] [111]
Position (Å) ED (e/Å ) Position (Å)) ED (e/Å ) Position (Å) ED (e/Å3)
3 3
The reliability indices are found to be very low. To analyse the MEM results,
one-dimensional profiles along different directions of the unit cell and two-
dimensional electron density contour maps have been constructed. The results of
the one-dimensional analysis have given in Table 3.32.
142 3 Results and Discussion on Metals and Alloys
Fig. 3.43 a 3-D representation of the electron density of pure Al super imposed with the
structural model. Isosurface level is 0.2. b 3-D representation of the electron density of
Al ? 0.215% Fe super imposed with the structural model. Isosurface level is 0.2. c 3-D
representation of the electron density of Al ? 0.304% Fe imposed with the structural super
imposed with the structural model. Isosurface level is 0.2
Fig. 3.44 a MEM electron density distribution of pure Al on the (100) plane. Contour range is
from 0.0 to 5.0 e/Å3. Contour interval is 0.03 e/Å3. b MEM electron density distribution of
Al ? 0.215% on the (100) plane. Contour range is from 0.0 to 5.0 e/Å3. Contour interval is
0.03 e/Å3. c MEM electron density distribution of pure Al ? 0.304% on the (100) plane. Contour
range is from 0.0 to 5.0 e/Å3. Contour interval is 0.03 e/Å3
Fig. 3.45 a MEM electron density distribution of pure Al on the (110) plane. Contour range is
from 0.0 to 5.0 e/Å3. Contour interval is 0.03 e/Å3. b MEM electron density distribution of pure
Al ? 0.215% Fe on the (110) plane. Contour range is from 0.0 to 5.0 e/Å3. Contour interval is
0.03 e/Å3. c MEM electron density distribution of pure Al ? 0.304% Fe on the (110) plane.
Contour range is from 0.0 to 5.0 e/Å3. Contour interval is 0.03 e/Å3
3.7 Aluminium Doped with Dilute Amounts of Iron Impurities 143
(a) (b)
140 0.25 140
0.5
Al Al+0.215 Fe Al
Al+0.215 Fe 0.20 Al+0.215 Fe 0.4
120 Al 120
Al+0.304 Fe Al+0.304 Fe 0.3 Al
Electron Density (e/Å 3)
3
100 100 0.2 Al+0.215 Fe
Al+0.304 Fe Al+0.304 Fe
1.0 1.3 1.6 1.9 2.2 2.5 2.8 0.7500 1.0875 1.4250 1.7625
80 80
60 60
40 40
20 20
0.00 0.41 0.82 1.23 1.64 2.05 2.46 2.87 3.28 3.69 0.0 0.6 1.2 1.8 2.4 3.0 3.6 4.2 4.8 5.4
Distance along [100] direction (Å) Distance along [110] direction (Å)
100 0.10
Al+0.304 Fe
0.08
0 1 2 3 4 5 6 7
80
60
40
20
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
Distance along [111] direction (Å)
Fig. 3.46 a One-dimensional MEM low electron density profiles of pure Al, Al ? 0.215% Fe
and Al ? 0.304% Fe along the [100] direction. The mid-bond variation of the electron density is
shown as insert in which the y-axis represents the electron density in e/Å3 and the x-axis, the
distance along the [100] direction in Å. b One-dimensional MEM low electron density profiles of
pure Al, Al ? 0.215% Fe and Al ? 0.304% Fe along the [110] direction. The mid-bond variation
of the electron density is shown as insert, in which the y-axis represents the electron density in e/
Å3 and the x-axis, the distance along the [110] direction in Å. c One-dimensional MEM low
electron density profiles of pure aluminium, Al ? 0.215% Fe and Al ? 0.304% Fe along the
[111] direction. The mid-bond variation of the electron density is shown as insert, in which the
y-axis represents the electron density in e/Å3 and the x-axis, the distance along the [111] direction
in Å
As expected from the substitution Al by Fe, the peak density decreases with the
increase of Fe impurity (Fig. 3.46). As more and more Fe atoms occupy the host
aluminium lattice substitutionally, the size of the charge cloud increases due to the
larger atomic radius of iron (1.72 Å—for aluminium, the atomic radius is 1.43 Å).
The electron density between atoms along [110] direction also decrease with
impurities (Table 3.32 and Fig. 3.45b) although it is almost a constant along the
other two directions, [100] and [111] (Table 3.32). The bonding direction shows
(which is [110] direction, because the first nearest neighbours are along this
direction), that the mid-bond electron density also decreases with the addition of
Fe impurities as shown in the inset of Fig. 3.46b. The mid-bond density along the
other two directions, [100] and [111] (insets of Fig. 3.46a, c), do not reveal any
trend with the Fe impurities, because the values of the densities involved are very
low along these non-bonding directions and the humps seen in the insets of
Fig. 3.46a, c are due to the interference of the electron density of the FCC atoms in
the unit cell (NNM—Non Nuclear Maxima, along the non-bonding directions).
The MEM electron density distribution of pure Al, Al ? 0.215% Fe and
Al ? 0.304% Fe on the (100) plane (Fig. 3.44a–c) shows clear densities with very
low noise levels. The large thermal vibration and smearing of charges is seen.
Hence the increased size of the electronic clouds can be attributed to the static and
dynamic thermal vibration parameters and also to the increased host lattice sub-
stitution of the Fe atoms.
3.7.5 Conclusion
The electron density distribution in the unit cell of Al with different iron impurity
has been studied precisely using the currently available most versatile technique,
the MEM analysis of the electron densities. Minute changes in the density details
and the incorporation of impurities in the host lattice have been revealed. The
electron density and the mid-bond electron density decrease with more Fe impurity
doping as evidenced from the 1D electron densities. Hence, the bonding strength in
Al decreases with Fe addition. In addition to the electron density details showing
up themselves as actual distributions of atomic electrons, the thermal parameters
also show a similar trend with respect to impurity addition. Moreover, an
expansion of the host lattice with the addition of more Fe impurities is observed
pictorially as revealed by the two-dimensional electron densities in Fig. 3.44.
Since, the crystallographic system presently studied is highly symmetric, there are
no appreciable changes in the mid-bond positions despite the increase in the
atomic sizes with more addition of Fe impurities.
References 145
References
Ak NF, Tekmen C, Ozedmir I, Soykan HS, Celik E (2003) Surf Coat Technol 174:1070–1073
Al-Quran FMF, Al-Itawi HI (2010) Eur J Sci Res 39:251–256
Banerjee NK, Mitra GB (1969) Indian J Pure Appl Phys 71:13
Bashir J, Butt NM, Khan QH (1988) Acta Cryst A44:638–639
Bowran DT, Finney JL (2003) J Chem Phys 118:8357–8372
Burley G (1967) Acta Cryst 1(23):1–5
Collins DM (1982) Nature 298:49–51
Coppens P, Volkov A (2004) Acta Cryst A60:357–364
Debay P, Menki H (1930) Physika Zeit 31:797–798
Egami T (1990) JIM 31(3):163–176
Enríquez JP, Mathew X (2003) J Cryst Growth 259:215–222
Eskandarany MS (2001) Mechanical alloying for fabrication of advanced engineering materials
(Chapter 1). Noyes publications, New York
Farrow CL, Juhas P, Liu JW, Bryndin D, Bozin ES, Bloch J, Proffen T, Billinge SJL (2007)
PDFfit2 and PDF gui: computer programs for studying nanostructures in crystals. J Phy Cond
Matt Phy 19:335219
Field DW, Medlin EH (1974) An X-ray sturdy of thermal vibrations in sodium. Acta Cryst
A30:234–238
Frils J, Jiang B, Spence J, Marthinsen K, Holmestad R (2004) Acta Cryst A60:402–408
Gull SF, Daniel GJ (1978) Nature 272:686–690
Hansen NK, Coppens P (1978) Acta Crys A34:909–921
Ikeda T, Takata M, Sakata M, Waliszewski J, Dobrzynski L, Porowski S, Jun J (1998) J Phys Soc
Jpn 67:4104–4109
Israel S, Saravanan R, Srinivasan N, Rajaram RK (2002) Phys Chem Solids 64:43–49
Israel S, Saravanan R, Srinivasan N, Rajaram RK (2003) J Phys Chem Solids 64:879–886
Israel S, Saravanan R, Rajaram RK (2004) Physica B 349:390–400
Izumi F, Dilanian RA (2002) Recent research developments in physics, vol 3, Part II. Transworld
Research Network, Trivandrum , pp 699–726
Jeng SP, Holloway H, Douglas A, Hoflund GB (1990) Surf Sci 235:175–185
Jeong IK, Thompson J, Proffen T, Perez A, Billinge SJL (2001) PDFGetX: a program for
obtaining the atomic pair distribution function from X-ray powder diffraction data. J Appl
Cryst 34:536
Kajitani T, Saravanan R, Ono Y, Ohno K, Isshiki M (2001) J Cry Growth 229:130–136
Kamal S, Jayaganthan R, Prakash S (2009) J Alloy Compd 472:378–389
Kazi IH, Wild PM, Moore TN, Sayer M (2003) Thin Solid Films 433:337–343
Kazi IH, Wild PM, Moore TN, Sayer M (2006) Thin Solid Films 515:2602–2606
Kim N, Kwon NH, Park Y, Choi J (2005) Electron Lett 41:982–984
Larson AC, Von Dreele RB (2004) General structure analysis system, Lansce, MS-H805, Los
Alamos national laboratory, Los Alamos, LAUR 86748
Linkoaho MV (1969) Acta Crystallographica A 25 III:450–455
Linkoaho MV (1972) Physica Scripta 5:271–272
Marabello D, Bianchi R, Gervasio G, Cargnoni F (2004) Acta Cryst A60:494–501
Mohanlal SK (1979) J Phys C Solid State Phys 12:651–653
Mohanlal SK, Pathinettam Padiyan D (1987) Cryst Latt Def Amorph Mat 14:651–653
Mohanlal SK, Pathinettam Padiyan D (1989) Radiat Eff Defects Solids 110:385–391
Momma K, Izumi F (2006) Commission on crystallography computing. IUCr Newslett 7:106–119
Murray GT (1997) Handbook of materials selection for engineering applications. Dekker, Inc.
New York
Na SH, Kim SH, Lee YW, Sohn DS (2002) J Korean Nucl Soc 34:60–67
Peng LM, Ren G, Dudarey SL, Whelan M (1996) J Acta Cryst A52:456–470
146 3 Results and Discussion on Metals and Alloys
The average and local structures of some elemental metals, alloys and dilute doped
alloys have been analysed from X-ray data using versatile techniques like Rietveld
method, MEM (maximum entropy method), multipole and PDF (pair distribution
function). The materials studied in this research work are listed below.
Metals
• Sodium (Na)
• Vanadium (V)
• Magnesium (Mg)
• Aluminium (Al)
• Titanium (Ti)
• Iron (Fe)
• Nickel (Ni)
• Copper (Cu)
• Zinc (Zn)
• Tin (Sn)
• Tellurium (Te)
Alloys
• Cobalt aluminium (CoAl)
• Nickel aluminium (NiAl)
• Nickel chromium (NiCr)
• Iron nickel (FeNi)
Dilute doped alloys
• Sodium chloride with iron impurities (Na1-xAgxCl)
• Aluminium, with iron impurities (0.215 and 0.304 wt% Fe)
The materials studied in this research work have been divided into several parts
for the exhaustive analytical study and the results of this research work are
presented below.
R. Saravanan and M. Prema Rani, Metal and Alloy Bonding: An Experimental Analysis, 147
DOI: 10.1007/978-1-4471-2204-3_4, Springer-Verlag London Limited 2012
148 4 Conclusion
• The average and local structure of simple metals Al, Ni and Cu have been
elucidated for the first time using MEM (Maximum Entropy Method), multipole
and PDF (Pair Distribution Function).
• The powder X-ray data sets were refined by Rietveld method using the software
program JANA 2000.
• The refined structure factors were used for the MEM and multipole analysis.
• The electron density distribution of Al, Ni and Cu has been mapped using the
MEM electron density values obtained through refinements.
• The mid-bond density in Ni is found to be larger among the three metal systems
along [110] direction. But along the [100] and [111] directions (non-bonding),
the element with the higher atomic number has a higher electron density.
• The thermal contribution to the charge density has been realized through the
dynamic and static multipole deformation maps.
• From multipole results there is appreciable contraction in Ni compared to the
other two systems.
4.2 Aluminium, Nickel and Copper 149
• The observed PDF of Al, Ni and Cu were obtained from the raw intensity data
using a software program PDFgetX.
• A comparison between the observed and calculated PDF has been carried out
using the software package PDFFIT.
• The PDF results give a clear picture of the local structure in these metals.
• The average and local structures of magnesium, titanium, iron, zinc, tin and
tellurium have been analysed using the MEM (maximum entropy method), and
PDF (pair distribution function).
• The structural parameters of the metals were refined from the powder X-ray data
with the well-known Rietveld powder profile fitting methodology using the
software package JANA 2006.
• The refined structure factors were used for the MEM analysis
• The one-, two- and three-dimensional electron density distributions of Mg, Ti,
Fe, Zn, Sn and Te have been mapped using the MEM electron density values
obtained through refinements.
• The variation of electron density distribution with respect to the thermal
vibration factor is clearly visualised.
• The electron density has been analysed qualitatively and quantitatively.
• The local structure has been analysed for the accurate nearest neighbour dis-
tances using PDF (pair distribution function).
• The precise electron density distribution and bonding in metal alloys CoAl and
NiAl are characterized using MEM (maximum entropy method) and multipole
method from reported single-crystal data sets.
• The mid-bond electron densities in these systems are found to be 0.358 and
0.251 e/Å3 respectively, for CoAl and NiAl in the MEM analysis.
• The two-dimensional maps and one-dimensional electron density profiles have
been constructed and analysed.
• The thermal vibration of the individual atoms Co, Ni and Al has also been
studied and reported.
• The thermal vibration of Al is more pronounced in CoAl than in NiAl.
• The contraction of atoms in CoAl and expansion of Ni and contraction of Al
atom in NiAl is found from the multipole analysis.
150 4 Conclusion
• The alloy Ni80Cr20 has been annealed and ball milled to study the effect of
thermal and mechanical treatments on the local structure and the electron
density distribution.
• The electron density between the atoms was studied by MEM (maximum
entropy method) and the local structure using PDF (pair distribution function).
• The electron density is found to be high for ball-milled sample along the
bonding direction.
• The particle sizes of the differently treated samples were realized by SEM and
through XRD.
• Clear evidence of the effect of ball milling is observed on the local structure and
electron densities.
• The alkali halide Na1-xAgxCl, with two different compositions (x = 0.03 and
0.10), was studied with regard to the Ag impurities in terms of the bonding and
electron density distribution using X-ray single-crystal data sets.
• Electron density distribution and hence the interaction between the atoms is
clearly revealed by maximum entropy method (MEM) and multipole analyses.
• The bonding in these systems has been studied using two-dimensional MEM
electron density maps on the (100) and (110) planes and one-dimensional
electron density profiles along the [100], [110] and [111] directions.
• The mid-bond electron densities between atoms in these systems are found to be
0.175 and 0.183 e/Å3, respectively, for Na0.97Ag0.03Cl and Na0.90Ag0.10Cl.
• Multipole analysis of the structure has been performed for these two systems,
with respect to the expansion/contraction of the ion involved.
• The electronic structure of pure and doped aluminium with dilute amounts of
iron impurities (0.215 and 0.304 wt% Fe) has been analysed using reported
X-ray data sets and the MEM (maximum entropy method) technique.
• Qualitative as well as quantitative assessment of the electron density distribution
in these samples has been made.
• An expansion of size of the host aluminium atom is observed with Fe impurities.
Minute changes in the density details and the incorporation of impurities in the
host lattice have been revealed.
4.7 Aluminium Doped with Dilute Amounts of Iron Impurities 151
• The electron density and the mid-bond electron density decrease with more Fe
impurity doping as evidenced from the one-dimensional electron densities.
• The bonding strength in Al decreases with Fe addition.
The broad applicability of structure determination and refinement with single
crystals and powders make the future of diffraction studies more exciting and
essential. The success of extraction of details mainly depends on three parameters,
choice of measurement device, pattern profile description and structure solving
algorithms. Modern instrumentation, sources yielding data of unprecedented
quality and modern analysis methods continue to increase our ability to harvest
useful information from the data.
In this research work the charge density of metals and alloys analysed using
versatile tools such as maximum entropy method, multipole and pair distribution
function have been done for the first time. These analyses give the fine details of
the electron environment in the unit cell. The changes in the properties of the
alloys due to the inclusion of dopant elements are clearly picturised. These local
and average analyses of the materials facilitate the better understanding of their
properties and enhance the engineering of new materials.