Contents. Physical Properties. Scalar, Vector. Second Rank Tensor. Transformation. 5 Representation Quadric. 6 Neumann s Principle

Similar documents
Physical Properties. Reading Assignment: 1. J. F. Nye, Physical Properties of Crystals -chapter 1

Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139

PHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review

Elastic Wave Theory. LeRoy Dorman Room 436 Ritter Hall Tel: Based on notes by C. R. Bentley. Version 1.

Tensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)

Cartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components

Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

Chapter 2 Introduction to Phenomenological Crystal Structure

4: birefringence and phase matching

Physics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know

Tensors, and differential forms - Lecture 2

Tensorial and physical properties of crystals

A short review of continuum mechanics

Wave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces

A.1 Appendix on Cartesian tensors

Sample Geometry. Edps/Soc 584, Psych 594. Carolyn J. Anderson

Chapter 8 Vectors and Scalars

Chapter 7. Kinematics. 7.1 Tensor fields

A Primer on Three Vectors

Part II Materials Science and Metallurgy TENSOR PROPERTIES SYNOPSIS

Mathematics for the Theory of Materials: Lectures 1 13

Vectors, metric and the connection

1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis.

Introduction to Polarization

Introduction to tensors and dyadics

Introduction to Matrix Algebra

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

Tensor Analysis in Euclidean Space

Elements of Rock Mechanics

3D and Planar Constitutive Relations

Some elements of vector and tensor analysis and the Dirac δ-function

Numerical methods to compute optical errors due to stress birefringence

Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur

PEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity

Summary of Fourier Optics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

Scalar & Vector tutorial

Semester University of Sheffield. School of Mathematics and Statistics

Math 234. What you should know on day one. August 28, You should be able to use general principles like. x = cos t, y = sin t, 0 t π.

Electro-optics. Chapter 7 Physics 208, Electro-optics Peter Beyersdorf. Document info

Functions of Several Variables

Chapter 9 Electro-Optics

ECE 185 ELECTRO-OPTIC MODULATION OF LIGHT

NIELINIOWA OPTYKA MOLEKULARNA

Rotational motion of a rigid body spinning around a rotational axis ˆn;

The structure of liquids and glasses. The lattice and unit cell in 1D. The structure of crystalline materials. Describing condensed phase structures

PHYSICAL PROPERTIES OF CRYSTALS

Math Review: Vectors and Tensors for Rheological Applications

MSE 201A Thermodynamics and Phase Transformations Fall, 2008 Problem Set No. 7

Symmetry. 2-D Symmetry. 2-D Symmetry. Symmetry. EESC 2100: Mineralogy 1. Symmetry Elements 1. Rotation. Symmetry Elements 1. Rotation.

Index Notation for Vector Calculus

VECTOR AND TENSOR ALGEBRA

Chap. 4. Electromagnetic Propagation in Anisotropic Media

Introduction to Tensor Notation

Phys 201. Matrices and Determinants

VECTORS, TENSORS AND INDEX NOTATION

Classical Mechanics Solutions 1.

example consider flow of water in a pipe. At each point in the pipe, the water molecule has a velocity

The tensors useful for geophysics and estimation of anisotropic polycrystalline physical properties

Cartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components

Linear Algebra Review

11.1 Three-Dimensional Coordinate System

Classical Mechanics in Hamiltonian Form

ASTR 320: Solutions to Problem Set 2

From the SelectedWorks of David Fraivert. David Fraivert. Spring May 8, Available at:

Mathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.

Mathematical Preliminaries

1 Index Gymnastics and Einstein Convention

Derivatives in General Relativity

Tensors - Lecture 4. cos(β) sin(β) sin(β) cos(β) 0

Scattering and Diffraction

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

The Matrix Representation of a Three-Dimensional Rotation Revisited

Lecture 21 Reminder/Introduction to Wave Optics

Name: Lab Partner: Section: In this experiment vector addition, resolution of vectors into components, force, and equilibrium will be explored.

Physics 6303 Lecture 2 August 22, 2018

FoMP: Vectors, Tensors and Fields

Physics I Keystone Institute Technology & Management Unit-II

Vector calculus. Appendix A. A.1 Definitions. We shall only consider the case of three-dimensional spaces.

Math review. Math review

Chap. 2. Polarization of Optical Waves

PLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur

Remark 3.2. The cross product only makes sense in R 3.

Group Representations

Mathematical Properties of Stiffness Matrices

11/29/2010. Propagation in Anisotropic Media 3. Introduction. Introduction. Gabriel Popescu

