Advanced Nanomaterials
Advanced Nanomaterials
Advanced Nanomaterials
1
Advanced nanomaterials
Cours support
This text is also partially used for the cours
Introduction to Nanomaterial
H.Hofmann
Powder Technology Laboratory
IMX
EPFL
Version 1 Sept 2009
NANOMATERIALS
2
Contents
1 INTRODUCTION 6
1.1 WHAT ARE NANOMATERIALS? 6
1.1.1 CLASSIFICATION OF NANOSTRUCTURED MATERIALS 6
1.1.2 WHY SO MUCH INTEREST IN NANOMATERIALS? 8
1.1.3 INFLUENCE ON PROPERTIES BY "NANO-STRUCTURE INDUCED EFFECTS" 9
1.2 SOME PRESENT AND FUTURE APPLICATIONS OF NANOMATERIALS 10
1.3 WHAT ARE THE FUNDAMENTAL ISSUES IN NANOMATERIALS? 13
2 ATOMS, CLUSTERS AND NANOMATERIALS 15
3 NANOCOMPOSITES SYNTHESIS AND PROCESSING 19
3.1 INTRODUCTION 19
3.2 INORGANIC NANOTUBES 20
3.2.1 INTRODUCTION 20
3.2.2 GENERAL SYNTHETIC STRATEGIES 26
3.3 FUNCTIONAL MATERIALS BASED ON SELF-ASSEMBLY OF POLYMERIC SUPRAMOLECULES 28
3.4 MOLECULAR BIOMIMETICS: NANOTECHNOLOGY THROUGH BIOLOGY 32
3.4.1 SELECTION OF INORGANIC-BINDING PROTEINS THROUGH DISPLAY TECHNOLOGIES 34
3.4.2 CHEMICAL SPECIFICITY OF INORGANIC-BINDING POLYPEPTIDES 36
3.4.3 PHYSICAL SPECIFICITY OF PEPTIDE BINDING 37
3.4.4 PEPTIDE-MEDIATED NANOPARTICLE ASSEMBLY 38
3.5 REFERENCES 40
4 MECHANICAL PROPERTIES 41
4.1 INTRODUCTION 41
4.2 METALS 41
4.2.1 GRAIN SIZE EFFECTS IN PLASTICITY AND CREEP 42
4.2.2 METAL PLASTIC DEFORMATION: A COMPARISON BETWEEN CU AND NI NANOPHASE SAMPLES
49
4.2.3 HARDNESS 53
NANOMATERIALS
3
4.3 CERAMICS / NANOCOMPOSITES 56
4.3.1 DENSITY 56
4.3.2 FRACTURE STRENGTH 56
4.3.3 STRENGTHENING AND TOUGHENING MECHANISMS 58
4.3.4 REDUCTION IN PROCESSING FLAW SIZE 60
4.3.5 CRACK HEALING (ANNEALING TREATMENT) 60
4.3.6 TOUGHENING (K-MECHANISMS) 61
4.3.7 GRAIN BOUNDARY STRENGTHENING MECHANISMS 63
4.3.8 THERMAL EXPANSION MISMATCH (SELSING MODEL) 63
4.3.9 AVERAGE INTERNAL STRESSES 64
4.3.10 LOCAL STRESS DISTRIBUTION 67
4.4 FINAL REMARKS ON STRENGTHENING AND TOUGHENING MECHANISMS 67
4.5 REFERENCES 69
5 THERMAL CONDUCTIVITY IN NANOSTRUCTURED MATERIAL 70
5.1 THERMAL CONDUCTIVITY OF THERMAL BARRIER COATINGS 70
5.1.1 THERMAL CONDUCTIVITY 70
5.1.2 LATTICE WAVES 71
5.1.3 INTERACTION PROCESSES 72
5.2 HIGH-TEMPERATURE THERMAL CONDUCTIVITY OF POROUS AL
2
O
3
NANOSTRUCTURES 76
5.2.1 THEORY 76
5.2.2 EXPERIMENT 82
5.2.3 RESULTS AND DISCUSSION 83
5.3 REFERENCES 88
6 THERMODYNAMIC 89
6.1 NANOTHERMODYNAMICS 89
6.1.1 HILLS THEORY 90
6.1.2 TSALLIS GENERALIZATION OF ORDINARY BOLTZMANN-GIBBS THERMOSTATICS 96
6.1.3 THERMODYNAMICS OF METASTABLE PHASE NUCLEATION ON NANOSCALE 100
6.1.4 NANOTHERMODYNAMIC ANALYSES OF CVD DIAMOND NUCLEATION 112
6.2 THERMODYNAMICS OF MELTING AND FREEZING IN SMALL PARTICLES 117
6.2.1 THEORY VANFLEET AND AL. MODEL 117
6.3 PHASE DIAGRAMS 128
6.3.1 GOVERNING EQUATIONS 128
NANOMATERIALS
4
6.3.2 MATHEMATICAL DESCRIPTION FOR NANO-PHASES OF SN-BI ALLOYS 130
6.3.3 PHASE DIAGRAM FOR ISOLATED NANO-PHASES OF SN-BI ALLOYS 133
6.4 CRYSTAL-LATTICE INHOMOGENEOUS STATE 135
6.5 CONCENTRATIONAL INHOMOGENEITY 137
6.6 REFERENCES 140
7 ELECTRONIC AND OPTICAL PROPERTIES OF NANOMATERIALS 141
7.1 INTRODUCTION 141
7.2 METALS 143
7.2.1 INTRODUCTION 143
7.2.2 ELECTRICAL CONDUCTIVITY 154
7.2.3 SURFACE PLASMONS 162
7.3 CARBON NANOTUBES 171
7.3.1 ELECTRONIC STRUCTURE 172
7.3.2 QUANTUM TRANSPORT PROPERTIES 174
7.3.3 NANOTUBE JUNCTIONS AND DEVICES 178
7.4 SEMICONDUCTOR 181
7.4.1 INTRODUCTION 181
7.4.2 BAND GAP MODIFICATION 185
7.4.3 ELECTRICAL PROPERTIES 190
7.4.4 OPTICAL PROPERTIES 212
7.5 REFERENCES 217
8 MAGNETISM 219
8.1 INTRODUCTION 219
8.1.1 CONCEPT 219
8.1.2 PHENOMENA 220
8.2 MAGNETIC PROPERTIES OF SMALL ATOMIC CLUSTERS 222
8.2.1 INTRODUCTION 222
8.2.2 SIZE DEPENDENCE 223
8.2.3 THERMAL BEHAVIOUR 225
8.2.4 RARE EARTH CLUSTERS 226
8.3 SMALL PARTICLE MAGNETISM 226
8.3.1 CLASSIFICATIONS OF MAGNETIC NANOMATERIAL 226
8.3.2 ANISOTROPY 230
NANOMATERIALS
5
8.3.3 SINGLE DOMAIN PARTICLES 232
8.3.4 SUPERPARAMAGNETISM 239
8.4 MAGNETOELECTRONICS SPINS 252
8.4.1 SPIN-POLARIZED TRANSPORT AND MAGNETORISISTIVE EFFECTS 252
8.4.2 SPIN INJECTION 256
8.4.3 SPIN POLARIZATION 257
8.5 GIANT MAGNETORESISTANCE (GMR) 259
8.6 STORAGE DEVICES 264
8.6.1 MAGNETIC DATA STORAGE : 264
8.6.2 SENSORS: 265
8.7 REFERENCES 266
9 APPLICATIONS 267
NANOMATERIALS
6
1 INTRODUCTION
1.1 What are nanomaterials?
1.1.1 Classification of Nanostructured Materials
Nanomaterials are materials which are characterised by an ultra fine grain size
(< 50 nm) or by a dimensionality limited to 50 nm. Nanomaterials can be created with
various modulation dimensionalities as defined by Richard W. Siegel: zero (atomic
clusters, filaments and cluster assemblies), one (multilayers), two (ultrafine-grained
overlayers or buried layers), and three (nanophase materials consisting of equiaxed
nanometer sized grains) as shown in the above Figure 1-1.
Figure 1-1 : Definition of nanomaterials following Siegel
i
Nanomaterials consisting of nanometer sized crystallites or grains and
interfaces may be classified according to their chemical composition and shape
(dimensionality), as discussed above. According to the shape of the crystallites or
grains we can broadly classify nanomaterials into four categories:
1. clusters or powders (MD=0)
2. Multilayers (MD=1)
3. ultrafine grained overyaers or buried layers (where the layer thickness or
the rod-diameters are <50 nm) (MD=2)
4. nanomaterials composed of equiaxed nanomter-sized grains (MD=3)
NANOMATERIALS
7
Figure 1-2 : Classification schema for nanomaterials according to their chemical composition and the
dimensionality (shape)of the crystallites (structural elements) forming the nanomaterial. The boundary
regions of the first and second family of nanomaterials are indicated in black to emphasize the
different atomic arrangements in the crystallites and in the boundaries. The chemical composition of
the (black) boundary regions and the crystallites is identical in the first family. In the second family, the
(black) boundaries are the regions where two crystals of different chemical composition are joined
together causing a steep concentration gradient.
The latter three categories can be further grouped into four families as shown in
Figure 1-2.
In the most simple case (first family in the Figure 1-2), all grains and interfacial
regions have the same chemical composition. Eg. Semicrystalline polymers
(consisting of stacked lamellae separated by non-crystalline region),
multilayers of thin film crystallites separated by an amorphous layer (a-
Si:N:H/nc-Si)
ii
etc.
As the second case, we classify materials with different chemical composition
of grains. Possibly quantum well structures are the best example of this family.
In the third family includes all materials that have a different chemical
composition of its forming matter (including different interfaces) eg. Ceramic of
alumina with Ga in its interface.
iii
NANOMATERIALS
8
The fourth family includes all nanomaterials formed by nanometer sized grains
(layers, rods or equiaxed crystallites) dispersed in a matrix of different
chemical composition. Precipitation hardened alloys typically belong to this
family. Eg. Nanometer sized Ni
3
Al precipitates dispersed in a nickel matrix-
generated by annealing a supersaturated Ni-Al solid solution- are an example
of such alloys. Most high-temperature materials used in modern jet engines
are based on precipitation-hardened Ni
3
Al/Ni alloys.
A large part of this definition has been described in an article by Gleiter.
iv
,
v
1.1.2 Why so much interest in nanomaterials?
These materials have created a high interest in recent years by virtue of their
unusual mechanical, electrical, optical and magnetic properties. Some examples
are given below:
Nanophase ceramics are of particular interest because they are more ductile
at elevated temperatures as compared to the coarse-grained ceramics.
Nanostructured semiconductors are known to show various non-linear optical
properties. Semiconductor Q-particles also show quantum confinement effects
which may lead to special properties, like the luminescence in silicon powders
and silicon germanium quantum dots as infrared optoelectronic devices.
Nanostructured semiconductors are used as window layers in solar cells.
Nanosized metallic powders have been used for the production of gas tight
materials, dense parts and porous coatings. Cold welding properties combined
with the ductility make them suitable for metal-metal bonding especially in the
electronic industry.
Single nanosized magnetic particles are mono-domains and one expects that
also in magnetic nanophase materials the grains correspond with domains,
while boundaries on the contrary to disordered walls. Very small particles have
special atomic structures with discrete electronic states, which give rise to
special properties in addition to the super-paramagnetism behaviour. Magnetic
nano-composites have been used for mechanical force transfer (ferrofluids),
for high density information storage and magnetic refrigeration.
NANOMATERIALS
9
Nanostructured metal clusters and colloids of mono- or plurimetallic
composition have a special impact in catalytic applications. They may serve as
precursors for new type of heterogeneous catalysts (Cortex-catalysts) and
have been shown to offer substantial advantages concerning activity,
selectivity and lifetime in chemical transformations and electrocatalysis (fuel
cells). Enantioselective catalysis were also achieved using chiral modifiers on
the surface of nanoscale metal particles.
Nanostructured metal-oxide thin films are receiving a growing attention for the
realisation of gas sensors (NO
x
, CO, CO
2
, CH
4
and aromatic hydrocarbons)
with enhanced sensitivity and selectivity. Nanostructured metal-oxide (MnO
2
)
find application for rechargeable batteries for cars or consumer goods. Nano-
crystalline silicon films for highly transparent contacts in thin film solar cell and
nano-structured titanium oxide porous films for its high transmission and
significant surface area enhancement leading to strong absorption in dye
sensitized solar cells.
Polymer based composites with a high content of inorganic particles leading to
a high dielectric constant are interesting materials for photonic band gap
structure produced by the LIGA.
1.1.3 Influence on properties by "nano-structure induced effects"
For the synthesis of nanosized particles and for the fabrication of
nanostructured materials, laser or plasma driven gas phase reactions, evaporation-
condensation mechanisms, sol-gel-methods or other wet chemical routes like inverse
micelle preparation of inorganic clusters have been used, that will be discussed later.
Most of these methods result in very fine particles which are more or less
agglomerated. The powders are amorphous, crystalline or show a metastable or an
unexpected phase, the reasons for which is far from being clear. Due to the small
sizes any surface coating of the nano-particles strongly influences the properties of
the particles as a whole. Studies have shown that the crystallisation behaviour of
nano-scaled silicon particles is quite different from micron-sized powders or thin films.
It was observed that tiny polycrystallites are formed in every nano-particle, even at
moderately high temperatures.
Roughly two kinds of "nano-structure induced effects" can be distinguished:
NANOMATERIALS
10
First the size effect, in particular the quantum size effects where the normal
bulk electronic structure is replaced by a series of discrete electronic levels,
and second the surface or interface induced effect, which is important because
of the enormously increased specific surface in particle systems.
While the size effect is mainly considered to describe physical properties, the
surface or interface induced effect, plays an eminent role for chemical processing, in
particular in connection with heterogeneous catalysis. Experimental evidence of the
quantum size effect in small particles has been provided by different methods, while
the surface induced effect could be evidenced by measurement of thermodynamic
properties like vapour pressure, specific heat, thermal conductivity and melting point
of small metallic particles. Both types of size effects have also been clearly separated
in the optical properties of metal cluster composites. Very small semiconductor (<10
nm), or metal particles in glass composites, and semiconductor/polymer composites
show interesting quantum effects and non-linear electrical and optical properties.
The numerous examples, which are not complete, by far, indicate that these
materials will most probably gain rapidly increasing importance in the near future. In
general, properties, production and characterisation methods and their inter-relations
are however not yet satisfactorily understood. Hence, efforts need to be made to
enable the directed tailoring of nanophase, nanoscopic and nanocomposite materials
needed for future technical and industrial applications.
1.2 Some present and future applications of nanomaterials
Here we list some of the present and future applications of nanomaterials that
has been reported in recent literature:
Table 1-1 : Some typical properties of nano-structured materials and possibilities of future applications
Property Application
NANOMATERIALS
11
Bulk
Single magnetic domain
Small mean free path of electrons in a solid
Size smaller than wavelength
High & selective optical absorption of metal particles
Formation of ultra fine pores due to
superfine agglomeration of particles
Uniform mixture of different kinds of superfine particles
Grain size too small for stable dislocation
Magnetic recording
Special conductors
Light or heat absorption, Scattering
Colours, filters, solar absorbers,
photovoltaics, photographic
material, phototropic material
Molecular Filters
R&D of New Materials
High strength and hardness of
metallic materials
Surface/ Interface
Large specific surface area
Catalysis, sensors
Large surface area, small heat capacity
Heat-exchange materials
Combustion Catalysts
Lower sintering temperature
Specific interface area, large boundary area
Superplastic behaviour of ceramics
Cluster coating and metallization
Multi-shell particles
Sintering accelerators
Nano-structured materials
ductile ceramics
Special resistors, temperature sensors
Chemical activity of catalysts
Tailored Optical elements
NANOMATERIALS
12
Table 1-2 : Some examples of present and potential applications with significant technological
impact:
vi
Technology Present Impact Potential Impact
Dispersions and
Coatings
Thermal barriers
Optical (visible and UV) barriers
Imaging enhancement
Ink-jet materials
Coated abrasive slurries
Information-recording layers
Enhanced thermal barriers
Multifunctional nanocoatings
Fine particle structure
Super absorbant materials (Ilford
paper)
Higher efficiency and lower
contamination
Higher density information storage
High Surface
Area
Molecular sieves
Drug delivery
Tailored catalysts
Absorption/desorption materials
Molecule-specific sensors
Particle induced delivery
Energy storage (fuel cells, batteries)
Grtzel-type solar cells, Gas sensors
Consolidated
Materials
Low-loss soft magnetic materials
High hardness, tough WC/Co
cutting tools
Nanocomposite cements
Superplastic forming of ceramics
Materials
Ultrahigh-strength, tough structural
materials
Magnetic refrigerants
Nanofilled polymer composites
Ductile cements
Bio-medical
aspects
Functionalised nanoparticles Cell labelling by fluorescent
nanoparticles
Local heating by magnetic
nanoparticles
Nanodevices
GMR read heads
Terabit memory and microprocessing
Single molecule DNA sizing and
sequencing
Biomedical sensors
Low noise, low threshold lasers
NANOMATERIALS
13
Nanotubes for high brightness displays
1.3 What are the fundamental issues in nanomaterials?
The fundamental issues in this domain of nanomaterials are:
(1) ability to control the scale (size) of the system,
(2) ability to obtain the required composition -
not just the average composition - but details such as defects, concentration
gradients, etc.,
(3) ability to control the modulation dimensionality,
(4) during the assembly of the nano-sized building blocks, one should be able to
control the extent of the interaction between the building blocks as well as the
architecture of the material itself.
More specifically the following issues have to be considered for the future
development of nanomaterials:
Development of synthesis and/or fabrication methods for raw materials
(powders) as well as for the nanostructured materials.
Better understanding of the influence of the size of building blocks in nano
structured materials as well as the influence of microstructure on the physical,
chemical and mechanical properties of this material.
Better understanding of the influence of interfaces on the properties of nano-
structured material.
Development of concepts for nanostructured materials and in particular their
elaboration.
Investigation of catalytic applications of mono- and plurimetallic nanomaterials
Transfer of developed technologies into industrial applications including the
development of the industrial scale of synthesis methods of nanomaterials and
nanostructured systems.
In the following chapters we will review the various developments that have been
revolutionising the application of nanomaterials. We will attempt to correlate the
NANOMATERIALS
14
improvements in the material properties that are achieved due to the fine
microstructures arising from the size of the grains and/or dimensionality.
NANOMATERIALS
15
2 ATOMS, CLUSTERS AND NANOMATERIALS
atoms
?
?
?
?
?
?
?
?
?
?
?
?
molecules cluster nanocrystallites
Figure 2-1 : Schematic representation of various states of matter
At the beginning of last century, increasing attention was focused on the physical
chemistry of colloidal suspensions. By referring to them as "the world of neglected
dimensions", Oswald was the first to realize that nanoscale particles should display
novel and interesting properties largely dependant on their size and shape.
vii
However, it is only in the last two decades that significant interest has been devoted
to inorganic particles consisting of a few hundred or a few dozen atoms, called
clusters. This interest has been extended to a large variety of metals and
semiconductors and is due to the special properties exhibited by these nanometer-
sized particles, which differ greatly from those of the corresponding macrocrystalline
material.
Matter that is constituted of atoms and molecules as such, has been widely classified
and satisfactorily explained. However, an ensemble of atoms, or molecules forming
the so-called Clusters are far from being properly understood. Elemental clusters
are held together by various forces depending on the nature of the constituting
atoms:
Inert gas clusters are weakly held together by van-der-waals interactions, eg. (He)
n
Semiconductor clusters are held with strong directional covalent bonds, eg. (Si)
n
Metallic clusters are fairly strongly held together by delocalised non-directional
bonding, eg. (Na)
n
NANOMATERIALS
16
No. Of Shells 1 2 3 4 5
No. of atoms M
13
M
55
M
147
M
309
M
561
Percentage of
atoms
92% 76% 63% 52% 45%
Figure 2-2 : Idealized representation of hexagonal close packed full-shell magic number clusters.
viii
Note that as the number of atoms increases, the percentage of surface atoms decreases.
1
10
100
1000
10000
100000
0 10 20 30
Cluster size (nm)
T
o
t
a
l
n
o
.
o
f
a
t
o
m
s
(1)
0
20
40
60
80
100
120
0 5 10 15 20 25
Size of cluster (nm)
S
u
r
f
a
c
e
a
t
o
m
s
(
%
)
(2)
Figure 2-3 : (1) Total Number of atoms with size of the cluster. (2) Number of surface atoms for a
hypothetical model sphere of diameter 0.5 nm and density 1000 Kgm
-3
with a mass of 6.5 10
-26
Kg
occupying a volume of about 6.5 x 10
-29
m
3
with a geometrical cross-section of 2 x 10
-19
m
2
(in terms of
atomic mass the sphere is considered to have a mass of 40 amu, where 1 amu = 1.67 x 10
-27
Kg).
Calculated by (a) considering dense structures (Square),
ix
and (b) method suggested by Preining (dark
circle).
x
Either elemental clusters or a mixture of clusters of different elements constitute the
vast expanding field of materials sciences called nanomaterials. One has to be clear
right at this stage that clusters are not a fifth state of matter, as sometimes believed,
but they are simply intermediate between atoms on one hand, and solid or liquid
state of matter on the other, with widely varying physical and chemical properties.
Depending on the number of atoms forming the cluster determines the percentage of
atoms that are exposed on the surface of the cluster. An example of such an
ensemble of metal atoms show the decreasing number of surface atoms with
NANOMATERIALS
17
increasing size of the cluster as shown in fig: . When an ensemble of atoms add up to
form a few nanometer sized clusters, they form what we call nanoparticles, since
only a few atoms forming clusters are called molecular clusters. Agglomeration of a
few atoms have been studied in great details by physicists working with molecular
beams. Today, the mystery related to larger ensemble of atoms (in other words
nanomaterial) are getting clearer due to active research being carried over across
the world over the last decade or two.
Table 2-1 : Idealized representation of the variation of cross section, total mass, number of molecules
and the effective surface atoms in clusters. Note (a) considering dense structures (Square),
ix
and (b)
method suggested by Preining (dark circle).
x
Size
Cross
section
Mass
No. Of
molecules
Fraction of molecules at surface
(%)
(nm) (10
-18
m
2
) (10
-25
Kg) a b
0,5 0,2 0,65 1
1 0,8 5,2 8 100 99
2 3,2 42 64 90 80
5 20 650 1000 50 40
10 80 5200 8000 25 20
20 320 42000 64000 12
In the table, we see that the smaller particles contain only a few atoms, practically all
at the surface. As the particle size increases from 1-10 nm, cross-section increases
by a factor of 100 and the mass number of molecules by a factor of 1000. Meanwhile,
the proportion of molecules at the surface falls from 100% to just 25%. For particles
of 20nm size, a little more than 10% of the atoms are on the surface.
Of course this is an idealized hypothetical case. If particles are formed by macro
molecules (that are larger than the present example), number of molecules per
particle will decrease and their surface fraction increase. The electronic properties of
these ensemble of atoms or molecules will be the result of their mutual interactions
so that the overall chemical behaviour of the particles will be entirely different from
the individual atoms or molecules that they are constituted of. They will also be
different from their macroscopic bulk state of the substance in question under the
same conditions of temperature and pressure.
Table 2-2 : Particle size, surface area and surface energy of CaCo
3
.
xi
(the surface energy of bulk
CacO3 (calcite) is 0.23 Jm
-2
.)
NANOMATERIALS
18
Size
(nm)
Surface area
(m
2
mol
-1
)
Surface energy
(J mol
-1
)
1 1.11 x 10
9
2.55 x 10
4
2 5.07 x 10
8
1.17 x 10
4
5 2.21 x 10
8
5.09 x 10
3
10 1.11 x 10
8
2.55 x 10
3
20 5.07 x 10
7
1.17 x 10
3
10
2
1.11 x 10
7
2.55 x 10
2
10
3
(1 m)
1.11 x 10
6
2.55 x 10
The idea of tailoring properly designed atoms into agglomerates has brought in new
fundamental work in the search for novel materials with uncharacteristic properties.
Among various types of nanomaterials, cluster assembled materials represent an
original class of nanostructured solids with specific structures and properties. In
terms of structure they could be classified in between amorphous and crystalline
materials. In fact, in such materials the short-range order is controlled by the grain
size and no long-range order exists due to the random stacking of nanograins
characteristic of cluster assembled materials. In terms of properties, they are
generally controlled by the intrinsic properties of the nanograins themselves and by
the interactions between adjacent grains. Cluster assembled films are formed by the
deposition of these clusters onto a solid substrate and are generally highly porous
with densities as low as about one half of the corresponding bulk materials densities
and both the characteristic nanostructured morphology and a possible memory effect
of the original free cluster structures are at the origin of their specific properties.
xii
From recent developments in the cluster source technologies (thermal, laser
vaporisation and sputtering),
xiii,xiv
it is now possible to produce intense cluster beams
of any materials, even the most refractory or complex systems (bimetallic,
xv
oxides
and so on), for a wide range of size from a few atoms to a few thousands of atoms.
NANOMATERIALS
19
3 NANOCOMPOSITES SYNTHESIS AND PROCESSING
3.1 Introduction
A major goal of material science is to produce hierarchical materials that are
ordered on all length scales, from the molecular (1100 ) via the nano (10100 nm)
to the meso (1100 mm) scale. In these materials, the larger-scale properties can be
controlled by choosing the appropriate molecular characteristics.
Nanocomposites are a special class of materials originating from suitable
combinations of two or more such nanoparticles or nanosized objects in some
suitable technique, resulting in materials having unique physical properties and wide
application potential in diverse areas that can be formed into a useful object which
can be subsequently used. Novel properties of nanocomposites can be derived from
the successful combination of the individual characteristics of parent constituents into
a single material. To exploit the full potential of the technological applications of the
nanomaterials, it is thus extremely important to endow them with good processability.
Nanocomposites are either prepared in a host matrix of inorganic materials (glass,
porous ceramics etc) or by using conventional polymers as one component of the
nanocomposites. The second type of nanocomposites which are a special class of
hybrid materials are termed polymeric nanocomposites. These materials are
intimate combinations (up to almost the molecular level) of one or more inorganic
materials (nanoparticles , eg.) with a polymer so that unique properties of the former
can be mixed with the existing qualities of the polymer to result in a totally new
material suitable for novel applications. In most of the cases such combinations
require blending or mixing of the components, taking the polymer in solution or in
melt form.
xvi
Resulting nanocomposites have found successful applications in
versatile fields viz. battery cathodes,
xvii
microelectronics,
xviii
nonlinear optics,
xix
sensors,
xx
etc.
NANOMATERIALS
20
3.2 Inorganic Nanotubes
xxi
3.2.1 Introduction
In 1991, Iijima observed some unusual structures of carbon under the transmission
electron microscope wherein the graphene sheets had rolled and folded onto
themselves to form hollow structures. Iijima called them nanotubes of carbon which
consisted of several concentric cylinders of grapheme sheets. Graphene sheets are
hexagonal networks of carbon and these layers get stacked one above the other in
the c-direction to form bulk graphite. Following the initial discovery, intense research
has been carried out on carbon nanotubes (CNTs). The nanotubes can be open-
ended or closed by caps containing five-membered rings. They can be multi-
(MWNTs) or singlewalled (SWNTs). We show a typical high-resolution electron
microscope (HREM) image of a multi-walled nanotube in Figure 3-1.
Figure 3-1 : A typical TEM image of a closed, multi-walled carbon nanotube. The separation
between the graphite layers is 0.34 nm. [C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans.,
2003, p.p. 1-24]
Depending on the way the graphene sheets fold, nanotubes are classified as
armchair, zigzag or chiral as shown in Figure 3-2. The electrical conductivity of the
nanotubes depends on the nature of folding.
NANOMATERIALS
21
Figure 3-2 : Schematic representation of the folding of a graphene sheet into (a) zigzag, (b)
armchair and (c) chiral nanotubes. [C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans., 2003,
p.p. 1-24]
Several layered inorganic compounds possess structures comparable to the structure
of graphite, the metal dichalcogenides being important examples. The metal
dichalcogenides, MX
2
(M = Mo, W, Nb, Hf; X = S, Se) contain a metal layer
sandwiched between two chalcogen layers with the metal in a trigonal pyramidal or
octahedral coordination mode. The MX
2
layers are stacked along the c-direction in
ABAB fashion. The MX
2
layers are analogous to the single graphene sheets in the
graphite structure (Figure 3-3).
NANOMATERIALS
22
Figure 3-3 : Comparison of the structures of (a) graphite and inorganic layered compounds such as
(b) NbS
2
/TaS
2
; (c) MoS
2
; (d) BN. In the layered dichalcogenides, the metal is in trigonal prismatic
(TaS
2
) or octahedral coordination (MoS
2
). [C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans.,
2003, p.p. 1-24]
When viewed parallel to the c-axis, the layers show the presence of dangling bonds
due to the absence of an X or M atom at the edges. Such unsaturated bonds at the
edges of the layers also occur in graphite. The dichalcogenide layers are unstable
towards bending and have a high propensity to roll into curved structures. Folding in
the layered transition metal chalcogenides (LTMCs) was recognized as early as
1979, well before the discovery of the carbon nanotubes. Rag-like and tubular
structures of MoS
2
were reported by Chianelli who studied their usefulness in
catalysis.
The folded sheets appear as crystalline needles in low magnification transmission
electron microscope (TEM) images, and were described as layers that fold onto
themselves (Figure 3-4).
NANOMATERIALS
23
Figure 3-4 : Low-magnification TEM images of (a) highly folded MoS2 needles and (b) a rolled sheet
of MoS2 folded back on itself. [C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans., 2003, p.p. 1-
24]
These structures indeed represent those of nanotubes. Tenne et al
xxii
first
demonstrated that Mo and W dichalcogenides are capable of forming nanotubes
(Figure 3-5 a).
Figure 3-5 : TEM images of (a) a multi-walled nanotube of WS
2
and (b) hollow particles (inorganic
fullerenes) of WS
2
. [C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans., 2003, p.p. 1-24]
Closed fullerene-type structures (inorganic fullerenes) also formed along with the
nanotubes (Figure 3-5 b). The dichalcogenide structures contain concentrically
nested fullerene cylinders, with a less regular structure than in the carbon nanotubes.
Accordingly, MX
2
nanotubes have varying wall thickness and contain some
amorphous material on the exterior of the tubes. Nearly defect-free MX
2
nanotubes
are rigid as a consequence of their structure and do not permit plastic deformation.
The folding of a MS
2
layer in the process of forming a nanotube is shown in the
schematic in Figure 3-6.
NANOMATERIALS
24
Figure 3-6 : Schematic illustration of the bending of a MoS
2
layer. [C.N.R. Rao, M. Nath, Inorganic
nanotubes, Dalton Trans., 2003, p.p. 1-24]
Considerable progress has been made in the synthesis of the nanotubes of Mo and
W dichalcogenides in the last few years (Table 3-1 and Table 3-2).
Table 3-1 : Synthetic strategies for various chalcogenide nanotubes [C.N.R. Rao, M. Nath,
Inorganic nanotubes, Dalton Trans., 2003, p.p. 1-24]
Table 3-2 : Synthetic strategies for various chalcogenide nanotubes [C.N.R. Rao, M. Nath,
Inorganic nanotubes, Dalton Trans., 2003, p.p. 1-24]
NANOMATERIALS
25
There has been some speculation on the cause of folding and curvature in the
LTMCs. Stoichiometric LTMC chains and layers such as those of TiS
2
possess an
inherent ability to bend and fold, as observed in intercalation reactions.
The existence of alternate coordination and therefore of stoichiometry in the LTMCs
may also cause folding. Lastly, a change in the stoichiometry within the material
would give rise to closed rings.
Transition metal chalcogenides possess a wide range of interesting physical
properties. They are widely used in catalysis and as lubricants. They have both
semiconducting and superconducting properties (see paragraph 7). With the
synthesis and characterization of the fullerenes and nanotubes of MoS
2
and WS
2
, a
wide field of research has opened up enabling the successful synthesis of nanotubes
of other metal chalcogenides. It may be recalled that the dichalcogenides of many of
the Group 4 and 5 metals have layered structures suitable for forming nanotubes.
Curved structures are not only limited to carbon and the dichalcogenides of Mo and
W. Perhaps the most well-known example of a tube-like structure with diameters in
the nm range is formed by the asbestos mineral (chrysotil) whose fibrous
characteristics are determined by the tubular structure of the fused tetrahedral and
octahedral layers. The synthesis of mesoporous silica with well-defined pores in the
220 nm range was reported by Beck and Kresge
xxiii
. The synthetic strategy involved
the self-assembly of liquid crystalline templates. The pore size in zeolitic and other
NANOMATERIALS
26
inorganic porous solids is varied by a suitable choice of the template. However, in
contrast to the synthesis of porous compounds, the synthesis of nanotubes is
somewhat more difficult.
Nanotubes of oxides of several transition metals as well as of other metals have been
synthesized employing different methodologies. Silica nanotubes were first produced
as a spin-off product during the synthesis of spherical silica particles by the
hydrolysis of tetraethylorthosilicate (TEOS) in a mixture of water, ammonia, ethanol
and D,L-tartaric acid. Since selfassembly reactions are not straightforward with
respect to the desired product, particularly its morphology, templated reactions have
been employed using carbon nanotubes to obtain nanotube structures of metal
oxides. Oxides such as V
2
O
5
have good catalytic activity in the bulk phase. Redox
catalytic activity is also retained in the nanotubular structure. There have been efforts
to prepare V
2
O
5
nanotubes by chemical methods as well.
Boron nitride (BN) crystallizes in a graphite-like structure and can be simply viewed
as replacing a CC pair in the graphene sheet with the iso-electronic BN pair. It can,
therefore, be considered as an ideal precursor for the formation of BN nanotubes.
Replacement of the CC pairs partly or entirely by the BN pairs in the hexagonal
network of graphite leads to the formation of a wide array of two-dimensional phases
that can form hollow cage structures and nanotubes. The possibility of replacing CC
pairs by BN pairs in the hollow cage structure of C
60
was predicted and verified
experimentally. BN-doped carbon nanotubes have been prepared. Pure BN
nanotubes have been generated by employing several procedures, yielding
nanotubes with varying wall thickness and morphology. It is therefore quite possible
that nanotube structures of other layered materials can be prepared as well. For
example, many metal halides (e.g., NiCl
2
), oxides (GeO
2
) and nitrides (GaN)
crystallize in layered structures. There is considerable interest at present to prepare
exotic nanotubes and to study their properties.
3.2.2 General synthetic strategies
Several strategies have been employed for the synthesis of carbon nanotubes. They
are generally made by the arc evaporation of graphite or by the pyrolysis of
hydrocarbons such as acetylene or benzene over metal nanoparticles in a reducing
atmosphere. Pyrolysis of organometallic precursors provides a one-step synthetic
NANOMATERIALS
27
method of making carbon nanotubes. In addition to the above methods, carbon
nanotubes have been prepared by laser ablation of graphite or electron-beam
evaporation. Electrochemical synthesis of nanotubes as well as growth inside the
pores of alumina membranes have also been reported. The above methods broadly
fall under two categories. Methods such as the arc evaporation of graphite employ
processes which are far from equilibrium. The chemical routes are generally closer to
equilibrium conditions. Nanotubes of metal chalcogenides and boron nitride are also
prepared by employing techniques similar to those of carbon nanotubes, although
there is an inherent difference in that the nanotubes of inorganic materials such as
MoS
2
or BN would require reactions involving the component elements or
compounds containing the elements. Decomposition of precursor compounds
containing the elements is another possible route.
Nanotubes of dichalcogenides such as MoS
2
, MoSe
2
and WS
2
are also obtained by
employing processes far from equilibrium such as arc discharge and laser ablation.
By far the most successful routes employ appropriate chemical reactions. Thus,
MoS
2
and WS
2
nanotubes are conveniently prepared starting with the stable oxides,
MoO
3
and WO
3
. The oxides are first heated at high temperatures in a reducing
atmosphere and then reacted with H
2
S. Reaction with H
2
Se is used to obtain the
selenides. Recognizing that the trisulfides MoS
3
and WS
3
are likely to be the
intermediates in the formation of the disulfide nanotubes, the trisulfides have been
directly decomposed to obtain the disulfide nanotubes. Diselenide nanotubes have
been obtained from the metal triselenides. The trisulfide route is indeed found to
provide a general route for the synthesis of the nanotubes of many metal disulfides
such as NbS
2
and HfS
2
. In the case of Mo and W dichalcogenides, it is possible to
use the decomposition of the precursor ammonium salt, such as (NH
4
)
2
MX
4
(X = S,
Se; M = Mo, W) as a means of preparing the nanotubes. Other methods employed
for the synthesis of dichalcogenide nanotubes include hydrothermal methods where
the organic amine is taken as one of the components in the reaction mixture (Table
3-1 and Table 3-2).
The hydrothermal route has been used for synthesizing nanotubes and related
structures of a variety of other inorganic materials as well. Thus, nanotubes of
several metal oxides (e.g., SiO
2
, V
2
O
5
, ZnO) have been produced hydrothermally.
Nanotubes of oxides such as V
2
O
5
are also conveniently prepared from a suitable
metal oxide precursor in the presence of an organic amine or a surfactant.
NANOMATERIALS
28
Surfactant-assisted synthesis of CdSe and CdS nanotubes has been reported. Here
the metal oxide reacts with the sulfidizing/selenidizing agent in the presence of a
surfactant such as TritonX.
Solgel chemistry is widely used in the synthesis of metal oxide nanotubes, a good
example being that of silica and TiO
2
. Oxide gels in the presence of surfactants or
suitable templates form nanotubes. For example, by coating carbon nanotubes
(CNTs) with oxide gels and then burning off the carbon, one obtains nanotubes and
nanowires of a variety of metal oxides including ZrO
2
, SiO
2
and MoO
3
. Solgel
synthesis of oxide nanotubes is also possible in the pores of alumina membranes. It
should be noted that MoS
2
nanotubes are also prepared by the decomposition of a
precursor in the pores of an alumina membrane.
Boron nitride nanotubes have been obtained by striking an electric arc between HfB
2
electrodes in a N
2
atmosphere. BCN and BC nanotubes are obtained by arcing
between B/C electrodes in an appropriate atmosphere. A greater effort has gone into
the synthesis of BN nanotubes starting with different precursor molecules containing
B and N. Decomposition of borazine in the presence of transition metal nanoparticles
and the decomposition of the 1 : 2 melamineboric acid addition compound yield BN
nanotubes. Reaction of boric acid or B
2
O
3
with N
2
or NH
3
at high temperature in the
presence of activated carbon, carbon nanotubes or catalytic metal particles has been
employed to synthesize BN nanotubes.
3.3 Functional Materials Based on Self-Assembly of
Polymeric Supramolecules
xxiv
Here, we describe some possibilities for preparing functional polymeric materials
using the "bottom-up" route, based on self-assembly of polymeric supramolecules.
Directed assembly leads to the control of structure at several length scales and
anisotropic properties. The physical bonds within the supramolecules allow controlled
cleavage of selected constituents. The techniques constitute a general platform for
constructing materials that combine several properties that can be tuned separately.
To achieve enhanced functionalities, the principal periodicity is at ~10 to 2000 .
There are established ways to accomplish this by using various architectures of block
copolymers, in which the structure formation is based on self-organization, that is, on
the repulsion between the chemically connected blocks. Depending on the
NANOMATERIALS
29
architecture, block length, and temperature, it is possible to obtain lamellar,
cylindrical, spherical, gyroid, or more complicated structures in the 100 to 2000
range. Also, rodlike moieties within the block copolymers can be used to further tailor
the structures in terms of shape persistency. However, self-organization renders only
the local structures. To fully realize the opportunities offered by the symmetry of the
self-organized structures to prepare materials with a strongly directional variation of
properties, additional mechanisms and interactions have to be invoked to obtain
macroscale order. This may be achieved by flow, by electric or magnetic fields, or by
using topographically patterned surfaces. One can further extend the structural
complexity by mixing block copolymers with additional polymers and inorganic
additives, thereby increasing the self-organization periods into the photonic band gap
regime. Block copolymers have also been used as templates for the synthesis of
inorganic materials, even allowing the creation of separate ceramic nano-objects.
To achieve even greater structural complexity and functionality, we can combine
recognition with self-organization. Lehn elaborated on the concept of recognition in
synthetic materials, whereby two molecules with molecularly matching
complementary interactions and shapes recognize each other and form a receptor-
substrate supramolecule. To achieve sufficient bonding, synergism of several
physical interactions is often required. Homopolymerlike supramolecules have been
constructed based on a combination of four hydrogen bonds and through
coordination. Supramolecules can spontaneously assemble or self-organize to form
larger structures.
A general framework for forming complex functional materials emerges. Molecules
are constructed that recognize each other in a designed way. The subsequent
supramolecules in turn form assemblies or self-organize, possibly even forming
hierarchies. The overall alignment of the local structures can be additionally improved
by electric or magnetic fields, by flow, or by patterned surfaces.
To illustrate recognition-driven supramolecule formation in polymers and the
subsequent self-organization and preparation of functional materials and nano-
objects, we focus on the comb-shaped architecture (Figure 3-7) encouraged by the
enhanced solubility of socalled hairy-rod polymers.
NANOMATERIALS
30
Figure 3-7 : Comb-shaped supramolecules and their hierarchical self-organization, showing primary
and secondary structures. Similar schemes can, in principle, be used both for flexible and rodlike
polymers. In the first case, simple hydrogen bonds can be sufficient, but in the latter case a synergistic
combination of bondings (recognition) is generally required to oppose macrophase separation
tendency. In (A through C), the self-organized structures allow enhanced processibility due to
plastization, and solid films can be obtained after the side chains are cleaved (D). Self-organization of
supramolecules obtained by connecting amphiphiles to one of the blocks of a diblock copolymer (E)
results in hierarchically structured materials. Functionalizable nanoporous materials (G) are obtained
by cleaving the side chains from a lamellae-within-cylinders structure (F). Disk-like objects (H) may be
prepared from the same structure by crosslinking slices within the cylinders, whereas nano rods (I)
NANOMATERIALS
31
result from cleaving the side chains from a cylinder-within-lamellae structure. Without loss of
generality, (A) is shown as a flexible polymer, whereas (B) and (C) are shown as rodlike chains. [O.
Ikkala, G.T. Brinke, Functional Materials Based on Self-assembly of polymeric supramolecules,
Science, New Series, Vol. 295, No. 5564 (Mar. 29, 2002), pp. 2407-2409]
The simplest case is a flexible polymer having bonding sites along its backbone
(Figure 3-7 A). Therefore, the backbone is typically polar, and repulsive nonpolar side
groups can be connected by complementary bonds, leading to comb-shaped
supramolecules, which in turn self-organize. We have extensively used hydrogen
bonding or coordination to bond side chains to the polymer backbone. Antonietti et
al.
xxv
have used ionic interactions in polyelectrolyte-surfactant complexes to form
comb-shaped polyelectrolyte surfactant complexes. The resulting self-organized
multidomain structures may be aligned, using, for example, flow, in order to approach
monodomains. One can also tune the properties by tailoring the nature of the side
chains. For example, if the side chains are partly fluorinated, low surface energy
results, which allows for applications that lead to reduced friction. In another case,
the backbone consists of the double helix of DNA, and self-organization is achieved
by ionically bonding cationic liposomes or cationic surfactants to the anionic
phosphate sites. This allows for materials design beyond the traditional scope of
biochemical applications. For example, dyes can be intercalated into the helices,
suppressing their aggregation tendency and leading to promising properties as
templates for photonic applications. In such a structure, the polymer backbone may
contain two or even more kinds of binding sites where different additives can be
bonded (Figure 3-7 B). Side chains can also have two separate functions. For
example, in addition to providing a repulsive side chain required for self-organization,
the side chains may contain an acidic group that acts as a dopant for a conjugated
polymer such as polyaniline, which leads to electronic conductivity. To introduce
further degrees of freedom in tailoring the self-organized phases and their
processing, polyaniline may first be doped by a substance such as camphor
sulphonic acid and subsequently connected to hexyl resorcinol molecules using their
two hydrogen bonds (Figure 3-7 C). The alkyl chains of the hydrogen-bonded hexyl
resorcinol molecules act as plasticizers, leading to thermoplastic processibility of the
otherwise infusible polymer. They enforce self-organization where camphor sulfonic
acid-doped polyaniline chains are confined in nanoscale conducting cylinders,
leading to increased conductivity. The concept can be applied even to rodlike
NANOMATERIALS
32
polymers, such as polypyridine, which consists of para-coupled heteroaromatic rings.
Its optical properties can be tuned based on camphor sulphonic acid. Subsequent
hydrogen bonding with alkyl resorcinol creates comb-shaped supramolecules, which
self-organize in lamellae in such a way that the material is fluid even without
additional solvents. Such a fluid state incorporating rigid polymeric rods is uncommon
and allows processing toward monodomains where the rods are aligned. Ultimately,
the plasticizing hydrogen-bonded alkyl resorcinol molecules can be removed by
evaporation in a vacuum oven, thus interlocking the chains in solid stable films
(Figure 3-7 D). In this way, efficient polarized luminance has been achieved.
To increase complexity, one can incorporate structural hierarchies. This can be
accomplished by applying within a single material different self-organization and
recognition mechanisms operating at different length scales. For example, block
copolymeric self-organization at the 100 to 2000 length scale and polymer-
amphiphile self-organization at the 10 to 60 length scale can be combined (Figure
3-7 E). After selective doping of one block, conductivity can be switched based on a
sequence of phase transitions.
3.4 Molecular biomimetics: Nanotechnology through
biology
xxvi
Molecular biomimetics. This is the marriage of materials science engineering and
molecular biology for development of functional hybrid systems, composed of
inorganics and inorganic-binding proteins. The new approach takes advantage of
DNA-based design, recognition,and self-assembly characteristics of biomolecules.
Traditional materials science engineering produces materials (for example, medium-
carbon steels depicted in the bright- and dark-field TEM images), that have been
successfully used over the last century. Molecular biology focuses on structure
function relations in biomacromolecules, for example, proteins.
NANOMATERIALS
33
Figure 3-8 : In molecular biomimetics, inorganic-binding proteins could potentially be used as (i)
linkers for nanoparticle immobilization; (ii) functional molecules assembled on specific substrates; and
(iii) heterobifunctional linkers involving two (or more) binding proteins linking several nanoinorganic
units. (I1: inorganic-1,I2: Inorganic-2, P1 and P2: inorganic specific proteins, LP:linker protein, FP:
fusion protein). [Sarikaya, C. Tamerler, A.K.Y. Jen, K. Schulten F. Baneyx, Molecular biomimetics:
nanotechnology through biology, Nature Materials, vol 2, 2003, p.p. 577-585]
In molecular biomimetics, a marriage of the physical and biological fields, hybrid
materials could potentially be assembled from the molecular level using the
recognition properties of proteins (Figure 3-8) under the premise that inorganic
surface-specific polypeptides could be used as binding agents to control the
organization and specific functions of materials.Molecular biomimetics simultaneously
offers three solutions to the development of heterofunctional nanostructures.
The first is that protein templates are designed at the molecular level through
genetics. This ensures complete control over the molecular structure of the
protein template (that is, DNA-based technology).
The second is that surface-specific proteins can be used as linkers to bind
synthetic entities, including nanoparticles, functional polymers, or other
nanostructures onto molecular templates (molecular and nanoscale
recognition).
The third solution harnesses the ability of biological molecules to self- and co-
assemble into ordered nanostructures. This ensures a robust assembly
NANOMATERIALS
34
process for achieving complex nano-, and possibly hierarchical structures,
similar to those found in nature (self-assembly).
The current knowledge of protein-folding predictions and surface-binding chemistries
does not provide sufficiently detailed information to perform rational design of
proteins. To circumvent this problem, massive libraries of randomly generated
peptides can be screened for binding activity to inorganic surfaces using phage and
cell-surface display techniques. It may ultimately be possible to construct a
molecular erector set, in which different types of proteins, each designed to bind to a
specific inorganic surface, could assemble into intricate, hybrid structures composed
of inorganics and proteins. This would be a significant leap towards realizing
molecularly designed, genetically engineered technological materials.
3.4.1 Selection of inorganic-binding proteins through display
technologies
There are several possible ways of obtaining polypeptide sequences with specific
affinity to inorganics. A number of proteins may fortuitously bind to inorganics,
although they are rarely tested for this purpose. Inorganic-binding peptides may be
designed using a theoretical molecular approach similar to that used for
pharmaceutical drugs. This is currently impractical because it is time consuming and
expensive. Another possibility would be to extract biomineralizing proteins from hard
tissues followed by their isolation, purification and cloning. Several such proteins
have been used as nucleators, growth modifiers, or enzymes in the synthesis of
certain inorganics. One of the major limitations of this approach is that a given hard
tissue usually contains many proteins, not just one, all differently active in
biomineralization and each distributed spatially and temporally in complex ways.
Furthermore, tissue-extracted proteins may only be used for the regeneration of the
inorganics that they are originally associated with, and would be of limited practical
use. The preferred route, therefore, is to use combinatorial biology techniques. Here,
a large random library of peptides with the same number of amino acids, but of
different sequences, is used to mine specific sequences that strongly bind to a
chosen inorganic surface.
Since their inception, well-established in vivo combinatorial biology protocols (for
example, phage display (PD) and cell-surface display (CSD)) have been used to
NANOMATERIALS
35
identify biological ligands and to map the epitope (molecular recognition site) of
antibodies. Libraries have also been screened for various biological activities, such
as catalytic properties or altered affinity and specificity to target molecules in many
applications including the design of new drugs, enzymes, antibodies, DNA-binding
proteins and diagnostic agents. The power of display technologies relies on the fact
that an a priori knowledge of the desired amino acid sequence is not necessary, as it
can simply be selected and enriched if a large enough population of random
sequences is available. In vitro methods, such as ribosomal and messenger RNA
display technologies, have been developed for increased library size (10
15
) compared
to those of in vivo systems (10
710
).
Combinatorial biology protocols can be followed in molecular biomimetics to select
polypeptide sequences that preferentially bind to the surfaces of inorganic
compounds chosen for their unique physical properties in nano- and biotechnology.
Libraries are generated by inserting randomized oligonucleotides within certain
genes encoded on phage genomes or on bacterial plasmids (step 1 in Figure 3-9).
Figure 3-9 : Phage display and cell-surface display. Principles of the protocols used for selecting
polypeptide sequences that have binding affinity to given inorganic substrates. [Sarikaya, C. Tamerler,
A.K.Y. Jen, K. Schulten F. Baneyx, Molecular biomimetics: nanotechnology through biology, Nature
Materials, vol 2, 2003, p.p. 577-585]
NANOMATERIALS
36
This leads to the incorporation of a random polypeptide sequence within a protein
residing on the surface of the organism (for example, the coat protein of a phage or
an outer membrane or flagellar protein of a cell; step 2). The eventual result is that
each phage or cell produces and displays a different, but random peptide (step 3). At
this stage, a heterogeneous mixture of recombinant cells or phages are contacted
with the inorganic substrate (step 4).Several washing cycles of the phages or the
cells eliminate non-binders by disrupting weak interactions with the substrate (step
5). Bound phages or cells are next eluted from the surfaces (step 6). In PD, the
eluted phages are amplified by reinfecting the host (step 7). Similarly in CSD, cells
are allowed to grow (steps 7, 8). This step completes a round of biopanning.
Generally, three to five cycles of biopanning are repeated to enrich for tight binders.
Finally, individual clones are sequenced (step 9) to obtain the amino acid sequence
of the polypeptides binding to the target substrate material.
3.4.2 Chemical specificity of inorganic-binding polypeptides
A genetically engineered polypeptide for inorganics (GEPI) defines a sequence of
amino acids that specifically and selectively binds to an inorganic surface. The
surface could be well defined, such as a single crystal or a nanostructure. It might
also be rough, or totally non-descriptive, such as a powder. Researchers have
focused on using materials that can be synthesized in aqueous environments under
physiological conditions (biocompatible) and that exhibit fairly stable surface
structures and compositions. These include noble metals (Pt and Pd) as well as
oxide semiconductors (Cu
2
O and ZnO) that were biopanned using either PD or
flagellar display (both studies unpublished). Some of the identified binders as well as
sequences selected by other researchers are listed in Table 3-3.
Table 3-3 : Examples of polypeptide sequences exhibiting affinity for various inorganics.
[Sarikaya, C. Tamerler, A.K.Y. Jen, K. Schulten F. Baneyx, Molecular biomimetics: nanotechnology
through biology, Nature Materials, vol 2, 2003, p.p. 577-585]
NANOMATERIALS
37
3.4.3 Physical specificity of peptide binding
Ideally, selection of sequences should be performed using an inorganic material of
specific morphology, size, crystallography or surface stereochemistry. In practice,
however, powders of various sizes and morphologies have been used for selection.
The sequence space should be largest for powders, as peptides can attach to
NANOMATERIALS
38
surfaces with various morphological features. On the other hand, because powders
are non-descriptive, the selected polypeptides may exhibit little or no homology and
solve the binding problem through different strategies. Binders selected for a given
size, morphology, crystallography or stereochemistry may share a higher degree of
homology. For example, a GEPI binding to a material of a certain size may also bind
to a smaller particle of the same material, but less strongly. Similarly, a GEPI binding
strongly to a specific crystallographic surface may bind with an altered affinity to
another surface of the same material. Finally, a GEPI strongly binding to a material of
composition A may bind less strongly to a material B with a different composition but
having similar structure (for example, perovskites). Therefore, if one seeks highly
specific binders, the physical and chemical characteristics of the material must be
known. An alternative approach is that, once a relatively large number of binders
have been identified by panning on powders, a subset specific for morphology, size
or surface could be identified on well-defined materials.
3.4.4 Peptide-mediated nanoparticle assembly
Organization and immobilization of inorganic nanoparticles in two- or three-
dimensional geometries are fundamental in the use of nanoscale effects.
a)
b)
c)
f)
NANOMATERIALS
39
d) e)
Figure 3-10 : Effect of GEPI on nanocrystal morphology. ac, One of the two mutants (RP1) from a
library of goldbinding GEPIs were tested in the formation of flat gold particles, shown in a, similar to
those formed under acidic (b) or boiling (c) conditions. Particles formed in the presence of vector-
encoded alkaline phosphatase and neutral conditions do not result in morphological change of gold
particles (not shown). d-e, The atomic force microscope images show quantum (GaInAs) dots
assembled on GaAs substrate; d, through high-vacuum (molecular beam epitaxy) strain-induced self-
assembly, and e, through 7-repeat GBP1. f, Schematic illustration of e. PS: polystyrene substrate, GA:
glutaraldehyde, GBP: 7- repeat GBP1, and gold: 12-nm-diameter colloidal gold particles. [Sarikaya, C.
Tamerler, A.K.Y. Jen, K. Schulten F. Baneyx, Molecular biomimetics: nanotechnology through biology,
Nature Materials, vol 2, 2003, p.p. 577-585]
For example, quantum dots can be produced using vacuum techniques, such as
molecular beam epitaxy, shown in Figure 3-10 d for the GaInAs/GaAs system.
However, this can only be accomplished under stringent conditions of high
temperature, very low pressures and a toxic environment. A desirable alternative
would be not only to synthesize inorganic nanodots under mild conditions, but also to
immobilize/self-assemble them. Inorganic particles have been functionalized with
synthetic molecules, including thiols and citrates, and with biological molecules, such
as lipids, amino acids, polypeptides and ligand-functionalized DNA. Using the
recognition properties of the coupling agents, novel materials have been generated
and controlled growth has been achieved. These molecules, however, do not exhibit
specificity for a given material. For example, thiols couple gold as well as silver
nanoparticles in similar ways. Likewise, citrate ions cap noble metals indiscriminately.
A desirable next step would be to use GEPIs that specifically recognize inorganics for
nanoparticle assembly. An advantage of this approach is that GEPI can be
genetically or synthetically fused to other functional biomolecular units or ligands to
produce heterobifunctional (or multifunctional) molecular entities. Figure 3-10 e and f
shows the assembly of nanogold particles on GBP1-coated flat polystyrene surfaces,
which resembles the distribution of quantum dots obtained by high-vacuum
deposition techniques (Figure 3-10 d). The homogenous decoration of the surface
with nanogold suggests that proteins may be useful in the production of tailored
nanostructures under ambient conditions and aqueous solutions. Furthermore, the
recognition activity of the protein could provide an ability to control the particle
distribution, and particle preparation conditions could allow size control. This
approach makes it possible to pattern inorganic-binding polypeptides into desirable
NANOMATERIALS
40
arrays to produce inorganic particles through templating using, for example, dip-pen
lithography.
3.5 References
Sarikaya, C. Tamerler, A.K.Y. Jen, K. Schulten F. Baneyx, Molecular biomimetics:
nanotechnology through biology, Nature Materials, vol 2, 2003, p.p. 577-585
O. Ikkala, G.T. Brinke, Functional Materials Based on Self-assembly of polymeric
supramolecules, Science, New Series, Vol. 295, No. 5564 (Mar. 29, 2002), pp. 2407-
2409
C.N.R. Rao, M. Nath, Inorganic nanotubes, Dalton Trans., 2003, p.p. 1-24
NANOMATERIALS
41
4 MECHANICAL PROPERTIES
4.1 Introduction
One of the very basic results of the physics and chemistry of solids is the insight that
most properties of solids depend on the microstructure, i.e. the chemical composition,
the arrangement of the atoms (the atomic structure) and the size of a solid in one,
two or three dimensions. In other words, if one changes one or several of these
parameters, the properties of a solid vary. The most well-known example of the
correlation between the atomic structure and the properties of a bulk material is
probably the spectacular variation in the hardness of carbon when it transforms from
diamond to graphite. The important aspects related to structure are:
atomic defects, dislocations and strains
grain boundaries and interfaces
porosity
connectivity and percolation
short range order
Defects are usually absent in either metallic or ceramic clusters of
nanoparticles because dislocations are basically unstable or mobile. The stress field
around a dislocation (or the electrostatic potential around charges and currents)
have to satisfy the Laplace equation:
2
= 0. This sets up an image dislocation
which pulls the defect out of the particle.
When these clusters are assembled under uniaxial pressure into a pellet, for
example, it is found that the individual clusters are packed very tightly into a
polycrystalline solid. Cluster-assembled materials often show close to 100% density.
A fully consolidated nanophase material looks very much like a normal, dense
polycrystalline aggregate, but at a far smaller scale.
4.2 Metals
The microstructure of a material is controlled by the processing steps chosen for its
fabrication. Such microstructural design affects the nature of the phases present,
their topology (i.e. geometrical distribution and interconnection) and their dispersion
NANOMATERIALS
42
(described by relevant size parameters). The full characterization of these
parameters is the domain of quantitative metallurgy.
Most of these size effects come about because of the microstructural constraint to
which a particular physical mechanism is subjected. Consider the classic case of
strengthening a metallic matrix by particles or grain boundaries: lattice dislocations
are forced, by the microstructural constraint, to bow out or pile up, which requires an
external stress characteristic of a microstructural parameter. The wall thickness
relative to the size of the microstructural inhomogeneity can control the macroscopic
behaviour.
In general, it is therefore the competition or coupling between two different size
dependencies that determines the properties of a material. One thus has to deal with
the interaction of two length scales: (1) is the dimension characteristic of the physical
phenomenon involved, called the characteristic length. (2) is some microstructural
dimension, denoted as the size parameter.
4.2.1 Grain size effects in plasticity and creep
4.2.1.1 Hall-Petch effects.
Strengthening of polycrystalline materials by grain size refinement is
technologically attractive because it generally does not adversely affect ductility
and toughness. The classical effect of grain size on yield stress ( ) can, among
other possibilities, be explained by a model invoking a pile-up of dislocations
against grain boundaries, which results in a dependence of the hardening
increment k
HP
on the square root of the grain size D
D
k
HP
= (4. 1)
Where k
HP
is a constant. This is the classical Hall-Petch effect.
In Figure 4-1, some of the available data have been plotted in a Hall-
Petch plot.
xxvii
NANOMATERIALS
43
Figure 4-1 : Compilation of yield stress data for several metallic systems. [R.A. Masumura, P.M.
Hazzledine, C.S. Pande, Yield stress of fine grained materials, 1998, Acta-Materialia]
It is seen that the yield stress-grain size exponent for relatively large grains
appears to be very close to -0.5 and generally this trend continues until the very fine
grain regime (~100 nm) is reached. The reported data show three different regions:
1. A region from single crystal to a grain size of about 1 mm where the classical
Hall-Petch description can be used.
2. A region for grain sizes ranging from about 1 mm to about 30 nm where the
Hall-Petch relation roughly holds, but deviates from the classical -0.5
exponent to a value near zero (to ascertain such behaviour, a wide range of
grain sizes extending into the ultra-fine grain size regime is required).
3. A region beyond a very small critical grain size where the Hall-Petch slope is
essentially zero, with no increase in strength on decreasing grain size or
where the strength actually decreases with decreasing grain size.
There is universal agreement regarding the first region, i.e. relatively large grain
sizes. Early hardness measurements had already established a distinct increase in
hardness as grain sizes decrease as compared to their annealed coarse grained
counterparts, and this increase follows the Hall-Petch relationship reasonably well.
The trend is less well established for finer grains (Region 2). Some of the deviation
from Hall-Petch strengthening could simply be due to pores in the material (as
evidenced by lower densities) leading also to a lower shear stress for the
NANOMATERIALS
44
deformation mode and lower shear modulus. Indeed, the lower modulus has been
ascribed to a decrease in bulk density. Additional complications arise due to
impurities at the grain boundaries such as oxides and impurities inside the grain
such as trapped or diffused gas. In spite of the above difficulties, once the totality of
the data is taken into consideration, it is fairly safe to conclude that the increase in
strength on grain refinement in the middle region is somewhat less than predicted by
the Hall-Petch relation.
The third region is much more controversial and is going to be discussed later.
4.2.1.2 Limits to Hall-Petch behaviour: dislocation curvature vs. grain size.
Whereas many metallic materials obey such a relationship over several
orders of magnitude in grain size, it is inevitable that the reasoning behind equation
D
k
HP
= (4. 1) must break down for very small grains. A clear limit for the
occurrence of dislocation plasticity in a poly-crystal is given by the condition that at
least one dislocation loop must fit into average grain (Figure 4-2).
Figure 4-2 : Grain size strengthening, as explained by pile-ups of dislocation loops against grain
boundaries (a). this mechanism must break down when the diameter d of the smallest loop no longer
fits into a grain of size D (b). [E. Arzt, Size effects in materials due to microstuctural and dimensional
constraints: A comparative review, 1998, Acta Meteriala]
NANOMATERIALS
45
The characteristic length, i.e. the loop diameter (
Gb
b
T
d
d
= =
2
) ( ), must now be
compared with the grain size D as the relevant size parameter
D d = ) ( (4. 2)
D
Gb
bD
T
d
= =
2
(4. 3)
Figure 4-3 illustrates schematically this limit on Hall-Petch behaviour:
conventional grain size strengthening can be expected only to the right of the heavy
line which signifies the limiting condition D d = ) ( (4. 2)
D
Gb
bD
T
d
= =
2
(4. 3).
Figure 4-3 : The limiting condition is shown as the heavy line where the shear strength is plotted
schemetically as a function of grain size D. Hall-Petch behavior can only be found to the right of this
line; abnormal or inverse behavior may result otherwise. The dotted line reflects schematically the
equation
0
ln
r
D
D
Gb
(4. 6). [E. Arzt, Size effects in materials due to microstuctural and
dimensional constraints: A comparative review, 1998, Acta Meteriala]
For Cu, as an example, the critical grain size estimated in this way is about
50nm; this value is in reasonable agreement with experimental results, as shown in
Figure 4-4.
NANOMATERIALS
46
Figure 4-4 : Inverse Hall-Petch behaviour in nanocrystalline Cu (H-H
0
denotes the hardness increment,
D the grain size): the classical behaviour breaks down at a grain size of about 50 nm, in agreement with
an estimate based on the loop diameter [equations
D d = ) (
(4. 2)
D
Gb
bD
T
d
= =
2
(4. 3)].
Replotted after Chokshi et al.
xxviii
[E. Arzt, Size effects in materials due to microstuctural and
dimensional constraints: A comparative review, 1998, Acta Meteriala]
The plastic behaviour of nanocrystalline materials with grain sizes below the
critical value is not fully clear. It has been argued that because of the viscous
behaviour of amorphous materials (which can be considered the limiting case for grain
refinement) the grain size strengthening effect will have to be reversed once the grain
size D starts to approach the grain-boundary thickness
b
. One possible explanation
for such a softening effect comes from a re-consideration of the line tension T
d
in
equation
D
Gb
bD
T
d
= =
2
(4. 3). The more refined expression
2
1
ln
4
r
r Gb
T
d
= (4. 4)
contains a lower (r
0
) and an upper (r
1
) cut-off distance for the stress field of the
dislocation. In conventional materials r
1
generally lies in the micrometer range and
therefore significantly exceeds r
0
(for which values between 2 and 10b are commonly
assumed); this justifies replacing the logarithmic term by a constant. However, in
nanocrystalline materials it is reasonable to equate r
1
to the grain size, which now
gives r
1
r
0
and makes T sensitive to the value of the grain size D. Therefore, we now
have a case in which the characteristic length (d) is a function of the size parameter
(D). The resulting strength increment is given by
NANOMATERIALS
47
0
2 r
D
n
D
Gb
= (4. 5)
This expression vanishes rapidly as the grain size D approaches the lower cut-
off distance r
0
. An even more refined expression has been obtained by Scattergood
and Koch
xxix
.
The dislocation density scales inversely with grain size D, the obstacle
spacing is L~1/ p~ D, which yields
0
ln
r
D
D
Gb
(4. 6)
This expression, which is schematically shown as dotted line in Figure 4-3,
reduces correctly to Hall-Petch behavior for D >>r
0
. It gives a possible interpretation of
grain-boundary softening behavior in nanocrystalline Cu and Pd.
4.2.1.3 Diffusional creep as a size effect
An alternative explanation of grain-boundary softening in very fine-grained materials
can be based on increasing contributions of diffusional creep. Diffusional processes in
a potential gradient [caused in this case by a normal stress gradient,] exhibit a natural
size effect because the length scale affects the magnitude of the gradient.
Figure 4-5 : Diffusional creep is driven by gradients in normal tractions on grain boundaries (a). Fine
arrows delineate the paths for transport of matter. This mechanism ceases to operate (b) once a grain
boundary dislocation loop no longer fits into a grain facet (d > D). Note the analogy with Figure 4-2 for
lattice dislocations [E. Arzt, Size effects in materials due to microstuctural and dimensional constraints:
A comparative review, 1998, Acta Meteriala]
For maintaining a constant strain rate
v
D C
kTD
1
2
(4. 7)
Here D
v
is the volume diffusivity, the atomic volume, and C
1
a dimensionless
constant of the order of 10. Accounting for grain-boundary diffusion (with diffusivity D
b
through a grain boundary with thickness
b
) gives
b b
D C
kTD
2
3
(4. 8)
In addition to this, the triple lines in nanocrystalline materials can also act as fast
diffusion paths. Equations
v
D C
kTD
1
2
(4. 7) and
b b
D C
kTD
2
3
(4. 8) reflect grain
size effects which are opposite in direction and far stronger than those of dislocation
plasticity (Hall-Petch effect). They are due to the increase, with finer grain size, in the
volume fraction of disordered' material which can act as short-circuit diffusion path,
and to the higher density of sinks and sources for matter.
It is still a matter of debate whether grain-boundary softening, which has occasionally
been reported for nanocrystalline materials, can be attributed to these effects at room
temperature.
One can note that in very small grains the rate of creep may no longer be controlled
by the diffusion step [as is tacitly assumed in equations
v
D C
kTD
1
2
(4. 7)
and
b b
D C
kTD
2
3
(4. 8)], but by the deposition and removal of atoms at the grain
boundaries. Ashby
xxx
have shown that for such interface-controlled diffusional creep
the grain size dependence is much weaker
2 / 1
2
1
4
D
D C
kTGb
eff
b
|
|
|
\
|
(4. 9)
This result was obtained by modeling the interface reaction as the climb motion of an
array of grain-boundary dislocations. Here D
eff
is an effective diffusivity, b
b
the Burgers
vector of a boundary dislocation and C
4
another numerical constant. The
NANOMATERIALS
49
D
1/2
-proportionality, which results from the assumption of a stress-dependent
dislocation density, is in better agreement with the data of Chokshi et al. (Figure 4-4).
However, because of the reduced grain size dependence, an even lower activation
energy (about 40 kJ/mol) for diffusion has to be assumed to predict realistic
deformation rates at room temperature.
Also, the motion of grain boundary dislocations is subject to a similar grain size limit as
for lattice dislocations: models based on their presence must break down once an
average grain facet of diameter D' can no longer accommodate a grain-boundary
dislocation loop [Figure 4-5(b)]. The corresponding limiting condition is, in analogy with
equation
D
Gb
bD
T
d
= =
2
(4. 3), given by
' D
Gb
b
= (4. 10)
The value of b
b
corresponds to the difference in Burgers vector between two lattice
dislocations and is therefore only a fraction of b. Hence, a stress window will exist in
which plasticity due to lattice dislocations is suppressed or slowed down [at stresses
below that given by equation
D
Gb
bD
T
d
= =
2
(4. 3)], but diffusion creep operates
because grain-boundary dislocations are still present and mobile.
4.2.2 Metal plastic deformation: A comparison between Cu and Ni
nanophase samples
In the previous paragraph we have studied the relation between yield stress and
grain size. Now two questions are still not answered:
1. Whether the grain boundaries in nanocrystalline materials are unusual or
whether they have the short main structure of most grain boundaries found in
conventional polycrystalline materials
2. What is the influence of the grain size on the plastic deformation mechanism
In this paragraph, we will focus on the question of dislocation activity in two different
materials Cu and Ni. Swygenhoven and al.
xxxi xxxii
worked on molecular dynamics
simulations on nanocrystalline Ni and Cu samples in the grain size range 5-12nm.
They studied the interfaces responsible for the plastic deformation, aiming at
NANOMATERIALS
50
providing a structural characterization of them. They also present evidence of two
competing mechanism for plastic deformation:
At the smallest grain sizes all the deformation is accommodated in the
grain boundaries.
At larger grain sizes lattice dislocation activity is observed.
For both materials, uniaxial deformation at the smallest loads reveals that Youngs
modulus equals the value for a polycrystalline material when the grain size is 10 nm
or higher. At smaller grains sizes a gradual reduction of the modulus is observed.
At high load, after a transient period following the application of the load, the strain
increases almost linearly with time for all the grain sizes.
Figure 4-6 : Strain rate as a function of mean grain size for Cu and Ni. [H. Van Swygenhoven, A.
Caro, D. Farkas, A molecular dynamics study of polycrystalline fcc metals at the nanoscale Grain
boundary structure and its influence on plastic deformation, 2001, Materials Science and Engineering
A.]
Figure 4-6 shows the strain rate versus the inverse of the grain size for the Ni and
Cu samples. The strain rates under these conditions are high compared with actual
experimental values, but for these small sample sizes (1025 nm) any relative
velocity is still four orders of magnitude smaller than the velocity of sound.
At the smallest grain sizes explored the strain rate for a given applied stress
increases with decreasing grain sizes. This behaviour indicates that the HallPetch
slope is negative at these very fine grain sizes. An energy balance indicates that at
these sizes the total amount of grain boundary remains constant during deformation.
These observations suggest an important characteristic of plasticity in nanophase
metals under the present conditions: there is no damage accumulation during
deformation, similar to the case of superplasticity. Careful examination of the
NANOMATERIALS
51
samples confirms the absence of intra-grain defects. Swygenhoven and al. [4]
discussed the deformation mechanism in terms of a model based on Grain Boundary
viscosity controlled by a self-diffusion mechanism at the disordered interface,
activated by thermal energy and stress.
When similar load is applied on samples with larger average grain sizes, the
deformation rate is much smaller. These observations indicate a transition to another
deformation mechanism. They analyzed the atomic structure of Ni and Cu samples
with larger grain sizes, deformed at those stress levels that give approximately the
same strain rate as in the sample with the very small grain sizes, and compare the
structures at similar values of plastic deformation. In this way they take into account
the different elastic contribution in very fine-grained samples due to the reduction of
the Youngs modulus.
Figure 4-7 : Slice of the 12 Ni sample, deformed until a plastic deformation level of 1.4%. Open gray
symbols are perfect f.c.c., full gray are good f.c.c., red are h.c.p., green and blue are other 12- and
non-12-coordinated atoms, respectively. [H. Van Swygenhoven, M. Spaczer, A. Caro, Microscopic
description of plasticity in computer generated metallic nanophase samples: A comparison between
Cu and Ni, 1999, Acta Materialia.]
Figure 4-7 shows a section of Ni 12 nm sample, deformed to a total deformation level
of 2.7% which means a plastic deformation level of 1.4% using a tensile load of 2.6
GPa. The stacking fault observed in this section is produced by motion of Shockley
partial dislocations generated and absorbed in opposite grain boundaries. Figure 4-8
shows the Shockley partial in Ni 12 nm just emitted from a triple point.
NANOMATERIALS
52
Figure 4-8 : Shockley partial dislocation traveling through a grain in the 12 nm Ni sample. Black atoms
are the hcp atoms forming the intrinsic stacking fault in a (1, 1, 1) plane. [H. Van Swygenhoven, A.
Caro, D. Farkas, A molecular dynamics study of polycrystalline fcc metals at the nanoscale Grain
boundary structure and its influence on plastic deformation, 2001, Materials Science and Engineering
A.]
Its glide plane is (1, 1, 1); the dislocation line at this time is mainly composed of two
segments, aligned approximately along 0 , 1 , 1 and 011 directions. The Burgers
vector being a/6 1 , 2 , 1 ; two possible additional partials could give perfect a/2 011
or a/2 0 , 1 , 1 dislocations. Both segments in the figure are of mixed character. The
dislocation velocity estimated from sections made at different deformation times is 4
A/psec, a tenth of the speed of sound. Due to the orientation of the particular grain
relative to the strain direction, the two slip directions for the perfect dislocations 011
and 0 , 1 , 1 have both a Schmid factor of 0.37, which is the largest possible value in
this grain.
Evidences of stacking faults inside the grains in Ni are observed at a slightly higher
grain size (11 nm) compared to Cu (8 nm), probably due to the higher stacking fault
energy. In Cu, Swygenhoven and al. have observed partial dislocations which glide
on slip systems that are not necessarily those favoured by the Schmid factor.
In Ni and Cu a change in deformation mechanism is observed
At the smallest grain sizes all deformation is accommodated in the grain
boundaries and grain boundary sliding, a process based on mechanical and
NANOMATERIALS
53
thermally activated single atomic jumps, dominates the contribution to
deformation.
At large grain sizes, a combination of sliding and intra-grain dislocation activity
is observed.
4.2.3 Hardness
Normally, one would expect that the introduction of a very high fraction of boundaries,
which is the case of nanostructured polycrystals, will result in a further increase in
hardness since the boundaries are thought to play the role of barriers to dislocation
motion, thus hindering deformation. The aforementioned experimental findings
prompted various researchers to propose models in terms of viscous flow or grain
boundary sliding, different dislocation mechanism taking place for small grain sizes,
or by rule-of-mixtures approaches.
In this following section, we are going to study a model
xxxiii
that obtains the total
hardness involving the grain size, the grain boundary width, the hardness
characteristics of coarse grain counterparts and an extra parameter describing the
arrangement of triple junctions.
In the present formulation the grain boundary phase is considered as having similar
mechanical behaviour to the grain interior phase, rather than being amorphous.
Moreover, the effect of triple grain boundary junctions and the associated
disclinations is considered. Such considerations amount, in effect, into assuming that
the strength of the nanostructured polycrystal does not depend only on their volume
fraction but also scales with the grain dimension in a way similar to the grain interior,
with an enhanced fraction of triple junctions providing a different scenario for the
available acting deformation mechanisms.
The total hardness may be assumed to obey a rule of mixtures relation of the form
H=H
G
(1-f) + H
GB
f (4. 11)
Where H, H
G
and H
GB
are the total, grain interior and grain boundary hardness
respectively, while f denotes the volume fraction of the grain boundaries region.
By assuming a tetrakaidecahedra grain configuration, H=H
G
(1-f) + H
GB
f (4. 11) can
be rewritten in the form
NANOMATERIALS
54
GB G
H
d
d d
H
d
d
H
3
3 3
3
3
) ( ) (
+
= (4. 12)
Where is the grain boundary width and d is the grain size.
The basic assumption of the present model is that a Hall-Petch relationship describes
the hardness of both phases, i.e.
H
G
=H
0G
+ k
G
d
-1/2
H
GB
=H
OGB
+ k
GB
d
-1/2
(4. 13)
Where the coefficients have their usual meaning though now assigned to each region
separately. For the dependence of hardness on the grain size d, a power of -1/2 is
used. It may be noted that H
G
is the hardness contribution of the bulk due to the
existence of the grain boundaries necessary to form dislocation pile-ups, while H
GB
is
an analogous contribution of the boundaries due to the existence of triple junctions
which act as sources and sinks of dislocations and may be thought to possess the
properties of disclination dipoles.
To find k
GB
the Hall-Petch approach (see paragraph 4.2.1.2) is used and we can find
|
|
|
|
\
|
=
0
1
0
1
ln
ln
r
r
r
r
k k
c
G GB
(4. 14)
Where r
1c
is the critical transition grain size. For grain sizes D<D
c
the movement of
the glide dislocation is no longer prohibited by the dislocation line, a situation known
as Orowan bypassing.
Thus, the relation predicting the total hardness of the nanocrystalline materials, when
taking into account the special contribution of the boundaries, is written in its final
form as
2 / 1
0
1
0
1
3
3 3
3
3
ln
ln
) ( ) (
|
|
|
|
|
\
|
|
|
|
|
\
|
+
+ = d
r
r
r
r
d
d d
d
d
k H H
c
G OG
(4. 15)
NANOMATERIALS
55
It is worth noting that equation
2 / 1
0
1
0
1
3
3 3
3
3
ln
ln
) ( ) (
|
|
|
|
|
\
|
|
|
|
|
\
|
+
+ = d
r
r
r
r
d
d d
d
d
k H H
c
G OG
(4. 15) bears some resemblance with a similar expression obtain for the yield
stress by using interphase disclination arguments.
|
|
\
|
H
l G
log
) 1 ( 2
(4. 16)
Where G is the shear modulus, the Poissons ratio, the disclination power,
H
the boundary length and a numerical factor estimated close to 0.1.
The hardness dependence on grain size for metals are given in Figure 4-9 a) and b)
and for intermetallics in Figure 4-9 c) and d).
Figure 4-9 : Hardness of coarse grained and nanocrystalline metals [(a),(b)] and intermetallics
NANOMATERIALS
56
[(c),(d)]. [D.A. Konstantinidis, E.C. Aifantis, On the anomalous hardness of nanocrystalline materials,
1998, Nanostructured Materials]
Further development of the model may be possible through the incorporation of a
different dislocation pile-up mechanism which may be active for very small grain
sizes. Another interesting issue is to study the influence of phase boundaries in
addition to grain boundaries so as to elaborate better the differences between pure
metals and intermetalics.
Finally, it is noted that the elastic properties of the grain boundaries could also be
considered different than those of the bulk due to the presence of the free boundary
volume. This free volume may be viewed as a kind of nanoporosity that could have
an effect on the properties of the grain boundaries and of the nanocrystalline material
as a whole.
4.3 Ceramics / Nanocomposites
4.3.1 Density
Nanocomposites with considerably improved mechanical properties have, to date,
only been achieved in hot-pressed materials. The strength depends on relative
density and decreased with increasing porosity. This might explain why Zhao et
al.
xxxiv
found almost no strength increase in presureless sintered Al
2
O
3
/SiC 5vol%
nanocomposites with a density of 98.3% whereas in hot-pressed materials with a
density of 99.9% considerable increase in strength was observed.
4.3.2 Fracture strength
Figure 4-10 shows a comparison of strength and fracture toughness of alumina-
based nanocomposites at room temperature as a function of SiC content as reported
by several groups. The strength of the monolithic alumina used as a reference varied
from 350 to 560 MPa. Therefore, a given relative improvement in strength of a
ceramic nanocomposite can be misleading. Following Niiharas original work
xxxv
-
, xxxvi , xxxvii , xxxviii , xxxix , xl , xli , xlii
the addition of only 5 vol% nanosize SiC increased the
strength to 1050 MPa. A further increase of SiC content lowers the strength to a
constant value of approximately 800 MPa. However, Niihara
explains the decrease of
strength value for higher SiC contents due to agglomeration problems.
NANOMATERIALS
57
3
1000
800
600
400
5
4
0
10 20 30
Si C content (vol %)
S
t
r
e
n
g
t
h
(
M
P
a
)
T
o
u
g
h
n
e
s
s
(
M
P
a
x
m
1
/
2
)
Figure 4-10 : Strength and toughness of Al
2
O
3
/SiC nanocomposites as a function of SiC
content (adapted from various works): () Niihara & Nakahira
xxxix
measured by three-point bend test
and Vickers indentation; ( ) Borsa et al.
xliii
by four-point bend test; () Zhao et al.
xxxiv
by four-point bend
test and indentation-strength method; (V) Davidge et al.
xliv
58
by three-point bend test and notched
beams
Figure 4-11 shows strength and fracture toughness for various Si
3
N
4
/SiC
nanocomposites as a function of SiC volume fraction. All data were obtained from
different publications of Niihara and his co-workers. Nanocomposites derived from
the classical powder route are represented by solid symbols and nanocomposites
fabricated from an amorphous Si-C-N powder by open symbols. With few exceptions,
all data point follows a defined trend. The toughness increased from approximately
5MPam for mopnolithic Si
3
N
4
to 6 MPam for manocomposites with 10 vol% SiC.
This increase is accompanied by a modest increase in strength from 900 MPa to
1100 MPa.
NANOMATERIALS
58
Figure 4-11 : Strength and fracture toughness for various Si
3
N
4
/SiC nanocomposites as a
function of SiC volume fraction
4.3.3 Strengthening and Toughening Mechanisms
The evidence is that Ai
2
O
3
/SiC nanocomposites show an explicit increase in strength
accompanied by a modest increase in toughness. Furthermore, grain boundaries are
strengthened in nanocomposites as manifested by the transcrystalline fracture mode
as well as by the increased resistance to wear and creep.
4.3.3.1 Critical flaw size reduction (c-mechanism): Zener grain size boundary
pinning
One of the main features of nanocomposites is that matrix becomes refined on
adding nanosize SiC. Following the Hall-Petch relation a refinement of the grain size
leads to higher strength. Furthermost, abnormal grain growth is reduced in
nanocomposites, thus leading to a narrower grain size distribution. The effect of grain
boundary pinning by small inclusions is described by Smith
xlv
107
after a semi
quantitative approach by Zener:
NANOMATERIALS
59
3
4
f
r
R
V
(4. 17)
Where R is the average matrix grain boundary radius of curvature, depending on the
radius (r) and the volume fraction (V
f
) of spherical inclusions. Equation
3
4
f
r
R
V
(4. 17) is
often given in the form
Dr/V
f
(4. 18)
Where D is the average matrix grain diameter.
Figure 4-12 : Experimental data where R is plotted as a function of 1/V
f
.
In Figure 4-12 the Zener model (equation
3
4
f
r
R
V
(4. 17) predicts larger matrix grains for larger SiC particle
sizes.
However, the reduction of the matrix grain size is obvious when taken into
account that monolithic alumina sintered under the same conditions possess a much
larger grain size of about 20 m.
4.3.4 Reduction in processing flaw size
A further explanation for the increased strength of the nanocomposites is a
reduction in the size of processing flaws. Fractographical studies on broken four-point
bend test beams have shown that the strength-determining processing flaws change
in size and in morphology from large-volume pores in alumina to crack-like flaws due
to SiC agglomerations in Al
2
O
3
/SiC nanocomposites.
xlvi
36
The different processing
flaw type results from the specific nanocomposite processing rather than from an
intrinsic nanocomposite effect. In Al
2
O
3
/SiC nanocomposite powder mixtures, hard
SiC agglomerates. The latter commonly cause large processing flaws such as voids
in alumina ceramics.
4.3.5 Crack healing (annealing treatment)
Zhao et al.
xxxiv
55
suggest that SiC particles only indirectly influence the strength by
enabling the compressive stress induced by the grinding process to be retained in the
surface region of the test specimens. Another theory is that cracks in nanocomposites
can heal during annealing. Alumina and Al
2
O
3
/SiC nanocomposites are indented with
NANOMATERIALS
61
a Vickers piramid to generate radial cracks. After annealing at 1300
C in Ar for 2 h
the material behave completely differently.
xlvii
110
Whereas cracks in alumina grow,
cracks in nanocomposites close, thus explaining the strength increase of annealing
nanocomposites.
Possible reasons for the crack healing have not yet found and more systematic work
is needed. However, by taking into account internal stresses, one has to distinguish
between the stresses introduced by the Vickers indentation (comparable to
machining-introduced stresses), stresses introduced by the thermal expansion
mismatch of Al
2
O
3
and SiC and, finally, stresses due to the thermal anisotropy of
alumina grains. Fang et al.
xlvii
110
observed that residual stresses introduced by
Vickers indentation fully relax in alumina after an annealing procedure whereas
compressive stresses are still present in nanocomposites.
4.3.6 Toughening (K-mechanisms)
4.3.6.1 R-curve effects
Concerning R curve effects one has to distinguish between mechanisms acting on
the crack wedge behind the crack tip and mechanism acting directly at or in front of
the crack tip. It is well known that manolithic alumina ceramics exhibits R curve
behaviour due to crack bridging mechanisms. The fracture mode is intercrystalline
with partially connected grains acting as ligaments.
xlviii
111
In nanocomposites any
toughening effects acting on the crack wedge behind the crack tip are unlikely
because of the lack of bridging elements. This is supported by the transgranular
fracture mode. Therefore, only mechanisms acting directly at or in front of the crack
tip can be assumed to be applicable to nanocomposites. Such micro-toughening
mechanisms do not necessarily lead to an increase in the toughness plateau value
but they can result in a steep rise of the R curve for very short crack lengths. This
would explain a higher strength, as schematically shown in Figure 4-13 where the
crack resistance is plotted as a function of crack length. The slope of the tangents on
the R curves represents the strength.
xlviii
111
By assuming the same initial flaw sizes
and plateau toughness values, a higher strength can be achieved for a sharply rising
R-curve.
NANOMATERIALS
62
Figure 4-13 : Schematical dependence of the toughness.
4.3.6.2 Crack Deflection
The interactions of a crack front with second-phase inclusions, such as spherical
particles in nanocomposites, depend on the differences in the thermoelastic
properties of the matrix and inclusions. If there are no differences, the planar crack
front will not be influenced.
xlix
112
Following Faber and Evans,
l
113
in the case of such
differences, the crack front will be deflected from planarity by a single particle or
twisted between two neighbouring particles. The toughening of the composite is a
result of diminishing the stress intensity directly at the tip of the deflected crack. The
extend of toughening increase can be obtained by calculating the local stress
intensities at the crack front. The toughness increase depends on the shape, volume
fraction and interparticle spacing of the reinforcement. As discussed later, crack
deflection can also explain the change in fracture mode.
Niihara has proposed a toughening effect due to crack deflection caused by
compressive residual stresses around the SiC particles.
xxxv
One requirement of the
crack deflection theory is strong interface between the SiC particles and the matrix.
Investigations of SiC/Al
2
O
3
interfaces revealed that the boundary is free of second
phases. 114.28 An alignment of the SiC particles with the Al
2
O
3
matrix can be
observed in some cases.
li
114
These examinations suggest strong adhesion between
the SiC particles and the matrix. Concerning SiC/Al
2
O
3
interfaces at SiC particles
located at Al
2
O
3
/Al
2
O
3
grain boundaries, Jiao et al.
li
114
have estimated that the
NANOMATERIALS
63
interfacial fracture energy between SiC and alumina is twice that of the alumina grain
boundary fracture energy.
4.3.7 Grain boundary strengthening mechanisms
Several authors
lii
35, 55, 59
have found that grain pullout is significantly reduced in
Al
2
O
3
/SiC nanocomposites during polishing, machining or abrasive wear.
Furthermore, the fracture mode changes to transcrystalline. All these observations
indicate that matrix grain boundaries are strengthened in nanocomposites. This effect
is probably the most important difference between nanocomposites and the pure
matrix. Conceivable reasons for the grain boundary strengthening are:
(1) deflection of a crack running along a grain boundary at an SiC particle into the
grain.
(2) strengthening of the grain boundaries due to local internal stresses.
Therefore, internal stresses are discussed in detail in the following sections.
4.3.8 Thermal expansion mismatch (Selsing model)
Since the microstructures of nanocomposite ceramics are formed during sintering at
high temperatures, differences in the thermal expansion coefficient of the matrix
(
matrix
) and of the nano-particles (
particle
) cause strains during cooling. These thermal
expansion misfit strains,, can be calculate by an integration over temperature.
The upper limit is taken as the temperature below which plastic deformation is
insignificant (T
plastic
) and the lower limit is the room temperature.
( )
0
* plastic
T
T particle matrix
dT =
(4. 19)
The thermal expansion misfit stress,
T
, inside a single spherical inclusion in infinite
matrix can be described by the following expression after Selsing:
liii
118
*
1 2
1
2
T
p
m
m p
v
v
E E
+
+
(4. 20)
NANOMATERIALS
64
E and are Youngs modulus and Poissons ratio of the matrix (m) and the particles
(p). The tangential,
T
, and the radial,
T
, stress distributions in the matrix around the
particle are given by:
3 3
2
T
T Tr T
r r
x x
( (
= =
( (
(4. 21)
Where r denotes the radius of the inclusion and x is the radial distance from the
inclusion surface. Assuming T
plastic
as 1500
o
C and the room temperature
thermoelastic data given in Table 4-1, equation (
*
1 2
1
2
T
p
m
m p
v
v
E E
+
+
(4. 20) leads
to a compressive hydrostatic stress inside the SiC particle of 2.0 GPa for the
Al
2
O
3
/SiC and 4.3 GPa for the MgO/SiC system. Tensile hydrostatic stresses of
500MPa and 600MPa can be calculated for the Si
3
N
4
/SiC and the Al
2
O
3
/TiN systems,
respectively. Evidently, a change in fracture mode occurs only in systems with high
compressive stresses within the nano-particles such as Al
2
O
3
/SiC and MgO/SiC.
Table 4-1 : Thermoelastical data for matrix and nanophase
E[Mpa] [10
-6
K
-1
]
Al
2
O
3
400 0.23 8.3
Si
3
N
4
300 0.27 3.2
MgO 300 0.18 14
TiN ~470 ~0.25 9.4
SiC 480 0.17 4.4
4.3.9 Average internal stresses
A model developed by Taya et al.
liv
119
incorporates the influence of the particle
volume fraction on the average residual microstress in the matrix for a composite
material with spherical inclusions showing a thermal expansion mismatch. Their
NANOMATERIALS
65
model is based on Eshelbys equivalent inclusions approach referring to the inelastic
eigenstrains. The average residual microstresses inside the matrix and the particles,
m
and
p
, are given by the following expressions:
*
*
2(1 )
(1 )( 2)(1 ) 3 (1 )
2
(1 )( 2)(1 ) 3 (1 )
f p
m f m f m
f
m
m f m f m
V
E V v V v
V
E V v V v
=
+ + +
=
+ + +
(4. 22)
where
1
.
1 2
p
m
p m
E
v
v E
+
=
In the case of
m
>
p
, the average thermal stresses are compressive inside the
particles and tensile in the matrix. For V
f
=0, equation
(
*
*
2(1 )
(1 )( 2)(1 ) 3 (1 )
2
(1 )( 2)(1 ) 3 (1 )
f p
m f m f m
f
m
m f m f m
V
E V v V v
V
E V v V v
=
+ + +
=
+ + +
(4. 22) provides the same value for the
stresses inside the particles as the Selsing model.
liii
118
Another boundary condition is
given by mechanical equilibrium with
m
(1V
f
)+
p
V
f
=0.
Figure 4-14 shows experimental results together with a prediction using equation
below. A stress-free temperature of 1500
o
C, below which plastic deformation is
insignificant, is assumed. (By considering that internal stresses at the surface of a
nanocomposite are not relieved after annealing at 1300
o
C, as observed by Fang et
al.,
xlvii
110
the stress-free temperature must be somewhere between 1300 and 1550
o
C
which is the minimum hot-pressing temperature to get fully dense nanocomposites).
NANOMATERIALS
66
1000
500
0
-500
-1000
-1500
-2000
-2500
0.0
0.1 0.2
0.3
Volume fraction V
f
A
v
e
r
a
g
e
r
e
s
i
d
u
a
l
s
t
r
e
s
s
[
M
P
a
]
<
m
>
<
p
>
Figure 4-14. Average residual microstresses in the matrix
m
and nano inclusions
p
for Al
2
P
3
/SiC
nanocomposites. ( , ) adapted from Todd et. al; ( , ) Ref. 120)
It has to be emphasised that all theoritical calculations concerning internal stresses
depend strongly on value used for the thermoelastic properties of the matrix and
inclusions.
Levin et al.
lv
123
have presented a model for the influence of the SiC particles on the
fracture toughness of nanocomposites. The average tensile stress field in the matrix
due to the thermal expansion mismatch of Al
3
O
2
and SiC reduces the fracture
toughness. As shown in Figure 4-14, the average microstress in the matrix is
approximately 100 MPa for an Al
2
O
3
/5 vol% SiC nanocomposites. SiC submicron
particles within the grains strengthen the grain boundaries because of compressive
radial stress components, thus increasing the fracture toughness via a change in the
fracture mode. It has been claimed that a net increase in toughness can only be
achieved for small SiC volume fractions because the strengthening effect of grain
boundaries is high and only in the case are the average internal stresses small. The
maximum increase of fracture toughness at 5 vol% SiC as shown by Niihara
xxxv
agrees with the model proposed by Levin et al.
lv
123
and is plausible when taking into
account the two opposite effects of grain boundary strengthening and average tensile
internal stresses.
NANOMATERIALS
67
4.3.10 Local stress distribution
For discussion of the grain boundary strengthening effect due to internal stresses it is
helpful to review the stress fields around the SiC particles. For the model presented in
this review, a superposition of stresses in an arrangement of several boundaries
including a grain boundary is assumed. Figure 4-15 shows the configuration assuming
nine particles in a cubic-body-centred arrangement with an average nearest particle
spacing equivalent to a SiC volume fraction of 2.5%. One particle is located in the
plane. The remaining particles are the same distance from that plane. Each particles
causes a stress field which can be calculated using the Selsing equation. The total
stress at each point in the x-y-plane highlighted in Figure 4-15 can easily be calculated
assuming a simple superposition. Obviously, the particle at the grain boundary
generates high tensile stresses immediately around it and the other eight particles
nearby to the plane generate compressive stresses up to 120 MPa. SiC particles
within the Al
2
O
3
grain boundaries strengthen the grain boundaries due to radial
compressive stresses. The influence of SiC particles within grain boundaries is still
unclear, but could be expected to generate countervailing tensions across the
boundary.
Figure 4-15 : Configuration for the stress distribution model
4.4 Final remarks on strengthening and toughening
mechanisms
Table 4-2 summarises all the mechanisms discussed in the modeling chapter. The
strength increase observed in nanocomposites can be explained by a decrease in
critical flaw size. As shown above, not only in the size of processing flaws decreased
NANOMATERIALS
68
but also their morphology changes completely. Furthermore, the matrix grain size is
reduced with a narrow size distribution due to the grain boundary pinning by inert SiC
particles. Dislocation network plays only a minor role.
A clear identifications of the toughening mechanism in nanocomposites remains
difficult because, first, the toughness increase is small or even absent and, second,
there is no single persuasive mechanism (Table 4-2). A balance exists between
mechanisms decreasing toughness, i.e. crack deflection or grain boundary
strengthening. All these mechanisms depend highly on specific processing and
microstructural details. However, it appears that grain boundary strengthening due to
local radial compressive stress components around SiC particles or due to a riveting
effect is the most likely nanocomposite effect. The mechanism explains the transition
in fracture mode as well as the improved resistance to wear and machining damage.
Table 4-2 : Summary of strengthening and toughening mechanisms
Mechanism Comment
Zener grain boundary
pinning mechanism (c-
mechanism)
Matrix grain sizes are drastically reduced
(typical for nanocomposites)
Reduction in processing
flaw size (c-mechanism)
Strength increase can be fully explained
by observed change in processing flaw
type (careful processing is very important)
R-curve effects (K-
mechanism)
Steep rising R-curve behaviour is
proposed (no experimental evidence)
Crack Deflection (K-
mechanism)
Cracks seem to be reflected at SiC
particles (importance for toughening is
unclear)
Thermal expansion
mismatch (grain
boundary strengthening)
Fracture mode is changed to
transcrystalline if
particle matrix
< (eg. For
Al
2
O
3
/SiC)
Average internal stresses
(grain boundary
Average tensile stresses in matrix if
particle matrix
< (toughness is reduced)
NANOMATERIALS
69
strengthening)
Local stress distribution
(grain boundary
strengthening)
Local compressive stresses in matrix
grain boundaries if
particle matrix
< (can
explain change in fracture mode)
4.5 References
R.A. Masumura, P.M. Hazzledine, C.S. Pande, Yield stress of fine grained materials,
1998, Acta Materialia
E. Arzt, Size effects in materials due to microstuctural and dimensional constraints: A
comparative review, 1998, Acta Meteriala
Chokshi, A. H., Rosen, A., Karch, J. and Gleiter, H., Scripta metall., 1989, 23, 1679
H. Van Swygenhoven, A. Caro, D. Farkas, A molecular dynamics study of
polycrystalline fcc metals at the nanoscale Grain boundary structure and its influence
on plastic deformation, 2001, Materials Science and Engineering A.
H. Van Swygenhoven, M. Spaczer, A. Caro, Microscopic description of plasticity in
computer generated metallic nanophase samples: A comparison between Cu and Ni,
1999, Acta Materialia.
D.A. Konstantinidis, E.C. Aifantis, On the anomalous hardness of nanocrystalline
materials, 1998, Nanostructured Materials
NANOMATERIALS
70
5 THERMAL CONDUCTIVITY IN NANOSTRUCTURED
MATERIAL
5.1 Thermal conductivity of thermal barrier coatings
lvi
5.1.1 Thermal conductivity
The theory of heat conduction by lattice waves in solids can be used to understand
how the conductivity is influenced by lattice defects, grain boundaries and extended
imperfections. At high temperatures one must also consider conduction by
electromagnetic radiation within the solid, which is governed by analogous
considerations. The thermal conductivity by mobile carriers, whether waves or
particles, can be expressed in general in the form
Cvl
3
1
= (5. 23)
where C is their specific heat per unit volume, v their speed and l their mean free
path. If the carriers are waves (lattice waves or electromagnetic waves) ranging over
a spectrum of frequencies f, this must be generalized to
( ) ( )df f vl f C
=
3
1
(5. 24)
where C( f )df is the contribution to the specific heat per unit volume from waves in
the frequency interval df, v is the group velocity of the waves and l(f) their
attenuation length, usually a function of frequency. Also, f
m
denotes an effective
upper limit to the spectrum because either v or l become very small for f\ f
m
. In the
present discussion f
m
is the upper frequency limit of the acoustic branch, as
discussed in the next section. The energy content of waves consist of quanta,
phonons or photons respectively, and these quanta can also be considered as
particles which carry heat. The present treatment emphasizes the wave nature of the
carriers of heat.
NANOMATERIALS
71
5.1.2 Lattice waves
Lattice waves are just elastic or ultrasonic waves, but their spectrum extends to very
high frequencies, where their wavelength is of atomic dimensions, and where the
relation between f and the wavelength is modified. At low frequencies =v/f, but at
high frequencies, due to the discrete atomic structure of the lattice, this relation holds
no longer. The displacement field of the wave has the form of a progressive wave, so
that the displacement u(r,t) varies with position r and time t as
( ) ( ) t i r iq t r u . exp , (5. 25)
where the wave-vector q has magnitude 2/ and points in the propagation direction,
and is the angular frequency = 2 f. At high frequencies is then no longer
proportional to q. The relation between frequency and wave-vector for a discrete
lattice is shown schematically in Figure 5-1.
Figure 5-1 : Dependence of frequency on the wave-vector q for a discrete lattice in the principal
direction. Curve A is the acoustic branch, the only type in simple lattices. Curve O is an optical branch,
present in complex crystals and describing vibration of atoms within molecular units. [P.G. Klemens,
M. Gell, Thermal conductivity of thermal barrier coatings, Materials Science and Engineering A245
(1998), p.p. 143-149]
The slope of the curve defines the group velocity v=d/dq, that is the speed by
which a wave transports energy.
The Debye theory of lattice vibrations makes the simplification that v is constant, but
that the lattice wave spectrum is terminated at the Debye frequency f
D
chosen so that
the number of normal modes agrees with the number of atoms. However, when the
basic structural unit contains N atoms of varying mass, the spectrum divides into
acoustic modes, of the progressive wave form, (curve A in Figure 5-1), and into
optical modes (curve O) vibrations of the N atoms in each structural unit relative to
NANOMATERIALS
72
each other. Only the former, with upper frequency f
m
= f
D
/ N
, have appreciable
group velocity and transport heat. This approximation is used here. In the limit
T > hf
m
/k, which is the case of interest here, each wave or normal mode has, in
thermal equilibrium, energy kT, and C(T) f. Here k and h are the Boltzmann and
Planck constants, respectively.
The difference between acoustic and optic vibrations is shown schematically in
Figure 5-2.
Figure 5-2 : In crystals composed of molecular units, such as ceramic oxides, the lattice waves are of
two types: acoustic and optical. [P.G. Klemens, M. Gell, Thermal conductivity of thermal barrier
coatings, Materials Science and Engineering A245 (1998), p.p. 143-149]
The optic vibrations are vibrations of the atoms in the lattice; the term optic refers to
their strong interaction with infrared electromagnetic waves in the case of ionic solids.
5.1.3 Interaction Processes
In an ideal solid, one that is perfectly harmonic (obeying linear elasticity), structurally
perfect and without external boundaries the lattice waves would be normal modes.
This means that the energy content of each wave would remain constant. Thermal
equilibrium between the waves could not be established. Real solids depart from this
ideal behavior. The lattice forces have anharmonic components, so that there are
departures from linear elasticity, or terms in the strain energy which are cubic in
strain. Also there are always lattice defects and, of course, boundaries. All these
cause some exchange of energy between the waves, which are no longer true
normal modes. The large number of normal modes results in sufficient randomness
in these interactions, so to establish thermal equilibrium. To treat thermal equilibrium
properties it suffices to know that these interactions exist. To study transport
NANOMATERIALS
73
properties and estimate the attenuation lengths one must consider these interactions
in detail.
The thermal conductivity is limited by various interaction processes, which transfer
energy between the waves and establish thermal equilibrium. These are
schematically shown in Figure 5-3 and correspond to the different types of departures
from ideal behavior.
Figure 5-3 : Interaction processes between lattice waves (a) due to anharmonic interactions, (b)
scattering of waves by lattice defects, and (c) reflection of waves by boundaries and interfaces. [P.G.
Klemens, M. Gell, Thermal conductivity of thermal barrier coatings, Materials Science and Engineering
A245 (1998), p.p. 143-149]
The intrinsic processes exchange energy between triplets of waves (three-phonon
interactions). These satisfy frequency conservation, as in all non-linear processes, as
well as wave-vector selection rules. The resulting intrinsic attenuation length is of the
form
( )
1 2
,
= T Bf T f l
i
(5. 26)
where
m
vf a B
3
, being the shear modulus and a
3
the volume per atom. It is thus
inversely proportional to the mean square thermal strain kT/a
3
. The intrinsic
conductivity becomes then
( ) T f v N
D i
/
4
3
2 3 / 2
= (5. 27)
and this describes, reasonably well, the thermal conductivity of structurally perfect
dielectric crystals near and above their Debye temperature hf
D
/k. Weak binding,
large atomic masses and structural complexity (large N) all tend to reduce the
intrinsic conductivity. Note that if one writes
|
|
\
|
|
|
\
|
=
|
|
\
|
|
|
\
|
=
0
0
arctan 1
arctan
f
f
f
f
f
f
f
f
m
m
i P
B
m
m
B
i B
Eq 5.2
Note that the equation for
P
also holds if f
0
> f
m
.
For a substitutional atom of mass M+M instead of M the change in velocity is
M
M
v
v
=
2
1
. Thus one finds
4 4 3 3
2
4
) (
f v a
M
M
c
f l
l
P
\
|
= Eq 5.3
where c is the defect concentration per atom. Effective values of M/M can be
obtained for other defects, including vacancies. If
i
is known, one can deduce l
i
(f, T)
and thus obtain
P
; this has been done for many systems, e.g. for cubic zirconia.
Cubic zirconia presents a further difficulty: it only exists if it contains stabilizing
solutes together with the accompanying vacancies needed for charge neutrality. Thus
the intrinsic thermal conductivity and l
i
must first be estimated in terms of its elastic
and anharmonic parameters. The resulting uncertainty of perhaps 30% in
i
translates into a 15% uncertainty in . There is rough agreement for various solutes
and concentrations at room temperatures. At high temperatures the theoretical
values are always less than the data. This was shown for several different solutes
and concentrations. This suggests the presence of a radiative component. Values of
B
can be calculated for various grain diameters, and the resulting lattice thermal
conductivities are shown for zirconia containing 7wt% yttria in Figure 5-6.
NANOMATERIALS
76
Figure 5-6 : Theoretical reductions in the thermal conductivity of ZrO
2
containing 7 wt% Y
2
O
3
due to
point defects and grain boundary scattering for various grain sizes. [P.G. Klemens, M. Gell, Thermal
conductivity of thermal barrier coatings, Materials Science and Engineering A245 (1998), p.p. 143-149]
The reduction is greatest at low temperatures. At high temperatures any radiative
component present would not depend markedly on grain size.
5.2 High-temperature thermal conductivity of porous Al
2
O
3
nanostructures
lvii
5.2.1 Theory
The two mechanisms important for the heat transport in nonmetallic solids at high (of
the order of Debye or higher) temperatures are:
1. the phonon relates to the movement of atoms in the crystalline lattice
2. radiative mechanisms is due to the electron transitions between the energy
levels in these atoms
If we neglect the relation between the energy levels and spacing between the atoms,
then both mechanisms can be considered independently. This assumption seems to
be reasonable, if the amplitude of the atomic oscillation is considerably lower than the
lattice constant. We can then write
r ph
+ = Eq 5.4
where is the thermal conductivity,
ph
and
r
are its phonon and radiative
components, and consider each of them separately.
NANOMATERIALS
77
5.2.1.1 Phonon contribution
To investigate the high temperature phonon transport in nanostructures we assume
that two mechanisms are dominant:
1. scattering at the grain boundaries
2. the phononphonon interaction.
At the same time, for the materials whose thermal conductivity is considerably
reduced due to nanostructuring the phonon scattering at the grain boundaries is of
most importance. This scattering determines the phonon mean free path that
becomes now comparable with the grain size and is independent of the phonon
wavelength. This means that the phonons whose wavelength is of about the lattice
constant determine the thermal conductivity at high temperatures. This distinguishes
nanostructured materials from single crystals or amorphous ones where the phonons
of different wavelength are of the same importance for TC (thermal conductivity).
We adopt the following model. We start with the phonon states that initially were
localized at the individual grains of the nanostructure and then consider their
transport as the phonon hopping between the grains. This allows us to reduce the
problem of the phonon transport to the problem of electrical conductivity of a net of
arbitrary resistors R
ij
, which correspond to relevant grain boundaries (between the
i-th and j-th grains). This yields the following expression for TC:
( )
/ 1
0
2
2
) (
3
dx x B
d a
S TZt k
B
ph
Eq 5.5
where
D
x
x
T
T
x
e
e x
x B = |
\
|
,
1
) 1 ( 2
9
) (
2
2
4
4
Eq 5.6
T is the temperature, T
D
is the Debye temperature, and k
B
are the Planck and
Boltzmann constants, a is the lattice constant, d is the grain size, S is the main area
of the intergrain boundary (the neck size), and Z is the coordination number.
Deviation of this parameters from its mean value is determined with the parameter of
disorder . It has been estimated as
1
2
2
1
|
|
\
|
+ =
Z
Eq 5.7
NANOMATERIALS
78
is the mean relative deviation of S .
( )
/ 1
0
2
2
) (
3
dx x B
d a
S TZt k
B
ph
Eq 5.5 results in
the temperaBture dependence of T
C
that increases at T<T
D
, and becomes constant
at T>T
D
. Thus, for T>T
D
( )
/ 1
0
2
2
) (
3
dx x B
d a
S TZt k
B
ph
Eq 5.5
yields
( )
d a
S Zt T k
D B
ph
2
2
20
= Eq 5.8
The main difference between
( )
/ 1
0
2
2
) (
3
dx x B
d a
S TZt k
B
ph
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
Eq 5.9
Here
i
=
i
(T) is the thermal conductivity of the single crystal; an influence of
impurities and other bulk mechanisms of phonon scattering can also be included in
this value. It is important that all parameters of
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
=
0
3
d l C
c
V r
Eq 5.10
where c, , and l() are the photon speed, frequency, and mean free path, is the
dielectric permittivity,
( )
d
e T c
d C
T k
V
B
|
|
\
|
=
1
/
3
3 2
Eq 5.11
is the heat capacity of the photons of the frequency .
Let us estimate the mean free path of the thermal photon ~k
B
T in a perfect
Al
2
O
3
single crystal.
( )
* / 2
~
2
2
m e N
c c
l
e
= Eq 5.12
where * /
2
m e N
e
= is the electric conductivity, e and m* are the charge and
effective mass of the electron, is the electron relaxation time; at high temperatures it
can be estimated as = /(k
B
T), so that
T k m
e N
B
e
*
2
= Eq 5.13
N
e
is the electron density. In a perfect crystal it can be estimated as
|
|
\
|
|
|
\
|
=
T k
E
T k m
N
B
g
B
e
2
exp
) 2 (
* 2
2
2 / 3
2
Eq 5.14
where for the band gap and effective mass values we have 7 . 8
g
E eV,
e
m m 35 . 0 *
(e and m
e
are the charge and mass of the free electron). Then
NANOMATERIALS
80
=
0
2
4
2 3 2
3 4
) 1 (
* / 2 3
x
e
B
r
e x
dx x
m e N c
T k
Eq 5.15
Using
T k m
e N
B
e
*
2
= Eq 5.13 and
=
0
2
4
2 3 2
3 4
) 1 (
* / 2 3
x
e
B
r
e x
dx x
m e N c
T k
Eq 5.15 for
T=2000 K yields ) /( 11 mK W
r
, and cm l 14 . This means that a perfect Al
2
O
3
crystal shows no resistance to the photon transport. The resistance could arise from
the impurities or intrinsic defects (the energy level of the defect closest to the top of
the valence band is E
i
=0.3 eV). The latter can be responsible for heat resistance if
their density is rather large, so that l becomes smaller than thickness of the specimen
L. The effect can be estimated by
=
0
2
4
2 3 2
3 4
) 1 (
* / 2 3
x
e
B
r
e x
dx x
m e N c
T k
Eq 5.15,
where N
e
~N
i
exp(E
i
/k
B
T) and m* is effective mass of the hole, N
i
is the defect
density. However, this yields
r
that is usually small, if l<<L.
It is easy to understand the reason. The electron density is low at low temperature.
This means long mean free paths both of electrons and photons [
T k m
e N
B
e
*
2
= Eq
5.13]. Increase of the temperature leads to increase of the photon part of the heat
capacity (which is proportional to T
3
), but also to rapid decrease of the photon mean
free path (which exponentially depends on temperature).
Appreciable contribution to the radiative component could originate from the electron
levels localized less than 0.1 eV from a band edge. They can be the electrons from
impurity or surface levels segregated from the conductive or valence bands
(electrons from dangling bonds). In the latter case, the density of the levels can be
roughly estimated as
D a
p
N
i
2
24
\
|
= + =
c l l l
SC eff
where
1
) (
=
SC
l Eq 5.16
Here,
2 2
) 1 )( 1 ( = p p , and p is the porosity
The effective mean free path ) (
eff
l substitutes ) ( l in ( ) ( )
=
0
3
d l C
c
V r
Eq
5.10. It is apparent this is approximately the smallest value between
2
1
) (
l and
4
) (
=
SC
l . In the perfect crystal ) ( ) (
SC
l l << , and the value of
T k
B
= is
effectively the upper limit for the integral ( ) ( )
=
0
3
d l C
c
V r
Eq 5.10. In the
real crystal the solution
0
of the equation ) ( ) (
0 0
SC
l l = can be either larger or
smaller than
T k
B
. In the first case the disorder has no appreciable influence on
r
, in
NANOMATERIALS
82
the second one the value of
T k
B
<
0
becomes the effective upper limit of
( ) ( )
=
0
3
d l C
c
V r
Eq 5.10; this decreases the radiative component
r.
5.2.2 Experiment
5.2.2.1 Sample preparation
Commercial high purity transition alumina powder (Baikowski CR125; samples
FD1/1FD1/3) as well as MgO doped (0.5 wt%) -alumina powder (Baikowski
CR11325; samples FD2/1FD2/3) were used. The powders were attrition milled for
3h with dispersant addition, freeze-dried and cold isostatically pressed. The samples
were consolidated by spark plasma sintering using a heating rate of 200 C/min up to
the sintering temperatures of 1400 C and holding t imes ranging from 0 to 5 min. The
pressure was applied from the beginning of the sintering cycle and released during
cooling. These different sintering conditions were chosen to produce different
microstructures (grain size, porosity) (Table 5-1).
Table 5-1 : Experimental data used in calculation. [L. Braginsky, V. Shklover, H. Hofmann, P.
Bowen, High-temperature thermal conductivity of porous Al2O3 nanostructures, Physical review B 70,
134201 (2004)]
5.2.2.2 Characterization of the samples
The porosity, grain size as well as the contact area between the grains were
determined using image analysis methods applied to SEM pictures of polished and
thermally etched samples. Additionally TEM studies of the samples FD1/1, FD1/2,
and FD1/3 sintered at different conditions (Figure 5-7) were performed with TEM
NANOMATERIALS
83
Philips CM300, operating at 300 kV. The samples were coated with carbon layer to
avoid charging.
Figure 5-7 : TEM images. (a) Smaller grain size, lower density sample FD1/2, scale bar 200 nm. (b)
Larger grain size, higher density sample FD1/3, scale bar 200 nm. (c) Two sintered grains in powder
FD1/2. The interface (indicated with an arrow) has a thickness, 1 nm. Both grains are crystalline, but
oriented differently and the crystalline fringes can be easily seen only on one of them. Scale bar 2 nm.
[L. Braginsky, V. Shklover, H. Hofmann, P. Bowen, High-temperature thermal conductivity of porous
Al2O3 nanostructures, Physical review B 70, 134201 (2004)]
5.2.2.3 Thermal conductivity measurements
The laser flash method was used for measurements of heat diffusivity in a TC-
3000H/L SINKU-RIKO device as described in Refs
lix,lx
.
5.2.3 Results and discussion
The results of the TEM study of our samples are presented in Figure 5-7, where the
structure of the interfaces between the Al
2
O
3
grains is shown. This is important,
because it influences the phonon propagation in nanopowder. The TEM study shows,
that
1. The interfaces between the aAl2O3 grains are very thin (~1 nm, Figure 5-7)
2. No epitaxial relation between the neighboring crystalline Al
2
O
3
grains was
observed.
Table 5-1 presents results of the structure measurements of our specimens together
with the values of the parameter t for each of them. These values of t ensure the best
fit of
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
= Eq 5.8. The results are close to that of Table 5-1.
Measured values of thermal conductivity as a function of temperature are presented
in Figure 5-9.
NANOMATERIALS
85
Figure 5-9 : (Color online) Temperature dependence of TC of the smaples: measured (dots) and
calculated (solid line). Different types of the dots indicate the different series of the measurements. [L.
Braginsky, V. Shklover, H. Hofmann, P. Bowen, High-temperature thermal conductivity of porous
Al2O3 nanostructures, Physical review B 70, 134201 (2004)]
We see that theoretical curves as a rule demonstrate more rapid decrease with
temperature than observed. This means that the expression
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
=
0
2
4
2 3 2
3 4
) 1 (
* / 2 3
x
e
B
r
e x
dx x
m e N c
T k
Eq 5.15, if N
e
is independent of
temperature. This is of particular importance for specimens of low porosity due to an
increase of the effective photon mean free path [
4
3 2
3
) (
1
) (
1 1
|
\
|
= + =
c l l l
SC eff
where
1
) (
=
SC
l Eq 5.16]. We have observed some increase of TC for the
specimen FD1/3; however, its value is comparable with the experimental errors. The
grain size dependence of TC following from
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
+
=
1
Eq 5.17
Where G
k
is the Kapitza conductance. Indeed,
d G
k
i
i
ph
+
=
1
Eq 5.17 results
from
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
/ 1
0
2 2
2
) (
3
dx x B
d a
S Zt T k
G
D B
k
Eq 5.19
Apparently G
k
is independent of the grain size, if
2
d S ; however, it depends on
temperature dependence of the factor
/ 1
0
) ( dx x B . This allows us to interpret our
NANOMATERIALS
87
phenomenological parameter
k
G t as the transparency of a grain boundary for the
short wavelength phonons.
Strictly speaking, the average value of x resulting from the integrand
(
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
+
=
1
Eq 5.17) could be impracticable. The parameter S determines
the influence of the small size pores. Indeed, this value should be equal to
2
d in the
absence of pores. If we assume a cubic shape of grains and suppose the pores lie
along the cube edges then for the porosity we obtain
(
\
|
|
\
|
=
d
x d
d
x d
p
9
4
1
4
3
2
Eq 5.20
Where x is the neck size (
2
x S = ). For the model of spherical grains such an
estimation has been carried out
lxii
. The main deviation of both estimations is due to
the large size pores, i.e., the pores of about or greater than the grain size. To
estimate their influence, we have to remove some resistors from our net, so that their
resistivity becomes equal to infinity; this means the absence of the neighbours for
some grains due to the large scale porosity. This factor can be related to the
cohesion C of each sample, which can be measured experimentally. Cohesion is an
average ratio of the number of boundaries grain/pore to the number boundaries
grain/grain along some straight line in the TEM pattern. Assuming a homogeneous
distribution of the pores, we can estimate the density of the nonremoved resistors in
the net as
C
C
=
2
Eq 5.21
Then the expression for TC should be multiplied by , if 1. In the case of large
porosity we have to take into account the problem of dead ends. In the more
general case the problem has to be considered with the percolation technique. The
resistivity of the net has an increase (and so the TC of our structure has a decrease)
NANOMATERIALS
88
at 3/(2Z), then the expression for TC should be multiplied by ( 3/(2Z))
1.6
. At
smaller values of the short-wavelength phonons becomes localized, and heat
transport would occur due to long-wavelength phonons only. Such a situation cannot
be realized in porous nanostructures, but it seems to be the case in composite
structures where the nanostructure is embedded into some amorphous media.
Our model
/ 1
0
2 1
) ( 3
) (
dx
S Zt x B T k d a k
S Zt x B
T k
D B i B
i
B ph
\
|
<
d
a
D
make the same contribution to
the TC. This is the case of phonon-phonon interaction when the larger number of the
phonons with larger frequencies compensates the smaller value of their mean free
path. For T>T
D
this leads to
~
i
(a/d). For the smaller temperatures we have to take
into account only the phonons with T k
B
< so that
~
i
(a/d)(T
D
/T). There are two
possibilities when this value might be important: at low temperatures and in the case
of the above-mentioned localization of the short-wavelength phonons.
5.3 References
P.G. Klemens, M. Gell, Thermal conductivity of thermal barrier coatings, Materials
Science and Engineering A245 (1998), p.p. 143-149
L. Braginsky, V. Shklover, H. Hofmann, P. Bowen, High-temperature thermal
conductivity of porous Al2O3 nanostructures, Physical review B 70, 134201 (2004)
NANOMATERIALS
89
6 THERMODYNAMIC
6.1 Nanothermodynamics
With the advancement of techniques of creating and characterizing materials, a huge
of small sizes grains (micrometers), nanosystems, molecular magnets, and atomic
clusters, has been formed and displays a variety of interesting physical and chemical
properties. It is important and timely to develop the new theoretical tools to address
experimental findings. On the other hand, the rapid progress in the synthesis and
processing of materials with the structures at the nanometer size has created a
demand for greater scientific understanding of the thermodynamics on nanoscale.
The issue of application of the thermodynamics on nanoscale has been continuously
attracted, since the nucleation reaction was discovered in the early 1930s. For
instance, a good example is the renowned publication of two books with the title
Thermodynamics of Small Systems by Hill in early 1960s, and recently, the
thermodynamics of small systems is renamed nanothermodynamics.
Traditionally, thermodynamics of large systems composing many particles has been
well developed. Classical thermodynamics describes the most likely macroscopic
behavior of large systems with the change of macroscopic parameters. The really
large systems of astrophysical objects as well as small systems containing a
relatively small number of constituents (at the nanometer scale) are excluded.
Therefore, there is a great deal of interest and activity in the present day to extend
the macroscopic thermodynamics and statistical mechanics to the nanometer scale
consisting of countable particles below the thermodynamic limit due to the recent
developments in nanoscience and nanotechnology. To generalize the
thermodynamics on scale, we need to well understand the unique properties of
nanosystems. We have seen previously, that one of the characteristic features of
nanosystems is their high surface-to-volume ratio. As results of surface effects
becoming increasingly important with decreasing size, and then, the Gibbs free
energy relatively increases for some thermodynamic equilibrium systems. Therefore,
the behavior of such nanoscopic clusters differs significantly from the usual
thermodynamic limit. On the other hand, it is clearly known that when the system size
NANOMATERIALS
90
decreases, one has to consider the fluctuations. Based on the nucleation reactions,
the first considerations are on the temperature fluctuations. The quantitative
measurements of temperature fluctuations were realized by superconducting
magnetometers. Interestingly, it is well explained in the following statement by the US
National Initiative on Nanotechnology that the fluctuations play an important role:
There are also many different types of time scales, ranging from 10
-15
s to several
seconds, so consideration must be given to the fact that the particles are actually
undergoing fluctuations in time and to the fact that there are uneven size
distributions. To provide reliable results, researchers must also consider the relative
accuracy appropriate for the space and time scales that are required; however, the
cost of accuracy can be high. The temporal scale goes linearly in the number of
particles N, the spatial scale goes as O(N log N), yet the accuracy scale can go as
high as N
7
to N! with a significant prefactor. Therefore, these valuable hints motivate
researchers to pursue the thermodynamic description at the nanometer size for the
nucleation of the metastable phase. Up to date, there are two kinds of fundamental
approaches to open out the thermodynamics on nanoscale based on the microscopic
and macroscopic viewpoints, respectively. One would go back to the fundamental
theorem of the macroscopic thermodynamics and establish the new formalism of the
nanothermodynamics by introducing the new function(s) presenting the fluctuations
or the surface effect of nanosystems. Another one could directly modify the equations
of the macroscopic thermodynamics and establish the new model of the
thermodynamics on nanoscale by incorporating the LaplaceYoung or Gibbs
Thomson relation presenting the density fluctuation of nanosystems in the
corresponding thermodynamic expressions. The fundamental outlines of these
approaches will be given in the following section.
6.1.1 Hills Theory
In the early 1960s, Hill addressed the subject of the thermodynamics of small
systems due to his interest in thermodynamics of polymers and macromolecules. In
order to clarify the relationship between the macroscopic thermodynamics and the
nanothermodynamics of Hill, first of all, go back to the fundamental theorem and
recapitulate the foundations of the thermodynamics of macroscopic systems.
NANOMATERIALS
91
In the case of the equilibrium thermodynamics of a macroscopic system, the
fundamental equation for the internal energy, U, in the absence of an external field is
expressed as
U(S, V, N) = TS - PV + N (6. 28)
Where S is the entropy (an extensive state function), and it is a function of the
extensive variables (U, V, N) in one-component system, T the absolute temperature,
P the pressure, V the volume, the chemical potential, and N is the number of
particles. The differential form of Eq. U(S, V, N) = TS - PV + N (6. 28) may be
represented as
dU = SdT + TdS VdP PdV + Nd + dN (6. 29)
On the other hand, the relationships among U, S, V, N, T, P, and can be expressed
as
V S
N
U
,
|
\
|
= (6. 30)
V N
S
U
T
,
|
\
|
= 6. 31)
N S
V
U
P
,
|
\
|
= (6. 32)
Eq. dU = SdT + TdS VdP PdV + Nd + dN (6. 29) will change into the
following form by employing one of the above three equations
V S
N
U
,
|
\
|
= (6.
30)-
N S
V
U
P
,
|
\
|
= (6. 32).
SdT VdP + Nd = 0 (6. 33)
This is the celebrated GibbsDuhem relation, and implies that the changes in the
intensive quantities (, T, P) are not independent. However, the usual choice (T, P) is
made in the literature, defining an equation of the state for the system. In particular,
the GibbsDuhem relation implies that
N
V
P
T
= |
\
|
(6. 34)
NANOMATERIALS
92
and
N
S
T
P
= |
\
|
(6. 35)
It is well known that three other functions, besides the internal energy, U, are very
useful in applications to the specific physical situations.
The enthalpy is
H (S,P,N) = U (S,V,N) + PV (6. 36)
The Helmholtz free energy
F (T,V,N) = U (S,V,N) TS (6. 37)
And the Gibbs free energy
G (T,P,N) = U (S,V,N) TS + PV (6. 38)
Eq. H (S,P,N) = U (S,V,N) + PV (6. 36) in eq. G (T,P,N) = U (S,V,N) TS + PV
(6. 38), you obtain
G (T,P,N) = H (S,P,N) TS (6. 39)
According to the dependence relationships of these functions (U, H, F, G) and their
continuity properties in their appropriate variables, four thermodynamic equations
called the Maxwell relations can be yielded as
V S
S
P
V
T
|
\
|
= |
\
|
(6. 40)
P S
S
V
P
T
|
\
|
= |
\
|
(6. 41)
V T
T
P
V
S
|
\
|
= |
\
|
(6. 42)
and
P T
T
V
P
S
|
\
|
= |
\
|
(6. 43)
If one combines eq. dU = SdT + TdS VdP PdV + Nd + dN (6. 29) with SdT
VdP + Nd = 0 (6. 33) we have the following Gibbs generalized relation
NANOMATERIALS
93
+ =
i
i i
dN PdV TdS dU (6. 44)
Equation
+ =
i
i i
dN PdV TdS dU (6. 44) applies to a microscopic system only.
For example, if the system is a small one-component spherical aggregate that has a
nonnegligible surface free energy term proportional to N
2/3
, equation
+ =
i
i i
dN PdV TdS dU (6. 44) is no longer applicable. This equation can be
generalized further to include small systems, by adding another term. But this
addition must be made at the ensemble level rather than at the single-system level.
The reason for this is that an ensemble of B equivalent and noninteracting small
systems is itself a macroscopic system, and hence we can use equation
+ =
i
i i
dN PdV TdS dU (6. 44), rewritten for an ensemble, as a reliable
starting point:
+ =
i
it i t t t
dN PdV TdS dU
U
t
BU
S
t
BS
V
t
BU
N
it
BN
i
(6. 45)
where t = total (i.e., the whole ensemble of small systems).
The term we add to eq N
it
BN
i
(6. 45) is EdB :
EdB dN PdV TdS dU
i
it i t t t
+ + =
(6. 46)
Where E( U
t
/ B)
St,Vt,Nit
. The motivation is that, for a small system, U
t
will also be
a function of B. E is a kind of system (rather than molecule) chemical potential, called
the sub-division potential. If the systems of the ensemble are macroscopic and are
subdivided in order to increase B in eq EdB dN PdV TdS dU
i
it i t t t
+ + =
(6.
46) (e.g., in an extreme case, each system is cut in half to double B ), this will not
have a noticeable effect on U
t
, S
t
, V
t
, N
it
are all held constant as in E( U
t
/ B).
That is, surface effects, edge effects, system rotation and translation, etc., are all
negligible for macroscopic systems, and hence E is essentially zero: the term E d B
NANOMATERIALS
94
in eq EdB dN PdV TdS dU
i
it i t t t
+ + =
(6. 46) does not contribute appreciably
to the equation. But the effects just mentioned are not negligible if an ensemble of
small systems is subdivided to increase B in eq EdB dN PdV TdS dU
i
it i t t t
+ + =
+ =
i
i i
dN PdV TdS dU (6. 49)
This is the same as eq
+ =
i
i i
dN PdV TdS dU (6. 44) for a macrosystem: a
small system has the same energy-heat-work relation as in a macrosystem (with the
N
i
constant) and also the same chemical potential terms. However, eq
+ =
i
i i
dN PdV TdS dU (6. 49) cannot be integrated to give eq
E N PV TS U
i
it i
+ + =
(6. 48) because U is no longer (for a small system) a linear
homogeneous function of S, V, and N
i
.
NANOMATERIALS
95
If we take differentials in eq E N PV TS U
i
it i
+ + =
(6. 48) and cancel the
terms in eq
+ =
i
i i
dN PdV TdS dU (6. 49), the remaining relation is
+ =
i
i i
N d VdP SdT dE (6. 50)
The left-hand side is zero in macroscopic thermodynamics. Note that T, P, and the
i
are all independent in eq
+ =
i
i i
N d VdP SdT dE (6. 50). A macrosystem has one
less degree of freedom because dE = 0.
Equation
+ =
i
i i
N d VdP SdT dE (6. 50) confirms that E is a function of
i
, T,
P. Note that derivatives of E, a composite intensive property, with respect to
intensive variables yield extensive variables:
i
i
i
p T
i
i
T
p
E
N
P
E
V
T
E
S
, ,
,
,
|
|
\
|
=
|
\
|
=
|
\
|
=
(6. 51)
This is a rather surprising result. Equations E N PV TS U
i
it i
+ + =
(6. 48),
+ =
i
i i
dN PdV TdS dU (6. 49) and
+ =
i
i i
N d VdP SdT dE (6. 50)
provide the fundamentals of the generalization of the thermodynamics of Gibbs to
include small systems. Although conventional macrothermodynamics is by far the
more widely known subject, it can be viewed as a limiting case of
nanothermodynamics. This is somewhat reminiscent of classical mechanics (large)
as a limiting case of quantum mechanics (small). In contrast to
macrothermodynamics, the thermodynamic properties of a small system will usually
be different in different environments. An example will clarify this statement.
Consider a rather rigid incompressible linear aggregate of N more or less spherical
molecules in an inert constant-temperature bath at T. N might be on the order of 10
2
or 10
3
. The environmental variables here are N and T (P can be ignored because of
the incompressibility). The aggregate will have a certain amount of entropy, S,
NANOMATERIALS
96
associated with translation, rotation, and the vibrational motion of and within the
molecules of the chain. Now suppose, in contrast, that the bath includes the
individual molecules at chemical potential and that these molecules can go on and
off the two ends of the aggregate. At equilibrium, there will be a mean number N of
molecules in the aggregate. The environmental variables here are and T. If is
chosen so that N in the , T case is the same as N in the N, T case above, the
entropy S will be larger in the , T case because of fluctuations in N that are not
present in the N, T case. However, this extra entropy becomes negligible as the
aggregate approaches macroscopic size (by adjusting ).
To study such effects using theoretical models, the applicable set or sets of
environmental variables must play a primary role. For given environmental variables,
the particular statistical mechanical partition function appropriate for these variables
must be chosen in order to deduce the thermodynamic properties associated with
these environmental variables. This is in contrast with the treatment of macroscopic
systems: the partition function used is optional and can be selected merely for
convenience.
It is interesting that the basic equations of nanothermodynamics are the same for all
sets of environmental variables but the thermodynamic properties are often not the
same.
It is worth noting that Chamberlin et al. have extended Hills idea by considering the
independent thermal fluctuations inside bulk materials. In detail, they adapt Hills
theory to obtain a mean-field model for the energies and size distribution of clusters
in condensed matter. Importantly, the model provides a common physical basis for
many empirical properties, including non-Debye relaxation, non-Arrhenius activation,
and non-classical critical scaling
6.1.2 Tsallis generalization of ordinary Boltzmann-Gibbs
thermostatics
The thermodynamic theory is on the basis of the Tsallis generalization of the
ordinary BoltzmannGibbs thermostatistics by relaxing the additive properties of the
thermodynamic quantities (the entropy, in particular) to include non-extensivity of
nanosystems. As described by Rajagopal et al. [], the nanothermodynamics differs
NANOMATERIALS
97
from Hills approach by considering that each of the nanosystems fluctuates around
the temperature of the reservoir, while nanosystems are coupled to the reservoir.
This means that the BoltzmannGibbs distribution has to be averaged over the
temperature fluctuations induced by the reservoir.
The formalism of nonextensive statistical mechanics has been developed over the
past 13 years as a powerful and beautiful generalization of ordinary statistical
mechanics. It is based on the extremization of the Tsallis entropies:
|
\
|
=
i
q
i q
p
q
S 1
1
1
(6. 52)
subject to suitable constraints. Here the p
i
are the probabilities of the physical
microstates, and q is the entropic index. The Tsallis entropies reduce to the
Boltzmann Gibbs (or Shannon) entropy
i
i
i
p p S
= log
1
for q 1, hence ordinary
statistical mechanics is contained as a special case in the generalized formalism.
There is growing experimental evidence that q 1 yields a correct description of
many complex physical phenomena, including hydrodynamic turbulence, scattering
processes in particle physics, and self-gravitating systems in astrophysics, to name
just a few. The Tsallis entropies are nonextensive. Given two independent
subsystems and with probabilities
i
p and
i
p , respectively, the entropy of the
composed system + (with probabilities
+
=
j i ij
p p p ) satisfies
+
+ =
q q q q q
S S q S S S ) 1 ( (6. 53)
Hence there is additivity for q = 1 only. This property has occasionally lead to
unjustified prejudices, using arguments of the type entropy must be extensive. In
the following, however, we will see that the nonadditivity property is not at all a bad
property, but rather a natural and physically consistent property for those types of
systems where nonextensive statistical mechanics is expected to work.
Broadly speaking, the formalism with q 1 has so far been observed to be relevant
for two different classes of systems. One class are systems with long-range
interactions, the other one are systems with fluctuations of temperature or of energy
dissipation. For the first class of systems, the concept of independent subsystems
does not make any sense, since all subsystems are interacting. Hence there is no
NANOMATERIALS
98
contradiction with the nonadditivity of entropy for independent subsystems since
independent subsystems do no exist for these systems (by definition of the long-
range interaction).
Let us thus concentrate on the other class, systems with fluctuations. Consider a
system of ordinary statistical mechanics with Hamiltonian H. Tsallis statistics with
q > 1 can arise from this ordinary Hamiltonian if one assumes that the temperature
1
is locally fluctuating. From the integral representation of the gamma-function one
can easily derive the formula
( ) ( )
= +
0
1
1
0
) ( 1 1
d f e H q
H
q (6. 54)
where
( ) ( )
(
)
`
|
|
\
|
0
1
1
1
1
1
0
1
exp
1
1
1
1
1
) (
q q
q
f
q
q
(6. 55)
is the probability density of the
2
distribution. The above formula is valid for arbitrary
Hamiltonians H and thus of great significance. The left-hand side of eq.
( ) ( )
= +
0
1
1
0
) ( 1 1
d f e H q
H
q (6. 54) is just the generalized Boltzmann
factor emerging out of nonextensive statistical mechanics. It can directly be obtained
by extremizing S
q
. The right-hand side is a weighted average over Boltzmann factors
of ordinary statistical mechanics. In other words, if we consider a nonequilibrim
system (formally described by a fluctuating ), then the generalized distribution
functions of nonextensive statistical mechanics are a consequence of integrating over
all possible fluctuating inverse temperatures , provided is
2
distributed.
Wilk et al. [] have also suggested a formula similar to eq.
( ) ( )
= +
0
1
1
0
) ( 1 1
d f e H q
H
q (6. 54) which is valid for
q < 1. However, in that case the distribution function f() is more complicated and
different from the simple form
( ) ( )
(
)
`
|
|
\
|
0
1
1
1
1
1
0
1
exp
1
1
1
1
1
) (
q q
q
f
q
q
NANOMATERIALS
99
(6. 55). In particular, for q < 1 one has f() = f(,H), i.e. the probability density
of the fluctuating temperature also depends on the Hamiltonian H. This is physically
less plausible, because usually one assumes that the heat bath is independent of the
Hamiltonian. Hence in the following we will restrict ourselves to the physically more
interesting case q > 1.
For q > 1, the
2
distribution obtained in eq.
( ) ( )
(
)
`
|
|
\
|
0
1
1
1
1
1
0
1
exp
1
1
1
1
1
) (
q q
q
f
q
q
(6. 55) is a universal distribution
that occurs in many very common circumstances. For example, it arises if is the
sum of squares of n Gaussian random variables, with n = 2/(q1). Hence one
expects Tsallis statistics with q > 1 to be relevant in many applications. For fully
developed turbulent hydrodynamic flows, where Tsallis statistics has been observed
to work very well,
1
is not the physical temperature of the flow but a formal
temperature related to the fluctuating energy dissipation rate times a time scale of the
order of the Kolmogorov time. In the application to scattering processes in collider
experiments,
1
is a fluctuating inverse temperature near the Hagedorn phase
transition.
The constant
0
in eq.
( ) ( )
(
)
`
|
|
\
|
0
1
1
1
1
1
0
1
exp
1
1
1
1
1
) (
q q
q
f
q
q
(6.
55) is the average of the fluctuating ,
= =
0
0
) ( ) ( d f E (6. 56)
(E denotes the expectation with respect to f()). The deviation of q from 1 can be
related to the relative variance of . One has
2
2 2
) (
) ( ) (
1
E
E E
q
= (6. 57)
We can now give physical meaning to eq.
+
+ =
q q q q q
S S q S S S ) 1 ( (6. 53) for q
> 1. Consider two independent subsystems and that are composed into one
system + . In system (as well as ) the fluctuations of temperature T are
NANOMATERIALS
100
expected to be larger than in + , due to the smaller size of (or ) as compared
to + . Remember that generally entropy is a measure of missing information on
the system. The entropy
+
q q
S S is larger than
+
q
S , since in the single systems the
probability distribution of T is broader, thus our missing information on these systems
is larger. Hence there must be a negative correction term to
+
q q
S S in eq.
+
+ =
q q q q q
S S q S S S ) 1 ( (6. 53) . It is physically plausible that this correction
term is proportional
to the relative strength of temperature fluctuations as given by q 1 and
to the entropies in the single
q
S and
q
S systems. Hence we end up with eq.
+
+ =
q q q q q
S S q S S S ) 1 ( (6. 53) .
For most physical applications the fluctuations of (or the relevant value of q) are
observed to be dependent on the spatial scale r. For example, in the turbulence
application detailed measurements of q(r) have been presented in []. q turns out to
be a strictly monotonously decreasing function of the distance r on which the velocity
differences are measured. Similarly, for the application to e
+
e
= (6. 59)
Further, for the spherical critical nucleus, Gibbs showed that with eq.
*
2
r
P P
T
v l
=
(6. 59), eq. ) ( *
V l T
P P V A G = (6. 58) becomes
NANOMATERIALS
103
2
3
) ( 3
16
*
v l
T
P P
G
=
(6. 60)
However, the quantity
T
could not be obtained by experiment, because it describes
such a surfacean imaginary physical object, i.e., the nucleus characterized by the
surface of tension. Therefore, the dependent relationships of
T
and pressure,
temperature, and composition of the parent phase, respectively, are not uncovered.
This limits the application of the nucleation theory to various cases of interest. On the
other hand, in order to describe the thermodynamic characterization of various
practical cases, we have to approximate
T
by a real physical quantity. Clearly, in
nearly all nucleation papers that followed Gibbs equation, e.g., one used the real
interface energy
0
between the bulk parent and new phases at the phase
equilibrium, i.e., at their coexistence, to replace the surface tension
T
of the
imaginary DS. However, to apply this famous formula of Gibbs, one has to know the
exact interface energy that is related to the radius of droplets and the droplet
reference pressure. Unfortunately, there are tremendous challenges because there is
no simple way to extract the interface energy from the force measurements in theory.
A given interface energy is a function of the many coordinates of the nanoparticles.
Lacking the knowledge of the exact interface energy, the first approximation is to use
the experimental interface energy of a flat interface, i.e.,
0
= . Actually, the surface
structure of droplets is different from that in the bulk, so that, strictly speaking, the
boundary surface never coincides with the equimolecular surface. Nevertheless, they
are usually close to each other and are often taken equally as the physical surface of
the droplet. By assuming
0
= , one can obtain the first form of the nucleation work
2
3
) ( 3
16
*
v l
P P
G
(6. 61)
In principle, by adopting the approximation, the validity of the Gibbss expression on
the basis of the LaplaceYoung equation should be limited to the sufficiently large
nuclei. Interestingly, the applications of the nucleation theorem in the analysis of
experimental data in various cases of nucleation implied that the LaplaceYoung
equation could predict well the size of the nuclei built up of less than a few tens of
molecules.
For instance, Tolbert and Alivisatos discussed the elevation of pressure in the solid
solid structural transformation as the crystallite size decreases in the high-pressure
NANOMATERIALS
104
system using the LaplaceYoung equation. Accordingly, it seems to be recognized
that the LaplaceYoung equation could be used to predict well the size of nuclei built
up of less than a few tens of molecules. On the other hand, based on the
thermodynamic identity, we have
=
l
v
P
P
m v l l l
dP V P P ) ( ) ( (6. 62)
Where V
m
is the molar volume of a new phase and ) (
l l
P and ) (
v l
P are the chemical
potential of matter in the new phase at the pressures P
l
and P
v
. When the critical
droplet and the metastable vapour locate the condition of the unstable equilibrium,
one can obtain
) ( ) (
l l v v
P P = (6. 63)
Furthermore, if we approximate P
l
by assuming that the droplet is incompressible,
and assume that V
m
is a constant. With Eqs.
=
l
v
P
P
m v l l l
dP V P P ) ( ) ( (6. 62) and
) ( ) (
l l v v
P P = (6. 63) becomes
m m
l l v v
v l
V V
P P
P P
=
=
) ( ) (
(6. 64)
Eq.
m m
l l v v
v l
V V
P P
P P
=
=
) ( ) (
(6. 64) turns into Eq.
2
3
) ( 3
16
*
v l
P P
G
(6. 61), one can obtain the second form of the nucleation work
2
2 3
) ( 3
16
*
=
m
V
G (6. 65)
As Obeidat et al. stated [], the form of the nucleation work is most useful if the
chemical potential difference between a new phase and its parent phase can be
obtained. However, the actual performance is quite complicated due to the lack of the
accurate potentials for most substances. In order to obtain the , we have to adopt
some necessary approximations. Generally, if we assume the supersaturated and
saturated vapors are the ideal gases and the droplet is an incompressible liquid, the
difference of the chemical potential is more commonly derived from the
approximate system. In detail, under the above assumptions, we have
NANOMATERIALS
105
ev l
P P = (6. 66)
Where P
ev
is the equilibrium-vapor pressure. With eq.
ev l
P P = (6. 66),
eq.
=
l
v
P
P
m v l l l
dP V P P ) ( ) ( (6. 62) becomes
) ( ) ( ) (
ev v m ev v v l
P P V P P + = (6. 67)
When the new bulk phase and the parent phase are in the state of thermodynamic
equilibrium, one can obtain
) ( ) (
ev v ev l
P P = (6. 68)
With eqs. ) ( ) ( ) (
ev v m ev v v l
P P V P P + = (6. 67) and ) ( ) (
ev v ev l
P P = (6. 68), eq.
m m
l l v v
v l
V V
P P
P P
=
=
) ( ) (
(6. 64) becomes
) ( ) ( ) (
ev v m ev v v v
P P V P P = (6. 69)
Under the ideal vapor condition, we can easily obtain
|
|
\
|
=
ev
v
ev v v v
P
P
kT P P ln ) ( ) ( (6. 70)
Where k is the Boltzmann constant, T the absolute temperature, and P
v
is the actual
pressure. With Eq.
|
|
\
|
=
ev
v
ev v v v
P
P
kT P P ln ) ( ) ( (6. 70),
Eq. ) ( ) ( ) (
ev v m ev v v v
P P V P P = (6. 69) becomes
) ( ln
ev v m
ev
v
P P V
P
P
kT
|
|
\
|
= (6. 71)
In Eq. ) ( ln
ev v m
ev
v
P P V
P
P
kT
|
|
\
|
= (6. 71), compared with the first term on the
right, the second term on the right is almost extremely small, and it is customary to
neglect it. Therefore, Eq.
2
2 3
) ( 3
16
*
=
m
V
G (6. 65) will become the third form of
the nucleation work
NANOMATERIALS
106
|
|
\
|
=
ev
v
m
P
P
kT
V
G
ln
3
16
*
2 3
(6. 72)
Applying the first two forms of the nucleation work requires the knowledge of the
droplet reference pressure or chemical potential. Usually, this information is
unavailable, and the experimental results are, instead, compared with the rates
predicted using the third form, because the supersaturation ratio is readily
determined from the experimental data. Naturally, the size of the critical nucleation,
the critical energy, the phase transition probability, and the nucleation rate would be
obtained by the determined nucleation work.
In summary, from the point of CNT above, one can see that there is an important
approximation, i.e., assuring
T
=
0
= . Namely, the surface tension (
T
) of a
specially chosen DS, the real interface energy (
0
) between the bulk parent and new
phases at the phase equilibrium, and the experimental interface energy of a flat
interface () are approximated to equal. Furthermore, CNT indicates that the
LaplaceYoung equation seems to be capable of predicting well the size of nuclei
built up of less than a few tens of molecules. However, on the other hand, it is well
known that the CNT describes that a stable new phase forms from a metastable
parent phase.
6.1.3.2 Thermodynamics of metastable phase nucleation in unstable of
thermodynamic equilibrium phase diagram
MPNUR, such as CVD (Chemical Vapor Deposition) diamond, CVD c-BN, HSRC
(Hydrothermal Synthesis and the Reduction of Carbide) diamond, and c-BN (cubic
boron nitride) nucleation, seems to be impossible from the standpoint of
thermodynamics, because the nucleation happens in the strongly unstable region of
the metastable structural state in terms of the thermodynamic equilibrium phase
diagram, and violates the second law of thermodynamics. For this issue, as Hwang
and Yoon stated regarding CVD diamond [], Something must be wrong either in
interpreting the experimental observation or in applying thermodynamics. In fact, in
the early 1965s, Garvie pointed out that MPNUR likely arises out of the capillary
pressure built up in the nuclei. Namely, the nanosize-induced additional pressure
could be so large that the high-pressure metastable phase tends to become more
stable than the low-pressure stable phase, as shown in Figure 6-1.
NANOMATERIALS
107
Figure 6-1 : A sketch map of MPNUR mechanism. A region is the metastable structural state of M
phase, and B region is the new stable state of M phase by the nanosize-induced additional pressure
driving. The inset shown in the sketch map displays the spherical nuclei nucleated on the hetero-
substrate. [C.X. Wang, G.W. Yang, Thermodynamics of metastable phase nucleation, Materials
Science and Engineering, 2005, p.p. 157-202]
Note that, in the following description, a phase is metastable or stable if it is stable or
metastable without the effect of the nanosize-induced additional pressure. In our
theoretical approach, it also is emphasized that the nanosize-induced additional
pressure is reasonably taken into account in the below analysis.
Generally, the Gibbs free energy is an adaptable measure of the energy of a state in
phase transformations among competing phases. At the given thermodynamic
condition, both stable and metastable phases can coexist, but only one of the two
phases is stable, with the minimal free energy, and the other must be metastable and
may transform into the stable state. Thermodynamically, the phase transformation is
promoted by the difference of the free energies. The Gibbs free energy of a phase
can be expressed as a function of the pressuretemperature condition, and
determined by a general coordinate or reactive coordinate. According to CNT, the
Gibbs free energy difference arises from the formation of spherical clusters in the
low-pressure gas phase is expressed as a function of radius r, pressure P, and
temperature T
( )
se se sn sn ne ne
m
S
A A A g
V
V
T P r G + + = ) , , ( (6. 73)
NANOMATERIALS
108
Where V
s
and V
m
are the volume of the spherical clusters with the metastable
structural phase and its molar volume, g the Gibbs free energy of molar volume
depending on the pressure P and temperature T in the phase transition, A
ne
and
ne
the interface area and the energy between the spherical clusters of the metastable
structural phase and the environment gas phase, A
sn
and
sn
the interface area and
the energy between the spherical clusters of the metastable phase and the hetero-
substrate, and A
se
and
se
are the interface area and the energy between the hetero-
substrate and the environment gas phase. The formation of the spherical clusters
with the metastable structural phase produces two interfaces, i.e., the interface A
ne
between the spherical clusters and the environment gas phase and the interface A
sn
between the spherical clusters and the hetero-substrate, and makes the original
interface A
se
(be equal to A
sn
) between the hetero-substrate and the environment
gas phase vanish. According to the geometry, the volume V
s
of the spherical clusters
of the metastable structural state, the interface area A
ne
between the spherical
clusters of the metastable state and the environment gas phase, and the interface
area A
sn
between the spherical clusters and the hetero-substrate are expressed as
( )( )
3
1 2
2 3
m m r
V
s
+
=
(6. 74)
( ) m r A
ne
= 1 2
2
(6. 75)
( )
2 2
1 m r A
sn
= (6. 76)
Where r is the curvature radius of spherical clusters of the metastable structural
phase and m is given by
ne
sn se
m
= = cos (6. 77)
Where is the contact angle between the spherical clusters of the metastable
structural state and the hetero-substrate, as shown in the inset of Figure 6-1. Here,
ne
is assumed to be approximately equal to the surface tension value of the
metastable structural phase (),
se
for the interface energy between the hetero-
substrate and the environment gas phase is taken to be equal to the surface tension
value of the hetero-substrate, and the interface between the spherical clusters and
the hetero-substrate is assumed to be incoherent interface; therefore
NANOMATERIALS
109
ne
se ne
sn
+
= (6. 78)
Thus, we can obtain
( )
( )( )
4
1 2
4
3
4
, , ,
2
2 3
m m
r
V
g
r T P T r G
m
+
|
|
\
|
+
= (6. 79)
Where the factor of
( )( )
4
1 2
) (
2
m m
f
+
= (6. 80)
is called as the heterogeneous factor, and its value is in the range of 01. Especially,
when the clusters nucleated on the homo-substrate, its value is 1.
According to thermodynamics, we have
( )
V
P
P T g
T
= |
\
|
,
(6. 81)
Then the difference of the Gibbs free energy per mole can be defined by
( ) ( )
= =
P
P
O
P V P P V VdP P T g P T g
0
) ( , ,
0
(6. 82)
Where V is the mole volume difference between the metastable and the stable
phase. When the conditions are near the equilibrium line, one can approximately
have g(T, P
0
) = 0. Thus, Eq. ( ) ( )
= =
P
P
O
P V P P V VdP P T g P T g
0
) ( , ,
0
(6. 82) would be defined as
( ) P V P T g = , (6. 83)
On the other hand, due to the nanosize-induced additional pressure P
n
, the clusters
enduring pressure will increase by the same amount. Under the assumptions of
spherical and isotropic clusters, the nanosize-induced additional pressure is denoted
by the LaplaceYoung equation, i.e.
r
P
n
2
= (6. 84)
Furthermore, as mentioned above, the nanosize-induced additional pressure can
drive the metastable phase regions into the stable phase region near the boundary
line of the high-pressure phase in the equilibrium phase diagram. Therefore, one can
NANOMATERIALS
110
obtain the size-dependent equilibrium phase boundary line between the metastable
and the stable phases, and it can be approximately defined as
r
P P
b
2
= (6. 85)
Where P
b
is the equilibrium phase boundary equation between the metastable and
the stable phases. From Figure 6-1, one can see that the equilibrium phase boundary
between the metastable and the stable phases can be expressed by
0 0
b T k P
b
+ = (6. 86)
Where k
0
and b
0
are the slope and intercept in the P coordinate axis of the
equilibrium phase boundary line between the metastable and the stable phases. With
Eq.
0 0
b T k P
b
+ = (6. 86), Eq.
r
P P
b
2
= (6. 85) can be defined as
r
b T k P
2
0 0
+ = (6. 87)
Therefore, P would change into
r
b T k P P
2
0 0
+ = (6. 88)
With Eq.
r
b T k P P
2
0 0
+ = (6. 88), Eq. ( ) P V P T g = , (6. 83) can be
denoted by
( )
|
\
|
+ =
r
b T k P V P T g
2
,
0 0
(6. 89)
With Eq.
r
b T k P P
2
0 0
+ = (6. 88), Eq.
( )
( )( )
4
1 2
4
3
4
, , ,
2
2 3
m m
r
V
g
r T P T r G
m
+
|
|
\
|
+
f r
V
r
b T k P V
r T P r G
m
|
|
|
|
\
|
+
|
\
|
+
= (6. 90)
When 0
) (
=
r
r G
, the critical size of the high-pressure phase is obtained as
NANOMATERIALS
111
P b T k
V
V
r
m
+
|
\
|
+
=
0 0
*
3
2
2
(6. 91)
Substituing Eq.
P b T k
V
V
r
m
+
|
\
|
+
=
0 0
*
3
2
2
(6. 91) into
Eq. ( ) ) ( 4
2
3
4
, ,
2
0 0
3
f r
V
r
b T k P V
r T P r G
m
|
|
|
|
\
|
+
|
\
|
+
= (6. 90), the critical
energy of the high-pressure phase nuclei is given by
( ) ) (
3
2
2
4
3
2
3
2
2
3
4
, ,
2
0 0
0 0
0 0
3
0 0
*
f
P b T k
V
V
V
V
P b T k
b T k P
V
V
P b T k
V
V
T P r G
m
m
m
m
|
|
|
|
|
\
|
|
|
|
|
|
\
|
+
|
\
|
+
+
|
|
|
|
\
|
+
+
+
|
|
|
|
|
\
|
+
|
\
|
+
=
6. 92)
On the other hand, it is well known that the phase transition is determined by the
probability. We have studied the nanosize effect on the probability of the phase
transformation based on the thermodynamic equilibrium phase diagram. The
probability of the phase transformation from the metastable phase to the stable
phase is related not only to the Gibbs free energy difference g(T, P), but also to an
activation energy (E
a
- g(T, P)), which is necessary for the phase transition, as
shown in Figure 6-2.
Figure 6-2 : The schematic diagram of Gibbs free energy vs. coordinate. [C.X. Wang, G.W. Yang,
Thermodynamics of metastable phase nucleation, Materials Science and Engineering, 2005, p.p. 157-
202]
NANOMATERIALS
112
When the two phases are at the equilibrium condition, i.e.,
g (T, P) = 0, E
a
is the maximum potential energy for both sides with respect to the
general coordinate. The general expression of the probability f of the phase
transformation from the initial states to final states is
|
\
|
|
\
|
=
RT
E
RT
P T g E
f
a a
exp
) , (
exp (6. 93)
Where R is the gas constant and g(T, P) is defined by Eq.
( )
|
\
|
+ =
r
b T k P V P T g
2
,
0 0
(6. 89). Accordingly, we have established the
thermodynamic approach at the nanometer size to quantitatively describe the
nucleation and the phase transition of the metastable phase in the strongly unstable
phase region of the metastable structural state in the thermodynamic equilibrium
phase diagram. In fact, the developed approach is a useful and effective theoretical
tool to address MPNUR, although it looks a little bit simple in thermodynamic.
Importantly, the validity of our thermodynamic theory has been substantively checked
by use in the nucleation of diamond and c-BN.
6.1.4 Nanothermodynamic analyses of CVD diamond nucleation
Generally, CVD (Chemical Vapor Deposition) diamond is usually a typical quasi-
equilibrium process, and the pressure is in the range of 102105 Pa and temperature
is in the range of 10001300 K. In the carbon phase diagram shown in Figure 6-3,
the general thermodynamic region of the diamond nucleation upon CVD is shown as
G region, which belongs to the strongly unstable or metastable region of the diamond
structural state, i.e., the stable region of graphite structural state.
NANOMATERIALS
113
Figure 6-3 : Carbon thermodynamic equilibrium phase diagram. The G region means a
metastable phase region of the diamond nucleation upon CVD; the D region means a new stable
phase region of the CVD diamond nucleation with respect to the effect of the nanosize-induced
additional pressure. The inset shows the relationship between the nanosize-induced additional
pressure and the nuclei size. [C.X. Wang, G.W. Yang, Thermodynamics of metastable phase
nucleation, Materials Science and Engineering, 2005, p.p. 157-202]
It is well known that the graphite nucleation would be prior to diamond nucleation in
the G region from the point of view of thermodynamics. Therefore, the diamond
nucleation would not happen unless the graphite nucleation is restrained or stopped.
For the issue, the most popular explanation is that the atomic hydrogen plays an
important role. Atomic hydrogen is an essential factor in CVD diamonds due to its
higher etching rate for the graphite phase and less etching rate for the diamond
phase. Unfortunately, some researchers have reported that diamond films are grown
upon CVD with a hydrogen-free environment. Further, Gruen [] concluded that CVD
diamonds do not require the reactant gas mixtures consisting primarily of hydrogen,
and the microstructure of diamond films can change continuously from micro- to
nanocrystalline when hydrogen is successively replaced by a noble gas such as
argon. Moreover, they pointed out that a chief function of the atomic hydrogen is to
reduce the secondary nucleation rates. Therefore, these experiments clearly indicate
that, besides enhancing the growth of the diamond nuclei, the atomic hydrogen may
have little function for the diamond primary nucleation. Our question is: would the
diamond nucleation be really in the strongly unstable region of the diamond structural
phase upon CVD?
NANOMATERIALS
114
The nucleation of CVD diamonds should happen in the D region of Figure 6-3 on the
basis of the nanosize-induced additional pressure. When we assume the surface
tension of diamond is 3.7 J/m
2
, the dependence relationships of the nanosize-
induced additional pressure based on the LaplaceYoung equation and the size of
diamond clusters can be obtained, as shown in the inset of Figure 6-3. From the
inset, one can see that the additional pressure increases with the crystal particles
size decreasing. Notably, when the radius is less than 4 nm, the additional pressure
goes up to above 2.0 GPa, which is above the phase equilibrium line shown as the D
region, i.e., the diamond stable region, in Figure 6-3. In other words, the nanosize-
induced additional pressure could drive the metastable region (G region) of the
diamond nucleation into the new stable region (D region) in the thermodynamic
equilibrium phase diagram of carbon. These deductions are supported by the
experimental cases from the CVD diamonds on non-diamond substrates. For
instance, Lee et al. reported that the size of the nuclei of CVD diamonds on Si
substrates is in the range of 26 nm []. Consequently, the nanosize-induced
additional pressure of 13 nm radius of the diamond nuclei would be enough to drive
the G region into the D region. Therefore, the nucleation of CVD diamonds should
happen in the D region based on our nanothermodynamic approach.
According to Eq. ( ) ) ( 4
2
3
4
, ,
2
0 0
3
f r
V
r
b T k P V
r T P r G
m
|
|
|
|
\
|
+
|
\
|
+
= (6.
90), = 3.7 J/m
2
, Vm = 3.417x10
-6
m
3
mol
-1
, V = 1.77x10
-6
m
3
mol
-1
, k
0
= 2.01x10
6
,
and b
0
= 2.02x10
9
Pa, one can obtain the relationship curves between the size of the
diamond critical nuclei and the pressure at the temperature of 1300 K upon the CVD
diamond case, and it is displayed in Figure 6-4, in which the inset shows the
dependent relations of the pressure and the critical radius at the given various
temperatures.
NANOMATERIALS
115
Figure 6-4 : The relationship curves between the radii of the critical nuclei and the pressure at the
temperature of 1300 K upon CVD systems. The inset shows the dependence relation of the pressure
and the critical radii at given various temperatures. One can see that, in faith, the radii of the critical
nuclei of diamond upon CVD for a broad range of pressures and temperatures are less than 5 nm.
Namely, the nucleation of CVD diamond could occur in the stable phase region of diamond. [C.X.
Wang, G.W. Yang, Thermodynamics of metastable phase nucleation, Materials Science and
Engineering, 2005, p.p. 157-202]
Clearly, we can see that the radii of the critical nuclei are less than 5 nm in a broad
range of the pressure and temperature. The diamond nucleation upon CVD seems to
be in the stable phase region of diamond due to the driving of the nanosize-induced
additional pressure of the diamond nuclei. In addition, from this figure, one also can
see the very weak dependence of the pressure on the critical radius. In other words,
at the given temperatures, the critical radii are hardly changed with the external
pressure change, because the external pressure is quite small compared with the
nanosize-induced additional pressure.
Based on Eq.
P b T k
V
V
r
m
+
|
\
|
+
=
0 0
*
3
2
2
(6. 91), the value of the surface tension of silicon
(1.24 J/m
2
), and the given parameters above, we display the relationship curves
between the pressure and the critical energy of CVD diamonds at the temperature of
1200 K, as shown in Figure 6-5, in which the inset displays the dependent relations
of the pressure on the critical energy at the given various temperature.
NANOMATERIALS
116
Figure 6-5 : The relationship curves between the pressure and the critical energy of the nucleation
upon CVD diamond at the temperature of 1200 K, and the inset showing the curves at various given
temperatures. [C.X. Wang, G.W. Yang, Thermodynamics of metastable phase nucleation, Materials
Science and Engineering, 2005, p.p. 157-202]
Obviously, one can see in Figure 6-5 that the critical energy of the diamond nuclei
slowly increases with the pressure increasing at a given temperature, and
approximately remains as unchanged. The case results from the too little external
pressure compared with the nanosize-induced additional pressure. These results
indicate that the critical energy of the diamond nucleation upon CVD is quite low
(10
-16
J), suggesting that the heterogeneous nucleation of CVD diamonds does not
require high forming energy.
Apparently, the low forming energy of the heterogeneous nucleation of CVD
diamonds implies that it is not difficult to nucleate diamond, and the diamond
nucleation could happen in the diamond stable region (G region) as shown in Figure
6-3. Based on the analyses above, the diamond nucleation upon CVD would happen
in the diamond stable region in the carbon phase diagram from the point of view of
thermodynamics. In fact, when the size of the crystalline particles is in the nanometer
scale, the additional pressure induced by the curvature of the nanometer-sized
particles is so high as to exceed the equilibrium pressure between diamond and
graphite, i.e., going up to break through the phase equilibrium boundary line, which
means that the additional pressure could drive the thermodynamic phase region of
the diamond nucleation from the metastable to the stable as shown in Figure 6-3.
Therefore, the diamond nucleation upon CVD is not actually the thermodynamic
paradox, and violating the second law of thermodynamics. These results indicate
that the presence of the atomic hydrogen is not a vital factor to grow diamonds upon
NANOMATERIALS
117
CVD from the viewpoint of the nanothermodynamics above. However, why did most
experimental studies all show that the atomic hydrogen plays a very important role in
CVD diamonds? In fact, it is recognized experimental evidence that the atomic
hydrogen etching the graphite phase (more etching rate to the graphite phase and
less etching rate to the diamond phase) and helping the sp
3
hybridization of carbon
atoms. Naturally, the diamond nucleation could be enhanced in the low-pressure gas,
only when the graphite phase forming is restrained or stopped by the atomic
hydrogen or other factors. Thus, the presence of the atomic hydrogen could increase
the rates of the diamond growth. Accordingly, the effect of the atomic hydrogen on
the diamond growth is much larger than that on the diamond nucleation upon CVD.
In other words, the influence of the atomic hydrogen on the diamond nucleation
would be small from the point of view of the experimental investigations involved in
how to enhance the diamond nucleation upon CVD.
In conclusion, aiming at a clear insight into the nucleation of CVD diamonds, we
studied the diamond nucleation from the point of the view of a nanoscaled
thermodynamics. Notably, these theoretical results show that the diamond nucleation
would happen in the stable phase region of diamond in the thermodynamic
equilibrium phase diagram of carbon, due to the nanosized effect induced by the
curvaturesurface tension of the diamond nuclei. In other words, at the nanometer
size, the diamond nucleation is prior to the graphite nucleation in competing growth
of diamond and graphite upon CVD.
6.2 Thermodynamics of melting and freezing in small
particles
6.2.1 Theory Vanfleet and al. model
Couchman [] developed a thermodynamic framework which Vanfleet and al. [] have
modified with Van Der Veens [] ideas to include surface melting. The combination of
their ideas results in a thermodynamic model that can help explain observed results.
Vanfleet and al. [] base their model on a spherical particle surrounded by another
material. This surrounding material will be referred to as the imbedding matrix. This
matrix can be as simple as the particles own vapor, or as complicated as a
crystalline or amorphous solid (although certain assumptions will be made shortly
NANOMATERIALS
118
about its effect on the particle). The particles temperature will be held constant, and
we will explore the Helmholtz free energy difference between the particle when it has
a liquid layer with a solid core and when it is fully solid. We assume that both the
solid and liquid phases are incompressible. Our interest is the shift of the melting and
freezing points due to finite size effects, so we must include both the surface energy
() and the surface stress (). We will assume that both of these surface quantities
are isotropic and independent of particle size and temperature. We will also assume
that the only effect of the matrix will be to vary the surface energies and surface
stresses of the particle. There can be confusion on the use of the terms surface
energy and surface stress. is the surface energy and has units of energy per area.
times the area is the energy attributed to the surface; it is the energy required to
make new surface. It can also be thought of as the energy of the unsatisfied bonds
on the surface, but as we will see later, it also includes the surface-induced order or
disorder in the bulk of the material. is the surface stress and also has units of
energy per area. This is the energy required to alter a given surface by elastically
deforming that surface, including the energy of the stretched or compressed bonds
within the surface. These two surface quantities are related. , the surface stress, is
the change in the total energy of the surface with a small change in its area.
A
A
dA
A d
+ = =
) (
(6. 94)
For the surface of a liquid the molecules are free to move and therefore cannot
sustain a bond stress. Thus, the second term in Eq.
A
A
dA
A d
+ = =
) (
(6. 94) is
zero and the surface stress is equal to the surface energy. The final surface will be
just like the original yet there will be more or less of it. Solids can sustain stresses.
Thus, the second term in Eq.
A
A
dA
A d
+ = =
) (
(6. 94) can be significant and is
of the same order of magnitude as the surface energy. It may be either positive or
negative. The choice of which surface quantity to use depends upon whether new
surface is being created () or an existing surface is being modified ().
The free energy of the fully solid particle is
F
solid
= (T, P) N
s
P
s
V
s
+
sm
A
s
(6. 95)
NANOMATERIALS
119
is used because the solid matrix surface must be created, and (T, P), a function
of the local temperature and pressure, is the chemical potential per particle. With a
spherical geometry assumed, the surface area (A
s
) and the volume (V
s
) are just
functions of the number of atoms in the particle (N
s
).
Here, and for the rest of the paragraph, a single subscript will refer to the volume
specified by the subscript, and a double subscript will refer to the surface specified by
the two subscripts. For example, in the above equation, the number of atoms in the
solid volume is N
s
the surface area of the solid volume is A
s
and the solid-matrix
surface energy is
sm
The local pressure (P) is the externally applied pressure (P
0
) plus the increase in
pressure due to the surface. The external pressure (P
0
) is all pressures except those
that arise from the particle itself. The increase in pressure, using the Laplace-Young
equation for a solid or liquid sphere of radius r, is
r
P P
2
0
= (6. 96)
The surface stress () is used because the surface is attempting to compress the
particle and change its area. The chemical potential is analyzed about the bulk
melting point (T
0
) and the externally applied pressure (P
0
). We will assume that the
externally applied pressure is small enough that there is no shift in the bulk melting
point from standard values. At this reference point, the chemical potential per particle
is the same for the solid and liquid phases. Only first-order terms are kept, to give
( ) ( ) ) ( ) ( , ,
0 0 0 0
P P
P
T T
T
P T P T
+ =
(6. 97)
That is
( ) ( ) ) ( ) ( , ,
0 0 0 0
P P T T s P T P T + + = (6. 98)
Where s is the entropy per particle and is the volume per particle. The pressure
due to the surface can be significant, but because of the assumed incompressibility
and small linear expansion coefficient, we can neglect second-order terms in
pressure. The second-order term in temperature is at most a few percent of the first-
order terms and will also be neglected. The final free energy of a crystalline particle is
( )
s sm s
s
sm
s
s
sm
s s s s solid
A V P
r
N
r
N T T s N P T F
+
|
|
\
|
+
|
|
\
|
+ + =
0 0 0 0
2 2
) ( , (6. 99)
NANOMATERIALS
120
Where the terms involving the surface stress will cancel.
For the liquid particle with the crystalline core we must split these free energy terms
between the phases. The subscript l will refer to the liquid layer and r
l
will be the
radius of the liquid-matrix surface. The subscript c will refer to the solid core that is
surrounded by the liquid layer and r
c
will be its radius.
F
mixed
= F
liquid
+ F
core
(6. 100)
( )
l lm l
l
lm
l
l
lm
l l l l liquid
A V P
r
N
r
N T T s N P T F
+
|
|
\
|
+
|
|
\
|
+ + =
0 0 0 0
2 2
) ( , (6. 101)
( )
c sl c
c
sl
l
lm
c
c
sl
l
lm
s c s c core
A V P
r r
N
r r
N T T s N P T F
+
|
|
\
|
+ +
|
|
\
|
+ + + =
0 0 0 0
2 2 2 2
) ( ,
(6. 102)
The volumes ( ) and entropies (s) are labeled for the phase (solid or liquid) to which
they apply. Again, surface stress terms cancel. This cancellation of the surface stress
is a result of the assumed incompressibility of the phases. An incompressible
material can have no net work done by changes in the pressure. Therefore, the
change in the chemical potential is countered by another term to leave no net change
in the free energy due to pressure.
The free energy difference between these two situations is now
F = F
mixed
- F
solid
(6. 103)
( )( ) ( )( ) ( )
s sm c sl l lm
s c l s s c c l l s c l
A A A
V V V P T T N s N s N s N N N P T F
+ +
+ + + =
0 0 0 0
,
(6. 104)
The total number of atoms is conserved (N
l
+ N
C
- N
s
) = 0. Entropy terms can be
written as the difference between the liquid and solid entropy per particle and the
number of liquid particles.
( )( ) ( )
|
|
\
|
= +
0
0 0
1
T
T
N s s T T T N s N s N s
l s l s s c c l l
(6. 105)
The latent heat (l
0
) per atom is (s
l
- s
s
)T
0
( )
s sm c sl l lm s c l l
A A A V V V P N
T
T
l F + + +
|
|
\
|
=
0
0
0
1 (6. 106)
NANOMATERIALS
121
There is a bulk term due to the temperature difference from the bulk value, a volume
change term, and surface terms that are proportional to the surface areas. From
geometrical considerations given below, the PV terms can be rewritten as a function
of N
l
.
( )
l
s l
s c l
N P V V V P
|
|
\
|
= +
1 1
0 0
(6. 107)
Where is the density (atoms per volume) of the labeled phase. This term is small
compared to the latent heat term and will be neglected. For example, at atmospheric
pressure this is about 10
-4
times the latent heat term, unless the temperature is very
near the bulk melting point.
We would like to take into account the ordering or disordering induced by the surface.
Vanfleet and al. took their cue from Van Der Veen and they derived a correction to
the surface energies, due to surface-induced ordering or disordering. A surface will
attempt to propagate its order into the volume of which it forms the surface. The order
will be related to a correction to the energy for each unit of volume. This effect will be
screened by intervening material and is analogous to a screened potential. Using the
screened potential decay of the interaction we will write that the energy, due to the
surface, at a distance r from an area dA will be
dA
r
r
dE
|
|
\
|
exp
(6. 108)
Where is the scale length for the screening. The total contribution, at a point, due to
the surface would then be the integral over the surface of the above. For a flat
surface this integration gives C exp (-t / ), where C is the constant of proportionality
(units of energy per unit volume) and t is the perpendicular distance from the surface.
Therefore, for a flat surface, the total contribution to the energy due to the surface is
the integral over the volume of C exp (-t / ), plus a surface term.
|
|
\
|
+
|
|
\
|
+ =
2 2
2
2
1 1
1
1
*
exp exp dV
z
C dV
z
C A E
surface
(6. 109)
The 1s and 2s designate which side of the surface is the volume to be integrated. If
we limit our interest to a thickness t on each side of the surface and integrate we
have
NANOMATERIALS
122
(
(
|
|
\
|
|
|
\
|
+
|
|
\
|
|
|
\
|
+ =
2
2
2 2
1
1
1 1
*
exp 1 exp 1
t
C
t
C A E
surface
(6. 110)
For ts much larger than the scale lengths, the energy per unit area due to the
surface is
2 2 1 1
*
C C + + = . This is the observed surface energy that we find in the
literature. The limit of t going to 0 gives *, which we will call the skin energy. If the
presence of the surface increases the order for a short distance into the bulk, the
energy will be lowered, and C will be negative. Increasing the disorder of the material
will raise the energy and C will be positive. If we now ask what is the energy per area
due to the surfaces, for non-interacting solid-liquid liquid-matrix surfaces a distance t
apart, where the solid and matrix are considered infinite, we have
s sl
l
l ls m ml
l
l lm sl lm
C
t
C C
t
C
A
E
+
|
|
\
|
|
|
\
|
+ +
|
|
\
|
|
|
\
|
+ + = exp 1 exp 1
* *
(6. 111)
Van Der Veen stated that there must be no difference between the energy of a solid
surface, and a solid surface with a liquid layer, if the liquid layer has zero thickness.
We will make that assumption but will see shortly that this imposes some conditions
under which this model can be applied. Requiring that the above equation give
sm
in
the limit of t = 0 results in
s sl m ml sl lm sm
C C + + + =
* *
. Letting t now be large
compared to the scale length, the mismatch between the observed values of the
surface energies is
( )
l ls lm s sl l ls m ml l lm sl lm sm lm sl sm
C C C C C C + = + = =
* *
(6. 112)
Because the presence of the surface will impose ordering on the liquid, both C
lm
and
C
ls
will be negative and therefore is positive. We know that there are materials
where this is not the case. When < 0 we cannot impose the constraint that
Eq.
s sl
l
l ls m ml
l
l lm sl lm
C
t
C C
t
C
A
E
+
|
|
\
|
|
|
\
|
+ +
|
|
\
|
|
|
\
|
+ + = exp 1 exp 1
* *
(6. 111)
give
sm
, when t = 0. < 0 is also the non-wetting condition, therefore, a flat liquid
layer would be unstable in that situation. Thus, for the moment, we will confine our
interest to > 0 and assume that Eq.
( )
l ls lm s sl l ls m ml l lm sl lm sm lm sl sm
C C C C C C + = + = =
* *
(6.
112) holds. The energy of these two surfaces separated by a thickness t can now be
written as
NANOMATERIALS
123
( )
( )
(
(
|
|
\
|
|
|
\
|
+ + =
(
(
|
|
\
|
|
|
\
|
+ + + =
|
|
\
|
l
sl
t
lm sm lm
l
lm sl sm sm
t
A E
t
A E
l
exp 1
exp 1
exp
(6. 113)
This last version is exactly the correction given by Van Der Veen.
To simplify the implementation with spherical geometry, we will make the
approximation that the corrections due to order or disorder propagating into the solid
or matrix are small compared to the contributions due to the order in the liquid. This
approximation will leave us with surface terms and the bulk liquid terms which will be
written as
=
=
) 1 (
l ls
l lm
C
C
(6. 114)
will be between 0 and 1 and will be approximated for the experiment in question.
divides between the order imposed on the liquid due to the surface with the matrix
and that due to the surface with the solid. If the solid and matrix were very similar, we
would expect to be near one half. The other extreme would be a very dissimilar
matrix or vapor, in which case would be small but non zero.
Even if the matrix has no order to impose on the liquid, the surface itself will impose
some order. Therefore, will never be zero.
Going back to integrating over the surface, but changing to a spherical geometry we
get
( )
|
|
\
|
|
|
\
| +
|
|
\
|
r R
r R
r
R
r E exp exp ) ( (6. 115)
R is the radius of the surface and r is the radius of the point in question. The total
energy due to the surface will again be found by integrating the above equation over
the volume. Using the approximation made above, we will limit our interest to
corrections to the bulk surface energies that involve . Using the geometry of a
liquid sphere with a solid core, the energy of the solid-liquid and liquid-matrix
surfaces is
sl lm sl sl lm lm
G G A A E + + + = (6. 116)
NANOMATERIALS
124
( ) ( ) ( )
(
|
|
\
|
+
|
|
\
| +
+ +
|
|
\
|
+ =
l
l
l l
l
c l
l c
l
c l
l c l l lm
r
r
r r
r
r r
r r G
2
exp
) (
exp
) (
exp 4
(6. 117)
( ) ( ) ( )
(
|
|
\
|
+ +
|
|
\
| +
+
|
|
\
|
+ =
l
c
l c
l
s l
l l
l
s l
l l l c sl
r
r
r r
r
r r
r r G
2
exp
) (
exp
) (
exp ) 1 ( 4
(6. 118)
This is more complicated than the simple form given by Van Der Veen, but it takes
into account the disappearance of the solid-liquid surfaces correction when the solid
core disappears.
Using the new correction to the solid-liquid liquid- matrix surface energies results in
the following free energy difference.
sl lm s sm c sl l lm l
G G A A A N
T
T
l F + + + +
|
|
\
|
=
0
0
1 (6. 119)
This can be written as a function of the number of atoms in the liquid layer by using
the total number of atoms, the solid and liquid densities, and the known geometry to
get
3
1
4
3
|
|
\
|
=
s
s
s
N
r
(6. 120)
3
1
4
3
|
|
\
|
=
s
l s
c
N N
r
(6. 121)
3
1
4
3
|
|
\
|
|
|
\
|
+
=
l
l
s
l s
l
N N N
r
(6. 122)
The thickness of the liquid layer, the volumes, and the areas of each surface follows
normally from these radii.
NANOMATERIALS
125
Figure 6-6 : Free energy difference from crystalline state as function of fraction of particles in melted
layer. Data are for a lead particle r = 5.0 nm at three different temperatures. (A) T=480 K,
(B) T=502 K and (C)T=540 K. [R.R. Vanfleet, J.M. Mochel, Thermodynamics of melting and freezing in
small particles, Surface Science, 1995, p.p. 40-50]
Figure 6-6 shows the free energy difference versus the fraction of the atoms in the
liquid layer for three different temperatures. The data are for an r = 5 nm lead particle
and the temperatures are as follows: (A) is T=480 K, (B) is 502 K, and (C) is 540 K.
For (A) and (B), there is a minimum in the energy for a liquid layer of finite thickness.
(B) corresponds to the layer being one monolayer in thickness. This minimum arises
from the propagation of the surface order into the liquid, and corresponds to surface
melting as predicted by Van Der Veen et al. []. (A) and (B) also show another local
minimum with all atoms melted, at the extreme right. This second minima is
separated from the surface melted state by an energy barrier. For (C), all atoms in
the liquid state is the minimum of energy. For a given material, both the presence or
the absence of these minima and the height of the barrier are functions of the size of
the particle, the temperature, the surface energies, the scale length, and the dividing
parameter .
NANOMATERIALS
126
Figure 6-7 : Melting point (M), freezing point (F), and temperature at which there is a monolayer of
surface melt (S) for lead, as a function of radius of the particle. G is the region below which a cluster
once disordered will remain disordered. Arrows indicate direction of temperature change when
transition occurs. [R.R. Vanfleet, J.M. Mochel, Thermodynamics of melting and freezing in small
particles, Surface Science, 1995, p.p. 40-50]
We have used bulk parameters in Figure 6-7 to calculate melting point, freezing
point, and temperature where the surface melting is equal to one monolayer in
thickness, as functions of particle radius for lead. Figure 6-7 also shows the region
where a disordered particle, once created, would remain disordered. This is the
frozen in region. The arrows indicate the direction of the temperature change for
which the transition would be expected to occur. The curves were calculated with a
simple approach of choosing a particle size and then adjusting the temperature until
the melting point, freezing point, surface melt and frozen in temperatures were
found. All these points were then strung together to form the curves in Figure 6-7 and
Figure 6-8.
NANOMATERIALS
127
Figure 6-8 : Melting point (M), freezing point (F), and temperature at which there is a monolayer of
surface melt (S) for platinum, as a function of radius of the particle. G is the region below which a
cluster once disordered will remain disordered. Arrows indicate direction of temperature change when
transition occurs. [R.R. Vanfleet, J.M. Mochel, Thermodynamics of melting and freezing in small
particles, Surface Science, 1995, p.p. 40-50]
For a given size and temperature, the free energy as a function of the number of
atoms in the liquid layer is calculated. The position and depth of the surface melt
minimum is found and the position is compared to the number of atoms in a single
monolayer. Based on this comparison, the temperature is either raised or lowered
and the calculation is repeated. This process is continued until the surface melting
temperature is found to within a step size. The melting point is based upon the
difference between the surface melt minimum and the maximum free energy.
The main features of Figure 6-7 are a difference between the melting and freezing
points, and the disappearance of this difference for small enough particles. In the
region of differing melting and freezing points, the dominant factor in determining
these critical temperatures is the height of the energy barrier. A lower temperature is
required for the liquid-to-solid transition to occur than for the solid-to-liquid transition.
The difference between melting and freezing disappears when the barrier becomes
so small that it can be easily overcome in either direction. At this point the absolute
stability becomes the important factor leading to identical melting and freezing points.
Another feature of Figure 6-7, is the lack of surface melting for the smallest particles.
This is where the melting barrier is surmountable before a surface melt layer will
form. The line labeled G in Figure 6-7, indicates the temperature below which a
particle once disordered would stay disordered due to a high free energy barrier.
NANOMATERIALS
128
6.3 Phase Diagrams
In order to determine the phase diagram of a small system whose thermodynamic
equilibrium is a metastable one determined in part by the size of the system, isolated
particles will be considered. The transition solid-liquid results in an equilibrium
between the solid and its melt. Further, in view of the experimental observation that
many liquids have a zero contact angle with their solids, this wetting feature of the
system will be incorporated in the model by having the solid in contact with its melt.
Even though the isolated particle will likely be crystalline and hence faceted, as a
simplification, spherical shapes will be assumed. With these initial assumptions and
conditions, the geometry of the system to be modeled is an isolated solid sphere
transforming on heating to a specific temperature to a liquid layer surrounding a
spherical solid core. The initial state of the particle before melting is a solid sphere.
Upon melting the solid sphere develops a sheath of liquid which will completely
surround a small solid core. This case is illustrated in Figure 6-9.
Figure 6-9 : Schematic drawing of isolated particle before melting (complete solid) and after melting
(liquid sheath exhibiting zero contact angle with the solid). [W.A. Jesser, G.J. Shiflet, G.L. Allen, J.L.
Crawford, Equilibrium phase diagrams of isolated nano-phases, Mat. Res. Innovat, 1999, p.p. 211-
216]
6.3.1 Governing equations
First, the thermodynamic potential that is a minimum at equilibrium is considered. For
an infinite reservoir in contact with a small, spherical particle through a wall that is
rigid (no change in volume of the reservoir), diathermal (heat exchanges may occur
with the reservoir), and impermeable (no exchange of matter between the system
and reservoir), it has been shown (see previous section) that the Helmholtz potential
is minimized at equilibrium. For the more realistic case of an infinite reservoir in
contact with the particle through a wall that is flexible, diathermal, and impermeable,
NANOMATERIALS
129
the Helmholtz potential is minimized at equilibrium in this case too, provided the
reservoir pressure is zero. Accordingly, thermodynamic equations for the Helmholtz
potential will be applied to the model of a sphere of solid isothermally melting to form
a spherical shell of liquid surrounding a solid core.
The approach is to express in analytical form the difference F in Helmholtz energy
between the final state, a liquid sheath of thickness t surrounding a solid core of
radius R
s
, and the initial state, a solid sphere of radius R
0
. The initial and final states
can be written in terms of the chemical potential at the actual pressure P,
temperature T, and composition X
i
in the small particle of volume V and area A and
the number of atoms in the spherical particle N
i
, where the subscript i refers to a
particular component, and the surface energies of the solid, liquid, and solid-liquid
interface are
s
,
l
,
sl
respectively.
The Helmholtz energy F
0
of the initial state can be written
s s s s
i
s
i
s
i
A V P N F
0 0 0 0 0 0
+ =
(6. 123)
Where the subscript 0 refers to the initial state of a solid sphere and the superscripts
s or l refer to the solid, or liquid phases. The Helmholtz energy for the final state of a
liquid sheath of thickness t surrounding a solid core of radius R
s
can be written as
s s s s
i
s
i
s
i
l l l l
i
l
i
l
i f
A V P N A V P N F + + + =
(6. 124)
The Helmholtz energy difference between the final and initial states, F is plotted
against the composition and fraction of liquid for constant materials parameters in
order to determine whether and where any minima exist on the energy surface. It will
be assumed that a regular solution model is adequate to represent the bulk chemical
potentials and that a linear relationship exists between the surface and interfacial
energies and composition.
Using a regular solution approach one can incorporate interaction parameters
l
for
the liquid and
s
for the solid. These parameters are taken as temperature, pressure
and composition independent. The final and simplified expression for the Helmholtz
free-energy difference F can be written as
( ) ( )
( ) ( )
( ) [ ]
sl s s l
B A
s s
B
s
A
s l l
B
l
A
l
B B A A
s
B
s
B
s
A
s
A
s l
B
l
B
l
A
l
A B
l
B A
l
A
l
A A
N X X N X X N X X X kT N X kT N
X X X X kT N X X X X C X C X kT N F
+ +
+ +
+ + + + + =
0 0 0 0 0 0 0
ln ln
ln ln ln ln ln ln
(6. 125)
NANOMATERIALS
130
Where
s s
B
s
B
s s
A
s
A
l l
B
l
B
l l
A
l
A
N N X N N X N N X N N X / ... , / ... , / ... , / ... (6. 126)
With
l
B
l
A
l
N N N + = and
s
B
s
A
s
N N N + = (6. 127)
and
s s
B
l l
B B
s s
A
l l
A A
N X N X N N X N X N + + ... , ...
0 0
(6. 128)
For convenience ( )
bulk
l l s s
B A
X X C / ...
,
represents the bulk compositions and activity
coefficients for components A and B respectively for the case of the bulk phase
diagram.
The relationship between interaction parameter and activity coefficient is
( )
l s
i
l s
i
kT X
, ,
ln 2 1 = (6. 129)
6.3.2 Mathematical description for nano-phases of Sn-Bi alloys
To find local minima in the Helmholtz free-energy difference between a spherical
solid particle and a liquid sheath surrounding a solid core, equation
( ) ( )
( ) ( )
( ) [ ]
sl s s l
B A
s s
B
s
A
s l l
B
l
A
l
B B A A
s
B
s
B
s
A
s
A
s l
B
l
B
l
A
l
A B
l
B A
l
A
l
A A
N X X N X X N X X X kT N X kT N
X X X X kT N X X X X C X C X kT N F
+ +
+ +
+ + + + + =
0 0 0 0 0 0 0
ln ln
ln ln ln ln ln ln
(6.
125) is plotted as a function of composition of the solid core and fraction of liquid
present in the two-phase state thereby generating an energy surface. The calculation
requires a number of terms to be specified. Typically one should have data for the
dependence of interface and surface energies as a function of composition and even
the interaction parameter used in the regular solution approximation can be limited
to a fixed composition range. In this section the calculations for a binary phase
diagram of tin and bismuth will be given. The solidus and liquidus lines of the Bi-rich
and Sn-rich sides of the bulk phase diagram provide the reference state of the bulk
phase diagram from which the small particle system of Bi-Sn is calculated. For the
Sn-rich side of the bulk phase diagram the equation of the solidus is closely
approximated by
NANOMATERIALS
131
883 . 0 10 62 . 6 10 73 . 5
3 6
+ =
T T X
s
Sn
(6. 130)
Table 6-1 : Values of parameters employed for the Bi-Sn binary phase diagram. [W.A. Jesser,
G.J. Shiflet, G.L. Allen, J.L. Crawford, Equilibrium phase diagrams of isolated nano-phases, Mat. Res.
Innovat, 1999, p.p. 211-216]
and that of the liquidus by
70 . 1 10 23 . 6 10 76 . 1
3 2 6
+ =
T T X
l
Sn
(6. 131)
while for the Bi-rich side of the phase diagram the equation of the solidus is taken as
639 . 1 10 01 . 3
3
+ =
T X
s
Bi
(6. 132)
and that of the liquidus as
06559 . 0 10 204 . 1
4
+ =
T X
l
Bi
(6. 133)
The values of the parameters used in equation 2 may be found in Table 6-1.
When the above parameters are substituted into equation
( ) ( )
( ) ( )
( ) [ ]
sl s s l
B A
s s
B
s
A
s l l
B
l
A
l
B B A A
s
B
s
B
s
A
s
A
s l
B
l
B
l
A
l
A B
l
B A
l
A
l
A A
N X X N X X N X X X kT N X kT N
X X X X kT N X X X X C X C X kT N F
+ +
+ +
+ + + + + =
0 0 0 0 0 0 0
ln ln
ln ln ln ln ln ln
(6.
125), a F surface versus composition and fraction liquid may be plotted and
analyzed for the presence of local minima. As an illustrative example of a typical F
surface Figure 6-10 is presented for a binary bismuth alloy with lead for two
temperatures.
NANOMATERIALS
132
Figure 6-10 : Helmholtz free energy surface plotted for the case of a binary bismuth alloy as a function
of composition and fractional amount of liquid for two temperatures. a) intermediate temperature
corresponding to a metastable solid-liquid equilibrium, b) elevated temperature showing a stable
equilibrium corresponding to complete liquid. [W.A. Jesser, G.J. Shiflet, G.L. Allen, J.L. Crawford,
Equilibrium phase diagrams of isolated nano-phases, Mat. Res. Innovat, 1999, p.p. 211-216]
Here one finds a relative minimum whose location is indicated by the broad dot and
an absolute minimum located at the point where only liquid exists. These two minima
are separated by a saddle-point maximum that prevents the system from jumping this
energy barrier. Accordingly, the temperature of this plot corresponds to a metastable
solid-liquid equilibrium. At significantly lower temperatures no local minimum exists
and complete solid is the stable state. At significantly higher temperatures, again no
local minimum exists and the saddle-point maximum disappears so that the free
energy surface is one that continuously decreases as one moves toward the point of
minimum free energy which corresponds to complete liquid. Figure 6-10 b)
represents this latter condition in which the minimum is indicated by the broad dot.
There is no progression of the local minimum continuously to the condition of
complete liquid. Rather there is a discontinuous, abrupt shift of the minimum from a
location away from the point of complete liquid to the point of complete liquid once
the saddle-point maximum disappears. Physically this corresponds to heating the
two-phase particle to a temperature at which a discontinuous transition occurs from
the state of a liquid sheath surrounding a solid core to the state of complete liquid.
Such an abrupt transition is experimentally observed for this binary alloy as well as
NANOMATERIALS
133
for the solid elements tin and bismuth both of which exhibit no stable solid-liquid
equilibrium.
6.3.3 Phase diagram for isolated nano-phases of Sn-Bi alloys
A temperature-composition phase diagram for isolated nano-spheres of Sn-Bi alloys
can be constructed as a function of size of the isolated spherical particle. This is
accomplished by tracking of the position of the local minimum in the Helmholtz free
energy surface as a function of temperature for several values of radius of the initial
solid sphere. The resultant phase diagram for the liquid-solid phase equilibria, i.e.,
the location of the liquidus and solidus lines are shown in Figure 6-11 for the case of
a spherical particle 40 nm in diameter.
Figure 6-11 : Phase diagram for the binary alloy Sn-Bi. The bulk phase diagram is shown
superimposed on the phase diagram for isolated spheres of alloy 40 nm in diameter. The nano-particle
phase diagram is the shaded region as given by the thermodynamic calculations and the data points
for each of eight compositions as determined from in-situ TEM measurements. [W.A. Jesser, G.J.
Shiflet, G.L. Allen, J.L. Crawford, Equilibrium phase diagrams of isolated nano-phases, Mat. Res.
Innovat, 1999, p.p. 211-216]
This Figure also includes the phase diagram for bulk Sn-Bi alloys and the
experimentally measured liquidus and solidus lines for the isolated spherical
particles. The experimental data reported here were measured inside a transmission
electron microscope modified for ultra high vacuum so that in-situ deposition of the
isolated nanospheres could be performed. These new data show that for the Sn-rich
side of the phase diagram, the liquid-solid two-phase field pinches off to a line
thereby indicating no stable equilibrium between the liquid and solid phases in this
NANOMATERIALS
134
region. There was no detectable pinching off observed on the Bi-rich side of the
diagram.
In addition to the pinching off of the two-phase field near the composition of pure tin,
there are two other noteworthy features in the phase diagram. First, one notices that
there is a substantial increase in solubility for the nano-particles as compared to that
in the bulk alloy. The solubility limit of tin in bismuth increased from 1.5% (bulk) to
about 43% in particles of radius 20 nm. A strong increase also occurred in the tin-rich
phase.
Second, because the initial melting event is observed to be discontinuous (the liquid
sheath forms with a minimum finite thickness of about 2 nm) and the final melting
event is also a discontinuous event with the isothermal jump from a thick liquid
sheath surrounding a solid core to a completely liquid sphere, the location on the
phase diagram of the locus of the temperatures at which these discontinuous
transitions occur do not represent liquidus and solidus lines that can be used to
construct the tie lines. The conventional tie lines that describe mass conservation and
are useful for lever-rule calculations extend isothermally beyond the two phase
regions. If a tie line is constructed between the liquidus and solidus lines shown on
the phase diagram, its ends do not correspond to either the solid composition or the
liquid sheath composition. In order to account for these compositions a longer tie line
is required. As the temperature is increased, the tie line changes its length. The two
ends of the tie line demark the compositions of the liquid and solid phases and
sweep out two areas on the phase diagram both of which are detached from the
drawn liquidus and solidus lines that demark the extent of the two phase field. This
feature is shown in Figure 6-12.
Figure 6-12 : Sn-rich segment of a Sn-Bi phase diagram for an isolated sphere of radius 20 nm
showing that the areas swept out by the ends of the tie lines are detached from the two-phase region.
NANOMATERIALS
135
Note that the two-phase regions and the ends of the tie lines pinch off at about 5% Bi. [W.A. Jesser,
G.J. Shiflet, G.L. Allen, J.L. Crawford, Equilibrium phase diagrams of isolated nano-phases, Mat. Res.
Innovat, 1999, p.p. 211-216]
One can incorporate the size-dependence of the phase diagram by making a three-
dimensional plot of liquidus and solidus temperatures as a function of both
composition and size. One can incorporate a size variable in the phase diagram by
recognizing that the Helmholtz free energy depends upon pressure which exhibits a
discontinuity P across a curved surface of
P = 2 / R (liquids) (6. 134)
Accordingly the equilibrium lines shift approximately linearly with reciprocal radius.
The axes of the three-dimensional phase diagram then become temperature,
reciprocal radius, and composition. Figure 6-13 shows a 3-D phase diagram for the
Sn-Bi alloy system.
Figure 6-13 : 3-D phase diagram for an isolated tin-bismuth alloy sphere showing the size
dependence of the liquidus and solidus lines as well as the eutectic isotherm over the range of
reciprocal radius 0 nm
1
to 0.05 nm
1
. Over this size range the eutectic temperature changes from 412
to 320 K and the solubility of Sn in Bi increased from 1.5% to 43%. The dashed lines on the left margin
are reference lines parallel to the 1/R axis. The solid lines on the left margin are projections of the
melting and eutectic temperatures. The solid lines on the bottom margin are projections of the
solubility limits and the eutectic composition. [W.A. Jesser, G.J. Shiflet, G.L. Allen, J.L. Crawford,
Equilibrium phase diagrams of isolated nano-phases, Mat. Res. Innovat, 1999, p.p. 211-216]
6.4 Crystal-lattice inhomogeneous state
In nanoparticles and nanocrystallities the fraction of surface (or
interface) atoms is large. Then, appreciable distinctions between nano-objects
NANOMATERIALS
136
and massive bodies in lattice parameters and the types of atomic structures are
plausible. These difference manifest themselves in normal and tangential
relaxations of nanoparticles.
(54)
An atom on the surface has fewer neighbours
than one in the volume, all the neighbours of the former being located on one
side only. This breaks the equilibrium of the interatomic forces and gives rise to
both the change of interatomic spacing (normal relaxation) and tensile
deformations which smooth verticles and edges (tangential relaxation).
A rough estimate of the average distortion of the interatomic distances a
in nanoparticles caused by surface tension is shown below:
(55)
'
2
r a
a
=
(6. 135)
where
2
is the numerical coefficient of the order of 1 (in the case of isotropic
spherical particles it may be estimated as
2
/
3
), is the volume compressibility
coefficient. According to the estimate, the lattice construction is about several
tenths of a percent at typical values of parameters ~ 10
-11
m
3
/J, 1 11 1 J/m
2
for
nanoparticles with radius r ~ 10 nm.
Here it should be mentioned that parameters entering this formula are, in
turn, size-dependent. The surface energy of nanoparticles, for example,
depends on r as follows:
(56)*
|
|
\
|
+ |
\
|
+ = ... 1 '
2
4 3
r
a
r
a
(6. 136)
where numerical coefficient 1 11 10 00 0
1 11 11 11 1
m
3
/J, 1 11 1 J/m
2
. According to this estimate
the surface energy of nanoparticles decreases as r decreases which is
intuitively clear, since both the number of bonds formed by the surface atoms
and their bond energy are reduced as size decreases. Size-dependent
corrections are, in a general case, functions of temperature.
(59)
The difference in
the lattice parameters of contacting nanoparticles may seriously influence the
process of interparticle sliding.
For several decades contradictory experimental results were
calling into question eq.
'
2
r a
a
=
5
(6. 137)
Where the numerical coefficient
5
is of the order of unity. For spherical
nanoparticles (nanocrystallites) it is close to 5. Using
s
3a as typical width of
the segregated layers the corresponding diameter of nanoparticles
(nanocrystallites) is d ~10nm. Thus, under some conditions the multicomponent
nanocrystals with grains of the size may be considered massive segregated
layers.
If the main contribution of driving force for segregation is the mismatch in
atomic radii of the matrix and foreign atoms,
s
and r
*
s
may be estimated for a
heterophase nanocrystal. The mismatch may be described by the relative
NANOMATERIALS
138
difference
s
= (a
f
/ a 1) where a
f
is the effective radius of a foreign atom. If the
nanocrystallites in a nanocomposite the components of which have shear
moduli G and G are elastically isotropic and possess non-slipping boundaries,
the characteristic length l
*
sl
is as follows:
(65)
( )
3 / 1
3 2
6
'
|
|
\
|
=
kT
G G a
a l
s
s
(6. 138)
Where the numerical coefficient
6
3.5, k is the Boltzman constant, and
T is temperature. This estimate was obtained by equation kT to the energy of
the elastic interaction of the foreign atom with the non-slipping interface.
Provided that other driving forces for segregation (heat of mixing, the difference
between surface energies of the component) are negligible or cancel each
other, when a nanocrystallite has a shear modulus G less than that G of the
neighbouring crystallites segregation in that nanocrystallite becomes impossible.
So, the concentration of foreign atoms in the vicinity of interfaces diminishes in
elastically softer nanocrystallites.
Segregation effects change the composition of the interior of
nanoparticles. It can be shown
(66)
that the average concentration c of the
second component for a two-component nanoparticle varies with its size,
following the relationship:
( )
r
c c c c
s
b s b
+ =
8
(6. 139)
where the numerical coefficient
8 88 8
3, c
s
and c
b
are concentrations in the
surface layer of width
s
and in bulk of nanoparticles, respectively. If
concentrations c
s
and c
b
were independent of r over the whole range of sizes,
this relationship would be of a general character. When size r is not too small,
one may estimate the change in the concentration of foreign atoms in the
smaller the size r, the stronger the concentrational inhomogeneity which may
reach as much as 10%. This conclusion agrees with the numerical simulations
of small-atomic clusters
(67)
and experiments
(68)
with CuNi nanoparticles.
Certain conclusions about the character of the independence of surface
concentration c
s
on nanoparticle (nanocrystallite) radius r can be drawn on the
basis of simple mass balance considerations. In the limit of small nanoparticles,
NANOMATERIALS
139
when all the atoms belong to the surface, the surface composition in an
ensemble of nanoparticles would be identical to the average composition in the
massive alloy.
(69)
Such a decrease of the surface segregation has been
observed experimentally.
(70)
In the opposite of very large particles they will
behave like semi-infinite crystals, and the surface composition would be similar
to that obtained in microscopic samples. In the intermediate range of sizes,
where the number of surface sites is not a negligibly small fraction of the total
number of sites available in the particle volume, there exists a regime where
segregation to the surface tends to decrease the bulk concentration of the
segregant.
The depletion of the segregated atoms, which would follow from
elementary mass-balance considerations, may be impossible if the physical
and/or chemical properties are size-dependent. For instance, the change in the
surface energy (see eq.
|
|
\
|
+ |
\
|
+ = ... 1 '
2
4 3
r
a
r
a
(6. 136)) may modify the
segregation effect. The lattice contraction (eq.
'
2
r a
a
=
(6. 135)):
( )
2 1
3
=
P
s
(6. 140)
where
1
and
2
are the volume compressibility of the segregant and the
solvent, respectively. The hydrostatic stress include by surface tension,
approximated as 2 / l, may, thus, change the relative size mismatch by
s
.
According to the well-known Hume-Rothery rule (if the relative size mismatch
exceeds the 15% limit, solid solubility will be restricted
(72)
) a critical size l
*
sol
should exist, above which the solid solubility remains qualitatively the same:
NANOMATERIALS
140
( )
s
sol
15 . 0
2 1
2
(6. 141)
6.6 References
C.X. Wang, G.W. Yang, Thermodynamics of metastable phase nucleation, Materials
Science and Engineering, 2005, p.p. 157-202
T. L. Hill, A Different Approach to Nanothermodynamics, Nano letters, 2001 Vol.1,
N5, p.p. 273-275
C. Beck, Non-additivity of Tsallis entropies and fluctuations of temperature,
Europhysics Letters, 2002, p.p. 329-333
R.R. Vanfleet, J.M. Mochel, Thermodynamics of melting and freezing in small
particles, Surface Science, 1995, p.p. 40-50
W.A. Jesser, G.J. Shiflet, G.L. Allen, J.L. Crawford, Equilibrium phase diagrams of
isolated nano-phases, Mat. Res. Innovat, 1999, p.p. 211-216
NANOMATERIALS
141
7 ELECTRONIC AND OPTICAL PROPERTIES OF
NANOMATERIALS
7.1 Introduction
Energy
Forbidden bandgap
HOMO
LUMO
Molecule
Nanoparticle
Semiconductor Metal
Crystalline
Atom
Cluster
Figure 7-1: Size quantization effect. Electronic state transition from bulk metal/semiconductor to small
cluster.
Nanocrystalline particles represent a state of matter in the transition region
between bulk solid and single molecule. As a consequence, their physical and
chemical properties gradually change from solid state to molecular behaviour with
decreasing particle size. The reasons for this behaviour can be summarised as two
basic phenomena.
First, as we have already discussed, owing to their small dimensions, the
surface-to-volume ratio increases, and the number of surface atoms may be similar
to or higher than those located in the crystalline lattice core, and the surface
properties are no longer negligible. When no other molecules are adsorbed onto the
nanocrystallites, the surface atoms are highly unsaturated and their electronic
contribution to the behaviour of the particles is totally different from that of the inner
atoms. These effects may be even more marked when the surface atoms are
NANOMATERIALS
142
ligated.
lxiii
This leads to different mechanical and electronic transport properties,
which account for the catalytic properties of the nanocrystalline particles.
lxiv
The second phenomenon, which only occurs in metal and semiconductor
nanoparticles, is totally an electronic effect. The band structure gradually evolves with
increasing particle size, i. e., molecular orbital convert into delocalised band states.
Figure above, shows the size quantization effect responsible for the transition
between a bulk metal or semiconductor, and cluster species. In a metal, the quasi-
continuous density of states in the valence and the conduction bands splits into
discrete electronic levels, the spacing between these levels and the band gap
increasing with decreasing particle size.
lxv,lxvi
In the case of semiconductors, the phenomenon is slightly different, since a
band gap already exists in the bulk state. However, this band gap also increases
when the particle size is decreased and the energy bands gradually convert into
discrete molecular electronic levels.
lxvii
If the particle size is less than the De Broglie
wavelength of the electrons, the charge carriers may be treated quantum
mechanically as "particles in a box", where the size of the box is given by the
dimensions of the crystallites.
lxviii
In semiconductors, the quantization effect that
enhances the optical gap is routinely observed for clusters ranging from 1 nm to
almost 10 nm. Metal particles consisting of 50 to 100 atoms with a diameter between
1 and 2 nm start to loose their metallic behaviour and tend to become
semiconductors.
lxix
Particles that show this size quantization effect are sometimes
called Q-particles or quantum dots that we will discuss later in this chapter.
Certain particle sizes, determined by the so-called magic numbers, are more
frequently observed than others. This is not only due to the thermodynamic stability
of certain structures but is also related to the kinetics of particle growth. In
semiconductor nanoparticles such as CdS, the growth of the initially formed smallest
particles with an agglomeration number k occurs by combination of the particles.
Thus, particles so formed would have the agglomeration number of 2k, 3k and so
on.
lxx
Metals have a cubic or hexagonal close-packed structure consisting of one
central atom, which is surrounded in the first shell by 12 atoms, in the second shell by
42 atoms, or in principle by 10n
2
+2 atoms in the nth shell.
lxxi
For example, one of the
NANOMATERIALS
143
most famous ligand stabilized metal clusters is a gold particle with 55 atoms (Au
55
)
first reported by G. Schmid in 1981.
lxxii
7.2 Metals
7.2.1 Introduction
In drive towards nanotechnology, there is tremendous interest in understanding the
electric and thermal transport through nanowires. The questions that need to be
answered pertain to how the electric and heat transport laws that govern the
macroscopic systems get modified at the level of a few molecules or atoms. Wires
and contacts having average radius R, in the range of the Fermi wavelength
F
, can
indeed be fabricated and have shown unusual mechanical and electronic properties.
The effects of reduced size and dimensionality, and their quantum properties have
been investigated actively in the last decade. Ultrathin nanowires and even
monatomic chains suspended between two metal electrodes have been produced.
The discovery of scanning tunnelling microscopy (STM) and later atomic force
microscopy (AFM) has sparked a tremendous recent research effort on nanoparticles
because many properties have become directly observable. The tunnelling
microscope has also made it possible to produce atomic size contacts and atomic
wires only several angstroms in length. Electrical conduction and force have been
measured concomitantly in the course of stretching of nanowires. Initially, much of
the work was carried out on GaAsAlGaAs heterostructures which were grown to
contain a thin conducting layer at the interface. The conducting layer is treated as a
two-dimensional (2D) electron gas in which a narrow constriction of desired width w,
and length l, can be created by applying a negative gate voltage.
The conductivity , which relates the electric current density to the electric field by
j = E, is expressible in terms of the areal charge density
S
for 2D electron gas of
effective mass m through the equation
*
2
m
e
S
= (7. 142)
NANOMATERIALS
144
The experimental quantity of interest, however, is the conductance G = I /V, which is
the ratio of the total current I to the voltage drop V across the sample of length l in
the direction of current flow. For 2D, since I = wj, one can also write
l
w
G = (7. 143)
For a 3D conductor, this relationship is valid provided that w is replaced by the cross
sectional area A orthogonal to the current flow direction. Similar expressions are also
valid for the thermal transport of energy. The thermal conductance related to the
energy (or heat) current J
x
through a sample between two reservoirs is given by
K
x
= J
x
T , where T is the temperature difference. Depending on whether the
energy is carried by electrons (x = e) or by phonons (x = p) the thermal conductance
is identified as electronic K
e
or phononic K
p
. Here our focus is on quantum transport
through materials of very small dimensions.
Novel size-dependent effects emerge as w and l are reduced towards atomic
dimensions in the nanometre range. The relationship expressed by equation
l
w
G =
(7. 143) holds in the diffusive transport regime where both w and l are greater
than the mean free path. As the width of the constriction decreases, there comes a
point where quantum mechanics makes its presence known. The quantum
confinement of a carrier in a strip of width w leads to the discretization of energy
levels given by
n
= n
2
h
2
/(8m w
2
). The number of these w-dependent transverse modes, which are
occupied, determines the conductance. Thus, rather than a simple linear
dependence of G on w, quantum mechanics forces this indirect dependence on w.
As w is altered, the energy spectrum changes and so does the number of occupied
modes below the Fermi energy and hence the conductance.
Simple physical considerations show that the number of transverse occupied states,
N 2w/
F
, increases with the width of the constriction. Since all these modes can
contribute to the conductance, one still expects the conductance to increase linearly
with w even in the nanodomain. This is almost the case, except with one important
distinction. The width w can change continuously, at least in principle. But the
NANOMATERIALS
145
number of modes, being an integer, can only vary in discrete steps. The concept
becomes physically more transparent if we rewrite
n
~
F
E
n N . The highest occupied
n
= E
F
gives N from a simple counting of the number of discrete eigenvalues. Since
it is a rare coincidence for an eigenvalue to exactly align with the Fermi energy, N is
taken to be the integer which corresponds to the highest occupied level just below
E
F
. Thus the effect of quantum mechanics due to the reduced dimension w is to
cause the conductance to change in discrete steps in a staircase fashion. We have
yet to find the step height. This simple view is necessarily modified when other
factors explained later start to play a crucial role. Obviously, one then requires a
more detailed analysis.
The characteristic length that comes into play is
F
, since only those electrons having
energy close to the Fermi energy carry current at low temperatures. If w >> F, the
number of conducting modes is very large and so is the conductance. This is typical
in metals where
F
0.2 nanometres is very small and consequently the observation
of discrete conductance variation requires w
F
. That is why the discrete
conductance behaviour was first observed in semiconductor heterojunctions having
very low electron density and
F
two orders of magnitude larger than in metals.
The effect of reduced length l on the conductance is even more striking. If the ohmic
regime were to hold in equation
l
w
G = (7. 143), one would expect G to
increase without limit (or resistance to reach zero) as l was reduced towards zero.
We shall see that there would always be finite residual resistance. Heisenbergs
uncertainty principle defines a natural characteristic length as the mean free path l
e
. If
l < l
e
, carriers can propagate without losing their initial momentum, and this domain is
referred to as the ballistic transport regime. Landauer pointed out that conduction is
transmission. If we follow the standard definition, that the conductance is a measure
of current through a sample divided by the voltage difference,
V
G
= (7. 144)
NANOMATERIALS
146
The quantum of electrical conductance occurs naturally as a consequence of
Heisenbergs uncertainty principle. The current is given by the rate of charge flow
I = Q/t
transit
. Now charge is quantized in units of elementary charge, e. Hence in
the extreme quantum limit, setting Q = e and recognizing that the transit time
should then be at least in the range of time implied by Heisenbergs uncertainty
principle, we get
t
e
= (7. 145)
Combining equation
V
G
(7. 146)
Next invoking Heisenbergs uncertainty principle, E t >> h, the expression for
ballistic conductance including the spin degeneracy becomes G = 2e
2
/h in the ideal
case. It has a maximum value of 810
5
1
. This is the step height, the conductance
per transverse mode.
The corresponding resistance h/2e
2
has a value of 12.9 k (it suffices to call it
12345 for ease of remembrance) and is attributed to the resistance at the contacts
where the conductor is attached to the electrodes or electron reservoirs. Thus
classically G w. Quantum mechanically it increases in discrete steps. It jumps by
2e
2
/h as w increases enough to permit one more transverse mode to be occupied
and hence available for conduction. More recently, the electronic and transport
properties of metallic point contacts and wires having average dimension in the range
of metallic
F
produced by STM and also by mechanical break junctions have
displayed various quantum effects. The two-terminal electrical conductance, G, of the
wire showed a stepwise variation with the stretch. The sudden jumps of conductance
in the course of the stretch have been taken as the realization of quantized
conductance at room temperature and thus have created much popular interest.
Metallic carbon nanotubes providing two conductance channels at the Fermi level are
considered almost perfect one-dimensional (1D) conductors. These tubes, which can
NANOMATERIALS
147
be grown to several micrometres in length with a variety of diameters and chirality,
display unique electromechanical properties and dimensionality effects. They can be
either semiconductor or metal depending on their chirality and diameter. The
nanotubes are flexible, and at the same time are very strong with high yield strength.
Their strength far exceeds that of any other fibre. These properties are being actively
explored with a view to novel applications in nanotechnology. It has been
demonstrated experimentally that an individual nanotube can sustain very high
current density; first the current increases with the bias voltage and it eventually
saturates at a value which is very high for a nanowire. Recent studies have shown
that the electronic properties of semiconducting carbon nanotubes can be modified
variably and reversibly by radial deformation, a feature that may lead to some device
applications.
The thermal conductance may exhibit quantum features depending on the size of the
nanostructure that transfers the phononic energy. Owing to the finite level spacing of
vibrational frequencies, the phononic energy transfer through an electrically non-
conducting nano-object (i.e. a molecule, atomic chain, or a single atom) between two
reservoirs shows behaviour similar to the ballistic electron transport. Recent
theoretical studies indicate that the thermal conductance of an acoustic branch at low
temperature is independent of any material property, and is linearly dependent on
temperature. At low temperature the discrete vibrational frequency spectrum of a
chain gives rise to a sudden increase of the thermal conductance.
Study of quantum effects in 1D conductors has seen a tremendous explosion
because of growing interest in nanoscience and the quest for novel nanodevices. In
this section we will discuss some theoretical methods used in the calculation of
quantum transport then about the electrical conductance through metallic nanowires.
After a brief description of the structure and electronic properties, paragraph 7.3
shows recent experimental and theoretical investigation on the electronic transport
through singlewall carbon nanotubes (SWNTs).
Quantum transport of electrons
In 1957, Landauer introduced a novel way of looking at conduction. He taught us to
view the conduction as transmission and gave the famous formula, which has been a
breakthrough in the conductance phenomena and in the physics of mesoscopic
systems. According to Landauer, the conduction is a scattering event, and the
NANOMATERIALS
148
transport is the consequence of the incident current flux. On the basis of the self-
consistency arguments for reflections and transmissions, he derived his famous
formula for a one-dimensional conductor yielding the conductance
G = (2e/h) /R (7. 147)
Where T and R = 1 T are transmission and reflection coefficients, respectively.
Later, Sharvin investigated the electron transport through a small contact between
two free-electron metals. Since the length of the contact is negligible, i.e. l 0, the
scattering in the contact was absent and hence T 1. By using a semiclassical
approach he found an expression for the conductance,
G
S
= (2e
2
/h)(Ak
2
F
/4) (7. 148)
Which has come to be known as Sharvins conductance.
The original Landauer formula, i.e. equation G = (2e/h) /R (7. 147), and
Sharvins conductance have seemed to be at variance, since the former yields G
as R 0 in the absence of scattering. The confusion in the literature has been
clarified by recognizing the fact that the original formula is just the conductance of a
barrier in a 1D conductor. Engquist and Anderson pointed out an interesting feature
by arguing that the conductance has a close relation with the type of measurement.
Landauers formula was extended to the conductance measured between the two
outside reservoirs within which the finite-length conductor is placed. The
corresponding two-terminal conductance formula, G = (2e
2
/h)T , incorporates the
resistance due to the contacts to the reservoirs. Accordingly, for perfect transmission,
T 1, the conductance is still finite and equal to 2e
2
/h. The corresponding resistance
R = G1 = h/2e
2
is attributed to the resistance arising from the reflections at the
contacts where the finite-length conductor is connected to the reservoirs. This leaves
unanswered the question of how the contact resistance can be independent of the
type of contact. Also the infinite conductance within the conductor at perfect
transmission requires an in-depth discussion. Whether the finite conductance of an
ideal channel, 2e
2
/h, is due to the contact or due to the capacity of the current-
transporting state seems to be a matter of interpretation.
Bttiker developed the multi-probe generalization of the theory. During the last
decade, important progress has been made in mesoscopic physics; the quantum
transport of electrons through constrictions has been a subject of many in-depth
NANOMATERIALS
149
studies. In particular, the variation of the current through a point contact created by
an STM tip and the measurements of the conductance G through a narrow
constriction between two reservoirs of 2D electron gas in a high-mobility GaAs
GaAlAs heterostructure showing step behaviour with a step height of 2e
2
/h have led
to basic confirmation of Landauers fundamental ideas.
The step behaviour of G has been treated in several studies. Here we provide a
simple and formal explanation by using a 1D idealized uniform constriction having a
width w. The electrons are confined in the transverse direction and have states with
quantized energy
i
. They propagate freely along the length of the constriction. The
propagation constant for an electron with
i
< E
F
is
( )
2 / 1
* 2
(
=
i F i
E
m
(7. 149)
Whenever the width of the constriction increases by
F
/2, a new subband with
energy
i
+
i
/2m dips below the Fermi level and contributes to the current under
the small bias voltage V. The current is
( ) ( ) [ ]
=
+ =
j
i
F i F i i
E D V e E D e n
i
1
2
(7. 150)
Here j is the index of the highest subband that lies below the Fermi level, i.e.
i
E
F
+ eV and
j+1
> E
F
+ eV, and n
i
is the degeneracy of the state i. By
assuming perfect transmission in the absence of any contact resistance and barrier
inside the constriction, i.e. T = 1, and by expressing the group velocity
i
and the
density of states D
i
() in terms of the subband energy =
i
+
i
/2m* and dividing
I by V we obtain
=
=
j
i
i
n
h
e
G
1
2
(7. 151)
Accordingly, each current-transporting state with energy in the range
E
F
< < E
F
+ eV contributes to G an amount 2en
i
/h.
For a uniform, infinite wall constriction, the states are non-degenerate, i.e. n
i
= 1, and
hence the increase of w by
F
/2 causes G to jump by 2e/h. As a result, the G(w)
curve exhibits a staircase structure. Since the level spacing in the transverse
quantization is rather small for the low electron density in the 2D electron gas system,
NANOMATERIALS
150
and
F
-2
, the sharp step structure is likely to be smeared out at T 10 K or at
finite bias voltage.
For a finite-length constriction the scattering at the contacts to the reservoirs (i.e. the
contact resistance) and at the non-uniformities or the potential barriers inside the
constriction affect the transmission. Then, the conductance of a single channel
expressed as
G = (2e/h)T (7. 152)
may deviate from the perfect quantized values. As a result, the effects, such as the
contact and potential barrier, surface roughness, impurity scattering, cause the sharp
step structure of G(w) to smear out. The local widening of the constriction or impurity
potential can give rise to quasi-bound (0D) states in the constriction and to resonant
tunnelling effects.
The length of the constriction, l, is another important parameter. In order to get sharp
step structure, l has to be greater than
F
, but smaller than the electron mean free
path l
e
; G(w) is smoothed out in a short constriction (l <
F
). Therefore, G(w) exhibits
sharp step structure if the constriction is uniform, and w
F
and
F
<< l < le. On the
other hand, the theoretical studies predict that the resonance structures occur on the
flat plateaus due to the interference of waves reflected from the abrupt connections
to the reservoirs. The stepwise variation of G with w or E
F
has been identified as the
quantization of conductance. This is, in fact, the reflection of the quantized
constriction states in the electrical conductance.
We now extend the above discussion to analyse the ballistic electron conductance
through a point contact or a nanowire, in which the electronic motion is confined in
two dimensions, but freely propagates in the third dimension. The point contact (or
quantum contact) created by the indentation of the STM, a nanowire (or a connective
neck) that is produced by retracting the tip from an indentation and also a metallic
SWNT are typical systems of interest. Nanowires created by STM are expected to be
round (though not perfectly cylindrical) and have radius R
F
at the neck. The neck
is connected to the electrodes by horn-like ends, and hence the radius increases as
one goes away from the neck. An extreme case for l 0 is Sharvins conductance,
G
S
= (2e/h)(R/
F
), with contact radius R, where the step structure is almost
smeared out and plateaus disappear. In the quantum regime, where the cross
NANOMATERIALS
151
section A
F
, G
S
should vanish when A is smaller than a critical cross section set
by the uncertainty principle. As l increases, a stepwise behaviour for G
S
develops.
The real point contacts and metallic nanowires differ from the perfect uniform
constriction in the way that the electrons are quantized in the neck. The size and form
of the neck, and also its electronic potential, influence the quantization of the
electronic states and the level spacing. The form of the neck may be irregular if the
cross section is large and comprises few atoms, but shows cylindrical symmetry for
an undeformed SWNT or metallic nanowire having a single atom at the neck. The
cylindrical symmetry becomes apparent with the degeneracy n
i
= 1, 2, 2, 1, . . for
i = 1, 2, 3, 4, . . .. For the most general constriction the expression for the two-
terminal conductance in equation
V
G
l
(z) with n = n
x
+ n
y
and energy [ ]
2 / 1
,
* / ) ( 2 ) 1 ( ) ( m z n z
l l n
+ = . The propagation
constant is given by
[ ] ) ( ) (
2
) (
,
z z E
m
z
l n l n
=
(7. 156)
The conductance of the wire is calculated by using the transfer-matrix method, where
V
l
(, z) is divided into discrete segments of equal widths, and in each segment
z
j
< z < z
j+1
, the average values are used. The coefficients A
nki
and B
nki
are
determined from the multiple boundary matching conditions. The total conductance is
calculated by integrating the expectation value of the current operator over the Fermi
surface.
A constriction having more general geometry and potential can be treated by
expressing the current-transporting state in the Laue expansion
=
n
j n
iG
j
z f e z
n
) ( ) , (
,
(7. 157)
in terms of transverse reciprocal-lattice vectors,G
i
. This method is convenient when
using the SCF potentials obtained from supercell calculations. Furthermore, the
conductance of rather long constrictions calculated by using the recursion scheme do
not diverge. The Schrdinger equation for
j
is discretized in the segments
z
j
< z < z
j+1
, and converted to a matrix equation. The details of the calculations can
be found elsewhere.
NANOMATERIALS
153
The Greens function method has been used extensively to study the electrical
conduction through nanowires and also SWNTs. In the self-consistent calculation,
where the wave functions of the bare electrodes,
0
, are evaluated from a realistic
Hamiltonian, the wave functions of the complete system including the constriction,
C
, are put into the LippmannSchwinger form:
, where x
is a
unit vector. To determine the electromagnetic field surrounding the particle, we solve
LaPlaces equation (the fundamental equation of electrostatics),
2
= 0, where is
the electric potential and the field E is related to by E = - . In developing this
solution, we apply two boundary conditions: (i) that is continuous at the sphere
surface and (ii) that the normal component of the electric displacement D is also
continuous, where D = E.
It is not difficult to show that the general solution to the LaPlace equation has angular
solutions which are just the spherical harmonics. In addition, the radial solutions are
of the form r
l
and r
-(l +1)
, where l is the familiar angular momentum label (l = 0, 1, 2,...)
NANOMATERIALS
165
of atomic orbitals. If we restrict our considerations for now to just the l = 1 solution
and if E
o
is in the x direction, the potential is simply = A r sin cos inside the
sphere (r < a) and = (-E
o
r + B/r) sin cos outside the sphere (r > a), where A
and B are constants to be determined. If these solutions are inserted into the
boundary conditions and the resulting is used to determine the field outside the
sphere, E
out
, we get
(
+ + = ) (
3
5 3
z z y y x x
r
x
r
x
E x E E
o o out
(7. 162)
Where R is the sphere polarizability and x
, y
, and z
+ + = ) (
3
5 3
z z y y x x
r
x
r
x
E x E E
o o out
(7. 162) is
the applied field and the second is the induced dipole field (induced dipole moment =
E
o
) that results from polarization of the conduction electron density.
For a sphere with the dielectric constants indicated above, the LaPlace equation
solution shows that the polarizability is
3
a g
d
= (7. 163)
with
0
0
2 +
=
i
i
d
g (7. 164)
Although the dipole field in eq
(
+ + = ) (
3
5 3
z z y y x x
r
x
r
x
E x E E
o o out
(7. 162) is
that for a static dipole, the more complete Maxwell equation solution shows that this
is actually a radiating dipole, and thus, it contributes to extinction and Rayleigh
scattering by the sphere. This leads to extinction and scattering efficiencies given by
) Im( 4
d ext
g x Q = (7. 165)
2
4
3
8
d sca
g x Q = (7. 166)
where x = 2a(
0
)
1/2
/. The efficiency is the ratio of the cross-section to the
geometrical cross-section a. Note that the factor g
d
from eq
0
0
2 +
=
i
i
d
g (7.
NANOMATERIALS
166
164) plays the key role in determining the wavelength dependence of these cross-
sections, as the metal dielectric constant
i
is strongly dependent on wavelength.
7.2.3.3 Quadrupole plasmon Resonances
For larger particles, higher multipoles, especially the quadrupole term (l=2) become
important to the extinction and scattering spectra. Using the same notation as above
and including the l=2 term in the LaPlace equation solution, the resulting field outside
the sphere, E
out
, now can be expressed as
( ) ( ) ( )
(
+ +
+
+ + + + = z xz y y x x
r
z
r
z z x x
E z z y y x x
r
x
r
x
E z z x x ikE x E E
o o o o out
5 3
7 5 5 3
(7. 167)
and the quadrupole polarizability is
5
a g
d
= (7. 168)
0
0
2 / 3 +
=
i
i
d
g (7. 169)
Note that the denominator of eq
0
0
2 / 3 +
=
i
i
d
g (7. 169) contains the factor 3/2
while in eq
0
0
2 +
=
i
i
d
g (7. 164) the corresponding number is 2. These factors
arise from the exponents in the radial solutions to LaPlaces equation, e.g., the
factors r
l
and r
-(l +1)
that were discussed above. For dipole excitation, we have l=1,
and the magnitude of the ratio of the exponents is (l + 1)/l = 2, while for quadrupole
excitation (l + 1)/l = 3/2. Higher partial waves work analogously.
Following the same derivation, we get the following quasistatic (dipole + quadrupole)
expressions for the extinction and Rayleigh scattering efficiencies:
( )
(
+ + = 1
30
12
Im 4
i q d ext
x
g
x
g x Q (7. 170)
)
`
+ + =
2
4
2
4
2
4
1
900 240 3
8
i q d ext
x
g
x
g x Q (7. 171)
NANOMATERIALS
167
7.2.3.4 Extinction for silver spheres
We now evaluate the extinction cross-section using the quasistatic expressions, eqs
4, 5 and 9,10 as well as the exact (Mie) theory. We take dielectric constants for silver
that are plotted in Figure 7-6 a and the external dielectric constant is assumed to be 1
(i.e., a particle in a vacuum). The resulting efficiencies for 30 and 60 nm spheres are
plotted in Figure 7-6 b,c, respectively.
Figure 7-6: (a) Real and imaginary part of silver dielectric constants as function of wavelength. (b)
Extinction efficiency, i.e., the ratio of the extinction cross-section to the area of the sphere, as obtained
from quasistatic theory for a silver sphere whose radius is 30 nm. (c) The corresponding efficiency for
a 60 nm particle, including for quadrupole effects, and correcting for finite wavelength effects. In b and
c, the exact Mie theory result is also plotted. [K. Lance Kelly, E. Coronado, L.L. Zhao, G.C. Schatz,
The Optical Properties of Metal Nanoparticles: The Influence of Size, Shape, and Dielectric
Environment, 2003, Journal of Physical Chemistry]
The cross-section in Figure 7-6 b shows a sharp peak at 367nm, with a good match
between quasistatic and Mie theory. This peak is the dipole surface plasmon
resonance, and it occurs when the real part of the denominator in eq
0
0
2 +
=
i
i
d
g
NANOMATERIALS
168
(7. 164) vanishes, corresponding to a metal dielectric constant whose real part
is -2. For particles that are not in a vacuum, the plasmon resonance condition
becomes
Re
i
/
0
= -2, and because the real part of the silver dielectric constant decreases with
increasing wavelength (Figure 7-6 a), the plasmon resonance wavelength for
0
> 1 is
longer than in a vacuum. The plasmon resonance wavelength also gets longer if the
particle size is increased above 30 nm. This is due to additional electromagnetic
effects that will be discussed later.
Figure 7-6 c presents the l = 2 quasistatic (dipole + quadrupole) cross-section as well
as the full Mie theory result for a radius 60 nm sphere. The quasistatic result also
includes a finite wavelength correction. We see that the dipole plasmon wavelength
has shifted to the red, and there is now a distinct quadrupole resonance peak at 357
nm. This quadrupole peak occurs when the real part of the denominator in eq 8
vanishes, corresponding to a metal dielectric constant whose real part is -3/2. For a
sphere of this size, there are clear differences between the quasistatic and the Mie
theory results; however, the important features are retained. Although Mie theory is
not a very expensive calculation, the quasistatic expressions are convenient to use
when only qualitative information is needed.
7.2.3.5 Electromanetic fields for spherical particles
So far, we have emphasized the calculation of extinction and Rayleigh scattering
cross-section; however, for certain properties, such as SERS and hyperRaman
scattering (HRS) intensities, it is the electromagnetic field at or near the particle
surfaces that determines the measured intensity. Thus, if E() is the local field for
frequency , then the SERS intensity is determined by
2 2
) ' ( ) ( E E where is
the Stokes-shifted frequency and the brackets are used to denote an average over
the particle surface. The HRS intensity is similarly (but approximately) determined by
2 4
) 2 ( ) ( E E . Also, when one makes an aggregate or array of metal
nanoparticles, the interaction between the particles is determined by the polarization
induced in each particle due to the fields E arising from all of the other particles.
At the dipole (dipole + quadrupole) level, the field outside a particle is given by eq
(
+ + = ) (
3
5 3
z z y y x x
r
x
r
x
E x E E
o o out
+ +
+
+ + + + = z xz y y x x
r
z
r
z z x x
E z z y y x x
r
x
r
x
E z z x x ikE x E E
o o o o out
5 3
7 5 5 3
(7. 167). These expressions determine the near-fields at the particle surfaces quite
accurately for small enough particles; however, the field beyond 100 nm from the
center of the particle exhibits radiative contributions that are not contained in these
equations. To describe these, we need to replace the dipole or quadrupole field by its
radiative counterpart. In the case of the dipole field, this is given by
( )
5 3
) ( 3
) 1 (
r
P r r P r
ikr e
r
p r r
e k E
ikr ikr
dipole
+
= (7. 172)
where P is the dipole moment. Note that this reduces to the static field in eq
(
+ + = ) (
3
5 3
z z y y x x
r
x
r
x
E x E E
o o out
+ =
f
ik fk photon i f
dE E E E E E w ) , ( ) (
2 2
(7. 177)
where denotes the joint density of valence and conduction band electronic states
with the same k. The matrix element is factored so that the dipole operator acts on
the electronic wave function within the unit cell, and the phase factors between unit
cells, exp(ik
v
r) and exp(ik
c
r), must overlap for the transition to be allowed (k
v
= k
c
).
In this case, the transition is dipole allowed within each unit cell. Further, the
transition dipole moment points in the same direction from unit cell to unit cell.
Thus, the electric field of light, which is very long wavelength compared to the size of
the unit cell, can drive all the transition dipoles in phase, and overall the transition is
allowed. In the language of solid state physics, this is a vertical transition, with k=0,
and hence CdSe is called a direct band gap material. At low temperature, the
electronic transition will be accompanied by vibrational excitation of totally symmetric
NANOMATERIALS
189
modes, in accordance with the Franck-Condon factors. Past the absorption threshold,
the unit cell matrix element for absorption is largely unchanged, and the absorption
efficiency rises steeply at higher incident photon energies because of the large
increase in the joint density of states for vertical transitions. Numerous experiments
show that quantum confinement dominantly affects the absorption spectrum by
changes in the envelope functions and the density of states as a function of the
energy, with relatively less effect on the intrinsic, unit cell based, matrix element for
absorption.
7.4.2.2 Quantization and energy level spacing
As the size is reduced, the electronic states may be viewed as superpositions of bulk
states. Hence, there is a shift to higher energy, the development of discrete features
in the spectra, and concentration of the oscillator strength into just a few transitions.
Qualitatively, all of these effects can be readily observed in the spectra of Figure
7-21, which show data for CdSe.
Figure 7-21: Optical absorption vs size for CdSe nanocrystals shows the shift to higher energy in
smaller sizes, as well as the development of discrete structure in the spectra and the concentration of
oscillator strength into just a few transitions. [A.P. Alivisatos, Perspectives on the physical chemistry of
semiconductor nanocrystals, 1996, Journal of Physical Chemistry]
The quantitative analysis of these spectra remains a difficult subject, for several
reasons:
The foregoing picture is a single-particle one and does not include the
substantial effects of correlation. In molecules this is analogous to trying to use
the highly approximate molecular orbital theory, instead of more advanced
quantum chemistry methods. Regrettably, the nanocrystals are too large to
NANOMATERIALS
190
describe using even moderately advanced methods that are routinely applied
to small molecules.
Further, in CdSe at least, the large atomic number of the Se ensures that the
coupling between spin and orbital momenta is very strong in the valence
bands (p bands). This coupling is in the j-j, and not the Russell-Saunders, L-S,
coupling regime. When translational symmetry is removed, the mixing of k
vectors can also result in different bands mixing together.
The shape of the crystallites, which is regular (tetrahedral, hexagonal prisms),
or spherical, or ellipsoidal, will determine the symmetry of the nanocrystals
and will influence the relative spacing of the levels.
Finally, surface energy levels are completely excluded from this simple
quantum confinement picture. Yet it seems apparent that surface states near
the gap can mix with interior levels to a substantial degree, and these effects
may also influence the spacing of the energy levels.
Theoretical approaches that directly include the influence of the surface, as well as
electronic correlation, are being developed rapidly.
7.4.3 Electrical properties
We talk about single-electronics whenever it is possible to control the
movement and position of a single or small number of electrons. To understand how
a single electron can be controlled, one must understand the movement of electric
charge through a conductor. An electric current can flow through the conductor
because some electrons are free to move through the lattice of atomic nuclei. The
current is determined by the charge transferred through the conductor. Surprisingly
this transferred charge can have practically any value, in particular, a fraction of the
charge of a single electron. Hence, it is not quantized.
This, at first glance counterintuitive fact, is a consequence of the displacement
of the electron cloud against the lattice of atoms. This shift can be changed
continuously and thus the transferred charge is a continuous quantity (see left side of
Figure 7-22).
NANOMATERIALS
191
+ + + + + +
- - - - - -
+ + + + + +
- - - - - -
+ + + + + +
- - - - - -
+ + +
+
+ + +
+ +
+ + +
+
+ + +
+ +
t
u
n
n
e
l
j
u
n
c
t
i
o
n
Figure 7-22: The left side shows, that the electron cloud shift against the lattice of atoms is not
quantized. The right side shows an accumulation of electrons at a tunnel junction.
If a tunnel junction interrupts an ordinary conductor, electric charge will move
through the system by both a continuous and discrete process. Since only discrete
electrons can tunnel through junctions, charge will accumulate at the surface of the
electrode against the isolating layer, until a high enough bias has built up across the
tunnel junction (see right side of Figure 7-22). Then one electron will be transferred.
Likharev has coined the term `dripping tap' as an analogy of this process. In other
words, if a single tunnel junction is biased with a constant current I, the so called
Coulomb oscillations will appear with frequency f = I/e, where e is the charge of an
electron (see Figure 7-23).
I
I
conduction current
Figure 7-23: Current biased tunnel junction showing Coulomb oscillations.
Charge continuously accumulates on the tunnel junction until it is energetically
favorable for an electron to tunnel. This discharges the tunnel junction by an
elementary charge e. Similar effects are observed in superconductors. There, charge
carriers are Cooper pairs, and the characteristic frequency becomes f = I/2e, related
to the so called Bloch oscillations.
NANOMATERIALS
192
Figure 7-24: The electron-box can be filled with a precise number of electrons.
The current biased tunnel junction is one very simple circuit, that shows the
controlled transfer of electrons. Another one is the electron-box (see Figure 7-24).
A particle is only on one side connected by a tunnel junction. On this side electrons
can tunnel in and out. Imagine for instance a metal particle embedded in oxide, as
shown in Figure 7-25.
oxide
metal
particle
Figure 7-25: Metal particle embedded in oxide. Tunneling is only possible through the thin top layer
of oxide.
The top oxide layer is thin enough for electrons to tunnel through. To transfer
one electron onto the particle, the Coulomb energy E
C
= e
2
/2C, where C is the
particles capacitance, is required. Neglecting thermal and other forms of energy, the
only energy source available is the bias voltage V
b
. As long as the bias voltage is
small enough, smaller than a threshold V
th
= e/C, no electron can tunnel, because
not enough energy is available to charge the island. This behavior is called the
Coulomb blockade. Raising the bias voltage will populate the particle with one, then
two and so on electrons, leading to a staircase-like characteristic.
NANOMATERIALS
193
It is easily understandable, that these single-electron phenomena, such as
Coulomb oscillations and Coulomb blockade, only matter, if the Coulomb energy
is bigger than the thermal energy. Otherwise thermal fluctuations will disturb the
motion of electrons and will wash out the quantization effects. The necessary
condition is,
2
2
c B
e
E k T
C
= > (7. 178)
Where k
B
is Boltzmann's constant and T is the absolute temperature. This means
that the capacitance C has to be smaller than 12 nF for the observation of charging
effects at the temperature of liquid nitrogen and smaller than 3 nF for charging effects
to appear at room temperature. A second condition for the observation of charging
effects is, that quantum fluctuations of the number of electrons on an island must be
negligible. Electrons need to be well localized on the islands. If electrons would not
be localized on islands one would not observe charging effects, since islands would
not be separate particles but rather one big uniform space. The charging of one
island with an integer number of the elementary charge would be impossible,
because one electron is shared by more than one island. The Coulomb blockade
would vanish, since no longer would a lower limit of the charge, an island could be
charged with, exist. This leads to the requirement that all tunnel junctions must be
opaque enough for electrons in order to confine them on islands. The `transparency'
of a tunnel junction is given by its tunnel resistance R
T
which must fulfill the following
condition for observing discrete charging effects, where h is Planck's constant. This
should be understood as an order-of-magnitude measure, rather than an exact
threshold.
2
25813
T
h
R
e
> = (7. 179)
Therefore, these effects are experimentally verifiable only for very small high-
resistance tunnel junctions, meaning small particles with small capacitances and/or
very low temperatures. Advanced fabrication techniques, such as the production of
granular films with particle sizes down to 1 nm, and deeper physical understanding
allow today the study of many charging effects at room temperature.
Based on the Coulomb blockade many interesting devices are possible, such as
precise current standard,
lxxiii
very sensitive electrometers,
lxxiv
logic gates
lxxv
and
NANOMATERIALS
194
memories
lxxvi
with ultra low power consumption, down-scalability to atomic
dimensions, and high speed of operation. Altogether, single-electronics will bring new
and novel devices and is a very promising candidate to partly replace MOS
technology in the near future.
7.4.3.1 Characteristic Energies
Charging an island with a small number of electrons is associated with two
kinds of energies:
1. One kind arises from electron-electron interaction
2. and the other is linked to the spatial confinement of electrons.
These energies are only significant for very small structures in the nanometer
regime. Otherwise they are masked by thermal fluctuations.
We discuss two other important characteristic energies, the work done by the voltage
sources and Helmholtz's free energy, and the role they play in single-electron device
behavior.
Electron Electron Interaction
An entirely classical model for electron-electron interaction is based on the
electrostatic capacitive charging energy. The interaction arises from the fact, that for
every additional charge dq which is transported to a conductor, work has to be done
against the field of already present charges residing on the conductor. Charging a
capacitor with a charge q requires energy, where C is the total capacitance of the
capacitor. A first simple model for an island is a conducting sphere which has a self-
capacitance of (see the next table)
2
4
c
q
E
d
= (7. 180)
Where is the permittivity of the surrounding dielectric material and d is the
diameter of the sphere. Inserting
2
2
c B
e
E k T
C
= > (7. 178) in
2
4
c
q
E
d
= (7.
180)) yields :
2 C d = (7. 181)
NANOMATERIALS
195
An indirect proportionality to the size of the particle. Figure 7-26 to Figure 7-28
compare this electrostatic energy to the other characteristic energies involved in
charging an island with a small number of electrons.
Change of E
F
for 1st additional electron in Si
Electrostatic charging energy e
2
/2c
thermal energy at 300 k
10
-1
10
0
10
1
10
2
10
-3
10
- 2
10
-1
10
1
10
0
Particle diameter [nm]
Change of E
F
for 2nd additional electron in Si
Change of E
F
for 3rd additional electron in Si
E
n
e
r
g
y
[
e
V
]
Figure 7-26: Comparison of the change in Fermi energy for the addition of one, two, and three
electrons in Si to the electrostatic charging energy.
10
1
Change of E
F
for 1st additional electron in Si
Electrostatic charging energy e
2
/2c
10
-1
10
0
10
1
10
2
10
-3
10
- 2
10
-1
10
0
Particle diameter [nm]
Change of E
F
for 2nd additional electron in Si
E
n
e
r
g
y
[
e
V
]
Figure 7-27: Comparison of the change in Fermi energy for the addition of one electron in Si and Al.
The values below 1 nm particle diameter are questionable (see later)
For the sake of simplicity, the relative permittivity was assumed to be 10. A more
accurate model for the capacitance of an island is the capacitance of two or more
spheres in a line. A calculation for characteristic spatial dimensions shows an
increase of about 15% of the total capacitance.
For a system of N conductors the charge on one conductor may be written as
NANOMATERIALS
196
1
N
i ij j
j
q C
=
=
( 1, 2, 3,......) i = (7. 182)
Where the
ij
C denote the elements of the capacitance matrix. The diagonal elements
of the capacitance matrix
ij
C are the total capacitances of the conductors i and the
off-diagonal elements
ij
C are the negative capacitances between conductor i and j.
The electrostatic charging energy is in analogy to
2
2
c B
e
E k T
C
= > (7. 178)
( )
1
1 1 1
1 1
2 2
N N N
c i i i j
ij
i i j
E q C q q
= = =
= =
(7. 183)
where C
-1
is the inverse of the capacitance matrix.
Systems with sufficiently small islands are not adequately described with this
model. They exhibit a second electron-electron interaction energy, namely the
change in Fermi energy, when charged with a single electron. One must distinguish
between metals and semiconductors, because of their greatly different free carrier
concentrations and the presence of a band gap in the case of semiconductors.
Metals have typically a carrier concentration of several 10
22
/cm
3
and undoped Si
about 10
10
/cm
3
at room temperature. The low intrinsic carrier concentration for
semiconductors is in reality not achievable, because of the ubiquitous impurities
which are ionized at room temperature. These impurities supply free charge carriers
and may lift the concentration, in the case of Si, to around 10
14
/cm
3
, still many orders
of magnitude smaller than for metals. Clearly, doping increases the carrier
concentration considerably, and can lead in the case of degenerated semiconductors
to metal like behavior. We are considering only undoped semiconductors to
emphasize the difference to metals.
The dependence of the Fermi energy E
F
on carrier concentration n can be expressed
as below (N. W. Ashcroft and N. D. Mermin, Solid State Physics, W. B. Saunders
Company, international edition, 1976). For metals it is given by
( )
2
2 3
2
*
3
2
F
E n
m
=
(7. 184)
and for semiconductors it is
NANOMATERIALS
197
,
2sinh
F F i net
i B
E E n
n k T
= (7. 185)
Where n
net
is the net carrier concentration defined as n
net
= n-p. Figure 7-26
compares the change in the Fermi energy for the addition of one, two, and three
electrons in Si to the electrostatic charging energy.
The first electron causes the biggest change in Fermi energy. The second and
third electron contribute a change smaller than the thermal energy for room
temperature. Hence, for single-electron devices one will try to build structures with
few free carriers, to achieve a high change in energy when charging an island.
Figure 7-27 compares Si with Al. It is clearly visible that the change is much
bigger in Si. This is caused by two reasons. First, the difference in free carrier
concentration. In metals much more carriers are available and therefore the addition
of one electron has not such a big impact anymore. Second in metals an additional
electron finds a place slightly above the Fermi level. In semiconductors due to the
energy gap the electron needs to be inserted considerably above the Fermi level.
Quantum Confinement Energies
With decreasing island size the energy level spacing of electron states
increases indirectly proportional to the square of the dot size. Taking an infinite
potential well as a simple model for a quantum dot, one calculates by solving
Schrdinger's equation the energy levels to
2
*
1
2
N
N
E
m d
| |
=
|
\
(7. 186)
Because of the lower effective mass in Si compared to Al, the quantum confinement
energy is bigger in Si (see Figure 7-28).
NANOMATERIALS
198
Confinement energy for Si
Confinement energy for metals
Electrostatic energy e
2
/2c
thermal energy at 300 k
thermal energy at 77 k
10
-1
10
0
10
1
10
2
10
-3
10
-2
10
-1
10
1
10
0
Particle diameter [nm]
Figure 7-28: Comparison of the quantum confinement energy in Si and Al.
For very small particle diameters (around 1 nm) the formulas
( )
2
2 3
2
*
3
2
F
E n
m
=
(7. 184),
,
2sinh
F F i net
i B
E E n
n k T
| |
=
|
\
(7. 186) cease to be valid. The concept of the effective mass is based on a
periodic lattice and starts to break down for small islands. One would have to apply
cluster theory to calculate the energy levels more precisely. In reality slightly smaller
confinement energies than predicted in
2
*
1
2
N
N
E
m d
| |
=
|
\
(7. 192)
The current I cannot exceed e/t since for moderate bias and temperature at most one
extra electron resides on the island at any time. The energy uncertainty of the
electron E is no larger than the applied voltage V
b
.
b
E eV < (7. 193)
Inserting t=e/I results in
2
b
T
V h
R
I e
= >> (7. 194)
NANOMATERIALS
201
In fact, more rigorous theoretical studies of this issue have supported this conclusion.
Experimental tests have also shown this to be a necessary condition for observing
single-electron charging effects.
lxxvii
7.4.3.4 The Double Tunnel Junction
Consider two tunnel junctions in series biased with an ideal voltage source as shown
in Figure 7-29. The charges on junction one, junction two, and on the whole island
can be written as
1 1 1 2 2 2
2 1 0 0
, , q CV q C V and
q q q q ne q
= =
= + = +
(7. 195)
respectively, with n
1
the number of electrons that tunneled through the first junction
entering the island, n
2
the number of electrons that tunneled through the second
junction exiting the island, and n=n
1
-n
2
the net number of electrons on the island.
A background charge q
0
produces generally a non-integer charge offset. The
background charge is induced by stray capacitances that are not shown in the circuit
diagram Figure 7-29 and impurities located near the island, which are practically
always present. Using
1 1 1 2 2 2
2 1 0 0
, , q CV q C V and
q q q q ne q
= =
= + = +
(7. 195) and V
b
=V
1
+V
2
gives
2 0 1 0
1 2
1 2
,
,
b b
C V ne q CV ne q
V V
C C
where C C C
+ +
= =
= +
(7. 196)
NANOMATERIALS
202
Figure 7-29: Two tunnel junctions in series biased with an ideal voltage source. The background
charge q
0
is non-integer, and n
1
and n
2
denote the number of tunneled electrons through junction one
and junction two, respectively.
With
2 0 1 0
1 2
1 2
,
,
b b
C V ne q CV ne q
V V
C C
where C C C
+ +
= =
= +
(7. 196) the electrostatic energy
stored in the double junction is
( )
2
2 2 2
1 2 0 1 2
1 2
2 2 2
b
C
C C V ne q q q
E
C C C
+
= + = (7. 197)
In addition, to get the free energy one must consider, the work done by the
voltage source. If one electron tunnels through the first junction the voltage source
has to replace this electron -e, plus the change in polarization charge caused by the
tunneling electron. V
1
changes according to
2 0 1 0
1 2
1 2
,
,
b b
C V ne q CV ne q
V V
C C
where C C C
+ +
= =
= +
(7.
196) by e C
The charge q
1
gets smaller, which means that the voltage source `receives'
polarization charge. The total charge that has to be replaced by the voltage source is
therefore
2
eC C
and the work done by the voltage source in case electrons tunnel
through junction one and junction two is accordingly
1 2
1
2 1
2
b
b
n eV C
W
C
and
n eV C
W
C
=
=
(7. 198)
The free energy of the complete circuit is
( ) ( ) ( )
{ }
2
2
1 2 1 2 0 1 2 2 1
1 1
,
2
c b b
F n n E W C C V ne q eV C n C n
C
(
= = + + +
(
(7. 199)
At zero temperature, the system has to evolve from a state of higher energy to
one of lower energy. At non-zero temperatures transitions to higher energy states are
NANOMATERIALS
203
possible, but have exponentially reduced probability. The change in free energy for
an electron tunneling through junction one and two is given by
( ) ( ) ( )
1 1 2 1 2 2 0
1, ,
2
b
e e
F F n n F n n V C ne q
C
(
= = +
(
(7. 200)
( ) ( ) ( )
2 1 2 1 2 1 0
, 1 ,
2
b
e e
F F n n F n n V C ne q
C
(
= = +
(
(7. 201)
The probability of a tunnel event will only be high, if the change in free energy is
negative - a transition to a lower energy state. The leading term in equations above
causes F to be positive until the magnitude of the bias voltage V
b
exceeds a
threshold which depends on the smaller of the two capacitances. This is the case for
all possible transitions starting from an uncharged island, n=0 and q
0
=0. For
symmetric junctions (C
1
= C
2
) the condition becomes
b
V e C
+ +
= =
= +
(7. 196), the new voltages across the
junctions are
( )
2
1
1
2
,
g b g g o
b g g o
C C V C V ne q
V
C
CV C V ne q
V
C
+ +
=
+ +
=
(7. 205)
with
1 2 g
C C C C
= + +
The electrostatic energy has to include also the energy stored in the gate capacitor,
and the work done by the gate voltage has to be accounted for in the free energy.
The change in free energy after a tunnel event in junctions one and two becomes
( ) { }
{ }
1 2 0
1 1 0
2
g b g g
b g g
e e
F C C V C V ne q
C
e
F V C V C ne q
C
(
= + +
(
= + +
(7. 206)
At zero temperature only transitions with a negative change in free energy,
1
0 F <
Or,
2
0 F <
are allowed. These conditions may be used to generate a stability plot in the V
b
-V
g
plane, as shown in Figure 7-35.
Figure 7-35: Stability plot for the SET transistor. The shaded regions are stable regions.
NANOMATERIALS
208
The shaded regions correspond to stable regions with an integer number of
excess electrons on the island, neglecting any non-zero background charge. If the
gate voltage is increased, and the bias voltage is kept constant below the Coulomb
blockade,
b
V e C
0
= 4 x 10
-7
H/m is the permeability of free space.
We can also define the magnetic field parameters by this way:
B =
0
H(1
+ ) =H (8. 208)
Where =
0
(1
+ ) is the permeability of the material and = M/H is the
susceptibility of the material.
NANOMATERIALS
220
Finally, we need to define the magnetic moment,
m
, through the following equation:
M =
m
/ V (8. 209)
It needs to be noted that in magnetic theory several unit systems are
commonly in use. Historically, workers in magnetic materials have used the cgs
(centimetre, gram, second) or Gaussian system. More recently attempts have been
made to change over the SI system. As you can see on the next figure, two SI
systems are used; we will focus on the Sommerfield convention.
Table 8-1 : Magnetic units. A is ampere, cm is centimeter, m is meter, emu is electromagnetic unit, B
is magnetic induction, H is magnetic field strength, M is magnetization of a substance per unit volume.
[K. J Klabunde, Nanoscale Materials in Chemistry, Wiley-Interscience, 2001]
8.1.2 Phenomena
8.1.2.1 Atomic origins
Magnetic field is relative motion of an electric charge and the observer. Thus
magnetism is a result of moving charges. From an atomic view of matter, there are
two electronic motions:
- the orbital motion of the electron (revolution of electron around the core)
- and the spin motion of the electron (revolve around its own axis).
These two electron motions are the source of macroscopic magnetic phenomena in
materials.
8.1.2.2 Classification of magnetic phenomena
We use the word magnetism very loosely when implying the mutual magnetic
attraction of pieces of iron. There are, however, several other classes of magnetic
materials that differ in kind and degree in their mutual interaction. We shall
distinguish in the following between:
NANOMATERIALS
221
- Ferromagnetism: Imagine the external field strength is increased (in a
ferromagnetic solid), then the magnetization rises at first slowly and then more
rapidly. Finally, M levels off and reaches a constant value, called the saturation
magnetization M
s
. When H is reduced to zero, the magnetization retains a positive
value, called the remanent magnetization or remanence M. It is this retained
magnetization which is utilized in permanent magnets. The remanent magnetization
can be removed by reversing the magnetic field strength to a value H
c
, called the
coercive field.
- Paramagnetism: Without an external magnetic field, orbiting electrons
magnetic moments are randomly oriented and thus they mutually cancel one another.
As a result, the net magnetization is zero. However, when an external field is applied
the individual magnetic vectors tend to turn into the field direction, which is
counteracted by thermal agitation alignment (see the Curie law).
- Diamagnetism: may be explained by postulating that the external
magnetic field induces a change in the magnitude of inner-atomic currents, that is,
the external field accelerates or decelerates the orbiting electrons, in order that their
magnetic moment is in the opposite direction from the external magnetic field.
- Antiferromagnetism: Antiferromagnetic materials exhibit, just as
ferromagnetics a spontaneous alignment of moments below a critical temperature.
However, the responsible neighbouring atoms in antiferromagnetics are aligned in an
antiparallel fashion. Antiferromagnetic materials are paramagnetic above the Neel
temperature T
N
- Ferrimagnetism: Ferrimagnets are similar to antiferromagnets in that two
sublattices exist that couple through a superexchange mechanism to create an
antiparallel alignment. However the magnetic moments on the ions of the sublattices
are not equal and hence they do not cancel; rather, a finite difference remains to
leave a net magnetization. This spontaneous magnetization is defeated by the
thermal energy above the Curie temperature, and then the system is paramagnetic
having a nonlinear 1/x versus T relationship.
NANOMATERIALS
222
8.2 Magnetic Properties of small atomic clusters
8.2.1 Introduction
Metal clusters are composed of a few atoms to hundreds of atoms which
demonstrate properties different from the bulk due to the large surface to volume
ratio compared to its bulk state.
lxxx
e.g. for a cluster of thousand atoms, one quarter of
the atoms lie on the surface. Even in the very compact icosahedral configuration, a
nickel cluster of 55 atoms has 32 atoms on the surface. This provides a unique
opportunity to study the evolution of the physical properties of bulk solid with the
variation of the number of atoms in the system under investigation. Hence, the
understanding of the properties of the clusters in isolation has become a subject of
considerable interest from the stand point of fundamental research.
Another important aspect of the cluster research is that since the physical
properties of the clusters are size specific, there is possibility that materials with
desired properties can be custom made by changing the size and composition of the
cluster aggregates. e.g. thin films with controlled nanostructure characteristics can be
obtained by low energy cluster beam deposition technique.
lxxxi
It is however essential
to understand the intrinsic cluster properties in order to control the final material
properties.
Molecular beam technique is ideal for measuring the properties of clusters in
isolation as the clusters are free from any external influences. Metal clusters can be
produced by making a supersonic beam of metal vapour from an oven, typically
known as cluster-source. In the molecular beam experiment a cluster source and one
cluster detector is required. With the development of the pulsed laser vaporization of
metal substrates inside a high pressure pulsed nozzle, it is possible to synthesize
clusters of almost any metal. Cluster detectors are typical ion detectors specially
designed for the particular mode of experiment.
One of the areas of cluster research currently being pursued is the
investigation of the magnetic properties of small atomic clusters. Since large numbers
of atoms of the clusters are on the surface, one expects the clusters of these
elements to exhibit larger magnetic moments. Being small particles, it is possible to
measure the magnetic moment directly by measuring the interaction with a magnetic
field in terms of their deflection from the original trajectories.
NANOMATERIALS
223
8.2.2 Size dependence
The magnetic moment of transition metals decreases in general with the increase of
number of atoms (size) in the clusters. However, the moment of different transition
metal clusters are found to depend differently on the number of atom per clusters (N).
0
0.5
1
1.5
2
2.5
3
0 100 200 300 400 500 600 700 800
Cluster Size (N) [Number of atoms]
Bulk Fe
Fe-clusters
Co-clusters
Bulk Ni
Ni-clusters
Rh-clusters
Bulk Co
M
a
g
n
e
t
i
c
m
o
m
e
n
t
p
e
r
a
t
o
m
Figure 8-1: Magnetic moment per atom in units of Bohr magnetron for iron, cobalt, nickel and rodium
clusters as a function of cluster size (number of atom). The enclosed data points show the sizes of the
respective clusters at which they attain bulk moments. Only the maximum error in the data for each
clusters are shown. Iron, cobalt and nickel data are from W. A. de Heer et al.
lxxxii
and Billas.
Erreur ! Signet
non dfini.
The rodhium data is from Cox et. al.
Erreur ! Signet non dfini.
Among the ferromagnetic 3d transition metals, Nickel clusters attain the bulk
moment at N=150, whereas Cobalt and Iron clusters reach at the corresponding bulk
limit at N=450 and N=550 respectively (Figure 8-1). In icosahedral configuration these
numbers correspond to 3, 4 and 5 shells respectively.
lxxxiii
However one needs to
confirm the bulk like behaviour only after measuring the temperature dependence of
NANOMATERIALS
224
magnetic moment for clusters. The spin imbalance at 3d-orbital has been correlated
to the measured moment as:
) (
2
) ( ) ( N
g
n n
x
x
= + (8. 210)
Where g
x
is the bulk g-factor for the respective metal. The results show that in the
lower size limit the clusters have high spin majority configuration and the behaviour is
more like an atom, but with the increase of size an overlap of the 3d-band with the
Fermi level occurs, which reduces the moment towards that of the bulk with slow
oscillations.
Erreur ! Signet non dfini.
In analogy with the layer-by-layer magnetic moments
in the studies of surface magnetism, the clusters are assumed to be composed of
spherical layers of atoms forming shells and the value of moment of an atom in a
particular shell is considered to be independent of the size of cluster. The magnetic
moment values of the different shells are optimized to fit with the total moment per
atom maintaining the overall trend of observed decrease with size. It is generally
observed that the per atom moment for the surface atoms in the 1st layer are close to
that of an atom and deep layers correspond the bulk. The intermediate layers show
an overall decreasing trend. For both nickel and cobalt, there are oscillations in the
values at least up to the 5th shell. For nickel, the second layer is found to have a
negative coupling with the surface atoms and the 3rd-4th layers in iron are
magnetically dead. The model gives more or less close agreement with the
experimental results for cobalt and nickel. But in a more realistic model one has to
account for both the geometric shell closing as well as the electronic structure.
However this model provides a good first order picture to understand the magnetic
behaviour of 3d transition metal clusters.
Though the 4d transition metals are non-magnetic, the small clusters of
rhodium are found to be magnetic following the calculations.
lxxxiv
Within the size
range N=10-20, the moment varies between 0.8 to 0.1
B
per atom but the clusters
become non-magnetic within N=100. Due to reduced size, the individual moments on
the atoms in clusters are aligned. First principal calculations show that the ground
state of Rh
6
is non-ferromagnetic, while Rh
9
, Rh
13
and Rh
19
clusters have non-zero
magnetic moments. Rh
43
is found to be non-magnetic as the bulk. Ruthenium,
palladium, Chromium, Vanadium and Aluminium clusters are found to be non-
magnetic.
Erreur ! Signet non dfini.
NANOMATERIALS
225
8.2.3 Thermal behaviour
Since ferromagnetic state requires the moments to remain mutually aligned
even at relatively high temperatures, it is important to study the thermal behaviour of
the magnetic clusters. Temperature dependent studies show that up to 300K the
moment remains constant for nickel clusters and then it decreases at higher
temperatures resembling closely the bulk behaviour. This indicates that the
interactions affecting the mutual alignment of moments with temperature are in the
same order of magnitude as in the bulk. Ni
550-600
is almost bulk-like. Cobalt clusters
also follow the same trend but the moment increases slightly, between 300K to 500K,
which might be due to a structural phase transition corresponding to the hcp-fcc
phase transition in bulk at 670K, where the moment increases by 1.5%. Clusters of
different size ranges of iron behave differently; Fe
50-60
shows a gradual reduction of
moments from 3
B
at 120K to 1.53
B
at 800K. Up to 400K, Fe
120-140
remains with a
constant moment of 3
B
and then decreases to finally level off at 700K to 0.7
B
. The
size range N=82-92 behaves like Fe
120-140
below 400K and like Fe
50-60
above 400K.
The moment for the higher size ranges decreases steadily with the increase of
temperature but Fe
250-290
levels off at 700K whereas Fe
500-600
levels off between 500-
600K to 0.4
B
. Comparing with the thermal behaviour of bulk iron, the level-off
temperature can be termed as the critical temperature (T
C
) for clusters of a particular
size and it decreases with the increase of size of the clusters. But, since the bulk
Curie temperature is 1043K, one realizes clearly that this trend must reverse at
higher sizes to reach the bulk. It is also interesting to note that for Fe
120-140
the
moment decreases more rapidly within 600-700K. The specific heat measurement for
both Fe
120-140
and Fe
250-290
also show a peak at 650K. This might be because iron
clusters undergo a structural phase transition in this temperature range similar to the
bulk iron which undergoes bcc-fcc-bcc phase transition, but only beyond the Curie
temperature. Since the structure of clusters cannot be detected in a molecular beam,
this remains an unsettled question. However, there is a competition between the
different structures for iron clusters,
Erreur ! Signet non dfini.
and it is more probable that a
structural transition occurs at low temperature for increasing cluster size but revealing
it needs a better characterization.
NANOMATERIALS
226
8.2.4 Rare earth clusters
The magnetic interaction in rare earth solids are RKKY type which is mediated
by the conduction electrons. Since in clusters, structure and filling up of the
conduction band differs from that in the bulk, rare-earth clusters are expected to
exhibit an altogether different magnetism compared to the bulk. Magnetic properties
of gadolinium, terbium and dysprosium clusters show a very size-specific
behaviour.
Erreur ! Signet non dfini.
Gd
22
behaves super-paramagnetically (8.3.4) but Gd
23
shows a behaviour like a rigid magnetic rotor as if the moment is locked to the cluster
lattice. Gd
17
shows super-paramagnetism near room temperature but at low
temperatures it shows locked-moment behaviour. It indicates that there are at least
two groups of clusters sharing the same number of atoms which might correspond to
two different isomers, either magnetic or structural. But it is not clear whether the two
isomers co-exist in the beam or there occurs a phase transition of any kind. The
behaviour of other clusters from N=11-30, is like Gd
22
or Gd
23
or Gd
17
. Similar is the
case for Tb and Dy clusters. Tb
17
, Tb
22
, Tb
27
, Dy
22
and Dy
27
show super-
paramagnetic behaviour but Tb
23
and Tb
26
show locked-moment behaviour. In fact
Tb
26
is temperature dependent like Gd
17
. Assuming the clusters to be spherical with a
moment of inertia and a body fixed moment which rotates with the cluster,
calculations show a very complex precession and nutation when the rotational kinetic
energy becomes comparable to the magnetic potential energy. With two parameters
viz. the rotational temperature and the magnetic moment, a fit of the experimental
data for Gd with several initial conditions and a rotational temperature of 5K, shows
an agreement at least for sizes lower than N=30.
8.3 Small particle magnetism
8.3.1 Classifications of magnetic nanomaterial
The magnetic behaviour of most experimental realizable systems of is result of
contributions from both interaction and size effects. The correlation between
nanostructure and magnetic properties suggests a classification of nanostructure
morphologies. The following classification is designed to emphasize the physical
mechanisms responsible for the magnetic behaviour. Figure 8-2 schematically
represents the classification. At one extreme are systems of isolated particles with
NANOMATERIALS
227
nanoscale diameters, which are denoted by type A. These non interacting systems
derive their unique magnetic properties strictly from the reduced size from the
components, with no contribution from interparticle interactions. At the other extreme
are bulk materials with nanoscale structure denoted by type D, in which a significant
fraction (up to 50%) of the sample volume is composed of grain boundaries and
interfaces. In contrast to type A systems, magnetic properties here are dominated by
interactions. The length of the interactions can span many grains and is critically
dependent on the character of interface. Due to this dominance of interaction and
grain boundaries, the magnetic behaviour of type D nanostructures cannot be
predicted simply by applying theories of polycrystalline materials with reduced length
scale. In type B particles, the presence of a shell can help prevent particle
interactions, but often at the cost of interaction between the core and the shell. In
many cases the shells are formed via oxidation and may themselves be magnetic. In
type C materials, the magnetic interactions are determined by the volume fraction of
the magnetic particles and the character of the matrix they are embedded in.
A
B
C
D
Figure 8-2 : Schematic representation of the different types of magnetic nanostructures. Type-A
materials include the ideal ultrafine particle system, with interparticle spacing large enough to
NANOMATERIALS
228
approximate the particles as noninteracting. Ferrofluids, in which magnetic particles are surrounded
by a surfactant preventing interactions, are a subgroup of Type-A materials. Type-B materials are
ultrafine particles with a core-shell morphology. Type-C nanocomposites are composed of small
magnetic particles embedded in a chemically dissimilar matrix. The matrix may or may not be
magnetoactive. Type-D materials consist of small crystallites dispersed in a noncrystalline matrix. The
nanostructure may be two-phase, in which nanocrystallites are distinct phases from the matrix, or the
ideal case in which both the crystallites and the matrix are made of the same material.
8.3.1.1 Ferrofluid
A ferrofluid is a synthetic liquid that holds small magnetic particles in a colloidal
suspension, with particles held aloft their thermal energy. The particles are
sufficiently small that the ferrofluid retains its liquid characteristics even in the
pressure of a magnetic field, and substantial magnetic forces can be induced which
results in fluid motion.
A ferrofluid has three primary components. The carrier is the liquid element in
which the magnetic particles are suspended. Most ferrofluids are either water based
or oil based. The suspended materials are small ferromagnetic particles such as iron
oxide, on the order of 10 20 nm in diameter. The small size is necessary to
maintain stability of the colloidal suspension, as particles significantly larger than this
well precipitate. A surfactant coats the ferrofluid particles to help maintain the
consistency of the colloidal suspension.
The magnetic properties of the ferrofluid are strongly dependent on particle
concentration and on the properties of the applied magnetic field. With an applied
field, the particles align in the direction of the field, magnetizing the fluid. The
tendency for the particles to agglomerate due to magnetic interaction between
particles is opposed by the thermal energy of the particles. Although, particles vary in
shape and size distribution, insight into fluid dynamics can be gained by considering
a simple spherical model of the suspended particles.
NANOMATERIALS
229
S
2R
, , , ,
Figure 8-3: Representative model of a typical ferrofluid.
The particles are free to move in the carrier fluid under the influence of an applied
magnetic field, but on average the particles maintain a spacing S to nearest
neighbours. In a low density fluid, the spacing S is much larger than the mean
particle radius 2R, and magnetic dipole dipole interactions are minimal.
Applications for ferrofluids exploit the ability to position and shape the fluid
magnetically. Some applications are:
rotary shaft seals
magnetic liquid seals, to form a seal between region of different pressures
cooling and resonance damping for loudspeaker coils
printing with magnetic inks
inertial damping, by adjusting the mixture of the ferrofluid, the fluid viscosity
may be change to critically damp resonances accelerometers,
level and attribute sensors
electromagnetically triggered drug delivery
Table 8-2: Material Parameters of a typical ferrofluid:
lxxxv
Sample Volume 1.7 x 10
-6
m
3
Electrical conductivity of particles
6 1
3 10 ( )
f
m
=
Electrical conductivity of ferrofluid
7 1
10 ( )
fluid
m
<
Volume percentage of particles 3% by volume
Initial magnetic permeability of particles
0
100
NANOMATERIALS
230
Particle mean radius R10
-8
m
Density of magnetic particles
FE
= 7.8 g/cm
3
Fluid density
f
= 1 g/cm
3
Fluid viscosity of carrier fluid
=
2
) / ( ln
exp
2
1
) (
0
V V
V
V f (8. 214)
Where V
0
is the most probable particle volume and is the standard deviation of
ln(V). The validity of this simple assumption has been experimentally verified in many
systems by transmission electron microscopy (TEM) and other indirect techniques
such as x-ray diffraction and magnetic granulometry.
One fundamental question related to the finite and small volume of the particles is
how such volume determines the internal domain structure. In magnetic bulk
materials, it is well known that there exists a multidomain structure constituted by
regions of uniform magnetization separated by Domain Walls (DWs). The size and
shape of these domains depends on the interplay between the exchange,
magnetostatic and anisotropy energies of the system. As the volume of the magnetic
system decreases, the size of the domains and the width of the walls are reduced,
modifying, at the same time, their inner structure. Below a certain critical volume, the
energy cost to produce a DW is greater than the corresponding reduction in the
magnetostatic energy. Consequently, the system no longer divides in smaller
domains, instead maintaining the magnetic structure of a single domain. This critical
value depends on the saturation magnetization of the particle, anisotropy energy and
exchange interactions between individual spins. For example, for spherical particles
the critical diameter is within 10800 nm. Typical values for Fe and Co metallic
particles are 15 and 35 nm, respectively, while for SmCo5 it is as large as 750 nm.
NANOMATERIALS
234
In many fine-particle systems, f(V) is narrowand centred at very small V
0
below the
critical value, so that all the particles constituting the system can be considered as a
single domain.
This is the situation that we are going to assume in what follows. In thermal
equilibrium, the resulting magnetization of a particle points in a direction tending to
minimize its total anisotropy energy, which for a single-domain particle can be
considered as proportional to its volume, at least, in first approximation (Figure 8-6).
The total anisotropy of the particle can often be assumed to have uniaxial character,
being determined by a single constant K. Therefore, the anisotropy energy can be
expressed as
E() = kV Sin
2
(8. 215)
where is the angle between anisotropy axis and magnetization, and KV is the
anisotropy energy barrier (8.3.2) separating both easy directions for magnetization,
which in zero magnetic field correspond to = 0 and = . Taking into account the
spread of volumes f(V) given by equation E( ) = kV Sin
2
(8. 215) and
(
=
2
) / ( ln
exp
2
1
) (
0
V V
V
V f (8. 214), it is a matter of fact that the magnetization
of the particles will behave as if it were driven by an energy barrier distribution
) (
0
B
E f , where
0
B
E = kV.
Under the application of a magnetic field H forming an angle with the anisotropy
axis, the energy of a particle is modified as
) cos sin sin cos (cos sin ) ( + + =
S
HM KV E (8. 216)
Where M
s
is the spontaneous magnetization of the particle, and the axis system is
defined in Figure 8-5.
NANOMATERIALS
235
Figure 8-5 : Definition of the axis system for a fine particle. The uniaxial-anisotropy axis is along the z-
axis. [X. Batle, A. Labarta, Finite-size effects in fine particles: magnetic and transport properties,
Journal of Physics, 2002, ppR15-R42].
Under zero field, the equilibrium direction for the magnetization vector M coincides
with the anisotropy axis, while when a magnetic field is applied, it rotates away from
the anisotropy axis towards the field direction, forming an angle with H, and the
stability problem can be reduced to a two-dimensional one ( = 0). In this case, and
if H is lower than a certain critical value known as switching field at zero temperature,
0
SW
H , E() has two minima at different heights which are separated by an energy
barrier. When the field is applied opposite to the magnetization, the unstable
minimum is separated from the stable one by an energy barrier E
B
, which can be
approximated by the following expression:
|
|
\
|
=
0
0
1 ) (
SW
B B
H
H
E H E (8. 217)
Where is a phenomenological exponent that depends on . Consequently,
0
SW
H is
the minimum value of the field at zero T, at which E
B
vanishes, and the magnetization
inverts its orientation. When H is applied along the anisotropy axis and in the
opposite direction to magnetization, equation
|
|
\
|
=
0
0
1 ) (
SW
B B
H
H
E H E (8. 217) is
exact with = 2. In this case, M irreversably rotates at zero T when H =
0
SW
H = H
a
=
2K/M
s
, where H
a
is the anisotropy field. If H is applied at an arbitrary angle , there is
NANOMATERIALS
236
not a simple analytical expression. However, equation
|
|
\
|
=
0
0
1 ) (
SW
B B
H
H
E H E (8.
217) works quite well with 1.5, provided that the applied field forms an angle of a
few degrees with the anisotropy axis. On the other hand, Pffeifer gave the following
phenomenological approximation for any angle
a
SW
H
H
0
14 . 1
86 . 0 + = (8. 218)
A method frequently used to characterize the nature and strength of the particle
anisotropy consists in the determination of
0
SW
H as a function of . The resulting
curve turns to be the so-called Stoner Wohlfarth (SW) astroid. In fact, the SW astroid
represents the set of fields (applied opposite to the particle magnetization) at which
irreversible jumps occur. For uniaxial anisotropy,
0
SW
H can be evaluated from
equation ) cos sin sin cos (cos sin ) ( + + =
S
HM KV E (8. 216) with the
conditions dE/d = dE/d = 0, since irreversible jumps occur when the energy
barrier vanishes and the maximum in E( ) is substituted by an inflection point. The
obtained expression is
2 / 3 3 / 2 3 / 2
0
) cos (sin +
=
a
SW
H
H (8. 219)
Many groups have attempted to study the reversal process in single particles but the
experimental precision did not yield quantitative information.
Single domain
Multi domain
Super
Paramagnetic
H
c
D
c
D
NANOMATERIALS
237
Figure 8-6 : Qualitative illustration of the behaviour of the coercivity in ultrafine particle systems are
the particles size changes.
As the particles size continues to decrease (towards some critical particle diameter,
D
c
) below the single domain value, the spins get increasing affected by thermal
fluctuations and the system becomes superparamagnetic. Particles with sufficient
shape anisotropy can remain single domain to much larger dimensions than their
spherical counterparts.
8.3.3.2 Single-Domain Characteristics:
In a granular magnetic solid with a low volume fraction, one has a collection of single
domain particles each with a magnetic axis along which all the moments are aligned.
In the absence of a magnetic field, parallel and anti parallel orientations along the
magnetic axis are energetically equivalent but separated by an energy barrier of CV,
where C is the total anisotropy per volume, and V is the particle volume. Since the
size of each single domain remains fixed, under an external field, only the magnetic
axes rotate. Thus the measured magnetization M of a granular magnetic field solid
with a collection of single domain particles is the global magnetization
.
cos
S
M H
M M
H
= =
(8. 220)
Where is the angle between the magnetic axis of a particle.
S
M , is the saturation
magnetization, H
is the external field, and the average <cos > is taken over many
ferromagnetic particles. The hysterisis loop of a granular solid is thus a signature of
the rotation of the magnetic area of the single-domain particles. This should be
contrasted with the hysterisis loop of a bulk ferromagnet, in which the sizes and the
direction of the domains are altered drastically under an external field.
NANOMATERIALS
238
M
H
M
S
M
R
M
C
(a)
T
B
T
H
C
(c)
M
S
M
M
S /2
T
B
T
e
c
d
f
b
a
T
B
(b)
(d)
T
Figure 8-7 : (a) Hysterisis loop at 5K; (b) temperature dependence of saturation magnetization (M
S
)
and remnant magnetization (M
R
=M
S/2
); (c) temperature dependence of field cooled (FC) and zero-field
cooled (ZFC) susceptibility. At the blocking temperature (T
B
), M
R
and H
C
vanish, whereas the ZFC
susceptibility shows a cusp like feature.
An example of a hysterisis loop of a granular magnetic solid at low temperature is
shown in Figure 8-7a. In the initial unmagnetized state with M = 0 at H = 0, the
magnetic axes of the particles are randomly oriented, each along its own magnetic
axis, which is determine by the total magnetic anisotropy of the particles. The
directions of the giant moments are random and static at low temperatures. A
saturation magnetization (M = M
S
) is realized under a large field when all the
magnetic axes are aligned. In the remnant state when H is reduced to H = 0, one
observes the remnant magnetization (M
R
), whose values at low temperature is M
R
=
M
S/2
. This is because the magnetic axes are oriented only in one hemisphere due to
the uniaxial anisotropy of the single-domain particle.
It should be noted that the initial M=0 state and the initial magnetization curve if in
Figure 8-7a does not reappear, whereas the field cycle part does. The simples way to
recover the initial M=0 state is to heat the sample above the blocking temperature
and cool the sample in zero field back to low temperatures.
Because of the single domain nature of the magnetic entities, the coercivity
(H
C
) of the ultra fine particles is much higher than that in bulk material.
NANOMATERIALS
239
8.3.4 Superparamagnetism
8.3.4.1 The SPM regime
Below the Curie temperature of a ferromagnet or ferrimagnet, all the spins are
coupled together and so cooperate to yield a large moment. This moment is bound
rigidly to the particle by one or more of the variety of anisotropies that we have
discussed (8.3.2). With decreasing particle size, the energy Barrier E
B
decreases
until the thermal energy kT can disrupt the bonding f the total moment to the particle.
Then this moment is free to move and respond to an applied field independent of the
particle. This moment is the particle moment and is equal to
p
=M
S
V (8. 221)
It can be quite large, thousands of Bohr magnetons. An applied field would tend to
align this giant (super) moment, but kT would fight he alignment just as it does in a
paramagnet. Thus, this phenomenon is called superparamagnetism.
It is good to notice that if the anisotropy is very weak (zero), one would expect that
the total moment
p
could point in any direction, hence the Langevin law would apply.
The phenomenon of superparamagnetism is, in fact, timescale-dependent due to the
stochastic nature of the thermal energy. The anisotropy energy E
B
represents an
energy barrier to the total spin reorientation; hence the probability for jumping this
barrier is given by Arrhenius law:
|
\
|
=
kT
E
B
exp
0
(8. 222)
Where
0
is the attempt time by supposing that the particle spins are rigidly coupled
and reverse as a whole. This pre-factor depends on a large variety of parameters:
temperature, gyromagnetic ratio, saturation magnetization, energy barrier, direction
of the applied field and damping constant.
0
is assumed to be constant and
equation |
\
|
=
kT
E
B
exp
0
(8. 222), together with this assumption, is known as the
NeelBrown model. If the characteristic time window of an experiment is much
shorter than at a fixed temperature, the particle moment remains blocked during
the observation period (blocked regime), while, in the opposite situation, the rapid
fluctuations of the particle moment mimics the atomic paramagnetism (SPM regime).
NANOMATERIALS
240
In the intermediate regime, the probability that the magnetization has not switched
after an observation time t is given by
|
\
|
=
t
t P exp ) ( (8. 223)
A method used to verify the NeelBrown model in real systems is the scaling of the
switching field distributions measured at different T and rates of field ramping.
Thermal activation leads to a switching field distribution from the SW model, the
mean switching field being given by
(
(
|
|
\
|
|
\
|
=
/ 1
1
0
ln 1
cT
E
T k
H H
B
B
SW SW
(8. 224)
Where ) /(
0
0
B SW B
E H k c = , ) / 1 (
0
SW
H H = and is the field sweeping rate. From
the equation above, it is evident that H
sw
measured at different T and can be
collapsed in a single master curve when they are plotted as a function of the scaling
variable (T ln(cT /
1
))
1/
. In fact, the collapse is obtained by choosing the
appropriate values of and c (obtaining = 1.5, as expected for measurements with
the field slightly deviated from the anisotropy axis, see paragraph 8.3.3.1). This
scaling procedure has been successfully applied to analyse data obtained above
about 0.5K from single-particle experiments carried out with different kind of
nanoparticles.
At very low temperatures, (below about 0.5 K) strong deviations from the NeelBrown
model have been reported in this kind of scaling when it was applied to single-particle
experiments. The deviations consisted in a saturation of H
sw
and (width of the
switching field distribution), and a dependence of H
sw
on faster than the prediction
of the NeelBrown model. All these anomalies have been attributed to the
prevalence of a nonthermal process by which the particle magnetization escapes
from the metastable state through the energy barrier by macroscopic quantum
tunnelling (MQT). The theoretical models propose that there is a crossover
temperature T
B
, below which quantum tunnelling dominates the reversal process,
whereas above T
B
, the escape rate from the metastable state is given by thermal
activation, and the NeelBrown model is accomplished. Consequently, T
B
is the
temperature at which the thermal and quantum rates coincide.
The main parameters determining the particular value of T
B
for a given system are:
NANOMATERIALS
241
the number of spins contained in each magnetic entity (which should be as
small as possible to get a high T
B
),
anisotropy energy (which should be as high as possible)
applied magnetic field.
In a sense, this quantum phenomenon is related to finite-size effects, considering that
T
B
can only be experimentally achievable in magnetic particles of small enough size
(at most, few nanometers).
Deduced values for are 0.5 and 0.035 for Co and BaFe122xCoxTixO19
nanoparticles, the former larger, as expected, being due to the effect of the
conduction electrons. Little information is available on other fine-particle systems.
Moreover, damping also plays an important role in MQT. In general, dissipation due
to, e.g. conduction electrons, strongly reduces quantum effects, in agreement with
the observation of MQT in BaFe122xCoxTixO19 for T < 0.4K, whereas this effect is
absent for Co particles down to 0.2K.
8.3.4.2 Magnetism relaxation
If the total magnetization of a fine-particle system is measured with an observation
time window which for some particle sizes is comparable to their reversal time, , the
magnetization will evolve during the experiment, a phenomenon known as magnetic
relaxation. The study of this non-equilibrium phenomenon is a common method to
experimentally characterize the energy spectrum of the system, since it is the
macroscopic signature of the distribution of energy barriers separating local minima
which correspond to different orientations of the particle moments. Magnetic
relaxation experiments are carried out by measuring the time-dependence of the
magnetization after the system has been driven out of thermal equilibrium by the
application or removal of a magnetic field. For instance, starting from the equilibrium
state at a given field and temperature, and then removing the field, the time evolution
of the magnetization for noninteracting particles can be written as follows
|
|
\
|
=
) (
) ( ) (
0
0
E
t
xp e E dEf M t M
(8. 225)
Where f (E)dE is the fraction of particles having energy barriers in between E and E
+ dE, and M
0
is the initial magnetization. In equation
|
|
\
|
=
) (
) ( ) (
0
0
E
t
xp e E dEf M t M
NANOMATERIALS
242
(8. 225) and for the sake of simplicity, the particle moment has been assumed
to be volume independent, although, in fact, it is proportional to it. This assumption
does not significantly affect the results, providing that for a logarithmic-linear
distribution, f (E) and Ef (E) (which is equivalent to Vf (V)) are similarly shaped
functions. Let us introduce the function p(t, T,E), which is defined by
(
\
|
|
|
\
|
=
kT
E t
E T t p exp exp ) , , (
0
(8. 226)
Taking into account that p(t, T,E) for a given t varies abruptly from 0 to 1, as the
energy barrier E increases, we can approximate p(t, T,E) by a step function whose
discontinuity E
c
(t, T ) moves towards higher energies as time elapses. Consequently,
the integral is cut off at the lower limit by the value of E
c
(t, T ) and can be written as
Ec
E dEf M t M ) ( ) (
0
(8. 227)
Where E
c
(t, T ) = kBT ln(t/
0
) is the only time-dependent parameter and signals the
position of the inflection point of p(t, T,E). In fact, the narrow energy range spanned
by the step at E
c
(t, T ) can be understood as the experimental window of the
measurements. The experimental window can be swept over f (E) by varying the
temperature or the observation time. From equation
Ec
E dEf M t M ) ( ) (
0
(8. 227), it
can be deduced that M(t) obtained after integration is a function of the parameter
E
c
(t, T ), which acts as a scaling variable. The existence of this scaling variable
implies that measuring the magnetization as a function of the temperature at a given
time is equivalent to measuring the magnetization as a function of ln(t/
0
) at a given
temperature. This timetemperature correspondence is characteristic of activated
processes governed by Arrhenius law. The validity of the scaling hypothesis is only
determined by the validity of the cut-off approximation. Therefore, the scaling will be
accomplished as long as the width of p(t, T,E), which is approximately given by kBT,
is small when compared to the width of the energy barrier distribution. Although this
condition seems to be very restrictive, it is the usual situation found in experimental
observations of the magnetic relaxation in fine-particle systems.
NANOMATERIALS
243
As an example of this scaling method, Figure 8-8 shows the scaling of the relaxation
curves obtained at different temperatures for a sample consisting on Fe
0.78
C
0.22
particles with a mean diameter of 3.6 nm, dispersed in stable dilution with a
hydrocarbon oil. In this sample, the dipolar interaction among particles was estimated
to be very small. Figure 8-8 shows all the relaxation curves corresponding to
different temperatures collapsed into a single master curve when plotted as a
function of the scaling variable T ln(t/ 0), being
0
= (3.5 5) 10
11
s, which is an
expected value for ferromagnetic fine particles.
Figure 8-8 : The resulting scaling for a sample composed of FeC nanoparticles dispersed in a
hydrocarbon oil is shown. Open and full circles correspond alternatively to adjoining temperatures,
which are indicated above the corresponding interval (K units). The solid line is the theoretical curve
calculated by fitting the scaling curve to equation
Ec
E dEf M t M ) ( ) (
0
(8. 227) with a
logarithmic-linear distribution of energies. The values of the fitted parameters are: = 0.44 0.05 and
E0 = 287 50K (position of the maximum of the energy barrier distribution). [X. Batle, A. Labarta,
Finite-size effects in fine particles: magnetic and transport properties, Journal of Physics, 2002,
ppR15-R42].
One of the most interesting aspects of these results is that, in fact, measuring the
relaxation at a given temperature is completely equivalent to measuring it at different
temperature but shifting the observation time window according to the law
NANOMATERIALS
244
T ln(t/
0
). In this sense, the method enables us to obtain the relaxation curve at a
given temperature, in a time range that is not experimentally accessible, by simply
dividing the T ln(t/
0
) axis by this temperature. For instance, in the case shown in
Figure 8-8, we can obtain the relaxation curve at the lowest measured temperature of
1.8K at times as large as 10
119
s, which is obviously an experimentally inaccessible
time. The validity of the T ln(t/
0
) scaling method to analyse the relaxation behaviour
has been proved in many systems of non or slightly interacting fine particles.
Moreover, from the T ln(t/
0
) scaling of relaxation data, insight can be gained into the
microscopic details of the energy barrier distribution producing the relaxation. If f (E)
is nearly constant in the range of energy barriers which can be overcome during the
observation time, then equation
Ec
E dEf M t M ) ( ) (
0
(8. 227) can be approximated
by
|
|
\
|
|
|
\
|
0
0
ln ) ( 1 ) (
t
E Tf k M t M
C B
(8. 228)
Where E
c
is the mean energy barrier within the experimental window. From equation
|
|
\
|
|
|
\
|
0
0
ln ) ( 1 ) (
t
E Tf k M t M
C B
(8. 228)) we can define the viscosity parameter as
T k E f
t
M
M
S
B C
) (
)) (ln(
1
0
=
= (8. 229)
Therefore, as E
c
is varied, the magnetic viscosity maps the energy barrier distribution
at low enough temperatures for which the width of p(t, T,E) (equation
(
\
|
|
|
\
|
=
kT
E t
E T t p exp exp ) , , (
0
Ec
E dEf M t M ) ( ) (
0
(8. 227). The differential of the thermoremanence relative to
the saturation magnetization versus the temperature is also shown in full circles for comparison. [X.
Batle, A. Labarta, Finite-size effects in fine particles: magnetic and transport properties, Journal of
Physics, 2002, ppR15-R42].
If f (E) is almost constant over a certain range of energies, then S will be linear in T
over the corresponding temperature range; consequently, relaxation as a function of
time will follow a time logarithmic decay at fixed T. This behaviour is only
accomplished in a narrow region around the maximum of a bell-shaped energy
barrier distribution (see for instance, the region around the inflection point of the
master curve shown in Figure 8-8).
T ln(t/
0
) scaling can be also applied to relaxation experiments measuring the
acquisition of magnetization of an initially demagnetized system under the application
of a magnetic field. In this kind of experiments, the field modifies the energy barriers
that are responsible for the time evolution, as well as the final equilibrium state
towards which the system relaxes. The main effect of the field over f (E) is the
splitting into two sub-distributions which are centred at lower and higher energy
values as the field increases, the former dominating the relaxation of the system. For
a noninteracting system, the relaxation curves collected at different temperatures and
under a fixed field collapse into a single master curve when they are plotted as a
function of T ln(t/
0
) variable and proper normalization of the magnetization data to
NANOMATERIALS
246
the final equilibrium values is performed. Moreover, for fields smaller than the
switching field corresponding to the smallest barriers, the relaxation curves at
different fields and fixed T can be collapsed into a single curve, in a way similar to T
ln(t/
0
) scaling for curves at fixed field, but in this case, the additional shifts in the T
ln(t/
0
) axis necessary to get field scaling give the field dependence of the mean
relaxing barriers present in the system.
8.3.4.3 Effects of the interparticle interaction
We have reviewed different aspects of the magnetic relaxation of fine-particle
systems with random orientations, volume distribution and negligible magnetic
interactions. In this situation, the system is in an SPM regime, for which the time
evolution of the total magnetization is simply governed by the thermal activation over
the individual energy barriers of each particle. When interparticle interactions (for
instance, dipolar interactions) are non-negligible the behaviour of the system is
substantially more complicated and the problem becomes non-trivial, even if the
spins of the particle are assumed to be coupled to yield a super-spin moment.
The main types of magnetic interactions that can be found in fine-particle assemblies
are:
dipoledipole interactions, which always exit
exchange interactions through the surface of the particles which are in close
contact
in granular solids, RKKY interactions through a metallic matrix when particles
are also metallic, and super exchange interactions when the matrix is
insulating.
Bearing in mind the anisotropic character of dipolar interactions, which may favour
ferromagnetic or antiferromagnetic alignments of the moments depending on
geometry, fine-particle systems have all the ingredients necessary to give rise to a
spin-glass state, namely, random distribution of easy axes and frustration of the
magnetic interactions. The complex interplay between both sources of magnetic
disorder determines the state of the assembly and its dynamical properties.
Magnetic interactions modify the energy barrier coming from the anisotropy
contributions of each particle and, in the limit of strong interactions, their effects
become dominant and individual energy barriers can no longer be considered, only
NANOMATERIALS
247
the total energy of the assembly being a relevant magnitude. In this limit, relaxation is
governed by the evolution of the system through an energy landscape with a
complex hierarchy of local minima similar to that of spin-glasses. It is worth noticing
that in contrast with the static energy barrier distribution arising only from the
anisotropy contribution, the reversal of one particle moment may change the energy
barriers of the assembly, even in the weak interaction limit. Therefore, the energy
barrier distribution may evolve as the magnetization relaxes.
Three models have been developed to introduce particle interaction:
- Shtrikmann and Wohlfarth
For a weak interaction limit, the relaxation time is
|
|
\
|
=
) (
exp
0
0
T T k
E
B B
B
(8. 230)
Where T
B
is the blocking temperature and T
0
is an effective temperature which
accounts for the interaction effects.
- Dormann et al
This model correctly reproduces the variation of the blocking temperature T
B
, as a
function of the observation time window of the experiment,
m
, at least in a range of
time covering eight decades. The increase of T
B
with the strength of the dipolar
interactions (e.g. increasing particle concentration or decreasing particle distances)
has been predicted by this model and also experimentally confirmed.
- Morup and Tronc
For the weak interaction limit, the opposite dependence of T
B
with the strength of the
interactions is predicted. Morup suggested that two magnetic regimes, governed by
opposite dependencies of T
B
, occur in interacting fine particles. At high temperatures
and/or for weak interactions, T
B
signals the onset of a blocked state and T
B
decreases as the interactions increase. In contrast, at high temperatures and/or for
strong interactions, a transition occurs from an SPM state to a collective state which
shows most of the features of typical glassy behaviour. In this case, T
B
is associated
with a freezing process and it increases with the interactions.
Non-equilibrium dynamics, showing ageing effects in the relaxation of the residual
magnetization, have been observed in interacting fine-particle assemblies, but some
important differences with spin-glasses have been established. In the last few years,
NANOMATERIALS
248
several works demonstrated that for concentrated frozen ferrofluids and assemblies
of nanoparticles dispersed in a polymer, the relaxation depends on the waiting time,
t
w
, spent at constant temperature before the magnetic field is changed (Figure 8-10).
Figure 8-10 : The relaxation rate S versus log10(t) at different wait times for the concentrated sample
(a) and for the most diluted sample (b). The measurements were performed on a ferrofluid consisting
of closely spherical particles of maghemite with a mean diameter of 7 nm. [X. Batle, A. Labarta, Finite-
size effects in fine particles: magnetic and transport properties, Journal of Physics, 2002, ppR15-R42].
This phenomenon is absent in the most diluted samples, confirming that it is due to
dipoledipole interactions (Figure 8-10). Ageing effects on the magnetic relaxation
are the fingerprint of the existence of many minima in the phase space. However, a
complex hierarchy of energy minima is not an exclusive feature of spin-glasses and
ageing has also been found in other cluster-glass systems not having a true phase
transition. In fact, the non-equilibrium dynamics of interacting fine particles largely
mimics the corresponding of spin-glasses, including memory experiments in which
ageing is studied at different heating or cooling rates and/or cycling the temperature.
The main differences between the dynamics of canonical spin-glasses and that of
fine particles are:
NANOMATERIALS
249
the dependence of the magnetic relaxation on t
w
is weaker than in spin-
glasses
in the collective state, the relaxation times are widely distributed and strongly
dependent on temperature
large particles are blocked in all time scales acting as a temporary random
field and not taking part in the collective state of the system.
Another feature that differentiates the low-temperature collective state from the spin-
glass one is its extreme sensitivity to the application of an external magnetic field. In
particular, it has been shown that the collective state of strong interacting particles
can be erased by a field-cooling process with an applied magnetic field of moderate
strength which yields an asperomagnetic state. It has also been demonstrated that
the dynamics of these systems are strongly affected by the initial magnetic moment
configuration, in such a way that the collective state determines the dynamic
behaviour only in low-cooling-field experiments, while at high cooling fields the
dynamics are mostly dominated by the intrinsic energy barriers of the individual
particles.
8.3.4.4 Surface effects
Surface effects result from the lack of translational symmetry at the boundaries of the
particle because of the lower coordination number there and the existence of broken
magnetic exchange bonds which lead to surface spin disorder and frustration.
Surface effects dominates the magnetic properties of the smallest particles since
decreasing the particle size increases the ratio of surface spins to the total number of
spins. For example, for fcc Co with a lattice constant of 0.355 nm, particles having
about 200 atoms will have diameters around 1.6 nm and 60% of the total spins will be
at the surface. This represents more than half of the spins. Consequently, the ideal
model of a superspin formed by all the spins of the particle pointing in the anisotropy
direction and coherently reversing due to thermal activation is no longer valid, since
misalignment of the surface spin yields strong deviations from the bulk behaviour.
This is true, even for particles with strong exchange interactions such as many
ferrimagnetic oxides with magnetically competing sublattices, since broken bonds
destabilize magnetic order giving rise to frustration which is enhanced with the
strength of the interactions. As a consequence of the combination of both finite-size
NANOMATERIALS
250
and surface effects, the profile of the magnetization is not uniform across the particle
and the magnetization of the surface layer is smaller than that corresponding to the
central spins.
Many experimental results for metallic and oxide particles indicate that the anisotropy
of fine particles increases as the volume is reduced because of the contribution of
what is known as surface anisotropy. For instance, the anisotropy per unit volume
increases by more than one order of magnitude for 1.8 nm fcc Co particles being
310
7
erg cm
3
compared with the bulk value of 2.710
6
erg cm
3
. Even an
anisotropy value one order of magnitude larger than the preceding case has been
reported for Co particles embedded in a Cu matrix. In fact, surface anisotropy has a
crystal-field nature and it comes from the symmetry breaking at the boundaries of the
particle. The structural relaxation yielding the contraction of surface layers and the
existence of some degree of atomic disorder and vacancies induce local crystal fields
with predominant axial character normal to the surface, which may produce easy-axis
or easy-plane anisotropies. This can be justified by noting that the axis of the local
crystal field,
^
n , may be evaluated from the dipole moment of the nearest-neighbour
atomic positions with respect to the position of a given surface atom as follows:
nn
j
i j
i P P n ) (
^
(8. 231)
Where P
i
is the position of the ith atom and the sum extends to the nearest
neighbours of this atom. Since at the surface some of the neighbours are missing, i n
^
is non-zero and directed approximately normal to the surface. The effect of these
local fields is obtained by adding a term of the form
2
KS . Both cases corresponding to K <0 (radial) and K >0 (tangential) are displayed. [X. Batle, A.
Labarta, Finite-size effects in fine particles: magnetic and transport properties, Journal of Physics,
2002, ppR15-R42].
The second contribution due to the existence of boundaries is a consequence of
strains related to lattice deformations occurring at the surface, which through
magnetostriction effects induce an additional surface anisotropy. Strain anisotropy
has been observed in thin films because of stresses induced at the interface between
substrate and film due to non-matching lattice constants. Depending on the structural
deformations and the nature of the thin films, strain anisotropy could yield
perpendicular anisotropy to the film plane, as has been observed in as-prepared thin
films of heterogeneous Co(Fe)Ag(Cu) alloys. However, for the majority of the
nanoparticle systems, the strain energy is weak, it being difficult to give a general
formulation, and, in particular, for free particles it is negligible.
An effective anisotropy energy per unit volume, K
eff
, could be obtained by adding the
core (i.e. bulk anisotropy) and surface contributions. For a spherical particle the
following phenomenological expression has been used to account for K
eff
:
S b eff
K
d
K K
6
+ = (8. 232)
Where K
b
is the bulk anisotropy energy per unit volume, K
s
is the surface density of
anisotropy energy and d is the diameter of the particle. It is worth noticing that for a
spherical particle and based on symmetry arguments, surface anisotropy (normal to
the surface) would average to zero. However, this is not true for a nanometric particle
NANOMATERIALS
252
with a few atomic layers. For instance, applying equation
S b eff
K
d
K K
6
+ = (8.
232) for a 2 nm particle of fcc Co with K
b
= 2.710
6
erg cm
3
and K
s
1 erg cm
2
,
the surface contribution to the total anisotropy, K
eff
, is about 3 10
7
erg cm
3
, which
is one order of magnitude larger than the bulk contribution, and K
eff
= 3.3 10
7
erg
cm
3
.
This example is representative of the major role of surface contribution to the total
anisotropy in fine-particle systems, for which the anisotropy energy is governed by
surface anisotropy.
8.4 Magnetoelectronics spins
Spintronics is a multidisciplinary field whose central theme is the active manipulation
of spin degrees of freedom in solid-state systems. The goal of spintronics is to
understand the interaction between the particle spin and its solid-state environments
and to make useful devices using the acquired knowledge.
In this chapter, we will focus on magnetoelectronics materials. Typically they cover
paramagnetic and ferromagnetic metals and insulators which utilize magnetorisistive
effects, realized, e.g., as magnetic read heads in computer hard drives, nonvolative
magnetic random access.
8.4.1 Spin-polarized transport and magnetorisistive effects
Some Scientifics sought an explanation for an unusual behaviour of resistance in
ferromagnetic metals. They realized that at sufficiently low temperatures, where
magnon scattering becomes vanishingly small, electrons of majority and minority
spin, with magnetic moment parallel and antiparallel to the magnetization of a
ferromagnet, respectively, do not mix in the scattering processes. The conductivity
can then be expressed as the sum of two independent and unequal parts for two
different spin projections, the current in ferromagnets is spin polarized.
8.4.1.1 Tunneling magnetoresistance (TMR)
Tunneling measurements play a key role in work on spin-polarized transport. It is
important to study N/F/N junctions, where N was a nonmagnetic metal and F was an
NANOMATERIALS
253
Eu-based ferromagnetic semiconductor. These investigations showed that when
unpolarized current is passed across a ferromagnetic semiconductor, the current
becomes spin-polarized.
Ferromagnet / insulator / superconductor (F/I/S) and Ferromagnet / insulator /
Ferromagnet (F/I/F) junctions were studies these last years. F/I/S have proved that
the tunneling current remains spin polarized even outside of the ferromagnetic region
and was used as a detector of spin polarization of conduction electrons. By analysing
the tunnelling conductance from F/I/S to the F/I/F junctions, Jullire (1975) formulated
a model for a change of conductance between the parallel () and antiparallel ()
magnetization in the two ferromagnetic regions F1 and F2, as depicted in Figure
8-12.
Figure 8-12 : Schematic illustration of electron tunneling in F/I/F tunnel junctions: (a) Parallel and (b)
antiparallel orientation of magnetizations with the corresponding spinresolved density of the d states in
ferromagnetic metals that have exchange spin splitting
ex
. Arrows in the two ferromagnetic regions
are determined by the majority-spin subband. Dashed lines depict spinconserved tunneling. [I. Zutic, J.
Fabian, S. Das Sarma, Spintronics: Fundamentals and applications, Reviews of modern physics,
volume 76, 2004, pp 323-410]
The corresponding tunneling magnetoresistance (TMR) in an F/I/F magnetic tunnel
junction (MTJ) is defined as
=
G
G G
R
R R
R
R
TMR (8. 233)
NANOMATERIALS
254
Where conductance G and resistance R=1/G are labeled by the relative orientations
of the magnetizations in F1 and F2 (it is possible to change the relative orientations,
between and , even at small applied magnetic fields ~10 G). TMR is a particular
manifestation of a magnetoresistance that yields a change of electrical resistance in
the presence of an external magnetic field. Due to spin-orbit interaction, electrical
resistivity changes with the relative direction of the charge current (for example,
parallel or perpendicular) with respect to the direction of magnetization. Assuming
constant tunnelling matrix elements and that electrons tunnel without spin flip,
equation
=
G
G G
R
R R
R
R
TMR (8. 233) yields
2 1
2 1
1
2
P P
P P
TMR
= (8. 234)
Where the polarization P
i
is expressed in terms of the spin-resolved density of states
for majority and minority spin in F
i
.
8.4.1.2 Giant magnetoresistance (GMR)
The prediction of Jullires model illustrates the spinvalve effect: the resistance of a
device can be changed by manipulating the relative orientation of the magnetizations
M
1
and M
2
, in F
1
and F
2
, respectively. Such orientation can be preserved even in the
absence of a power supply, and the spin-valve effect displays the giant
magnetoresistance (GMR) effect can be used for nonvolatile memory applications.
GMR structures are often classified according to whether the current flows parallel
(CIP, current in plane) or perpendicular (CPP, current perpendicular to the plane) to
the interfaces between the different layers, as depicted in Figure 8-13.
Figure 8-13 : Schematic illustration of (a) the current in plane (CIP), (b) the current perpendicular to
the plane (CPP) giant magnetoresistance geometry. [I. Zutic, J. Fabian, S. Das Sarma, Spintronics:
Fundamentals and applications, Reviews of modern physics, volume 76, 2004, pp 323-410]
Most of the GMR applications use the CIP geometry, while the CPP version is easier
to analyze and relates to the physics of the tunneling magnetoresistance effect. The
NANOMATERIALS
255
size of magnetoresistance in the GMR structures can be expressed analogously to
Equation
=
G
G G
R
R R
R
R
TMR (8. 233), where parallel and
antiparallel orientations of the magnetizations in the two ferromagnetic regions are
often denoted by P and AP, respectively (instead of and ). Realization of a
large roomtemperature GMR enabled a quick transition from basic physics to
commercial applications in magnetic recording.
One of the keys to the success of magnetoresistance-based applications is their
ability to control the relative orientation of M
1
and M
2
. The flow of spin-polarized
current can transfer angular momentum from carriers to ferromagnet and alter the
orientation of the corresponding magnetization, even in the absence of an applied
magnetic field. This phenomenon is known as spin-transfer torque. It was also shown
that the magnetic field generated by passing the current through a CPP giant
magnetoresonance device could produce roomtemperature magnetization reversal.
In the context of ferromagnetic semiconductors additional control of magnetization
was demonstrated optically and electrically to perform switching between the
ferromagnetic and paramagnetic states. More information about GMR are given in
8.5 paragraph.
The challenge remains to preserve such spin polarization above room temperature
and in junctions with other materials, since the surface (interface) and bulk magnetic
properties can be significantly different.
8.4.1.3 Applications
While many existing spintronic applications are based on the GMR effects, the
discovery of large room-temperature TMR has renewed interest in the study of
magnetic tunnel junctions, which are now the basis for the several magnetic random-
access memory prototypes. Future generations of magnetic read heads are expected
to use MTJs instead of CIP giant magnetoresonance. To improve the switching
performance of related devices it is important to reduce the junction resistance, which
determines the RC time constant of the MTJ cell. Consequently, semiconductors,
which would provide a lower tunneling barrier than the usually employed oxides, are
being investigated both as the nonferromagnetic region in MTJs and as the basis for
an all-semiconductor junction that would demonstrate large TMR at low
temperatures.
NANOMATERIALS
256
8.4.2 Spin injection
Many materials in their ferromagnetic state can have a substantial degree of
equilibrium carrier spin polarization. However, this alone is usually not sufficient for
spintronic applications, which typically require current flow and/or manipulation of the
nonequilibrium spin (polarization). The importance of generating nonequilibrium spin
is not limited to device applications; it can also be used as a sensitive spectroscopic
tool to study a wide variety of fundamental properties ranging from spin-orbit and
hyperfine interactions to the pairing symmetry of high-temperature superconductors
and the creation of spin-polarized beams to measure parity violation in high-energy
physics.
Nonequilibrium spin is the result of some source of pumping arising from transport,
optical, or resonance methods. Once the pumping is turned off the spin will return to
its equilibrium value. While for most applications it is desirable to have long spin
relaxation times, it has been demonstrated that short spin relaxation times are useful
in the implementation of fast switching.
The Feher effect is based on the principle that the hyperfine coupling between the
electron and nuclear spins, together with different temperatures representing electron
velocity and electron spin populations, is responsible for the dynamical nuclear
polarization.
Figure 8-14 : Pedagogical illustration of the concept of electrical spin injection from a ferromagnet (F)
into a normal metal (N). Electrons flow from F to N: (a) schematic device geometry; (b) magnetization
M as a function of positionnonequilibrium magnetization M (spin accumulation) is injected into a
normal metal; (c) contribution of different spin-resolved densities of states to both charge and spin
transport across the F/N interface. Unequal filled levels in the density of states depict spin-resolved
electrochemical potentials different from the equilibrium value
0
. [I. Zutic, J. Fabian, S. Das Sarma,
Spintronics: Fundamentals and applications, Reviews of modern physics, volume 76, 2004, pp 323-
410]
NANOMATERIALS
257
When a charge current flowed across the F/N junction (Figure 8-14), spin-polarized
carriers in a ferromagnet would contribute to the net current of magnetization entering
the nonmagnetic region and would lead to nonequilibrium magnetization M,
depicted in Figure 8-14(b), with the spatial extent given by the spin diffusion length.
Such a M, which is also equivalent to a nonequilibrium spin accumulation, was first
measured in metals by Johnson and Silsbee (1985, 1988d). In the steady state M is
realized as the balance between spins added by the magnetization current and spins
removed by spin relaxation.
8.4.3 Spin Polarization
Spin polarization not only of electrons, but also of holes, nuclei, and excitations can
be defined as
P
x
= X
s
/X (8. 235)
The ratio of the difference X
s
= X
X
-
, and the sum X
s
= X
+ X
-
of the spin-
resolved components for a particular quantity X. To avoid ambiguity as to what
precisely is meant by spin polarization, both the choice of the spin-resolved
components and the relevant physical quantity X need to be specified.
Conventionally, is taken to be or + (numerical value +1) for spin up, and or -
(numerical value -1) for spin down, with respect to the chosen axis of quantization. In
ferromagnetic etals it is customary to refer to () as carriers with magnetic moment
parallel (antiparallel) to the magnetization or, equivalently, as carriers with majority or
minority spin.
The spin polarization of electrical current or carrier density, generated in a
nonmagnetic region, is typically used to describe the efficiency of electrical spin
injection. Silsbee (1980) suggested that the nonequilibrium density polarization in the
N region, or equivalently the nonequilibrium magnetization, acts as the source of spin
electromotive force (emf) and produces a measurable spin-coupled voltage
M V
S
. Using this concept, also referred to as spin-charge coupling, Silsbee
proposed a detection technique consisting of two ferromagnets F1 and F2 (Figure
8-15) separated by a nonmagnetic region.
NANOMATERIALS
258
Figure 8-15 : Spin injection, spin accumulation, and spin detection: (a) two idealized completely
polarized ferromagnets F1 and F2 with parallel magnetizations are separated by the nonmagnetic
region N; (b) density-of-states diagrams for spin injection from F1 into N, accompanied by the spin
accumulation-generation of nonequilibrium magnetization M. At F2 in the limit of low impedance
(Z=0) electrical spin is detected by measuring the spin-polarized current across the N/F2 interface. In
the limit of high impedance (Z=) spin is detected by measuring the voltage Vs~M developed across
the N/F2 interface; (c) spin accumulation in a device in which a superconductor is occupying the
region between F1 and F2. [I. Zutic, J. Fabian, S. Das Sarma, Spintronics: Fundamentals and
applications, Reviews of modern physics, volume 76, 2004, pp 323-410]
F1 serves as the spin injector (spin aligner) and F2 as the spin detector. This could
be called the polarizer-analyzer method, the optical counterpart of the transmission of
light through two optical linear polarizers.
From Figure 8-15 it follows that the reversal of the magnetization direction in one of
the ferromagnets would lead either to Vs -Vs , in an open circuit (in the limit of
large impedance Z), or to the reversal of charge current j -j, in a short circuit (at
small Z). Spin injection could be detected through the spin accumulation signal as
either a voltage or a resistance change when the magnetizations in F1 and F2 are
changed from parallel to antiparallel alignment.
The generation of nonequilibrium spin polarization has a long tradition in magnetic
resonance methods. However, transport methods to generate carrier spin
polarization are not limited to electrical spin injection. For example, they also include
scattering of unpolarized electrons in the presence of spin-orbit coupling and in
materials that lack inversion symmetry, adiabatic and nonadiabatic quantum spin
pumping, and proximity effects.
NANOMATERIALS
259
8.5 Giant Magnetoresistance (GMR)
8.5.1.1 Origin of the GMR effect
For a macroscopic electrical conductor the resistance R is given by R=P/A, which is
independent of shape, cross section and above all the scattering process involved
therein. Especially, electron-electron and electron-proton scattering process bring the
out-of-equilibrium distribution of accelerated electrons back to the equilibrium
distribution thereby contribution to resistance. This usual picture loses its validity
when the dimension of the material becomes of the order of relevant collision-mean-
free path ( ). Now for an elastic scattering process the conductivity due to spin up
and spin down electron is given by
1 2
= + (8. 236)
For an elastic scattering process
lxxxvi
For the occurrence of the GMR effect it is crucial that the angles between the
magnetization directions of the ferromagnetic layer can be modified by the application
of an applied magnetic field, the electron scattering probability at the interfaces or
within the bulk of the layers is spin-dependent, and in the case of parallel
magnetizations, the layer averaged electron mean-free path for (at least) one spin
direction is larger than the thickness of the non-magnetic spacer layer.
Now, each layer in a GMR- structure acts as a spin selective valve. In case of
parallel-alignment-layer the contribution towards conductance due to spin up electron
is highly leading to higher valve of total conductance. However, for antiparallel
alignments the magnetic layers results in appreciable scattering for electrons in both
spin direction and hence a lower total conductance. If the F layers are much thinner
than the spin up conductance in the parallel magnetization state is limited by
strong spin-dependent scattering in the buffer layer and in the AF layer. The contrast
with spindown conductance can then be enhanced by increasing the F-layer
thickness. It is to be noted that only a part of F-layer which is of the order of
/2 from
the interface responds to the above situation.
NANOMATERIALS
260
R
/
R
s
a
t
(
%
)
free F-layer thickness (nm)
0 20 40 60
8
6
4
2
0
12
10
8
6
4
2
0
(a) 293 K
(b) 5 K
R
/
R
s
a
t
(
%
)
Figure 8-16 : Dependence of the MR ratio on the free magnetic layer thickness, measured at 293 K
(a) and at 5K (b), for (X/Cu/X/Fe
50
Mn
50
) spin valves, grown on 3nm Ta buffer layer on Si(100)
substracts, with X=Co (squares), X=Ni
66
Fe
16
Co
18
(traingles) and (X= Ni
80
Fe
20
) (+), measured at 5K.
Layer thickness: t
Cu
= 3nm (but 2.5 for F=(X= Ni
80
Fe
20
); t
F
(pinned layer)=5nm (but 6nm for F=
Ni
66
Fe
16
Co
18
)
lxxxvii
.
Requirements for F/N combinations with a high GMR ratio are:
the F and N materials possess, for one spin direction, very similar electronic
structures, whereas for the other spin direction the electronic structures are
very dissimilar,
the electron transmission probability (transmission without diffusive scattering)
through the F/N interfaces is large for the type of electrons for which the
electronic structure in the F and N layers are very similar,
and the crystal structures of the F and N materials match very well.
8.5.1.2 Thermal stability
A number of factors control the use and processing of exchange-biased spin-
valve layered structures at elevated temperatures:
NANOMATERIALS
261
Exchanging biasing field decreases with increasing temperature. For some
exchange biasing materials this effect determines the temperature range
within which sensors can be applied.
The induced magnetocrystalline anisotropy of the free layer may be affected
by heating the material.
As multilayers are thermodynamically metastable system the enhanced
diffusion rates during annealing will inversibly degrade the layered structure.
Diffusion may also affect the effective coupling between the layers. There are
three contributions in interlayer coupling:
Coupling via ferromagnetic bridges between the ferromagnetic layers.
The mangetostatic Neel-type coupling.
Interlayer exchange coupling.
8.5.1.3 Factors determining the switching field interval
8.5.1.3.1 Induced magnetic anisotropy
This is believed to be due to the result of a very small degree of divertional pair order
of atom a in the otherwise random solid solution. The anisotropy field depends upon
the composition, the layer thickness, and the material between which the layers are
sandwiched. In case of a GMR spin valve the anisotropy effect is relevant for the
thickness below 15 nm. The effect is most likely due to a region with a decreased
degree of pair order situated in the ferromagnetic layer close to the interface with the
buffer layer.
The result shows that addition of Co to permalloy is unfavourable as far as switching
field interval is concerned. But one of the advantages is that the range of linear
operation of GMR sensor elements is improved. Annealing the sample can affect the
anisotropy field. It is often necessary to anneal the sample after decomposition in a
magnetic field that is directly perpendicular to field direction during growth, in order to
rotate the exchange direction over an angle of 90
o
. This leads to cross anisotropy
configuration. Therefore it is suggestive to anneal with the field parallel to the induced
anisotropy before annealing briefly with perpendicular field. On the other hand strain
results in an additional magnetocrystalline anisotropy due to magnetostriction. So for
sensor element coefficient of magnetostriction should be small. Ni
66
Co
18
Fe
16
alloy is
a good candidate for application in GMR sensor elements.
NANOMATERIALS
262
8.5.1.3.2 Superimposed AMR effect:
lxxxviii
If the current direction is parallel to the easy axis of the free F layer, the
magnetoresistance is given by
2
R(GMR)(1-cos )
R(H)= + R(AMR) sin
2
(8. 237)
The AMR-term represents the resistance change due to the change of angle (
2
)
between the magnetization of the free layer and the current direction. The presence
of AMR effect increases the sensitivity
R
H
R