Index Notation for Vector Calculus

Similar documents
VECTORS, TENSORS AND INDEX NOTATION

Vector analysis and vector identities by means of cartesian tensors

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

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

A Primer on Three Vectors

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

Vectors and Matrices Notes.

The Matrix Representation of a Three-Dimensional Rotation Revisited

NIELINIOWA OPTYKA MOLEKULARNA

Basic mathematics for nano-engineers (II)

1.4 LECTURE 4. Tensors and Vector Identities

PHY481: Electromagnetism

PHY481: Electromagnetism

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

Module 4M12: Partial Differential Equations and Variational Methods IndexNotationandVariationalMethods

1.3 LECTURE 3. Vector Product

Homework 7-8 Solutions. Problems

This document is stored in Documents/4C/vectoralgebra.tex Compile it with LaTex. VECTOR ALGEBRA

Vector analysis. 1 Scalars and vectors. Fields. Coordinate systems 1. 2 The operator The gradient, divergence, curl, and Laplacian...

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

Classical Mechanics Solutions 1.

Physics 6303 Lecture 2 August 22, 2018

Introduction to relativistic quantum mechanics

Tensors, and differential forms - Lecture 2

Introduction to Tensor Notation

MP2A: Vectors, Tensors and Fields

1 Index Gymnastics and Einstein Convention

FoMP: Vectors, Tensors and Fields

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

Mathematical Preliminaries

1. Rotations in 3D, so(3), and su(2). * version 2.0 *

Supplementary Notes on Mathematics. Part I: Linear Algebra

Tensor Calculus (à la Speedy Gonzales)

Quantum Physics II (8.05) Fall 2004 Assignment 3

Classical Mechanics in Hamiltonian Form

Math review. Math review

Vector and Tensor Calculus

Introduction to Vector Spaces

Valiantly Validating Vexing Vector Verities

VECTOR AND TENSOR ALGEBRA

1 Vectors and Tensors

EN221 - Fall HW # 1 Solutions

Physics 221A Fall 2017 Appendix E Introduction to Tensor Analysis

Rigid Bodies. November 5, 2018

Dr. Tess J. Moon, (office), (fax),

2.20 Fall 2018 Math Review

Tensor Calculus. arxiv: v1 [math.ho] 14 Oct Taha Sochi. October 17, 2016

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

A.1 Appendix on Cartesian tensors

Electrodynamics. 1 st Edition. University of Cincinnati. Cenalo Vaz

Gravity theory on Poisson manifold with R-flux

1.13 The Levi-Civita Tensor and Hodge Dualisation

Vectors. September 2, 2015

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

Mathematical Preliminary

Introduction to Tensor Calculus and Continuum Mechanics. by J.H. Heinbockel Department of Mathematics and Statistics Old Dominion University

Professor Terje Haukaas University of British Columbia, Vancouver Notation

Mathematics for the Theory of Materials: Lectures 1 13

ECS 178 Course Notes QUATERNIONS

( ) by D n ( y) 1.1 Mathematical Considerations. π e nx Dirac s delta function (distribution) a) Dirac s delta function is defined such that

Chapter 7. Kinematics. 7.1 Tensor fields

Computation of the mtx-vec product based on storage scheme on vector CPUs

We are already familiar with the concept of a scalar and vector. are unit vectors in the x and y directions respectively with

The groups SO(3) and SU(2) and their representations

Lecture 19 (Nov. 15, 2017)

Physics 342 Lecture 26. Angular Momentum. Lecture 26. Physics 342 Quantum Mechanics I

INTEGRALSATSER VEKTORANALYS. (indexräkning) Kursvecka 4. and CARTESIAN TENSORS. Kapitel 8 9. Sidor NABLA OPERATOR,

Tensor Analysis in Euclidean Space

Faculty of Engineering, Mathematics and Science School of Mathematics

BLAS: Basic Linear Algebra Subroutines Analysis of the Matrix-Vector-Product Analysis of Matrix-Matrix Product

2 Tensor Notation. 2.1 Cartesian Tensors

Chapter 2 The Continuum Equations

Linear Algebra (Review) Volker Tresp 2017

Discrete Applied Mathematics