Fundamentals of Linear Elasticity

PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)

AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES

Mathematics for Graphics and Vision

12. Rigid Body Dynamics I

16.20 Techniques of Structural Analysis and Design Spring Instructor: Raúl Radovitzky Aeronautics & Astronautics M.I.T

Texture, Microstructure & Anisotropy A.D. (Tony) Rollett

A matrix over a field F is a rectangular array of elements from F. The symbol

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

1. Background. is usually significantly lower than it is in uniaxial tension

Course Name : Physics I Course # PHY 107

Transcription:

Physical Properties

Contents 1 Physical Properties 2 Scalar, Vector 3 Second Rank Tensor 4 Transformation 5 Representation Quadric 6 Neumann s Principle

Physical Properties of Crystals - crystalline- translational symmetry, long range order - amorphous- no long range order ex) glass - physical properties amorphous- isotropic crystalline- anisotropic magnitude of physical properties depends on direction

Crystalline vs. Non-crystalline Crystalline materials... atoms pack in periodic, 3D arrays typical of: -metals -many ceramics -some polymers crystalline SiO2 Non-crystalline materials... atoms have no periodic packing occurs for: -complex structures -rapid cooling amorphous" = non-crystalline noncrystalline SiO2

ex) gypsum (CaSO 4 2H 2 O) Thermal Conductivity Monoclinic heated metal tip wax gypsum *ellipsoidal rather than circular

Physical Properties - scalar (zero rank tensor)- non-directional physical quantities, a single number ex) density, temperature - vector (first rank tensor)- magnitude and direction an arrow of definite length and direction ex) mechanical force, electric field, temperature gradient three mutually perpendicular axes Ox 1,Ox 2,Ox 3 components E = [ E1, E2, E3]

Physical Properties - second rank tensor - mechanical analogy central ring-2 pairs of springs at right angle springs on opposite sides are identical but have a different spring constant to perpendicular pair force (cause vector) displacement (effect vector) if a force is applied in a general direction, the displacement will not be in the same direction as the applied force (depends on relative stiffness)

Second Rank Tensor - problem solving 1. find components of the force F in the direction of each of the two springs 2. work out the displacement which each force component would produce parallel to each spring 3. combine two orthogonal displacement to find the resultant displacement

Second Rank Tensor F sinθ y kf 2 sinθ F cosθ θ kf 1 cosθ φ resultant displacement F x 1. force F = [ Fcos θ, Fsin θ] 2. spring constant along x and y are k and k, respectively 3. displacement [ kfcos θ, kfsin θ] 1 2 resultant displacement tanφ = k k 2 1 1 2 tanθ

Second Rank Tensor - consequences 1. in an anisotropic system, the effect vector is not, in general, parallel to the applied cause vector. 2. in two-dimensional example, there are two orthogonal directions along which the effect is parallel to the cause. 3. an anisotropic system can be analyzed in terms of components along these orthogonal principal directions, termed principal axes. along these principal axes, the values of the physical property are termed the principal values.

Second Rank Tensor -in 3-D, general direction- direction cosines, l,m,n -a force F is applied in a general direction resulting in a displacement D at some angle ϕ to F -component of in the direction of D F = Dcosϕ D Dcosϕ D K = = F F K= Kk ( 1, k2, k3) -component of F F F along principal axes Fx = lf, Fy = mf, Fz = nf - component of D along principal axes D = klf, D = k mf, D = k nf x 1 y 2 z 3

Second Rank Tensor D = Dl+ Dm+ Dn K F x y = k Fl + k Fm + k Fn 2 2 2 1 2 3 2 2 2 = ( kl 1 + k2m + k3n ) F DF = = kl + km + kn F 2 2 2 1 2 3 z - variation of a property K with direction - representation surface direction cosine lmn,, x y z l =, m=, n= r r r x y z K= kl + k m + kn = k( ) + k ( ) + k ( ) r r r 2 2 2 2 2 2 1 2 3 1 2 3

2 let 1, 1/ kx + k y + kz = 1 2 2 2 1 2 3 Second Rank Tensor rk= r= K if k, k, k are positive, kx + k y + kz = 1 (ellipsoid) 2 2 2 1 2 3 1 2 3 normal form of the equation of an ellipsoid 2 2 2 x y z + + = 1 ( a, bc. : semiaxes) 2 2 2 a b c representation surface semiaxes: 1 k 1 1 1,, k k 2 3 In any general direction, the radius is equal to the value of 1/ K in that direction.

