Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Chapter 8 – Symmetry in Crystal Physics – p. 1 9. Symmetry in Crystal Physics 9.1. Description Physical Properties of Crystals by Tensors r r v v Isotropic material: P = χE (el. polarization, el. susceptibility, el. field); P || E physical properties, which describe the relation between vectors are generally described by scalars, implying that both vectors are parallel. r Anisotropic materials: P; E are not necessarily parallel. Suitable description? Example: Electric susceptibility χ r r v v Pv = χ vv Ev and Ph = χ hv Ev In this 2D example, two scalars are required to describe the effect of the vertical field strength, same for the horizontal field strength. General case in three dimensional space: Chapter 8 – Symmetry in Crystal Physics – p. 2 χ r 11 P = χ 21 χ 31 χ12 χ 22 χ32 χ13 v χ 23 E χ33 or P1 = χ11 E1 + χ12 E2 + χ13 E3 P2 = χ 21 E1 + χ 22 E2 + χ 23 E3 P3 = χ 31 E1 + χ 32 E2 + χ33 E3 Abbreviated version: 3 Pi = ∑ χ ij E j ; i = 1,2,3 j =1 Einstein’s Notation: Pi = χ ij E j ; i, j = 1, 2,3 Ultimate abbreviation for people, who don’t like summation-symbols: Summation over all indices, which occur twice in the same term. The dielectric susceptibility is said to be a second rank tensor. General definition of a tensor: A tensor is defined as a set of 3r components which describes a physical quantity. In addition a tensor has to have special transformation properties upon a change of the basis of the coordinate system. These transformation properties will be discussed in section 9.4. r is called the rank of a tensor: • Scalar: 1 = 30 component, zero-rank tensor • Vector: 3 = 31 component, 1st-rank tensor • Matrix: 9 = 32 components, 2nd-rank tensor • general: 3r components, r-rank tensor General relation between two physical properties A and B: Bijk ... n = aijk ... npqr...u Apqr...u ; i, j , k ,..., n, p, q, r ,..., u = 1, 2,3 with • A: f-rank tensor • B: g-rank tensor Chapter 8 – Symmetry in Crystal Physics – p. 3 • a: (f+g)-rank tensor Another example: stress, strain, elastic modulus and Hooke’s law The stress tensor: Body acted on by external forces: state of stress. We consider a simple case: (1) case of homogeneous stress: stress is independent of position in the body, (2) static equilibrium, (3) no body forces or torques. Stress is described by 9 stress components σ = dF / dA (stress, area of body element, force acting on area), i.e. second rank tensor σ ij : The fact that the stress components represent a tensor is to be proved. Proves can be fount in the textbooks on crystal physics (see literature) . Obviously, the stress has 3 normal components σ11 , σ 22 , σ33 and 6 sheer components σ12 , σ13 , σ 23 , σ 21 , σ 31 , σ 32 . Condition (2) (static equilibrium) immediately requires, that there are certain relationships between the sheer components (otherwise there would be a torque acting on the body): σ12 = σ 21 , σ 23 = σ32 , σ 31 = σ13 , the stress tensor is a symmetric second rank tensor: σ ij = σ ji . (Latter relationship still holds in the case of inhomogeneous stress or body forces, see textbooks). Chapter 8 – Symmetry in Crystal Physics – p. 4 The strain tensor: Strain in 1 dimension: Strain: ε = dx dl In a body, along one direction, of the material in all three directions might change. Therefore: Extension: eij = ∂xi ∂li The extension contains a symmetrical part and a antisymmetrical part. The antisymmetrical part describes a rotation of the body, the symmetrical part describes the actual deformation: Strain tensor εij (symmetrical part of the extension): ε11 ε12 ε 13 ε12 ε22 ε23 e11 ε13 1 ε23 = (e12 + e21 ) 2 ε33 1 (e13 + e31 ) 2 1 (e12 + e21 ) 2 e22 1 (e23 + e32 ) 2 1 (e13 + e31 ) 2 1 (e23 + e32 ) 2 e33 Thus, the strain tensor is also a symmetric second rank tensor. With the previous definitions we can reformulate Hook’s law for a general solid: εij = sijkl σ kl The elastic modulus sijkl is a 4th rank tensor, in principle containing 81 components. Chapter 8 – Symmetry in Crystal Physics – p. 5 9.2 Examples of Tensors Representing Physical Properties Similar as in the case of the discussed examples many physical properties can be described in tern of tensors of different rank. Here some examples (from E. Hartmann, Introduction to Crystal Physics): 9.3. Polar and Axial Vectors and Tensors We consider a basis transformation from an old coordinate system êi to a new coordinate system eˆi ' , described by a transformation matrix a ji . Both sets of basis vectors are chosen orthonormal, yielding a −1 ji = a ij . Chapter 8 – Symmetry in Crystal Physics – p. 6 Normally, a vector pi in the old coordinate system would be described in the new coordinate system as: p'i = aij p j (note the summation over j according to the Einstein notation). Such a (normal) vector is also called a polar vector. In physics, there are some vectors, however, which have slightly different transformation properties. These vectors are usually connected to a definition involving the vector product: r e1 r r r c = a × b = a1 r e2 r e3 a2 a3 b1 b2 b3 r r (the vector product describes a vector of length c = ab sin ∠ab , is perpendicular to a a and b , r r r and which forms a right handed system a, b , c ). We consider a basis transformation, which changes the coordinate system from right-handed to left handed or vice versa, such as the inversion i: Chapter 8 – Symmetry in Crystal Physics – p. 7 In general, a basis transformation of a polar vector generates a change of sign, if the transformation changes the hand. A change of hands is deduced from the determinant aij of the transformation matrix • aij = 1 : Transformation leaves hand of the axes unchanged (rotation) • aij = −1 : Transformation changes hand of the axes (inversion, reflection) Transformation of a axial vector: p' i = aij a ij p j (note the summation over j according to the Einstein notation). Physical examples of axial vector (1st rank axial tensor) are the angular momentum r r r r r r ( L = mr × v ) or the magnetic flux density ( F = Qv × B ). In general a tensor a connecting two properties A and B via Bijk ... n = aijk ... npqr...u Apqr...u is axial if either A or B are axial. In every other case it is polar. 9.4. Transformation Properties of Tensors Again, we consider a basis transformation to a new coordinate system. A scalar (0-rank tensor) does not change upon this operation: T '= T (0-rank tensor) A (polar) vector pi is described in the new coordinate system as: p'i = aij p j (or T 'i = aijT j 1st-rank tensor) For a (polar) 2nd-rank tensor pk = Tkl ql we obtain with p'i = aik pk and q'l = alj q j or ql = a jl q ' j : p'i = aikTkl a jl q' j = T 'ij q' j or Chapter 8 – Symmetry in Crystal Physics – p. 8 T 'ij = aik a jlTkl (2nd-rank tensor) General transformation laws for tensors upon basis transformation: Tensor Rank Polar Tensor Axial Tensor 0 T 'ij = Tkl 1 T 'i = aijT j T 'i = aij aijT j 2 T 'ij = aik a jlTkl T 'ij = aij aik a jlTkl 3 T 'ijk = ail a jm aknTlmn T 'ijk = aij ail a jm aknTlmn n T 'ijk ... n = aip a jq akr ...anuTpqr...u T 'ijk ... n = aij aipa jq akr ...anuTpqr...u 9.5. Intrinsic Symmetry of Physical Properties A k-rank tensor has up to 3n different components. However, the number of independent components is much smaller in most cases, either due to intrinsic symmetries of the physical property described (this section) or due to the crystal symmetry (section 9.6). 9.5.1 Symmetry by Definition Some properties are defined such that the corresponding tensors exhibit an inner symmetry. Examples: Strain tensor εij : 2nd-rank, symmetric second rank tensor, 6 independent components Stress tensor σ ij : 2nd-rank, symmetric second rank tensor, 6 independent components Elastic modulus sijkl : 4th rank tensor with sijkl = s jikl = sijlk = s jilk , 36 independent components. 9.5.2 Equilibrium Properties and Thermodynamic Arguments Chapter 8 – Symmetry in Crystal Physics – p. 9 For tensors describing equilibrium properties, thermodynamic relations significantly reduce the number of independent components. Example: We consider elastic, electric and magnetic work plus heat exchange in a crystal (stress, strain, el. field, el. polarization, magn. field, magn. polarization, temperature, entropy): ∂ε dεij = ij dσ kl + ∂σ kl ∂P dPi = i dσ kl + ∂σ kl ∂J dJ i = i dσ kl + ∂σ kl ∂S dσ kl + dS = ∂σ kl ∂εij ∂Ek ∂Pi ∂Ek ∂J i ∂Ek ∂S ∂Ek dEk dEk dEk dEk ∂ε + ij dH k ∂H k ∂P + i dH k ∂H k ∂J + i dH k ∂H k ∂S dH k + ∂H k ∂ε + ij dT ∂T ∂P + i dT ∂T ∂J + i dT ∂T ∂S + dT ∂T According to the 1st and 2nd law of thermodynamics we obtain (reversible process): dU = dw + dq = σ ij dεij + E k dPk + H l dJ l + TdS We change the set of independent variables by introducing a Gibb’s free energy: dG = U − σ ij εij − Ek Pk − H l J l − TS yielding dG = − εij dσ ij − Pk dEk + J l dH l + SdT From comparison with the total differential of ∂G dσ ij + ∂G dG = ∂E k ∂σ ij ∂G ∂G dEk + dH l + dT ∂ H ∂ T l we obtain ∂G ∂σ = −εij ; ij ∂G = − Pk ; ∂E k ∂G ∂H k = − J k ; ∂G = −S ∂T From the commutability of the second derivatives (Schwartz theorem) it follows for the dielectric susceptibility χ jk : Chapter 8 – Symmetry in Crystal Physics – p. 10 - ∂P ∂ 2G 2 ∂P = − ∂ G = j = χ jk χ kj = k = − ∂E ∂E ∂E ∂E j ∂E j ∂Ek k j k Therefore, the dielectric susceptibility tensor is symmetric. Similar arguments hold for the magnetic susceptibility ψ lk = (∂J l / ∂H k ) and the elastic modulus sijkl = (∂εij / ∂σ kl ) . Moreover it follows that the tensors describing direct and reciprocal effects are identical. A an example we consider the piezoelectrical (stress -> el. polarization) and the reverse piezoelectrical effect (el. field -> strain): ∂Pk 2 2 ∂ε = − ∂ G = − ∂ G = ij = d kij ∂σ ∂σ ∂E ∂E ∂σ ∂E ij ij k k ij k Similar relations are found for the piezomagnetic effect qlij = (∂J l / ∂σ ij ) = (∂εij / ∂H l ) , the λlk = (∂J l / ∂Ek ) = (∂Pk / ∂H l ) , magneto-electrical polarization piezocaloric-effect α lk = (∂ε ij / ∂T ) = (∂S / ∂σ ij ) , pyroelectric and electrocaloric effect pk = (∂Pk / ∂T ) = (∂S / ∂Ek ) and pyromagnetic and thermal dilatation magneto-caloric and effect ml = (∂J l / ∂T ) = (∂S / ∂H l ) . Thus the above set of equations can be simplified significantly: dεij = sijkl dσ kl + d kij dEk + qlij dH l + α ij dT dPk = d kij dσij + χ kl dEk + λlk dH k + pk dT dJ l = qlij dσ ij + λlk dEk + ψ lm dH k + ml dT dS = α ij dσ ij + pk dEk + ml dH k + cd ln T 9.5.3. Transport Properties and Onsager’s Principle In irreversible thermodynamics, transport processes are described by sets of corresponding thermodynamic forces X i and fluxes ji , chosen such that σ& = ji X i corresponds to the rate of entropy production. For the corresponding linear system of transport equations ji = Lij X j , Onsager’s reciprocity relation states that Chapter 8 – Symmetry in Crystal Physics – p. 11 Lij = L ji . Example: Electric and heat transport 1 ∂Φ (el. potential); flux: electrical flux density j q k T ∂xk • Electric transport: force X q k = • Linear transport equations: j q i = Lq ik 1 ∂Φ ∂Φ (classical definition: j q i = −σ ik ) T ∂xk ∂xk è Onsager relation: Lik = Lki or σ ik = σ ki (conductivity tensor is symmetric) ∂(1 / T ) ; flux: heat flux density j Q k ∂xk • Heat transport: force X Q k = • Definition of Peltier and Seebeck effect: j Q i = −kik dT dΦ + πik j q k with j q k = σ kl dxk dxl j q i = σ ik E β k + σ ik with dΦ dT with E β k = − βkl dxk dxl d (1 / T ) d (1 / T ) dT dT = = T −2 it follows: dx dT dx dx j Q i = −kik (−T 2 ) d (1 / T ) 1 dΦ + πik σ klT dxk T dxl j q i = −σ ik βkl (−T 2 ) d (1 / T ) 1 dΦ + σ ikT dxl T dxk è Onsager relation: βklT = π kl (relation between Seebeck and Peltier effect) 9.6 Neumann’s Principle: Crystal Symmetry and Tensor Symmetry GT ⊇ GC with ( ⊇ : subgroup or equal) GT : symmetry group of tensor T GC : symmetry group of crystal C Chapter 8 – Symmetry in Crystal Physics – p. 12 The group of symmetry elements of a physical property includes all symmetry elements of the crystal, i.e. • Polar tensor: Tijk ... n = a ip a jq a kr ...anuT pqr... u • Axial tensor: Tijk ... n = aij aip a jq a kr ...a nuT pqr... u Must be fulfilled for all transformation matrices aij corresponding to symmetry operations of the point group of the crystal (crystal class) (remark: in general, it is not necessary to test all symmetry operations, but only the set of generating operations of the point group). Neumann’s principle may reduce the number of independent components of a tensor or may require that some tensor elements must vanish. Example 1: We consider the el. susceptibility tensor χ ij defined via Pi = χ ij E j (polar symmetrical tensor of rank 2). We consider a crystal which belongs to crystal class 4 (crystallographic) or C4 (Schoenfliess). The generating element is C4. aC 4 0 − 1 0 = 1 0 0 . 0 0 1 We obtain: Chapter 8 – Symmetry in Crystal Physics – p. 13 χ11 χ12 χ 13 χ13 0 − 1 0 χ11 χ12 χ13 0 1 0 χ 23 = 1 0 0 χ12 χ 22 χ 23 − 1 0 0 χ 33 0 0 1 χ13 χ 23 χ 33 0 0 1 0 − 1 0 − χ12 χ11 χ13 = 1 0 0 − χ 22 χ12 χ 23 0 0 1 − χ 23 χ13 χ 33 χ 22 − χ12 − χ 23 = − χ12 χ11 χ13 − χ χ13 χ 33 23 χ12 χ 22 χ 23 Thus: (a) χ11 = χ 22 (b) χ12 = − χ 21 ⇒ χ12 = 0 (c) χ13 = χ 23 = − χ 23 ⇒ χ13 = χ 23 = 0 Therefore, the electric polarizability tensor of a crystal belonging to class 4 has only two independent components: χ11 χ12 χ 13 χ12 χ 22 χ 23 χ13 χ11 χ 23 = 0 χ 33 0 0 χ11 0 0 0 χ 33 Example 1: We consider the process of second harmonic / sum frequency generation, which is described by the second order hyperpolarizability χ ijk , a polar tensor of rank 3 ( Pi = χ ij E j + χ ijk E j E k + ... ; for a fields Eω , Eω ' oscillating with frequency ω , ω ' , χ ijk generates frequency components with frequency 2 ω , 2 ω ' , ω + ω ' , ω − ω ' ): Chapter 8 – Symmetry in Crystal Physics – p. 14 We consider a crystal class which contains the inversion i described by the transformation matrix 0 −1 0 − 1; i = j a = 0 − 1 0 or a i ij = −δ ij = 0; i ≠ j 0 0 − 1 i Neumann’s principle requires that χ ijk = a i il a i jm a i kn χ lmn = ∑ ( −δ il )( −δ jm )( −δ kn )χ lmn lmn = ( −1)( −1)(−1) χ lmn = − χ lmn Therefore, χ ijk = 0 , i.e. second harmonic / sum frequency generation is forbidden in crystals / media with inversion symmetry. Note: In media with inversion symmetry, all processes described by polar tensors of odd order and axial tensors of even order are forbidden! 9.8 Contracted Matrix Notation In literature, a special matrix notation is often used to simplify the representation of higher rank tensors. Example: Hook’s law: εij = sijkl σ kl The elastic modulus is a 4th rank tensor with a maximum of 81 components, but symmetry with respect to suffixes i,j and k,l reduces the number of free components to 36. We contract the index pairs i,j to a single index m and k,l to n according to the following rule: (1,1) → 1; (2,2) → 2; (3,3) → 3; (2,3) → 4; (3,1) → 5; (1,2) → 6 Using the definition εij ; m = 1,2,3 εm = 2 εij ; m = 4,5,6 Chapter 8 – Symmetry in Crystal Physics – p. 15 σ n = σ ij ; n = 1,2,3,4,5,6 s mn sij ; 2 s ; ij = 2 sij ; 4 sij ; m, n = 1,2,3 m = 1,2,3; n = 4,5,6 m = 1,2,3; n = 4,5,6 m, n = 4,5,6 Hook’s law is expressed as: ε 1 c11 ε 2 c21 ε c 3 = 31 ε 4 c41 ε c 5 51 ε c 6 61 c12 c13 c14 c15 c22 c23 c24 c25 c32 c33 c34 c35 c42 c43 c44 c45 c52 c53 c54 c55 c62 c63 c64 c65 c16 σ 1 c26 σ 2 c36 σ 3 c46 σ 4 c56 σ 5 c66 σ 6 or ε n = s nmσ m Remarks: (1) The coefficients in the definition of the contracted notations are necessary to take into account the reduced number of terms in the summations. (2) This is only a notation trick! The transformation properties remain unchanged. 9.9 Value of a Physical Property in a Given Direction The value T of a physical property described by a 2nd rank tensor via pi = Tij q j in the r r r direction of q is defined as T = p|| / q , where p|| is the component of p parallel to q . Example: Chapter 8 – Symmetry in Crystal Physics – p. 16 - r r r Electrical conductivity σ , electrical field: E = En with n : unit vector in field direction. r r r r E E j Ei E j|| = j ⋅ n = j ⋅ = ji i = σ ij E E E E E With σ nr = j|| / E it follows: σ nr = σ ij ni n j This equation can be used to derive tensor components from a physical measurement or vice versa. 9.10 Geometrical Representation: The Representation Quadratic For the important group of symmetric 2nd-rank tensors (as an example we again consider the conductivity ji = σ ij E j ), there is a simple geometrical way of representation, the so called representation quadratic defined as: σ ij xi x j = 1 or σ 11 x12 + σ 22 x 22 + σ 33 x32 + 2σ 23 x2 x3 + 2σ 31 x3 x1 + 2σ 12 x1 x2 = 1 This is a second degree surface, in most cases it corresponds to an ellipsoid (tensor ellipsoid): It can be shown that upon basis transformation the representation quadratic behaves like a symmetrical 2nd-rank tensor. Thus the transformation properties of the tensor can be derived from a (graphical) inspection of the transformation properties of the representation quadratic. Chapter 8 – Symmetry in Crystal Physics – p. 17 We can choose a basis transformation to a coordinate system, in which σ ij is diagonal: 0 σ '1 0 σ ' = 0 σ '2 0 0 0 σ ' 3 The directions of this special set of basis vectors are referred to as the principle axes of the tensor. In the new basis, the tensor ellipsoid points along the coordinate axes: the representation quadratic takes a simple form σ '11 x '12 +σ ' 22 x ' 22 +σ ' 33 x ' 32 = 1 and physical equations involving the tensor become particularly simple: j '1 = σ '1 E '1 , j'2 = σ '2 E '2 , j'3 = σ '3 E '3 . The representation quadratic has two important geometrical properties (example: el. conductivity ji = σ ij E j ): r (a) The radius r in a given direction n is related to the physical property σ nr in this r r r direction via σ nr = 1 / r 2 . The resulting component of j parallel to E = En is j|| = σ nr E . (Proof: From σ ij xi x j = 1 with xi = rni , we obtain r 2σ ij ni n j = 1 . With the result from section 9.9 ( σ nr = σ ij ni n j ), we obtain σ nr = 1 / r 2 .) Chapter 8 – Symmetry in Crystal Physics – p. 18 r (b) The direction of j is along the normal of the representation quadratic at the endpoint of the radius (without proof, see textbooks). Example 2: Optical properties of crystals Optical properties of an isotropic medium: • r r D = ε 0εE with ε 0ε dielectrical permittivity or r r η0ηD = E with η0η dielectrical impermeability • Maxwell equations ( µ = 1 ): velocity of electromagnetic wave v = refractive index: n = c =c η ε c 1 = ε = v η Optical properties of an anisotropic medium: • Di = ε 0ε ij E j with ε 0ε ij dielectrical permittivity tensor or η0η ij D j = Ei with η0η ij dielectrical impermeability tensor (both 2nd-rank symmetric) • Maxwell equations ( µ = 1 ): (for proof see textbooks, e.g. Nye) In general, two plane polarized waves with different velocity may be propagated along one direction (double refraction). Chapter 8 – Symmetry in Crystal Physics – p. 19 Graphical representation: We consider the representation quadratic of the relative dielectric impermeability tensor ηij , the so called indicatrix (note: (a) principal axes are chosen; (b) ni are called the principal refractive indices, but the refractive index is not a tensor!): x12 x 22 x 32 + + =1 η x + η 2 x + η x = 1 or n12 n22 n32 2 1 1 2 2 2 3 3 The indicatrix has the following important property (for lengthy proofs see textbooks): We consider wave propagation along 0P. The central section through the indicatrix, perpendicular to the propagation direction is an ellipse. The axes of this ellipse represent the r two polarisation of D and the semi-axes 0A and 0B are identical to the refractive indices n A and n A for the two waves. From these properties of the indicatrix and Neumann’s principle we can immediately classify all crystal classes with respect to their optical properties: (a) Optical anaxial crystals: cubic (classes 23 , m3 , 432 , 432 , m3m ): • indicatrix is a sphere (several Cn / Sn axes with n>2) Chapter 8 – Symmetry in Crystal Physics – p. 20 - ⇒ no double refraction in any direction v (b) Optical uniaxial crystals: teragonal ( 4 , 4 , 4 m , 422 , 4mm , 42m , 4 m mm ), trigonal ( 3 , 3 , 32 , 3m , 3m ), hexagonal ( 6 , 6 , 6 m , 622 , 6mm , 6m 2 , 6 m mm ) • indicatrix is an ellipsoid of revolution along principal symmetry axes (one Cn / Sn axes with n>2) ⇒ no double refraction along principal symmetry axis (one optical axis) (a) Optical biaxial crystals: triclinic ( 1 , 1 ), monoclinic ( 2 , m , 2 m ), orthorhombic ( 222 , mm 2 , mmm ): • indicatrix is a triaxial ellipsoid (no Cn / Sn axes with n>2) Chapter 8 – Symmetry in Crystal Physics – p. 21 - ⇒ no double refraction along two axes (optical axes) 9.11 Curie’s Principle Often, crystal properties are considered under some external influence (electrical field, strain, etc.). Here, Curie’s principle states: GC~ = GC ∩ G E with ( ∩ greatest common subgroup) GC~ : symmetry group of crystal C under external influence of E GC : symmetry group of crystal C GE : symmetry group of external influence E Example: electro-optical and photoelastic effects We consider the change of the relative dielectric impermeability tensor under the influence of an electric or stress and expand it in term of a power series: ηij = ηij0 + rijk E k + k ijkl E k E l + π ijkl σ kl + ... • rijk : linear electro-optical tensor → linear electro-optical effect (Pockels effect) • k ijkl : quadratic electro-optical tensor → quadratic electro-optical effect (Kerr effect) • π ijkl : piezoptical tensor → photoelastic effect As an example, we consider the Pockels effect in ADP (ammonium-dihydrogen-phosphate). The crystal class of ADP is 42m , belonging to the tetragonal system. The crystal is optical Chapter 8 – Symmetry in Crystal Physics – p. 22 uniaxial. For wave propagation along the principal symmetry axis, no double diffraction occurs. We apply an electric field along the principal symmetry axis. The group of the electric field vector contains an ∞–fold rotation axis and ∞ mirror planes containing the axis (group ∞m) . As a set of common symmetry elements two mirror planes and a C2 axis survive. We identify the common subgroup mm2 , which belongs to the orthorhombic crystal system, i.e. an optical biaxial group. → The optical uniaxial ADP crystal becomes biaxial if an electric field along z (principal symmetry axis) is applied and double refraction occurs along z. Electro-optical and photoelastic effects are usually employed in optical elements which change the direction of polarization or modulate the intensity.