u z u y u x ChE 342 Vectors 1 VECTORS Figure 1 Basic Definitions Vectors have magnitude and direction: u = i u x + j u y + k u z (1)

GROUP THEORY PRIMER. New terms: tensor, rank k tensor, Young tableau, Young diagram, hook, hook length, factors over hooks rule

A Brief Introduction to Tensor

MATRICES. a m,1 a m,n A =

Electromagnetism HW 1 math review

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

Conservation Laws. Chapter Conservation of Energy

Lecture notes on introduction to tensors. K. M. Udayanandan Associate Professor Department of Physics Nehru Arts and Science College, Kanhangad

Classical Mechanics. (2 nd Edition) Cenalo Vaz. University of Cincinnati

NOTES ON DIFFERENTIAL FORMS. PART 3: TENSORS

Vectors. (Dated: August ) I. PROPERTIES OF UNIT ANSTISYMMETRIC TENSOR

02 - tensor calculus - tensor algebra tensor calculus

On Expected Gaussian Random Determinants

AMS526: Numerical Analysis I (Numerical Linear Algebra)

Formalism of the Tersoff potential

Citation for published version (APA): Halbersma, R. S. (2002). Geometry of strings and branes. Groningen: s.n.

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.

Lecture Notes on The Mechanics of Elastic Solids

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

Properties of Matrix Arithmetic

%#, for x = a & ' 0, for x $ a. but with the restriction that the area it encloses equals 1. Therefore, whenever it is used it is implied that

Appendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk

TENSOR ALGEBRA. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

611: Electromagnetic Theory II

A Brief Introduction to Tensors

Week 15-16: Combinatorial Design

Transcription:

Index Notation for Vector Calculus by Ilan Ben-Yaacov and Francesc Roig Copyright c 2006 Index notation, also commonly known as subscript notation or tensor notation, is an extremely useful tool for performing vector algebra. Consider the coordinate system illustrated in Figure 1. Instead of using the typical axis labels x, y, and z, we use x 1, x 2, and x 3,or x i i =1, 2, 3 The corresponding unit basis vectors are then ê 1, ê 2, and ê 3,or ê i i =1, 2, 3 The basis vectors ê 1, ê 2, and ê 3 have the following properties: ê 1 ê 1 =ê 2 ê 2 =ê 3 ê 3 =1 (1) ê 1 ê 2 =ê 1 ê 3 =ê 2 ê 3 =0 (2) x 3 a 3 a e 3 e 1 e2 a 2 x 2 a 1 x 1 Figure 1: Reference coordinate system. 1