Second Rank Tensor - electric field E current density i) if conductor is isotropic and obeys Ohm's law j = σ E E = [ E, E, E ] j = [ j, j, j ] 1 2 3 1 2 3 j = σe, j = σe, j = σe 1 1 2 2 3 3 ii) if conductor is anisotropic j = σ E + σ E + σ E 1 11 1 12 2 13 3 j = σ E + σ E + σ E 2 21 1 22 2 23 3 j = σ E + σ E + σ E 3 31 1 32 2 33 3 j isotropic anisotropic

Second Rank Tensor - physical meaning of σ ij if field is applied along x, E = [ E,0,0] 1 1 j = σ E j = σ E j = σ E 1 11 1 2 21 1 3 31 1 - conductivity - nine components speified in a square array σ11 σ12 σ13 σ21 σ22 σ 23 σ31 σ32 σ 33 second rank tensor, components, leading diagonal * the number of subscripts equals the rank of tensor

in general p= [ p, p, p ] q = [ q, q, q ] 1 2 3 1 2 3 p= Tq+ Tq+ Tq 1 11 1 12 2 13 3 p = T q + T q + T q 2 21 1 22 2 23 3 p = T q + T q + T q 3 31 1 32 2 33 3 Second Rank Tensor T T T T T T T T T 11 21 31 21 22 23 31 32 33

Second Rank Tensor p= Tq+ Tq+ Tq= Tq 1 11 1 12 2 13 3 1 j j= 1 3 3 3 p = T q+ T q + T q = T q p = Tq ( i= 1, 2,3) 2 21 1 22 2 23 3 2 j j i ij j j= 1 j= 1 3 31 1 32 2 33 3 3 j j i ij j j= 1 3 p = T q+ T q + T q = T q p = Tq( i= 1, 2,3) -Einstein summation convention: when a letter suffix occurs twice in the same term, summation with respect to that suffix is to be automatically understood. j dummy suffix, p = Tq = Tq i ij j ik k j i free suffix

Second Rank Tensor -in an equation written in this notation, the free suffixs must be the same in all the terms on both sides of the equation: while the dummy suffixs must occur as pairs in each term. ex) A+ BCD = EF ij ik kl lj ik kj i, jfree suffixs kl, dummy suffixs ( C B D = B C D ) kl ik lj ik kl lj -in this book, the range of values of all letter suffixs is 1,2,3 unless some other things is specified.

p= Tq+ Tq+ Tq 1 11 1 12 2 13 3 p = T q + T q + T q 2 21 1 22 2 23 3 p = T q + T q + T q 3 31 1 32 2 33 3 Transformation - q p ( T determine), arbitrary axes chosen j i ij - different set of axes different set of coefficients - both sets of coefficents equally well represent the same physical quantity - there must be some relation between them - when we change the axes of reference, it is only our method of representing the property that changes; the property itself remains the same. T ij

- transformation of axes a change from one set of mutually perpendicular axes to another set with same origin first set: x, x, x, second set: x', x', x' 1 2 3 1 2 3 angular relationship old new 1 2 3 x' a a a 1 11 12 13 x' a a a 2 21 22 23 x' a a a 3 31 32 33 a : cosine of the angle between x ' and x ( a ):matrix ij Transformation x x x i j ij

Direction Cosines, a ij - ( )-nine component- not independent a ij - only three independent quantities are needed to define the transformation. - six independent relation between nine coefficients a + a + a = 1 2 2 2 11 12 13 a a 11 21 a a ik jk ij + a a + a a = = δ 12 22 13 23 0 ( orthogonality relation) Kronecker delta δ = 1 ( i = j) ij 0 ( i j)

30 3 3 60 2 2 σ11 σ12 σ13 cos30 cos 60 cos90 σ σ σ = cos120 cos 30 cos 90 21 22 23 σ31 σ32 σ 33 cos90 cos90 cos 0 3/2 1/2 0 = 1/2 3/2 0 0 0 1 1 1 2 2 2 ( 3 / 2) (1/ 2) 0 1 + + = ( 3/2) ( 1/2) + 1/2 ( 3/2) + 0 0= 0

Transformation - transformation of vector components p p, p, p with respect to x, x, x 1 2 3 1 2 3 p', p', p' with respect to x', x', x' 1 2 3 1 2 3 p' = p cos x x' + p cos x x' + p cos x x' 1 1 1 1 2 2 1 3 3 1 = a p + a p + a p 11 1 12 2 13 p' = a p + a p + a p 2 21 1 22 2 23 3 p' = a p + a p + 3 31 1 32 2 33 3 in dummy suffix notation new in terms of old: p' old in terms of new: a p p 3 = a p i ij j = a p ' i ji j

Transformation - transformation of components of second rank tensor p = Tq with respect to x, x, x i ij j 1 2 3 p' = T' q' with respect to x', x', x' i ij j 1 2 3 p' p q q' ( : in terms of) p' = a p p = Tq q= aq' i ik k k kl l l jl p' = a p= atq= ataq' p' = T' q' T T ' i ik k ik kl l ik kl jl j i ij j aat ij ik jl kl = = a at ij ki lj kl ' j

Transformation T aat aat aat aat ' ij = ik jl kl = ik j 1 k 1+ ik j 2 k 2 + ik j 3 k 3 = aat + aa T + aa T i1 j1 11 i1 j2 12 i1 j3 13 + a a T + a a T + a a T i2 j1 21 i2 j2 22 i2 j3 23 + a a T + a a T + a a T i3 j1 31 i3 j2 32 i3 j3 33

Definition of a Tensor -a physical quantity which, with respect to a set of axes xi, has nine components Tij that transform according to equations ' ij ik jl kl -a second rank tensor- physical quantity existing in its own right, and quite independent of the particular choice of axes -when we change the axes, the physical quantity does not change, but only our method of representing it. ( ) a ij - : array of coefficient relating two set of axes T = aat T = - symmetric ij ji anti-symmetric (skew-symmetric) T T ij = T ji

Representation Quadric - geometrical representation of a second rank tensor - consider the equation - if S xx = 1 S :coefficients ij i j ij S x + S xx + S xx 2 11 1 12 1 2 13 1 3 + S xx+ S x + S xx 2 21 2 1 22 2 23 2 3 + S xx+ S xx + S x = 1 2 31 3 1 32 3 2 33 3 S ij = S ji S x + S x + S x + 2S x x + 2S x x + 2S x x = 1 2 2 2 11 1 22 2 33 3 23 2 3 31 3 1 12 1 2 - general equation of a second-degree surface (quadric) referred to its center as origin

Representation Quadric - transformed to new axes Ox ' x= ax' x = ax' i ki k j lj l Saax' x' = 1 ij ki lj k l S' x' x' = 1 where S' = aas kl k l kl ki lj ij - compared with second rank tensor transformation law T' ij = aat ik jl kl (identical) if S = S ij ji coefficient S ij i of the quadric transform like the components of a symmetrical tensor of the second tank.

Representation Quadric - a representation quadric can be used to describe any symmetrical second-rank tensor, and in particular, it can be used to describe any crystal property which is given by such a tensor. - principal axes principal axes- three directions at right angles such that Sijxx i j = 1 takes the simpler form S x + S x + S x = 1 2 2 2 1 1 2 2 3 3

Representation Quadric S11 S21 S31 S1 0 0 S = S S S 0 S 0 ij 21 22 23 2 S31 S32 S 33 0 0 S 3 S, S, S : principal components 1 2 3 2 2 2 x y z + + = 1 2 2 2 a b c 1 1 1 representation quadric- semi axes,, S S S 1 2 3

Representation Quadric

Representation Quadric - in a symmetric tensor refered to arbitrary axes, the number of independent components is six. - if the tensor is refered to its principal axes, the number of independent components is reduced to three. - the number of degree of freedom is nevertheless still six, for three independent quantities are needed to specify the directions of the axes, and three to fix the magnitudes of the principal components.

Representation Quadric - simplification of equations when referred to pricipal axes p= Sq ( T replaced by symmetric S) i ij j ij ij p = Sq, p = Sq, p = Sq 1 1 1 2 2 2 3 3 3 - for example, consider electrical conductivity j1 = σ1e1, j2 = σ2e2, j3 = σ 3E3 ( σ1, σ2, σ3: principal conductivities) - if E is parallel to Ox1, so E2 = E3 = 0 j2 = j3 = 0 j is parallel to Ox1 - if E = [ E, E,0], 1 2 j1 = σ1e1, j2 = σ2e2, j3 = 0 E and j not parallel

Effect of Crystal Symmetry on Crystal Properties - Neumann's Principle the symmetry elements of any physical properties of a crystal must include the symmetry elements of the point group of the crystal - physical properties may, and often do, possess more symmetry than the point group. - ex1) cubic crystals - optically isotropic physical property (isotropic) possesses the symmetry elements of all the cubic point groups.