2 Index Notation We now introduce the Kronecker delta symbol δ ij. δ ij has the following properties: { 0 i j δ ij = i, j =1, 2, 3 (3) 1 i = j Using Eqn 3, Eqns 1 and 2 may be written in index notation as follows: ê i ê j = δ ij i, j =1, 2, 3 (4) In standard vector notation, a vector A may be written in component form as A =A x î + A y ĵ + A zˆk (5) Using index notation, we can express the vector A as A = A 1 ê 1 + A 2 ê 2 + A 3 ê 3 (6) = A i ê i i=1 Notice that in the expression within the summation, the index i is repeated. Repeated indices are always contained within summations, or phrased differently a repeated index implies a summation. Therefore, the summation symbol is typically dropped, so that A can be expressed as A =A i ê i A i ê i (7) i=1 This repeated index notation is known as Einstein s convention. Any repeated index is called a dummy index. Since a repeated index implies a summation over all possible values of the index, one can always relabel a dummy index, i.e. A = A i ê i = A j ê j = A k ê k etc. A 1 ê 1 + A 2 ê 2 + A 3 ê 3 (8)

Index Notation 3 The Scalar Product in Index Notation We now show how to express scalar products (also known as inner products or dot products) using index notation. Consider the vectors a and b, which can be expressed using index notation as a = a 1 ê 1 + a 2 ê 2 + a 3 ê 3 = a i ê i (9) b = b1 ê 1 + b 2 ê 2 + b 3 ê 3 = b j ê j Note that we use different indices (i and j) for the two vectors to indicate that the index for b is completely independent of that used for a. We will first write out the scalar product a b in long-hand form, and then express it more compactly using some of the properties of index notation. ( ) a b = a i ê i b j ê j i=1 j=1 = [(a i ê i ) (b j ê j )] i=1 j=1 = [a i b j (ê i ê j )] (commutative property) i=1 j=1 = (a i b j δ ij ) (from Eqn3) i=1 j=1 Summing over all values of i and j, we get a b = a 1 b 1 δ 11 + a 1 b 2 δ 12 + a 1 b 3 δ 13 + a 2 b 1 δ 21 + a 2 b 2 δ 22 + a 2 b 3 δ 23 + a 3 b 1 δ 31 + a 3 b 2 δ 32 + a 3 b 3 δ 33 = a 1 b 1 δ 11 + a 2 b 2 δ 22 + a 3 b 3 δ 33 = a 1 b 1 + a 2 b 2 + a 3 b 3 = a i b i i=1 = a i b i = a j b j = a k b k

4 Index Notation Doing this in a more compact notation gives us a b =(a i ê i ) (b j ê j ) = a i b j ê i ê j = a i b j δ ij = a i b i = a 1 b 1 + a 2 b 2 + a 3 b 3 Notice that when we have an expression containing δ ij, we simply get rid of the δ ij and set i = j everywhere in the expression. Example 1: Kronecker delta reduction Reduce δ ij δ jk δ ki : Here we can see that δ ij δ jk δ ki = δ ik δ ki (remove δ ij, set j = i everywhere) = δii (remove δ ik, set k = i everywhere) = δ ii = 1=1+1+1=3 i=1 i=1 δ ii =3 (Einstein convention implied) (10) Note also that δ ij δ jk = δ ik (11) Example 2: r and ˆr in index notation (a) Express r using index notation. r = x 1 ê 1 + x 2 ê 2 + x 3 ê 3 = x i ê i

Index Notation 5 (b) Express ˆr using index notation. ˆr = r r = r ( r r) 1/2 = x i ê i (x j x j ) 1/2 (c) Express a ˆr using index notation. a ˆr = a r r = a i x i (x j x j ) 1/2 The Cross Product in Index Notation Consider again the coordinate system in Figure 1. Using the conventional righthand rule for cross products, we have ê 1 ê 1 =ê 2 ê 2 =ê 3 ê 3 =0 ê 1 ê 2 =ê 3 ê 2 ê 1 = ê 3 ê 2 ê 3 =ê 1 ê 3 ê 2 = ê 1 (12) ê 3 ê 1 =ê 2 ê 1 ê 3 = ê 2 To write the expressions in Eqn 12 using index notation, we must introduce the symbol ɛ ijk, which is commonly known as the Levi-Civita tensor, the alternating unit tensor, or the permutation symbol (in this text it will be referred to as the permutation symbol). ɛ ijk has the following properties: ɛ ijk =1 ɛ ijk = 1 ɛ ijk =0 if (ijk) is an even (cyclic) permutation of (123), i.e. ɛ 123 = ɛ 231 = ɛ 312 =1 if (ijk) is an odd (noncyclic) permutation of (123), i.e. ɛ 213 = ɛ 321 = ɛ 132 = 1 if two or more subscripts are the same, i.e. ɛ 111 = ɛ 112 = ɛ 313 =0 etc.

6 Index Notation Hence, we may rewrite the expressions in Eqn 12 as follows: ê 1 ê 2 = ɛ 123 ê 3 ê 2 ê 1 = ɛ 213 ê 3 ê 2 ê 3 = ɛ 231 ê 1 ê 3 ê 2 = ɛ 321 ê 1 (13) ê 3 ê 1 = ɛ 312 ê 2 ê 1 ê 3 = ɛ 132 ê 2 Now, we may write a single generalized expression for all the terms in Eqn 13: ê i ê j = ɛ ijk ê k (14) Here ɛ ijk ê k ɛ ijk ê k (k is a dummy index). That is, this works because k=1 ê 1 ê 2 = ɛ 12k ê k = ɛ 12k ê k = ɛ 121 ê 1 + ɛ 122 ê 2 + ɛ 123 ê 3 =ê 3 The same is true for all of the other expressions in Eqn 13. Note that ê i ê i = ɛ iik ê k =0, since ɛ iik for all values of k. ɛ ijk is also given by the following formula. k=1 ɛ ijk = 1 (i j)(j k)(k i) i, j, k =1, 2, 3 (15) 2 This is a remarkable formula that works for ɛ ijk if you do not want to calculate the parity of the permutation (ijk). Also note the following property of ɛ ijk : ɛ ijk = ɛ jik = ɛ kji i.e. switching any two subscripts reverses the sign of the permutation symbol (or in other words ɛ ijk is anti-symmetric). Also, ɛ ijk = ɛ kij = ɛ jki i.e. cyclic permutations of the subscripts do not change the sign of ɛ ijk. These

Index Notation 7 properties also follow from the formula in Eqn 15. Now, let s consider the cross product of two vectors a and b, where a = a i ê i b = bj ê j Then a b =(a i ê i ) (b j ê j )=a i b j ê i ê j = a i b j ɛ ijk ê k Thus we write for the cross product: a b = ɛ ijk a i b j ê k (16) All indices in Eqn 16 are dummy indices (and are therefore summed over) since they are repeated. We can always relabel dummy indices, so Eqn 16 may be written equivalently as a b = ɛ pqr a p b q ê r Returning to Eqn 16, the kth component of a b is ( ) a b = ɛ ijka i b j k where now only i and j are dummy indices. Note that the cross product may also be written in determinant form as follows: ê 1 ê 2 ê 3 a b = a 1 a 2 a 3 (17) b 1 b 2 b 3 The follwoing is a very important identity involoving the product of two permutation symbols. δ il δ im δ in ɛ ijk ɛ lmn = δ jl δ jm δ jn (18) δ kl δ km δ kn

8 Index Notation The proof of this identity is as follows: If any two of the indices i, j, k or l, m, n are the same, then clearly the lefthand side of Eqn 18 must be zero. This condition would also result in two of the rows or two of the columns in the determinant being the same, so therefore the right-hand side must also equal zero. If (i, j, k) and (l, m, n) both equal (1,2,3), then both sides of Eqn 18 are equal to one. The left-hand side will be 1 1, and the right-hand side will be the determinant of the identity matrix. If any two of the indices i, j, k or l, m, n are interchanged, the corresponding permutation symbol on the left-hand side will change signs, thus reversing the sign of the left-hand side. On the right-hand side, an interchange of two indices results in an interchange of two rows or two columns in the determinant, thus reversing its sign. Therefore, all possible combinations of indices result in the two sides of Eqn 18 being equal. Now consider the special case of Eqn 18 where n = k. In this case, the repeated index k implies a summation over all values of k. The product of the two permutation symbols is now δ il δ im δ ik ɛ ijk ɛ lmk = δ jl δ jm δ jk δ kl δ km δ kk (note δ kk =3) =3δ il δ jm 3δ im δ jl + δ im δ jk δ kl (19) δ ik δ jm δ kl + δ ik δ jl δ km δ il δ jk δ km =3δ il δ jm 3δ im δ jl + δ im δ jl δ il δ jm + δ im δ jl δ il δ jm (from Eqn 11) Or finally ɛ ijk ɛ lmk = δ il δ jm δ im δ jl (20)

Index Notation 9 Eqn 20 is an extremely useful property in vector algebra and vector calculus applications. It can also be expressed compactly in determinant form as ɛ ijk ɛ lmk = δ il δ im δ jl δ jm (21) The cyclic property of the permutation symbol allows us to write also ɛ ijk ɛ klm = δ il δ jm δ im δ jl To recap: ê i ê j = δ ij and a b = a i b i ê i ê j = ɛ ijk ê k and a b = ɛ ijk a i b j ê k These relationships, along with Eqn 20, allow us to prove any vector identity. Example 3: The scalar triple product Show that a ( b c) = c ( a b)= b ( c a) a ( b c) = (a i ê i ) (ɛ jkm b j c k ê m ) = ɛ jkm a i b j c k (ê i ê m ) = ɛ jkm a i b j c k δ im or = ɛ jki a i b j c k a ( b c) =ɛ ijk a i b j c k From our permutation rules, it follows that a ( b c) = ɛ ijk a i b j c k = ɛ kij c k a i b j = c ( a b) = ɛ jki b j c k a i = b ( c a)