Effect of Crystal Symmetry on Crystal Properties - ex2) trigonal system (tourmaline, 3 m) - optical properties (variation of refractive index with direction - indicatrix) indicatrix for 3m- ellipsoid of revolution about triad axis (optic axis) ellipsoid of revolution- vertical triad axis three vertical planes of symmetry (extra- center of symmetry, other symmetry elements) - the symmetry of a physical property a relation between certain measurable quantities associated with the crystal

Effect of Crystal Symmetry on Crystal Properties - all second-rank tensor properties are centrosymmetric. p = Tq i ij j - p = T ( q ) T : unchanged i ij j ij - symmetric second-rank tensor- 6 independent components - symmetry of crystal reduces the number of independent components - consider representation quadric for symmetric second rank tensor

Fick's first law Ji = Dij x j o at 1150 C Anisotropic Diffusion of Ni in Olivine c ex) Ni diffusion in olivine((mg,fe) SiO, orthorhombic) 2 4 D =4.40x10 cm /s, D =3.35x10 cm /s, D =124x10 cm /s -14 2-14 2-14 2 x y z 1 1 1 a:b:c= : : = 0.48: 0.55: 0.09 D D D x y z

Magnitude of a Property in a Given Direction - defintion in general, if pi = Sq ij j, the magnitude Sof the property [ Sij ] in a certain direction is obtained by applying q in that direction and measuring p / q, where p is the componet of p parallel to q - ex) electrical conductivity the conductivity σ in the direction of E is defined to be the component of j parallel to E divided by E, that is, j / E

Magnitude of a Property in a Given Direction - analytical expression (i) referred to principal axes direction cosine: l1, l2, l3 E = [ le 1, le 2, le 3 ] j = [ σ11 le, σ2le 2, σ33 le] component of j parallel to E j = l 2 1 σ E+ l σ E+ l σ E 2 2 1 2 2 3 3 magnitude of conductivity in the direction l i σ = l σ + l σ + l σ 2 2 2 1 1 2 2 3 3

Magnitude of a Property in a Given Direction - analytical expression (ii) referred to general axes l : direction cosine of E referred to general axes i Ei = Eli component of j parallel to E je / E in suffix notation je / E conductivity in the direction je E σ= σ l l i i ij j σ = = 2 2 ij i j σ EE E i l i i i

Geometrical Properties of Representation Quadric - length of the radius vector let Pbe a general point on the ellipsoid: σ xx = 1 direction cosines of OP: l x = rl where OP = r σ 1 ( σ σ ) 2 r ijll i j = = ijll i j σ = r r = 2 1/ 1/ σ i i i ij i j special cases- radius vectors in the directions of semi-axes of lengths 1/ σ, 1/ σ,1/ σ 1 2 3

Geometrical Properties of Representation Quadric - in general, any symmetric second-rank tensor property S ij 2 1/ 1/ S = r r = S - the length r of any radius vector of representation quadric is equal to the reciprocal of square root of magnitude S of the property in that direction

Geometrical Properties of Representation Quadric - radius-normal property Ox i principal axes of σ ij E= [ le 1, le 2, le 3 ] j= [ σ11 le, σ2le 2, σ33 le] direction cosines of j are proportional to σ l, σ l, σ l 11 2 2 33 2 if P is a point on σ1x1+ σ x such that OP is parallel to E P = ( rl, rl, rl ) where OP = r 1 2 3 tangent plane at P rl σ x + rl σ x + rl σ x = 1 1 1 1 2 2 2 3 3 3 + σ x = 1 2 2 2 2 3 3

Geometrical Properties of Representation Quadric - radius-normal property normal at P has direction cosines proportional to lσ, l σ, l σ 1 1 2 2 3 3 hence normal at P is parallel to if pi = Sq ij j, the direction of pfor a given q may be found by first drawing, parallel to q a radius vector OP of the representation quadric, and then taking the normal to the quadric at P. j

Heat Flow in a Crystal - a point source of heat on the face of a crystal of a tetragonal mineral (uniaxial) isothermal surface: (001) plane- circle, (100) plane- ellipse heat flow-radially away from P, thernal gradient-normal to isothermal surface heat flow thermal gradient (in general, not parallel) -the resistivity along the c-axis is less than that normal to it

T - resistivity = -rh ij x i Heat Flow in a Crystal dt / dr cosθ r= (long rod experiment) h T - conductivity hi = kij x h cosθ k= (thin flat plate) dt / dr j j

Heat Flow in a Crystal